LCDG4McPartManager Class Reference

#include <LCDG4McPartManager.hh>

Collaboration diagram for LCDG4McPartManager:

[legend]List of all members.

Public Member Functions
	LCDG4McPartManager ()
	LCDG4McPartManager (std::map< G4LogicalVolume *, G4int > *m, int dbflg=0)
virtual	~LCDG4McPartManager ()
void	SetDetTypeMap (std::map< G4LogicalVolume *, G4int > *m)
void	SetDebugFlag (int dbflg)
void	PreTracking (const G4Track *aTrack)
void	PostTracking (const G4Track *aTrack, const G4TrackingManager *aTrackMgr)
LCDMcPart *	CheckPrimaryParticle (const G4Track *aTrack)
void	MakeLCDMcPart (const G4Track *aTrack)
const G4int &	GetShowerFlag (G4int parentTrkID) const
void	Clear ()
int	GetDetType (G4LogicalVolume *lvol) const
int	GetNumMCparticles () const
void	PrintMcPart ()
void	PrintTrkid2MctagMap ()
G4int	trkidToMctag (G4int trkid) const
LCDMcPart *	GetMcPartFromTag (G4int mctag) const
LCDMcPart *	GetMcPartFromTrkid (G4int trkid) const
bool	CheckLastProcess (G4String proc) const
bool	IsShowerProcess () const
bool	IsInOutput (G4int trkid) const
bool	IsInOutput () const
bool	IsShowered ()
bool	IsLeft () const
bool	IsInteracted () const
bool	IsDecayed () const
bool	IsStopped () const
void	GetDaughterCandidates (const G4int dau1, const G4int dau2, const G4int pdgidref, std::vector< LCDMcPart * > &cands)
void	SetCalorimeterEntrancePoint (const G4ThreeVector &pos, const G4ThreeVector &mom)
Static Public Member Functions
LCDG4McPartManager *	GetPointer ()
Public Attributes
std::vector< LCDMcPart * >	MCparticles
std::vector< G4int >	trkids
Private Member Functions
void	PreAssignedDecay (G4TrackVector *SecTrkV)
void	EraseAllPreassignedDescendants (G4int m_mctag)
void	AssignTagToTrack (const G4Track *trk, const G4int tag)
void	AssignTagToTrack (const G4Track *trk, const LCDMcPart *mcp)
void	AssociateFromComplexLists (std::list< G4Track * > &trks, std::list< LCDMcPart * > &mcps)
void	FindSingleBestAssociation (std::vector< G4Track * > &trks2, std::vector< LCDMcPart * > &mcps2)
void	DumpLists (std::list< G4Track * > &trks, std::list< LCDMcPart * > &mcps)
void	DumpVectors (std::vector< G4Track * > &trks, std::vector< LCDMcPart * > &mcps)
Private Attributes
int	m_nmcprim
int	m_debug
G4int	m_trkid
int	m_mctag
G4Track *	curTrack
const G4ParticleDefinition *	dpart
const G4StepPoint *	postStep
const G4VProcess *	lastProc
std::map< G4LogicalVolume *, G4int > *	m_pdettype
std::map< G4int, G4int >	m_trkid2McTag
std::map< G4int, G4int >	m_trkid2dettype
std::map< G4int, G4int >	m_showerFlag
Static Private Attributes
LCDG4McPartManager *	_me = NULL

---

Constructor & Destructor Documentation

LCDG4McPartManager::LCDG4McPartManager (   )

Definition at line 41 of file LCDG4McPartManager.cc. References _me, m_trkid2dettype, and m_trkid2McTag. : m_nmcprim(0), m_debug(0), m_pdettype(0) { m_trkid2McTag.clear(); m_trkid2dettype.clear(); if ( _me == NULL ) _me = this; }

LCDG4McPartManager::LCDG4McPartManager (  std::map< G4LogicalVolume *, G4int > *  m, int  dbflg = 0 )

Definition at line 49 of file LCDG4McPartManager.cc. References _me, m_trkid2dettype, and m_trkid2McTag. : m_nmcprim(0), m_debug(dbflg), m_pdettype(m) { m_trkid2McTag.clear(); m_trkid2dettype.clear(); if ( _me == NULL ) _me = this; }

virtual LCDG4McPartManager::~LCDG4McPartManager (   )  [inline, virtual]

Definition at line 19 of file LCDG4McPartManager.hh. {};

---

Member Function Documentation

void LCDG4McPartManager::AssignTagToTrack (  const G4Track *  trk, const LCDMcPart *  mcp )  [inline, private]

Definition at line 93 of file LCDG4McPartManager.hh. References AssignTagToTrack(), and LCDMcPart::GetParticleTag(). { AssignTagToTrack( trk, mcp->GetParticleTag() ); }

void LCDG4McPartManager::AssignTagToTrack (  const G4Track *  trk, const G4int  tag )  [private]

Definition at line 822 of file LCDG4McPartManager.cc. References LCDG4Trk2McPart::GetMcPartTag(), LCDMcPart::getTrkID(), m_debug, m_trkid2McTag, MCparticles, LCDG4Trk2McPart::SetMcPartTag(), and LCDMcPart::setTrkID(). Referenced by AssignTagToTrack(), and PreTracking(). { assert(tag>0); // set the tag in UserTrkInfo LCDG4Trk2McPart* trkInfo = (LCDG4Trk2McPart*)trk->GetUserInformation(); assert(trkInfo); // trk->trackID may be zero here, so better check it G4int trkid = trk->GetTrackID(); if(trkid==0) { // if trkid=0, set trkinfo tag to -tag (incomplete connection) trkInfo->SetMcPartTag(-tag); if(m_debug & 0x10) { G4cout<< "Negate tag (incomplete association): tag=" << trkInfo->GetMcPartTag() << G4endl; } } else { // connection can be made complete trkInfo->SetMcPartTag(tag); // set trkid for this LCDMcPart -- tags start at 1, not zero LCDMcPart* mcp = MCparticles[tag-1]; G4int trkid2 = mcp->getTrkID(); if(trkid2) G4cout <<"AssignTagToTrack trkid2!=0: trkid(trac)="<< trkid <<", trkid(mcpart)="<< trkid2 << G4endl; assert( trkid2==0 || trkid2==trkid ); mcp->setTrkID(trkid); // save tag in trkid->tag map if(m_debug & 0x10) { G4cout<< "AssignTag2Trk: Set trkid2McTag["<<trkid<<"]="<<tag<< G4endl; } m_trkid2McTag[trkid] = tag; } }

void LCDG4McPartManager::AssociateFromComplexLists (  std::list< G4Track * > &  trks, std::list< LCDMcPart * > &  mcps )  [private]

Definition at line 863 of file LCDG4McPartManager.cc. References LCDMcPart::GetParticleID(). Referenced by PreAssignedDecay(). { list<G4Track*>::const_iterator itrk; list<LCDMcPart*>::const_iterator imcp; // debugging printout ========== // G4cout<<"\n*** ComplexList: sizes(trk,mcp)="<< trks.size() <<' '<< mcps.size() << G4endl; // DumpLists(trks,mcps); //=========== // loop over mcps to build the vector of distinct PDGIDs vector<int> pdgids; for(imcp=mcps.begin(); imcp!=mcps.end(); imcp++) { pdgids.push_back( (*imcp)->GetParticleID() ); } // sort sort(pdgids.begin(), pdgids.end()); // eliminate repeated ids vector<int>::iterator newend = unique(pdgids.begin(), pdgids.end()); int newsize = newend-pdgids.begin(); // loop over unique pdgids vector<LCDMcPart*> mcps2; vector<G4Track*> trks2; for( int iid=0; iid<newsize; iid++ ) { // for each pdgid int idref = pdgids[iid]; // build subvector of LCDMcPart* for( imcp=mcps.begin(); imcp!=mcps.end(); imcp++) { LCDMcPart* pmcp = *imcp; if( pmcp->GetParticleID() == idref ) { mcps2.push_back(pmcp); } } assert( mcps2.size()>0 ); // build subvector of G4Track* for( itrk=trks.begin(); itrk!=trks.end(); itrk++) { G4Track* ptrk = *itrk; if( ptrk->GetDefinition()->GetPDGEncoding() == idref ) trks2.push_back(ptrk); } // find all possible matches while( trks2.size()>0 && mcps2.size()>0 ) { // G4cout<< "\nSubvectors:"; // DumpVectors(trks2,mcps2); // non trivial case if(mcps2.size()>1 || trks2.size()>1) FindSingleBestAssociation(trks2, mcps2); AssignTagToTrack( trks2[0], mcps2[0] ); trks.erase( find(trks.begin(), trks.end(), trks2[0]) ); mcps.erase( find(mcps.begin(), mcps.end(), mcps2[0]) ); trks2.erase( trks2.begin() ); mcps2.erase( mcps2.begin() ); } // All possible matches have been found, there might be some trks or // mcps left. Reset both before going to next PDGID mcps2.clear(); trks2.clear(); } }

bool LCDG4McPartManager::CheckLastProcess (  G4String  proc  )  const

Definition at line 609 of file LCDG4McPartManager.cc. References lastProc. Referenced by IsDecayed(), IsLeft(), IsShowerProcess(), IsStopped(), and PostTracking(). { return lastProc->GetProcessName().contains(proc); }

LCDMcPart * LCDG4McPartManager::CheckPrimaryParticle (  const G4Track *  aTrack  )

Definition at line 297 of file LCDG4McPartManager.cc. References LCDMcPart::Get4Momentum(), LCDMcPart::GetParticleID(), LCDMcPart::GetParticleTag(), GeV, m_debug, m_nmcprim, and MCparticles. Referenced by PreTracking(). { LCDMcPart* result = NULL; // trkids for track and parent G4int trkid = aTrack->GetTrackID(); G4int pdgidref = aTrack->GetDefinition()->GetPDGEncoding(); // Look for this particle in the vector<LCDMcPart*> G4LorentzVector emomref = aTrack->GetDynamicParticle()->Get4Momentum(); if(m_debug & 0x20) { G4int ptrkid = aTrack->GetParentID(); // always zero for input particles G4cout << " CheckPrimaryParticle called, MCparticles.size() = " << MCparticles.size()<<", trkID,parID=" << trkid <<',' << ptrkid <<", idref ="<< pdgidref <<", Emomref="<< emomref/GeV << G4endl; } // loop over primary MCparticles backwards vector<LCDMcPart*>::iterator iter = MCparticles.begin(); LCDMcPart* pmcp; G4int imcp = 0; G4int pdgid = 0; G4LorentzVector emomvec; vector<G4double> diffs; vector<LCDMcPart*> mcps; for( ; imcp<m_nmcprim; imcp++, iter++) { pmcp = *iter; if(!pmcp) continue; // require same PDGID if( pmcp->GetParticleID() != pdgidref ) continue; // build a list of differences on vector 3-momenta G4LorentzVector delta = pmcp->Get4Momentum() - emomref; diffs.push_back( delta.vect().mag2() ); mcps.push_back( pmcp ); } // Return false if no matching LCDMcPart was found if(diffs.size()==0) return result=NULL; // sort vector vector<G4double> sorted(diffs); sort(sorted.begin(), sorted.end()); G4double dmin = sorted[0]; unsigned int isave= find( diffs.begin(), diffs.end(), dmin ) - diffs.begin(); if(m_debug & 0x20) { G4double dmin2 = sorted.size()>1 ? sorted[1] : 9e10; G4cout <<"sorted: d(Min,Min2) =" << dmin <<' '<< dmin2 <<", size="<< sorted.size() <<", id="<< mcps[isave]->GetParticleID() << G4endl; } // criteria for finding a match if( dmin < 1e-6 ) { result = mcps[isave]; pdgid = result->GetParticleID(); assert( pdgid==pdgidref ); } if(m_debug & 0x20) { G4cout<<"LCDG4McPartManager::CheckPrimaryParticle:" <<" Primary trkid="<< trkid <<", tag="<< result->GetParticleTag() <<", pdgid="<< pdgid << G4endl; } // // set particle starting position // G4cout<< "Set starting position: "<< aTrack->GetVertexPosition() << G4endl; // result->SetPosition( aTrack->GetVertexPosition() ); // // set also parent ending position // LCDMcPart* parent = GetMcPartFromTag( result->GetParentTag() ); // parent->SetTermPosition( aTrack->GetVertexPosition() ); return result; }

void LCDG4McPartManager::Clear (   )  [inline]

Definition at line 34 of file LCDG4McPartManager.hh. References m_nmcprim, m_trkid2dettype, m_trkid2McTag, MCparticles, and trkids. Referenced by LCDG4EventAction::EndOfEventAction(), and LCDG4readStdFile::MakeMcParts(). { m_nmcprim = 0; m_trkid2McTag.clear(); m_trkid2dettype.clear(); for(unsigned int i=0; i<MCparticles.size(); ++i) { if(MCparticles[i]) { // G4cout<< "deleting MCparticles[" << i << "]\n"; delete MCparticles[i]; } } MCparticles.clear(); trkids.clear(); }

void LCDG4McPartManager::DumpLists (  std::list< G4Track * > &  trks, std::list< LCDMcPart * > &  mcps )  [private]

Definition at line 1000 of file LCDG4McPartManager.cc. { list<G4Track*>::const_iterator itrk; list<LCDMcPart*>::const_iterator imcp; G4cout << "\n List of mcps:"; for(imcp=mcps.begin(); imcp!=mcps.end(); imcp++) { G4cout<<" ("<< (*imcp) <<')'; if( *imcp ) G4cout<<' '<< (*imcp)->GetParticleID() <<','; } G4cout <<"\n List of trks:"; for(itrk=trks.begin(); itrk!=trks.end(); itrk++) { G4cout<<' '<< (*itrk)->GetDefinition()->GetParticleName(); } G4cout << G4endl; }

void LCDG4McPartManager::DumpVectors (  std::vector< G4Track * > &  trks, std::vector< LCDMcPart * > &  mcps )  [private]

Definition at line 1018 of file LCDG4McPartManager.cc. { vector<G4Track*>::const_iterator itrk; vector<LCDMcPart*>::const_iterator imcp; G4cout << "\n List of mcps:"; for(imcp=mcps.begin(); imcp!=mcps.end(); imcp++) { G4cout<<' '<< (*imcp)->GetParticleID(); } G4cout <<"\n List of trks:"; for(itrk=trks.begin(); itrk!=trks.end(); itrk++) { G4cout<<' '<< (*itrk)->GetDefinition()->GetParticleName(); } G4cout << G4endl; }

void LCDG4McPartManager::EraseAllPreassignedDescendants (  G4int  m_mctag  )  [private]

Definition at line 1035 of file LCDG4McPartManager.cc. References LCDMcPart::GetDaughterTags(), MCparticles, and LCDMcPart::SetGenStatus(). Referenced by PostTracking(). { LCDMcPart* parmcp = MCparticles[m_mctag-1]; if(!parmcp) { G4cout<<"EraseAllPreassDescends: parmcp=0: m_mctag="<< m_mctag<<G4endl; return; } // reset parent generator status code to FINAL STATE int FinalState = 1; parmcp->SetGenStatus(FinalState); // recursively erase all descendants G4int dau1,dau2; parmcp->GetDaughterTags(dau1,dau2); if(dau1>0) { for(int itag=dau1; itag<=dau2; itag++) { // erase grandchildren EraseAllPreassignedDescendants(itag); // erase this child delete MCparticles[itag-1]; MCparticles[itag-1] = NULL; G4cout<< "*** Erasing tag="<< itag << G4endl; } } }

void LCDG4McPartManager::FindSingleBestAssociation (  std::vector< G4Track * > &  trks2, std::vector< LCDMcPart * > &  mcps2 )  [private]

Definition at line 937 of file LCDG4McPartManager.cc. References LCDMcPart::Get4Momentum(). { // loop over McParts vector<G4Track*> trks3; vector<LCDMcPart*> mcps3; vector<G4double> diffs; for( unsigned int i=0; i<mcps2.size(); i++) { LCDMcPart* pmcp = mcps2[i]; G4LorentzVector emomref = pmcp->Get4Momentum(); // loop over G4Tracks for( unsigned int j=0; j<trks2.size(); j++) { G4Track* ptrk = trks2[j]; G4LorentzVector deltaPvec = ptrk->GetDynamicParticle()->Get4Momentum() - emomref; // store deltaPs, mcps and trks diffs.push_back( deltaPvec.vect().mag2() ); mcps3.push_back(mcps2[i]); trks3.push_back(trks2[j]); } } // Smallest element in diffs is the best match in terms of mom components vector<G4double> sorted(diffs); sort(sorted.begin(), sorted.end()); G4double dmin = sorted[0]; unsigned int isave = find( diffs.begin(), diffs.end(), dmin ) - diffs.begin(); if(sorted[0]>1.e+3) { G4cout<<"SimpleList problem? dmin="<< sorted[0] << G4endl; // Debugging printout for(unsigned int i=0; i<diffs.size(); i++) { G4cout<<"SimpleList: mcps,trks,diff="<< mcps3[i] <<' '<< trks3[i] <<' '<< diffs[i] << G4endl; } G4cout<< "Best match:"<< mcps3[isave]<<' '<< trks3[isave] <<' ' << diffs[isave] << G4endl; } // swap positions of first and best-matching elements // mcps swap unsigned int bestpos = find(mcps2.begin(), mcps2.end(), mcps3[isave]) - mcps2.begin(); if(bestpos!=0) { mcps2[bestpos] = mcps2[0]; mcps2[0] = mcps3[isave]; } // trks swap bestpos = find(trks2.begin(), trks2.end(), trks3[isave]) - trks2.begin(); if(bestpos!=0) { trks2[bestpos] = trks2[0]; trks2[0] = trks3[isave]; } }

void LCDG4McPartManager::GetDaughterCandidates (  const G4int  dau1, const G4int  dau2, const G4int  pdgidref, std::vector< LCDMcPart * > &  cands )

Definition at line 756 of file LCDG4McPartManager.cc. References LCDMcPart::GetParticleID(). { for(int i=dau1; i<=dau2; i++) { LCDMcPart* test = MCparticles[i-1]; if(test->GetParticleID()==pdgidref) { cands.push_back(test); } } }

int LCDG4McPartManager::GetDetType (  G4LogicalVolume *  lvol  )  const [inline]

Definition at line 48 of file LCDG4McPartManager.hh. Referenced by IsInteracted(), and LCDG4SteppingAction::UserSteppingAction(). { return (*m_pdettype)[lvol]; }

LCDMcPart* LCDG4McPartManager::GetMcPartFromTag (  G4int  mctag  )  const [inline]

Definition at line 57 of file LCDG4McPartManager.hh. References MCparticles. Referenced by GetMcPartFromTrkid(), MakeLCDMcPart(), and PostTracking(). { return MCparticles[mctag-1]; }

LCDMcPart * LCDG4McPartManager::GetMcPartFromTrkid (  G4int  trkid  )  const

Definition at line 585 of file LCDG4McPartManager.cc. References GetMcPartFromTag(), IsInOutput(), and trkidToMctag(). Referenced by MakeLCDMcPart(), and PreTracking(). { if( IsInOutput(trkid) ) return GetMcPartFromTag( trkidToMctag( trkid ) ); else return NULL; }

int LCDG4McPartManager::GetNumMCparticles (   )  const [inline]

Definition at line 49 of file LCDG4McPartManager.hh. References MCparticles. { return MCparticles.size(); }

LCDG4McPartManager* LCDG4McPartManager::GetPointer (   )  [inline, static]

Definition at line 16 of file LCDG4McPartManager.hh. References _me. Referenced by LCDG4PrimaryGeneratorAction::GeneratePrimaries(), LCDG4LcioHelper::getCalHitLcio(), and LCDG4LcioHelper::getTrkHitLcio(). { return _me; }

const G4int & LCDG4McPartManager::GetShowerFlag (  G4int  parentTrkID  )  const [inline]

Definition at line 152 of file LCDG4McPartManager.hh. References m_showerFlag. Referenced by LCDG4EventAction::EndOfEventAction(), and LCDG4SteppingAction::UserSteppingAction(). { return m_showerFlag.find(trkid)->second; }

bool LCDG4McPartManager::IsDecayed (   )  const

Definition at line 642 of file LCDG4McPartManager.cc. References CheckLastProcess(). Referenced by PostTracking(). { return CheckLastProcess( G4String("Decay") ); }

bool LCDG4McPartManager::IsInOutput (   )  const

Definition at line 605 of file LCDG4McPartManager.cc. References m_trkid. Referenced by GetMcPartFromTrkid(), PostTracking(), PreTracking(), and SetCalorimeterEntrancePoint(). { return IsInOutput(m_trkid); }

bool LCDG4McPartManager::IsInOutput (  G4int  trkid  )  const

Definition at line 595 of file LCDG4McPartManager.cc. References LCDMcPart::getTrkID(), MCparticles, and trkidToMctag(). Referenced by LCDG4SteppingAction::UserSteppingAction(). { LCDMcPart* pMcPart = 0; G4int tag = trkidToMctag(trkid); if(tag<=0) G4cout<< "IsInOutput: trkid="<< trkid<< ", tag="<< tag << G4endl; if(tag>0) pMcPart = MCparticles[tag-1]; if( pMcPart ) return ( pMcPart->getTrkID() == trkid ); else return false; }

bool LCDG4McPartManager::IsInteracted (   )  const

Definition at line 662 of file LCDG4McPartManager.cc. References GetDetType(), IsShowerProcess(), and postStep. Referenced by IsShowered(), and PostTracking(). { // always check for a showering process if( !IsShowerProcess() ) return false; G4VPhysicalVolume* physVol = postStep->GetPhysicalVolume(); G4LogicalVolume* lvol=0; if(physVol) lvol = physVol->GetLogicalVolume(); G4int dtype = GetDetType(lvol); if( dtype >= 10 && dtype <= 20 ) { // inside calorimeters or muon system // to be flagged as showered, not interacted return false; } else { // showering process outside calorimeters, keep particle as INTERACTED return true; } }

bool LCDG4McPartManager::IsLeft (   )  const

Definition at line 657 of file LCDG4McPartManager.cc. References CheckLastProcess(), and curTrack. Referenced by PostTracking(). { double kinE = curTrack->GetKineticEnergy(); return ( CheckLastProcess( G4String("Transportation") ) && kinE>0 ); }

bool LCDG4McPartManager::IsShowered (   )

Definition at line 681 of file LCDG4McPartManager.cc. References curTrack, IsInteracted(), IsShowerProcess(), m_debug, m_mctag, m_showerFlag, m_trkid, and MCparticles. Referenced by PostTracking(). { G4int parTrkId = curTrack->GetParentID(); if(m_debug&0x80) { G4cout << "IsShowered(): trkid,parTrkId="<< m_trkid <<' '<< parTrkId <<", parent ShowerFlag="<< m_showerFlag[parTrkId] << G4endl; } // catch most obvious cases first if( m_showerFlag[parTrkId]>1 ) { // definitely a shower! // keep track of particle generations inside showers m_showerFlag[m_trkid] = m_showerFlag[parTrkId] + 1; if(m_debug&0x80) G4cout << "Showered, showerFlag set to "<< m_showerFlag[m_trkid] << G4endl; return true; } /* Now for most tricky cases: Shower initiators (showerFlag=1) are those who enter the calorimeters, as opposed to shower secondaries (showerFlag>1), which are produced inside the calorimeters. Some initiators may have non-shower children, which are produced before the calorimeters and may become shower initiators themselves. So the trick here is to separate and identify all these possibilities. */ m_showerFlag[m_trkid]=0; bool originator = IsShowerProcess() && !IsInteracted(); // G4cout <<"More checks: IsShowerProc,IsInteracted = " << IsShowerProcess() // <<' '<< IsInteracted() << G4endl; if( m_showerFlag[parTrkId]==0 ) { // parent was not part of a shower, so return true only if // child is an originator if( originator ) { // child has showered in calorimeters m_showerFlag[m_trkid] = 1; return true; } } else if(m_showerFlag[parTrkId]==1) { // parent started a shower, but child may have been produced // before the calorimeters. More checks are needed if( MCparticles[m_mctag-1]->getTrkID() != m_trkid ) { // no LCDMcPart, so child is just a shower secondary m_showerFlag[m_trkid] = m_showerFlag[parTrkId] + 1; if(m_debug&0x80) G4cout <<"Just a shower secondary, not saved"<< G4endl; return true; } else { // it has its own LCDMcPart, but it may still be an originator if( originator ) { // child has showered in calorimeters m_showerFlag[m_trkid] = 1; if(m_debug&0x80) G4cout <<"Both parent and child are originators."<< G4endl; return true; } } } else { G4cout << "***** LCDG4McPartManager error, please investigate!"<< G4endl; exit(1); } if(m_debug&0x80) G4cout<<" Not a shower particle!" << G4endl; return false; }

bool LCDG4McPartManager::IsShowerProcess (   )  const

Definition at line 613 of file LCDG4McPartManager.cc. References CheckLastProcess(), and dpart. Referenced by IsInteracted(), and IsShowered(). { if( CheckLastProcess( G4String("Inelastic") ) ) return true; // When pions stop, they are absorbed by a nucleus producing // several byproducts like neutrons, protons, deuterons, tritons, etc if( CheckLastProcess( G4String("AbsorptionAtRest") ) ) return true; // EM interaction before cal // gamma G4int pdgid = dpart->GetPDGEncoding(); if( pdgid==22 ) { // for now all gamma processes are considered showering return true; } // electron else if( abs(pdgid)==11 ) { // for now all e+/e- processes are considered showering return true; } else if( abs(pdgid)==111 ) { // pizero // pizero decays are considered showering if( CheckLastProcess("Decay") ) return true; } return false; }

bool LCDG4McPartManager::IsStopped (   )  const

Definition at line 646 of file LCDG4McPartManager.cc. References CheckLastProcess(), and curTrack. Referenced by PostTracking(). { double kinE = curTrack->GetKineticEnergy(); if(kinE==0) { // G4cout << "E=0: lastproc="<< lastProc->GetProcessName() << G4endl; // inelastics are included in interacted or showered if( CheckLastProcess( G4String("Inelastic") ) ) return false; return true; } return false; }

void LCDG4McPartManager::MakeLCDMcPart (  const G4Track *  aTrack  )

Definition at line 375 of file LCDG4McPartManager.cc. References dpart, LCDMcPart::Get4Momentum(), LCDMcPart::GetGenStatus(), GetMcPartFromTag(), GetMcPartFromTrkid(), LCDG4Trk2McPart::GetParent4Momentum(), LCDMcPart::GetParentTag(), LCDMcPart::GetParticleID(), LCDMcPart::GetParticleTag(), LCDMcPart::GetTermPosition(), LCDMcPart::getTrkID(), m_debug, m_mctag, m_trkid, m_trkid2McTag, MCparticles, LCDMcPart::Set4Momentum(), LCDMcPart::SetChargeFromPid(), LCDMcPart::SetG4Status(), LCDMcPart::SetGenStatus(), LCDG4Trk2McPart::SetMcPartTag(), LCDMcPart::SetParentTag(), LCDMcPart::SetParticleID(), LCDMcPart::SetParticleTag(), LCDMcPart::SetPosition(), LCDMcPart::setTrkID(), and trkidToMctag(). Referenced by LCDG4SteppingAction::UserSteppingAction(). { G4int trkid = aTrack->GetTrackID(); assert(trkid==m_trkid); // if parent does not have its own LCDMcPart, bail out G4int ptrkid = aTrack->GetParentID(); LCDMcPart* parMcPart = GetMcPartFromTrkid(ptrkid); if( !parMcPart ) { if(m_debug & 0x40) { G4cout<<"MakeLCDMcPart: denial: No LCDMcPart for parent" << G4endl; } return; } // If parent decay was preassigned by generator, the daughter's // LCDMcPart already exists, and daughter tag was already found // at parent's PostTracking() G4int parStat = parMcPart->GetGenStatus(); if( parStat==2 || parStat==3 ) { // make sure tag is already assigned here... G4int tag = trkidToMctag(trkid); if(tag>0) { // some printout LCDMcPart* mcp = MCparticles[tag-1]; if(m_debug & 0x40) { G4cout<<"Using existing LCDMcPart: tag="<< tag <<", pdgid="<< mcp->GetParticleID() <<", parTag="<< parMcPart->GetParticleTag() <<", genStat="<< mcp->GetGenStatus() << G4endl; } // Update momentum of preassigned decay products HepLorentzVector& mom = mcp->Get4Momentum(); mom.setVect( aTrack->GetStep()->GetPreStepPoint()->GetMomentum() ); mom.setE( aTrack->GetStep()->GetPreStepPoint()->GetTotalEnergy() ); return; } else { // This should never happen, as neg tags become positive in PreTracking assert(false); } } //****** Split parent G4Track into two McParts ******** // If produced by an interaction (not a decay) const G4VProcess* proc = aTrack->GetCreatorProcess(); if( (!proc->GetProcessName().contains("Decay")) ) { /* At this point, the particle being processed (DAU) has been produced * from an non-destructive interaction of its parent (PAR) in the * tracker volume. * In order to save the production point of DAU, one has to split PAR * into two McParts (PAR->PAR2+DAU+X) at that point, so that the * production point of DAU and PAR2 is the endpoint of PAR. Other * daughter particles (X) may also be produced in the interaction. */ // Check endpoint of PAR. If it is the same as starting point of // current track, that means PAR2 was already created previously, // so no need to create another one... if( parMcPart->GetTermPosition() == aTrack->GetVertexPosition() ) { if(m_debug & 0x40) { G4cout<< " PAR was already split at this very same point: " << parMcPart->GetTermPosition() <<", no need to split it again."<< G4endl; } } else { // create a new McPart for PAR2, so the starting point is saved int par2Tag = MCparticles.size()+1; LCDMcPart* par2Mcp = new LCDMcPart( par2Tag, parMcPart, aTrack->GetVertexPosition() ); // PAR McPart's children to be inherited by PAR2 McPart G4int oldParentTag = parMcPart->GetParticleTag(); for(int i=1; i<par2Tag; i++) { // par2Tag excluded, of course! LCDMcPart* mcp = GetMcPartFromTag(i); // make sure mcp is valid, as some McParts may have been erased if( mcp && mcp->GetParentTag() == oldParentTag ) { if(m_debug & 0x40) { G4cout << "Resetting parent tags: "<< oldParentTag <<" -> "<< par2Tag <<" for tag="<< i << G4endl; } mcp->SetParentTag( par2Tag ); } } // 4-momentum of PAR2 McPart after interaction // Makoto: save momentum in daughter's TrackInfo, and use it here LCDG4Trk2McPart* secInfo = (LCDG4Trk2McPart*)aTrack->GetUserInformation(); const HepLorentzVector& emom = secInfo->GetParent4Momentum(); par2Mcp->Set4Momentum( emom ); // save McPart* into list if(m_debug & 0x40) { G4cout << "Creating new parent McPart: tag="<< par2Tag <<", emom="<< emom << G4endl; } MCparticles.push_back( par2Mcp ); } // else } // Next, resume processing DAU //******************************** // create a new LCDMcPart LCDMcPart* mc = new LCDMcPart(); mc->SetParticleTag( MCparticles.size() + 1 ); // particles created by Geant4, GenStatus=4 mc->SetGenStatus(4); if(mc->getTrkID()>0) { G4cout<<"problem: TrkID already set!! tag="<< mc->GetParticleTag() <<", old=" << mc->getTrkID() <<", new="<< trkid << G4endl; } mc->setTrkID(trkid); G4int newTag = mc->GetParticleTag(); // G4cout<< "MakeLCDMcPart: Set trkid2McTag["<<trkid<<"]="<<newTag<< G4endl; m_trkid2McTag[trkid] = newTag; m_mctag = newTag; // Set G4VUserTrackInfo so that SensitiveDetectors can see it LCDG4Trk2McPart* TrkInfo = (LCDG4Trk2McPart*)(aTrack->GetUserInformation()); TrkInfo->SetMcPartTag(newTag); if(m_debug & 0x40) { G4cout << " Setting UserTrkInfo: trkid,newTag = "<< aTrack->GetTrackID() <<' '<< newTag << G4endl; } //((LCDG4Trk2McPart*)aTrack->GetUserInformation())->SetMcPartTag(newTag); G4ThreeVector pos = aTrack->GetVertexPosition(); G4ThreeVector mom = aTrack->GetVertexMomentumDirection(); G4double Ekin = aTrack->GetVertexKineticEnergy(); // sets momentum magnitude if(abs(mom.mag() - 1.0) < 1e-4) { G4double mass = aTrack->GetDefinition()->GetPDGMass(); G4double pmag = sqrt(Ekin*Ekin + 2*mass*Ekin); // mom.setMag(Ekin); mom.setMag(pmag); } else { G4cout<< "Problem: mom.mag() != 1?!?"<<G4endl; } Hep3Vector vpos(pos.x(),pos.y(),pos.z()); HepLorentzVector vmom4(mom.x(),mom.y(),mom.z(), aTrack->GetTotalEnergy()); G4ParticleDefinition* dpart = aTrack->GetDefinition(); G4int pdgid = dpart->GetPDGEncoding(); mc->SetG4Status(0); mc->SetParticleID(pdgid); mc->SetChargeFromPid(pdgid); mc->SetPosition(vpos); mc->Set4Momentum(vmom4); // debugging printout if(m_debug & 0x40) { G4cout<<"Creating new LCDMcPart: tag="<< newTag <<", pdgid="<< mc->GetParticleID() <<", parTag="<< parMcPart->GetParticleTag() <<", genStat="<< mc->GetGenStatus() << G4endl; } // add to list if( newTag==(int)MCparticles.size()+1 ) { // note: parentage is only kept at daughters, as in the current scheme, // parents can only keep track of adjacent-tagged daughters mc->SetParentTag( parMcPart->GetParticleTag() ); MCparticles.push_back(mc); } }

void LCDG4McPartManager::PostTracking (  const G4Track *  aTrack, const G4TrackingManager *  aTrackMgr )

Definition at line 176 of file LCDG4McPartManager.cc. References CheckLastProcess(), dpart, EraseAllPreassignedDescendants(), GetMcPartFromTag(), LCDG4Trk2McPart::GetMcPartTag(), LCDMcPart::GetParentTag(), LCDMcPart::GetPosition(), LCDMcPart::GetTermPosition(), IsDecayed(), IsInOutput(), IsInteracted(), IsLeft(), IsShowered(), IsStopped(), lastProc, m_debug, m_mctag, m_showerFlag, m_trkid, MCparticles, postStep, PreAssignedDecay(), LCDMcPart::SetG4Status(), LCDMcPart::SetPosition(), and LCDMcPart::SetTermPosition(). Referenced by LCDG4TrackingAction::PostUserTrackingAction(). { // get tag and trkid for this track LCDG4Trk2McPart* trkInfo = (LCDG4Trk2McPart*)aTrack->GetUserInformation(); assert( m_mctag == trkInfo->GetMcPartTag() ); // Find out what happened at very last step dpart = aTrack->GetDefinition(); postStep = aTrack->GetStep()->GetPostStepPoint(); lastProc = postStep->GetProcessDefinedStep(); // Deal with all possibilities // set shower particle tag in IsShowered() m_showerFlag[m_trkid]=0; IsShowered(); // Save startPoint,endPoint if this particle has its own LCDMcPart const G4ThreeVector& startPoint = aTrack->GetVertexPosition(); const G4ThreeVector& endPoint = postStep->GetPosition(); if( IsInOutput() ) { LCDMcPart* pMcPart = MCparticles[m_mctag-1]; pMcPart->SetPosition(startPoint); pMcPart->SetTermPosition(endPoint); // Make sure parent endpoint is the same as daughter starting point int parentTag = pMcPart->GetParentTag(); if(parentTag>0) { LCDMcPart* parentMcPart = GetMcPartFromTag( parentTag ); if( parentMcPart->GetTermPosition() != pMcPart->GetPosition() ) { if(m_debug&0x80) { G4cout<< "Resetting parent endPoint: tag="<< pMcPart->GetParentTag() << ", from " << parentMcPart->GetTermPosition() <<" to "<< pMcPart->GetPosition() << G4endl; } parentMcPart->SetTermPosition(pMcPart->GetPosition()); } } // classify track and assign Geant4 status if( IsShowered() ) { if(m_debug&0x80) G4cout<<" PostTrk: Showered: "<< G4endl; pMcPart->SetG4Status(9); } else if( IsDecayed() ) { if(m_debug&0x80) G4cout<<" PostTrk: Decay: "<< G4endl; pMcPart->SetG4Status(1); } else if( IsInteracted() ) { if(m_debug&0x80) G4cout<<" PostTrk: Interacted: "<< G4endl; pMcPart->SetG4Status(2); } else if( IsLeft() ) { if(m_debug&0x80) G4cout<<" PostTrk: Left: "<< G4endl; pMcPart->SetG4Status(3); } else if( IsStopped() ) { if(m_debug&0x80) G4cout<<" PostTrk: Stopped: "<< G4endl; pMcPart->SetG4Status(4); } else { // should never reach this point assert(false); } } // loop over secondaries G4TrackVector* secTrkV = aTrackMgr->GimmeSecondaries(); G4Track *trk = 0; if(m_debug&0x80) { G4cout<<"PostTrk: # secondaries = "<< secTrkV->size() <<G4endl; } for(unsigned int isectrk = 0; isectrk < secTrkV->size(); isectrk++) { trk = (*secTrkV)[isectrk]; // set m_trkid2McTag map so that shower secondaries' tags point to // the same tag as its parent tag LCDG4Trk2McPart* trkInfo = (LCDG4Trk2McPart*)trk->GetUserInformation(); if(trkInfo == 0) { trk->SetUserInformation(new LCDG4Trk2McPart(m_mctag)); } if(m_debug&0x80) { // trk->GetTrackID() returns zero at this point if(isectrk==0) G4cout<<"PostTrk: secondaries= "; G4cout<< trk->GetDefinition()->GetParticleName(); if( isectrk == secTrkV->size()-1) G4cout << G4endl; else G4cout <<", "; } } //printf("LCDG4MCP::Postracking--m_mctag=%d\n",m_mctag); // Check preassigned decay -- if so, tags in daughter's trkInfos // will be reset const G4DecayProducts* pKids = aTrack->GetDynamicParticle()->GetPreAssignedDecayProducts(); if( pKids ) { if(m_debug&0x100) { G4cout<< "*************** Check preassigned decays ********" << G4endl; pKids->DumpInfo(); } if( CheckLastProcess( G4String("Decay") ) ) { // try to associate secondaries to LCDMcParts PreAssignedDecay(secTrkV); } else if( CheckLastProcess( G4String("Inelastic") ) ) { // an inelastic interaction happens before the preassigned decay // erase the LCDMcParts for the children of the preassigned decay EraseAllPreassignedDescendants(m_mctag); } } // finished processing for this track m_trkid = 0; }

void LCDG4McPartManager::PreAssignedDecay (  G4TrackVector *  SecTrkV  )  [private]

Definition at line 776 of file LCDG4McPartManager.cc. References AssociateFromComplexLists(), m_mctag, and MCparticles. Referenced by PostTracking(). { G4Track *trk = 0; // build list of G4Track secondaries list<G4Track*> secTrks; for(unsigned int i=0; i<secTrkV->size(); i++) { trk = (*secTrkV)[i]; // make sure we don't pick any other secondaries than // those from preassigned decays, which have associated G4PrimaryParticles bool isPrimary=(trk->GetDynamicParticle()->GetPrimaryParticle() != NULL); if( isPrimary ) secTrks.push_back(trk); } // remember that m_trkid,m_mctag refer to the parent! list<LCDMcPart*> mcParts; mcParts.push_back( MCparticles[m_mctag-1] ); while(mcParts.size()>0) { // replace particles in the list by its immediate daughters list<LCDMcPart*> moms(mcParts); mcParts.clear(); for(list<LCDMcPart*>::iterator itmcp=moms.begin(); itmcp!=moms.end(); itmcp++) { G4int dau1,dau2; (*itmcp)->GetDaughterTags( dau1, dau2 ); if(dau1>0) { for(int i=dau1; i<=dau2; i++) { LCDMcPart* mcp = MCparticles[i-1]; mcParts.push_back( mcp ); } } } // Make associations between G4Tracks and LCDMcParts AssociateFromComplexLists( secTrks, mcParts ); } //printf("LCDG4MCP::Postracking--m_mctag=%d\n",m_mctag); }

void LCDG4McPartManager::PreTracking (  const G4Track *  aTrack  )

Definition at line 70 of file LCDG4McPartManager.cc. References AssignTagToTrack(), CheckPrimaryParticle(), curTrack, LCDMcPart::GetGenStatus(), GetMcPartFromTrkid(), LCDG4Trk2McPart::GetMcPartTag(), LCDMcPart::GetParticleID(), LCDMcPart::GetParticleTag(), GeV, IsInOutput(), m_debug, m_mctag, m_nmcprim, m_trkid, m_trkid2McTag, and MCparticles. Referenced by LCDG4TrackingAction::PreUserTrackingAction(). { // set member data for track id being processed m_trkid = aTrack->GetTrackID(); if(m_debug & 0x10) { G4ThreeVector pvec = aTrack->GetDynamicParticle()->GetMomentum(); G4cout<< "PreTrking: trkid=" << m_trkid << ", id="<< aTrack->GetDefinition()->GetParticleName() << ", parID="<< aTrack->GetParentID() <<", E="<< aTrack->GetDynamicParticle()->GetTotalEnergy()/GeV // <<", pvec=("<< pvec.x()/GeV <<", "<< pvec.y()/GeV <<", "<< pvec.z()/GeV <<')' << G4endl; } // m_nmcprim is cleared at EndOfEventAction if(m_nmcprim == 0) { // initialize m_nmcprim for each event m_nmcprim = MCparticles.size(); } // creates TrkInfo with tag=-1 if non-existing curTrack = const_cast<G4Track*>(aTrack); if(aTrack->GetUserInformation() == 0) { if(m_debug&0x10) { G4cout<<" Create invalid TrkInfo (tag=-1) for this track!"<< G4endl; } curTrack->SetUserInformation( new LCDG4Trk2McPart(-1) ); } // If no G4parent, it certainly comes from input file. // Find its LCDMcPart and save its tag LCDMcPart* mcp = NULL; G4int parTrkid = aTrack->GetParentID(); if(parTrkid==0) { mcp = CheckPrimaryParticle(aTrack); assert( mcp ); // connect LCDMcPart <--> G4Track if(mcp) AssignTagToTrack( aTrack, mcp ); // debugging printout if(m_debug & 0x10) { G4cout<<"primary: <"<< m_trkid <<','<< mcp->GetParticleTag() <<'>' <<", pdgid="<< mcp->GetParticleID() <<", trkid2mctag.size="<< m_trkid2McTag.size() << G4endl; } } // if( parentTrackID == 0 ) else { // It might come from input file (preassigned decays) or not // does immediate parent have an LCDMcPart? if( IsInOutput(parTrkid) ) { // yes it does // Is parent decay preassigned? LCDMcPart* mcpPar = GetMcPartFromTrkid(parTrkid); G4int parStat = mcpPar->GetGenStatus(); if( parStat==2 || parStat==3 ) { // yes // Check if tag assigned from PostTracking of parent // (negative tag stored in trkInfo) LCDG4Trk2McPart* trkInfo = (LCDG4Trk2McPart*)aTrack->GetUserInformation(); G4int tag = trkInfo->GetMcPartTag(); if( tag<0 ) { // set to positive tag and complete G4Track <--> LCDMcPart connection tag *= -1; // Some invalid trackinfo have been created elsewhere with tag=-1, // we don't want to get those here... assert(tag!=1); AssignTagToTrack( aTrack, tag ); // debugging printout if(m_debug & 0x10) { G4cout<<"reusing: <"<< m_trkid <<','<< tag <<'>' <<", pdgid="<< MCparticles[tag-1]->GetParticleID() <<", trkid2mctag.size="<< m_trkid2McTag.size() << G4endl; } } else { // Tracks showing up here are secondaries from a GenStatus=2 parent // (preassigned decays), but other than its preassigned decay // products. // These secondaries are usually from EM interactions before the // calorimeters, so they don't get associated to any preexisting // McPart here, and a new one might be created at its first step. } } // endif(genStatus==2) } // endif(parent has its own LCDMcPart) } // endif(G4ParentID()!=0) LCDG4Trk2McPart* trkInfo = (LCDG4Trk2McPart*)aTrack->GetUserInformation(); m_mctag = trkInfo->GetMcPartTag(); // This tag may be invalid at this point (-1) // set trkid to tag map for this track, if not done yet if(m_trkid2McTag[m_trkid] == 0) { // G4cout<<"PreTrk: Setting trkid2McTag["<<m_trkid<<"]="<<m_mctag<< G4endl; m_trkid2McTag[m_trkid] = m_mctag; } }

void LCDG4McPartManager::PrintMcPart (   )

Definition at line 550 of file LCDG4McPartManager.cc. References LCDMcPart::Get4Momentum(), LCDMcPart::GetChargeFromPid(), LCDMcPart::GetG4Status(), LCDMcPart::GetGenStatus(), LCDMcPart::GetParentTag(), LCDMcPart::GetParticleID(), LCDMcPart::GetPosition(), and MCparticles. { // vector<LCDMcPart*>::const_iterator part; // for(part = MCparticles.begin(); part != MCparticles.end(); part++) { for(unsigned int i=0; i<MCparticles.size(); i++) { LCDMcPart *mcp = MCparticles[i]; G4cout << "GenStatus:" << mcp->GetGenStatus() << G4endl; G4cout << "G4Status:" << mcp->GetG4Status() << G4endl; G4cout << "PDGID:" << mcp->GetParticleID() << G4endl; G4cout << "Charge:" << mcp->GetChargeFromPid() << G4endl; G4cout << "Position:" << mcp->GetPosition() << G4endl; G4cout << "Momentum:" << mcp->Get4Momentum() << G4endl; G4cout << "Parent Index->" << mcp->GetParentTag() << G4endl; G4cout << G4endl; } }

void LCDG4McPartManager::PrintTrkid2MctagMap (   )

Definition at line 567 of file LCDG4McPartManager.cc. References m_trkid2McTag. { std::map<G4int,G4int>::const_iterator iter = m_trkid2McTag.begin(); for (; iter != m_trkid2McTag.end(); iter++) { G4cout<<"trkid, mcid = <"<< iter->first <<','<< iter->second <<">"<<G4endl; } }

void LCDG4McPartManager::SetCalorimeterEntrancePoint (  const G4ThreeVector &  pos, const G4ThreeVector &  mom )

Definition at line 58 of file LCDG4McPartManager.cc. References IsInOutput(), m_mctag, m_trkid, MCparticles, LCDMcPart::SetCalorimeterEntranceMomentum(), and LCDMcPart::SetCalorimeterEntrancePoint(). Referenced by LCDG4SteppingAction::UserSteppingAction(). { if( IsInOutput(m_trkid) ) { LCDMcPart* mcp = MCparticles[ m_mctag-1 ]; mcp->SetCalorimeterEntrancePoint(pos); mcp->SetCalorimeterEntranceMomentum(mom); } }

void LCDG4McPartManager::SetDebugFlag (  int  dbflg  )  [inline]

Definition at line 22 of file LCDG4McPartManager.hh. References m_debug. Referenced by LCDG4EventAction::BeginOfEventAction(). { G4cout << "***** Setting debug flag to "<< dbflg << G4endl; m_debug = dbflg; }

void LCDG4McPartManager::SetDetTypeMap (  std::map< G4LogicalVolume *, G4int > *  m  )  [inline]

Definition at line 21 of file LCDG4McPartManager.hh. References m_pdettype. { m_pdettype = m; }

G4int LCDG4McPartManager::trkidToMctag (  G4int  trkid  )  const [inline]

Definition at line 53 of file LCDG4McPartManager.hh. References m_trkid2McTag. Referenced by LCDG4EventAction::EndOfEventAction(), LCDG4LcioHelper::getCalHitLcio(), GetMcPartFromTrkid(), LCDG4LcioHelper::getTrkHitLcio(), IsInOutput(), and MakeLCDMcPart(). { return m_trkid2McTag.find(trkid)->second; }

---

Member Data Documentation

LCDG4McPartManager * LCDG4McPartManager::_me = NULL [static, private]

Definition at line 38 of file LCDG4McPartManager.cc. Referenced by GetPointer(), and LCDG4McPartManager().

G4Track* LCDG4McPartManager::curTrack [private]

Definition at line 128 of file LCDG4McPartManager.hh. Referenced by IsLeft(), IsShowered(), IsStopped(), and PreTracking().

const G4ParticleDefinition* LCDG4McPartManager::dpart [private]

Definition at line 129 of file LCDG4McPartManager.hh. Referenced by IsShowerProcess(), MakeLCDMcPart(), and PostTracking().

const G4VProcess* LCDG4McPartManager::lastProc [private]

Definition at line 131 of file LCDG4McPartManager.hh. Referenced by CheckLastProcess(), and PostTracking().

int LCDG4McPartManager::m_debug [private]

Definition at line 123 of file LCDG4McPartManager.hh. Referenced by AssignTagToTrack(), CheckPrimaryParticle(), IsShowered(), MakeLCDMcPart(), PostTracking(), PreTracking(), and SetDebugFlag().

int LCDG4McPartManager::m_mctag [private]

Definition at line 127 of file LCDG4McPartManager.hh. Referenced by IsShowered(), MakeLCDMcPart(), PostTracking(), PreAssignedDecay(), PreTracking(), and SetCalorimeterEntrancePoint().

int LCDG4McPartManager::m_nmcprim [private]

Definition at line 122 of file LCDG4McPartManager.hh. Referenced by CheckPrimaryParticle(), Clear(), and PreTracking().

std::map<G4LogicalVolume*,G4int>* LCDG4McPartManager::m_pdettype [private]

Definition at line 144 of file LCDG4McPartManager.hh. Referenced by SetDetTypeMap().

std::map<G4int,G4int> LCDG4McPartManager::m_showerFlag [private]

Definition at line 147 of file LCDG4McPartManager.hh. Referenced by GetShowerFlag(), IsShowered(), and PostTracking().

G4int LCDG4McPartManager::m_trkid [private]

Definition at line 125 of file LCDG4McPartManager.hh. Referenced by IsInOutput(), IsShowered(), MakeLCDMcPart(), PostTracking(), PreTracking(), and SetCalorimeterEntrancePoint().

std::map<G4int,G4int> LCDG4McPartManager::m_trkid2dettype [private]

Definition at line 146 of file LCDG4McPartManager.hh. Referenced by Clear(), and LCDG4McPartManager().

std::map<G4int,G4int> LCDG4McPartManager::m_trkid2McTag [private]

Definition at line 145 of file LCDG4McPartManager.hh. Referenced by AssignTagToTrack(), Clear(), LCDG4McPartManager(), MakeLCDMcPart(), PreTracking(), PrintTrkid2MctagMap(), and trkidToMctag().

std::vector<LCDMcPart*> LCDG4McPartManager::MCparticles

Definition at line 84 of file LCDG4McPartManager.hh. Referenced by AssignTagToTrack(), CheckPrimaryParticle(), Clear(), LCDG4EventAction::EndOfEventAction(), EraseAllPreassignedDescendants(), LCDG4PrimaryGeneratorAction::GeneratePrimaries(), LCDG4readStdFile::GeneratePrimaryVertex(), GetMcPartFromTag(), GetNumMCparticles(), IsInOutput(), IsShowered(), MakeLCDMcPart(), LCDG4readStdFile::MakeMcParts(), PostTracking(), PreAssignedDecay(), PreTracking(), PrintMcPart(), and SetCalorimeterEntrancePoint().

const G4StepPoint* LCDG4McPartManager::postStep [private]

Definition at line 130 of file LCDG4McPartManager.hh. Referenced by IsInteracted(), and PostTracking().

std::vector<G4int> LCDG4McPartManager::trkids

Definition at line 85 of file LCDG4McPartManager.hh. Referenced by Clear().

---

The documentation for this class was generated from the following files:

• LCDG4McPartManager.hh
• LCDG4McPartManager.cc

---