LCDG4: LCDG4SteppingAction class Reference

00049                                                                    {
00050 
00051   // this should be uncommented for geantino scanning
00052 //   UserSteppingDebug(theStep);
00053 
00054   G4Track * theTrack = theStep->GetTrack();
00055 
00056   // Check whether one needs cal entrance point for this track
00057   if( theTrack->GetCurrentStepNumber() == 1) {
00058     // reset flag
00059     calEntryNeeded=false;
00060     pvol = theTrack->GetVolume();
00061     lvol = pvol->GetLogicalVolume();
00062 //     G4cout<< "Step1: trkid="<< theTrack->GetTrackID()
00063 //        <<", detType="<< m_mcmgr->GetDetType(lvol)
00064 //        <<", detName="<< pvol->GetName() << G4endl;
00065     G4int dtype = m_mcmgr->GetDetType(lvol);
00066     if( dtype < 10 || dtype > 20 ) {
00067       // if before calorimeters, we might need entrance point
00068 //       G4cout<< "*** calEntryNeeded = true"<< G4endl;
00069       calEntryNeeded=true;
00070     }
00071   }
00072 
00073   if( calEntryNeeded ) {
00074     pvol = theTrack->GetNextVolume();
00075     lvol = 0;
00076     if(pvol) lvol = pvol->GetLogicalVolume();
00077     // look for cal entrance point
00078     G4int detType = m_mcmgr->GetDetType(lvol);
00079     // All calorimeter components go from 10 to 15
00080     if( detType >= 10 && detType <= 15 ) {
00081       // save it
00082       const G4StepPoint* stepPoint = theStep->GetPostStepPoint();
00083       const G4ThreeVector pos = stepPoint->GetPosition();
00084       const Hep3Vector mom = stepPoint->GetMomentum();
00085       m_mcmgr->SetCalorimeterEntrancePoint(pos,mom);
00086       calEntryNeeded = false;
00087     }
00088   }
00089 
00090   // Look at secondaries
00091   G4TrackVector* secVec = fpSteppingManager->GetSecondary();
00092 
00093 //   G4cout<< "# of Secondaries: "<< secVec->size()
00094 //      <<": along,rest,post="<< fpSteppingManager->GetfN2ndariesAtRestDoIt()
00095 //      <<' '<< fpSteppingManager->GetfN2ndariesAlongStepDoIt()
00096 //      <<' '<< fpSteppingManager->GetfN2ndariesPostStepDoIt()
00097 //      << G4endl;
00098 
00099   G4int nSecsInStep
00100     = fpSteppingManager->GetfN2ndariesAtRestDoIt()
00101     + fpSteppingManager->GetfN2ndariesAlongStepDoIt()
00102     + fpSteppingManager->GetfN2ndariesPostStepDoIt();
00103   G4int nSecStart = secVec->size() - nSecsInStep;
00104 
00105   // loop over secondaries generated at this very last step
00106   for(int i=0; i<nSecsInStep; i++) {
00107     if( i+nSecStart<0 ) continue;
00108     G4Track* psec =(*secVec)[i+nSecStart];
00109     G4double energy = 0;
00110     if(psec) {
00111       energy = psec->GetDynamicParticle()->GetKineticEnergy();
00112     }
00113     else {
00114       G4cout <<" psec = NULL, i,nSec(Start,InStep,Rest,Along,Post,size) = "
00115              << i <<' '<< nSecStart 
00116              <<' '<< nSecsInStep 
00117              <<' '<< fpSteppingManager->GetfN2ndariesAtRestDoIt()
00118              <<' '<< fpSteppingManager->GetfN2ndariesAlongStepDoIt()
00119              <<' '<< fpSteppingManager->GetfN2ndariesPostStepDoIt()
00120              <<' '<< secVec->size()
00121              << G4endl;
00122     }
00123     if( energy>1.*MeV ) {
00124       // Get trkInfo from secondary
00125       LCDG4Trk2McPart* secInfo = (LCDG4Trk2McPart*)psec->GetUserInformation();
00126       if(secInfo==0) {
00127         LCDG4Trk2McPart* parInfo
00128           = (LCDG4Trk2McPart*)theTrack->GetUserInformation();
00129         secInfo = new LCDG4Trk2McPart( parInfo->GetMcPartTag() );
00130         psec->SetUserInformation( secInfo );
00131       }
00132       const HepLorentzVector& emom = theTrack->GetDynamicParticle()->Get4Momentum();
00133       secInfo->SetParent4Momentum( emom );
00134     }
00135   }
00136 
00137   // check particles in their first step only
00138   if( theTrack->GetCurrentStepNumber() > 1) return;
00139 
00140 //      G4cout << "Step1: "<< theTrack->GetTrackID() <<' '
00141 //         << theTrack->GetDefinition()->GetParticleName()
00142 //         << ", PreStepPt: "<< theStep->GetPreStepPoint()->GetPosition()
00143 //         << ", PostStepPt: "<< theStep->GetPostStepPoint()->GetPosition()
00144 //         << ", PreVolume: "<< theStep->GetPreStepPoint()->GetPhysicalVolume()->GetName()
00145 //         << ", parent= "<< theTrack->GetParentID()
00146 //         << G4endl;
00147 
00148   if( theTrack->GetParentID() == 0) return;
00149 
00150   G4int parTrkID = theTrack->GetParentID();
00151 
00152   // If decay product, suggest saving it into output file
00153   // Final decision will be made by McPartManager
00154   if( theTrack->GetCreatorProcess()->GetProcessName().contains("Decay") ) {
00155 
00156 //     G4ParticleDefinition* partDef = theTrack->GetDefinition();
00157 //     G4cout << "Decay: suggesting new LCDMcPart: pvol,detType,trkid,name = "
00158 //         <<' '<< theTrack->GetVolume()->GetName()
00159 // //      <<' '<< m_mcmgr->GetDetType(lvol)
00160 //         <<' '<< theTrack->GetTrackID() <<' '<< partDef->GetParticleName()
00161 //         << G4endl;
00162 
00163     m_mcmgr->MakeLCDMcPart(theTrack);
00164     return;
00165   }
00166 
00167   // Require minimum Ekin for interaction prods, but not for decay prods
00168   if( theTrack->GetKineticEnergy() < 1.0*MeV ) return;
00169 
00170   // Require starting point to be before the calorimeters
00171   pvol = theTrack->GetVolume();
00172   lvol = 0;
00173   if(pvol) lvol = pvol->GetLogicalVolume();
00174   if( m_mcmgr->GetDetType(lvol) > 9 )  return;
00175 
00176 //   G4ParticleDefinition* partDef = theTrack->GetDefinition();
00177 //   G4cout << "Create new LCDMcPart? pvol,detType,trkid,name = "
00178 //       <<' '<< pvol->GetName() <<' '<< m_mcmgr->GetDetType(lvol)
00179 //       <<' '<< theTrack->GetTrackID()
00180 //       <<' '<< partDef->GetParticleName() << G4endl;
00181 
00182   // In order to be saved, its parent must also be saved
00183   // Check parent's showerFlag
00184   G4int showerFlag = m_mcmgr->GetShowerFlag(parTrkID);
00185   if( showerFlag==0 ) {
00186     // parent did not shower
00187 //     G4cout <<"Yes, parent has showerFlag=0" << G4endl;
00188     m_mcmgr->MakeLCDMcPart(theTrack);
00189   }
00190 
00191   else if( showerFlag==1 ) {
00192     // parent originated shower, but maybe child was produced
00193     // before parent entered the calorimeters
00194 //     G4cout <<"Maybe, parent has showerFlag==1" << G4endl;
00195     m_mcmgr->MakeLCDMcPart(theTrack);
00196   }
00197 
00198   else {
00199     // child from several generations down the showering, probably
00200     // a backscatter from the calorimeters.  Don't save into output
00201 //     G4cout <<"No, parent has showerFlag = "<< showerFlag << G4endl;
00202   }
00203 
00204   // at this point we know if the track will be saved or not
00205   // if not saved, do not bother with its calorimeter entrance point
00206   if( !m_mcmgr->IsInOutput(theTrack->GetTrackID()) ) {
00207 //     G4cout<< "Track not saved into output, resetting calEntryNeeded=false"
00208 //        << G4endl;
00209     calEntryNeeded=false;
00210   }
00211 }

Public Member Functions
	LCDG4SteppingAction ()
	LCDG4SteppingAction (LCDG4McPartManager *mcmgr)
virtual	~LCDG4SteppingAction ()
virtual void	UserSteppingAction (const G4Step *)
void	UserSteppingDebug (const G4Step *)
G4double	GetTotalRadlen ()
G4double	GetTotalNuclearInterLength ()
void	ResetTotals ()
Private Attributes
LCDG4McPartManager *	m_mcmgr
double	totStepLen
double	totRadLen
double	totIntLen
G4VPhysicalVolume *	pvol
G4LogicalVolume *	lvol
bool	calEntryNeeded
LCDG4SteppingAction Class Reference

Public Member Functions

Private Attributes

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation
