FCALSteppingAction Class Reference

#include <FCALSteppingAction.hh>

Inheritance diagram for FCALSteppingAction:

Collaboration diagram for FCALSteppingAction:

Public Member Functions
	FCALSteppingAction ()
virtual	~FCALSteppingAction ()
virtual void	UserSteppingAction (const G4Step *)
void	initialize (G4int)
G4double	GetOutOfWorldTracks (G4int, G4int)
G4double	GetSecondaries (G4int, G4int)
G4double	GetEdepFCAL (G4String)
Public Member Functions inherited from G4UserSteppingAction
	G4UserSteppingAction ()
virtual	~G4UserSteppingAction ()
void	SetSteppingManagerPointer (G4SteppingManager *pValue)

Private Attributes
G4int	EventNo
G4int	IDnow
G4int	IDold
G4int	IDout
G4int	NTracks
G4int	NSecondaries
G4double	OutOfWorldTracksData [6000][11]
G4double	Secondaries [6000][11]
G4double	EdepFCALEm
G4double	EdepFCALHad
G4ThreeVector	PrimaryVertex
G4ThreeVector	PrimaryDirection
G4ThreeVector	SecondaryVertex
G4ThreeVector	SecondaryDirection
G4ThreeVector	Distance
G4ThreeVector	VectorProduct
G4ThreeVector	VectorProductNorm
G4double	VectorProductMagnitude
G4double	DistOfClosestApproach

Additional Inherited Members
Protected Attributes inherited from G4UserSteppingAction
G4SteppingManager *	fpSteppingManager

Detailed Description

Definition at line 42 of file FCALSteppingAction.hh.

Constructor & Destructor Documentation

FCALSteppingAction()

FCALSteppingAction::FCALSteppingAction

(

)

Definition at line 57 of file FCALSteppingAction.cc.

                                      :IDold(-1),IDout(-1)
{;}

~FCALSteppingAction()

FCALSteppingAction::~FCALSteppingAction ( )

virtual

Definition at line 62 of file FCALSteppingAction.cc.

{;}

Member Function Documentation

GetEdepFCAL()

G4double FCALSteppingAction::GetEdepFCAL

(

G4String

FCAL

)

Definition at line 239 of file FCALSteppingAction.cc.

                                                      {
  if(strcmp(FCAL,"FCALEm") == 0) {
    return EdepFCALEm;
  } else {
    if(strcmp(FCAL,"FCALHad") == 0) {
      return EdepFCALHad;}
  }
  return 0.0; 
} 

Here is the caller graph for this function:

GetOutOfWorldTracks()

G4double FCALSteppingAction::GetOutOfWorldTracks	(	G4int	i,
		G4int	j
	)

Definition at line 231 of file FCALSteppingAction.cc.

                                                                {
  return OutOfWorldTracksData[i][j];
}

Here is the caller graph for this function:

GetSecondaries()

G4double FCALSteppingAction::GetSecondaries	(	G4int	i,
		G4int	j
	)

Definition at line 235 of file FCALSteppingAction.cc.

                                                           {
  return Secondaries[i][j];
}

Here is the caller graph for this function:

initialize()

void FCALSteppingAction::initialize

(

G4int

Nev

)

Definition at line 215 of file FCALSteppingAction.cc.

                                             {
  EventNo = Nev;
  NTracks = 0;
  NSecondaries = 0;
  EdepFCALEm = EdepFCALHad = 0.;

  for(G4int i=0; i<6000; i++)
    {
      for(G4int j=0; j<11; j++) 
    { 
      OutOfWorldTracksData[i][j] = 0.;
      Secondaries[i][j] = 0.; 
    }
    };
}

Here is the caller graph for this function:

UserSteppingAction()

void FCALSteppingAction::UserSteppingAction ( const G4Step *  astep )

virtual

Reimplemented from G4UserSteppingAction.

Definition at line 67 of file FCALSteppingAction.cc.

{ 
  // Get Edep
   G4double Edep = astep->GetTotalEnergyDeposit();

   // Get Track
     G4Track* aTrack = astep->GetTrack();

   // Get Touchable History
  G4TouchableHistory* theTouchable =  (G4TouchableHistory*)(aTrack->GetTouchable());

  // Energy deposit in FCAL1 and FCAL2
  if(Edep != 0.) 
    {
      G4VPhysicalVolume* physVol = theTouchable->GetVolume();
      
      if(strcmp(physVol->GetName(),"FCALEmModulePhysical")== 0 ||
     strcmp(physVol->GetName(),"F1LArGapPhysical") == 0) 
    {
      EdepFCALEm = EdepFCALEm + Edep;
    };
      
      if( (strcmp(physVol->GetName(), "FCALHadModulePhysical") == 0) ||
      (strcmp(physVol->GetName(), "CuPlateAPhysical") == 0) ||
      (strcmp(physVol->GetName(), "CuPlateBPhysical") == 0) ||
      (strcmp(physVol->GetName(), "WAbsorberPhysical") == 0) ||
      (strcmp(physVol->GetName(), "F2RodPhysical") == 0) ||
      (strcmp(physVol->GetName(), "F2LArGapPhysical") == 0) ) 
    {
      EdepFCALHad = EdepFCALHad + Edep;
    };
    };

  // Get Tracks properties
  G4int TrackID = aTrack->GetTrackID();
  G4int ParentID = aTrack->GetParentID();
  // Get Associated particle
  const G4DynamicParticle * aDynamicParticle = aTrack->GetDynamicParticle();
  G4ParticleDefinition * aParticle = aTrack->GetDefinition();
  G4String ParticleName = aParticle->GetParticleName();
  
  IDnow = EventNo + 10000*TrackID+ 100000000*ParentID;
  
  if(IDnow != IDold)
    {
      IDold = IDnow;
      
      // Get the primary particle
      if(TrackID==1 && ParentID==0 && (aTrack->GetCurrentStepNumber()) == 1)
    {
      PrimaryVertex    = aTrack->GetVertexPosition(); 
      PrimaryDirection = aTrack->GetVertexMomentumDirection();
      
      NSecondaries = 1;   
      Secondaries[NSecondaries][1] = aParticle->GetPDGEncoding();
      Secondaries[NSecondaries][2] = PrimaryVertex.x();
      Secondaries[NSecondaries][3] = PrimaryVertex.y();
      Secondaries[NSecondaries][4] = PrimaryVertex.z();
      Secondaries[NSecondaries][5] = (aDynamicParticle->GetMomentum()).x(); 
      Secondaries[NSecondaries][6] = (aDynamicParticle->GetMomentum()).y();
      Secondaries[NSecondaries][7] = (aDynamicParticle->GetMomentum()).z();
      Secondaries[NSecondaries][8] = aDynamicParticle->GetTotalMomentum();
      Secondaries[NSecondaries][9] = aDynamicParticle->GetTotalEnergy();
      Secondaries[NSecondaries][10] = aDynamicParticle->GetKineticEnergy();
      
      G4cout << " ****  Primary : " << EventNo << G4endl;
      G4cout << " Vertex : " << PrimaryVertex << G4endl;
    }
      
      
      // Get secondaries in air close to the primary tracks (DCA < 2.mm) 
      G4double DCACut = 2.*mm;
      G4String Material = aTrack->GetMaterial()->GetName();
      G4ThreeVector TrackPos = aTrack->GetVertexPosition();
      
      if(TrackID != 1 && ParentID == 1 && (strcmp(Material,"Air")==0) && (TrackPos.z() > 135.*cm)) 
    {
      SecondaryVertex = aTrack->GetVertexPosition();
      SecondaryDirection = aTrack->GetVertexMomentumDirection();
      
      // calculate DCA of secondries to primary particle
      Distance = PrimaryVertex - SecondaryVertex ;
      VectorProduct = PrimaryDirection.cross(SecondaryDirection);
      if(VectorProduct == G4ThreeVector() &&  
         PrimaryDirection != G4ThreeVector() && SecondaryDirection != G4ThreeVector()) 
        {
          G4ThreeVector Temp = Distance.cross(PrimaryDirection);
          VectorProduct = Temp.cross(PrimaryDirection);
        };    
          
VectorProductMagnitude = VectorProduct.mag();
      if(VectorProductMagnitude == 0.) 
        {
          VectorProductNorm = G4ThreeVector();
        } else {
          VectorProductNorm = (1./VectorProduct.mag()) * VectorProduct ;
        };    
      DistOfClosestApproach = Distance * VectorProductNorm ;
      
      if(std::abs(DistOfClosestApproach) < DCACut) 
        {
          NSecondaries++;         
          Secondaries[0][0] = NSecondaries;
          Secondaries[NSecondaries][1] = aParticle->GetPDGEncoding();
          Secondaries[NSecondaries][2] = (aTrack->GetVertexPosition()).x();
          Secondaries[NSecondaries][3] = (aTrack->GetVertexPosition()).y();
          Secondaries[NSecondaries][4] = (aTrack->GetVertexPosition()).z();
          Secondaries[NSecondaries][5] =(aDynamicParticle->GetMomentum()).x(); 
          Secondaries[NSecondaries][6] = (aDynamicParticle->GetMomentum()).y();
          Secondaries[NSecondaries][7] = (aDynamicParticle->GetMomentum()).z();
          Secondaries[NSecondaries][8] = aDynamicParticle->GetTotalMomentum();
          Secondaries[NSecondaries][9] = aDynamicParticle->GetTotalEnergy();
          Secondaries[NSecondaries][10] =aDynamicParticle->GetKineticEnergy();
        };  
    };
    };


  // Get the World leaving particle
  if(aTrack->GetNextVolume() == 0) {
    if(IDnow != IDout) {
      IDout = IDnow;

      NTracks++;

      OutOfWorldTracksData[0][0] = NTracks;

      OutOfWorldTracksData[NTracks][1] = aParticle->GetPDGEncoding();

      OutOfWorldTracksData[NTracks][2] = (aTrack->GetVertexPosition()).x();
      OutOfWorldTracksData[NTracks][3] = (aTrack->GetVertexPosition()).y();
      OutOfWorldTracksData[NTracks][4] = (aTrack->GetVertexPosition()).z();

      OutOfWorldTracksData[NTracks][5] = (aDynamicParticle->GetMomentum()).x();
      OutOfWorldTracksData[NTracks][6] = (aDynamicParticle->GetMomentum()).y();
      OutOfWorldTracksData[NTracks][7] = (aDynamicParticle->GetMomentum()).z();
      
      OutOfWorldTracksData[NTracks][8] = aDynamicParticle->GetTotalMomentum();

      OutOfWorldTracksData[NTracks][9] = aDynamicParticle->GetTotalEnergy();

      OutOfWorldTracksData[NTracks][10] = aDynamicParticle->GetKineticEnergy();      
    };
  };
  
 
}

Here is the call graph for this function:

Member Data Documentation

Distance

G4ThreeVector FCALSteppingAction::Distance

private

Definition at line 59 of file FCALSteppingAction.hh.

DistOfClosestApproach

G4double FCALSteppingAction::DistOfClosestApproach

private

Definition at line 60 of file FCALSteppingAction.hh.

EdepFCALEm

G4double FCALSteppingAction::EdepFCALEm

private

Definition at line 55 of file FCALSteppingAction.hh.

EdepFCALHad

G4double FCALSteppingAction::EdepFCALHad

private

Definition at line 55 of file FCALSteppingAction.hh.

EventNo

G4int FCALSteppingAction::EventNo

private

Definition at line 51 of file FCALSteppingAction.hh.

IDnow

G4int FCALSteppingAction::IDnow

private

Definition at line 52 of file FCALSteppingAction.hh.

IDold

G4int FCALSteppingAction::IDold

private

Definition at line 52 of file FCALSteppingAction.hh.

IDout

G4int FCALSteppingAction::IDout

private

Definition at line 52 of file FCALSteppingAction.hh.

NSecondaries

G4int FCALSteppingAction::NSecondaries

private

Definition at line 53 of file FCALSteppingAction.hh.

NTracks

G4int FCALSteppingAction::NTracks

private

Definition at line 53 of file FCALSteppingAction.hh.

OutOfWorldTracksData

G4double FCALSteppingAction::OutOfWorldTracksData[6000][11]

private

Definition at line 54 of file FCALSteppingAction.hh.

PrimaryDirection

G4ThreeVector FCALSteppingAction::PrimaryDirection

private

Definition at line 57 of file FCALSteppingAction.hh.

PrimaryVertex

G4ThreeVector FCALSteppingAction::PrimaryVertex

private

Definition at line 57 of file FCALSteppingAction.hh.

Secondaries

G4double FCALSteppingAction::Secondaries[6000][11]

private

Definition at line 54 of file FCALSteppingAction.hh.

SecondaryDirection

G4ThreeVector FCALSteppingAction::SecondaryDirection

private

Definition at line 58 of file FCALSteppingAction.hh.

SecondaryVertex

G4ThreeVector FCALSteppingAction::SecondaryVertex

private

Definition at line 58 of file FCALSteppingAction.hh.

VectorProduct

G4ThreeVector FCALSteppingAction::VectorProduct

private

Definition at line 59 of file FCALSteppingAction.hh.

VectorProductMagnitude

G4double FCALSteppingAction::VectorProductMagnitude

private

Definition at line 60 of file FCALSteppingAction.hh.

VectorProductNorm

G4ThreeVector FCALSteppingAction::VectorProductNorm

private

Definition at line 59 of file FCALSteppingAction.hh.

---

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

• FCALSteppingAction.hh
• FCALSteppingAction.cc