BesEmcHit.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//      BOOST --- BESIII Object_Oreiented Simulation Tool                    //
//---------------------------------------------------------------------------//
// Descpirtion: EMC detector
// Author: Fu Chengdong
// Created: Sep 4, 2003
// Comment:
//---------------------------------------------------------------------------//
//
#include "EmcSim/BesEmcHit.hh"
 
#include "G4UnitsTable.hh"
#include "G4ios.hh"
#include <iomanip>
 
G4Allocator<BesEmcHit>      BesEmcHitAllocator;
G4Allocator<BesEmcTruthHit> BesEmcTruthHitAllocator;
 
using namespace std;
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
BesEmcHit::BesEmcHit() {
  TotalEdepCrystal        = 0.;
  TotalTrackLengthCrystal = 0.;
  EdepCrystal             = 0.;
  TrackLengthCrystal      = 0.;
  EdepCasing              = 0.;
  PositionCrystal         = G4ThreeVector( 0, 0, 0 );
  TimeCrystal             = 0.;
  PartId                  = 0;
  NumTheta                = 0;
  NumPhi                  = 0;
  trackIndex              = 0;
  g4Index                 = 0;
  momentum                = G4ThreeVector( 0, 0, 0 );
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
BesEmcHit::~BesEmcHit() {}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
BesEmcHit::BesEmcHit( const BesEmcHit& right ) : G4VHit() {
  EdepCrystal        = right.EdepCrystal;
  TrackLengthCrystal = right.TrackLengthCrystal;
  EdepCasing         = right.EdepCasing;
  PositionCrystal    = right.PositionCrystal;
  TimeCrystal        = right.TimeCrystal;
  PartId             = right.PartId;
  NumTheta           = right.NumTheta;
  NumPhi             = right.NumPhi;
  trackIndex         = right.trackIndex;
  g4Index            = right.g4Index;
  momentum           = right.momentum;
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
const BesEmcHit& BesEmcHit::operator=( const BesEmcHit& right ) {
  EdepCrystal        = right.EdepCrystal;
  TrackLengthCrystal = right.TrackLengthCrystal;
  EdepCasing         = right.EdepCasing;
  PositionCrystal    = right.PositionCrystal;
  TimeCrystal        = right.TimeCrystal;
  PartId             = right.PartId;
  NumTheta           = right.NumTheta;
  NumPhi             = right.NumPhi;
  trackIndex         = right.trackIndex;
  g4Index            = right.g4Index;
  momentum           = right.momentum;
  return *this;
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
int BesEmcHit::operator==( const BesEmcHit& right ) const {
  return ( this == &right ) ? 1 : 0;
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
void BesEmcHit::Draw() {}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
void BesEmcHit::Print( G4int verboseLevel ) {
  if ( verboseLevel > 0 )
    G4cout << "Hit in crystal:" << NumTheta << "," << NumPhi << G4endl
           << "Energy deposited:" << G4BestUnit( EdepCrystal, "Energy" ) << G4endl;
  if ( verboseLevel > 1 )
    G4cout << "Hit time        :" << G4BestUnit( TimeCrystal, "Time" ) << G4endl
           << "    position    :" << G4BestUnit( PositionCrystal, "Length" ) << G4endl;
  if ( verboseLevel > 2 )
    G4cout << "Track length    :" << G4BestUnit( TrackLengthCrystal, "Length" ) << G4endl;
}
 
void BesEmcHit::Print() {
  G4cout << "time: " << TimeCrystal << " edep: " << EdepCrystal << G4endl;
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
BesEmcTruthHit::BesEmcTruthHit() {
  m_trackIndex   = -1;
  m_g4TrackId    = -1;
  m_hitEmc       = -1;
  m_PDGCode      = 0;
  m_PDGCharge    = 0;
  m_particleName = G4String();
  m_edep         = 0;
  m_time         = 0;
  m_momentum     = G4ThreeVector( 0, 0, 0 );
  m_position     = G4ThreeVector( 0, 0, 0 );
  m_hitMap.clear();
}
 
BesEmcTruthHit::~BesEmcTruthHit() {}
 
BesEmcTruthHit::BesEmcTruthHit( const BesEmcTruthHit& right ) : G4VHit() { *this = right; }
 
const BesEmcTruthHit& BesEmcTruthHit::operator=( const BesEmcTruthHit& right ) {
  m_hitMap.clear();
  std::map<Identifier, G4double>::const_iterator iHitMap;
  if ( this != &right )
  {
    for ( iHitMap = right.Begin(); iHitMap != right.End(); iHitMap++ )
    { Insert( iHitMap->first, iHitMap->second ); }
 
    m_identify     = right.m_identify;
    m_trackIndex   = right.m_trackIndex;
    m_g4TrackId    = right.m_g4TrackId;
    m_hitEmc       = right.m_hitEmc;
    m_PDGCode      = right.m_PDGCode;
    m_PDGCharge    = right.m_PDGCharge;
    m_particleName = right.m_particleName;
    m_edep         = right.m_edep;
    m_time         = right.m_time;
    m_momentum     = right.m_momentum;
    m_position     = right.m_position;
  }
 
  return *this;
}
 
void BesEmcTruthHit::Print() {
  G4cout << "Id: " << m_identify << "\tTrack Index: " << m_trackIndex
         << "\tG4 Track Id: " << m_g4TrackId << "\tHit Emc: " << m_hitEmc
         << "\tTotal Energy: " << m_edep << "\nPDGCode: " << m_PDGCode
         << "\tCharge: " << m_PDGCharge << "\tParticle Name: " << m_particleName
         << "\nGloble Time: " << m_time << "\tMomentum: " << m_momentum.mag()
         << "\tPosition: " << m_position << G4endl;
 
  std::map<Identifier, G4double>::iterator iHitMap;
  for ( iHitMap = m_hitMap.begin(); iHitMap != m_hitMap.end(); iHitMap++ )
  { G4cout << iHitMap->first << "\t" << iHitMap->second << G4endl; }
}
 
std::map<Identifier, G4double>::const_iterator BesEmcTruthHit::Begin() const {
  return m_hitMap.begin();
}
 
std::map<Identifier, G4double>::const_iterator BesEmcTruthHit::End() const {
  return m_hitMap.end();
}
 
std::map<Identifier, G4double>::const_iterator BesEmcTruthHit::Find( Identifier id ) const {
  return m_hitMap.find( id );
}
 
G4double BesEmcTruthHit::GetEHit( Identifier id ) { return m_hitMap[id]; }
 
void BesEmcTruthHit::AddEHit( Identifier id, G4double energy ) {
  if ( energy > 0 ) m_hitMap[id] += energy;
}
 
void BesEmcTruthHit::Insert( Identifier id, G4double energy ) {
  if ( energy > 0 ) m_hitMap[id] = energy;
}
 
G4int BesEmcTruthHit::Size() const { return m_hitMap.size(); }