---

SteppingAction Class Reference

#include <SimG4Core/Application/interface/SteppingAction.h>

List of all members.

Public Member Functions
	SteppingAction (EventAction *ea, const edm::ParameterSet &ps)
void	UserSteppingAction (const G4Step *aStep)
	~SteppingAction ()
Public Attributes
SimActivityRegistry::G4StepSignal	m_g4StepSignal
Private Member Functions
void	catchLowEnergyInVacuumHere (const G4Step *aStep)
void	catchLowEnergyInVacuumNext (const G4Step *aStep)
bool	initPointer ()
bool	isThisVolume (const G4VTouchable *touch, G4VPhysicalVolume *pv)
Private Attributes
G4VPhysicalVolume *	calo
EventAction *	eventAction_
bool	initialized
bool	killBeamPipe
double	theCriticalDensity
double	theCriticalEnergyForVacuum
G4VPhysicalVolume *	tracker
int	verbose

---

Detailed Description

Definition at line 13 of file SteppingAction.h.

---

Constructor & Destructor Documentation

SteppingAction::SteppingAction	(	EventAction *	ea,
		const edm::ParameterSet &	ps
	)

Definition at line 13 of file SteppingAction.cc.

  : eventAction_(e),
    killBeamPipe(p.getParameter<bool>("KillBeamPipe")),
    theCriticalEnergyForVacuum(p.getParameter<double>("CriticalEnergyForVacuum")*MeV),
    theCriticalDensity(p.getParameter<double>("CriticalDensity")*g/cm3),
    verbose(p.getUntrackedParameter<int>("Verbosity",0)), 
    initialized(false), tracker(0), calo(0) {}

SteppingAction::~SteppingAction

(

)

Definition at line 21 of file SteppingAction.cc.

{}

---

Member Function Documentation

void SteppingAction::catchLowEnergyInVacuumHere

(

const G4Step *

aStep

)

[private]

Definition at line 52 of file SteppingAction.cc.

References g.

Referenced by UserSteppingAction().

                                                                    {
  G4Track * theTrack = aStep->GetTrack();
  double theKenergy = theTrack->GetKineticEnergy();
  if (theTrack->GetVolume()!=0) {
    double density = theTrack->GetVolume()->GetLogicalVolume()->GetMaterial()->GetDensity();
    if (theKenergy <= theCriticalEnergyForVacuum && theKenergy > 0.0 &&
        density <= theCriticalDensity && theTrack->GetDefinition()->GetPDGCharge() != 0 &&
        theTrack->GetTrackStatus() != fStopAndKill) {
      if (verbose>1)
        edm::LogInfo("SimG4CoreApplication") 
          <<   " SteppingAction: LoopCatchSteppingAction:catchLowEnergyInVacuumHere: "
          << " Track from " << theTrack->GetDefinition()->GetParticleName()
          << " of kinetic energy " << theKenergy/MeV << " MeV "
          << " killed in " << theTrack->GetVolume()->GetLogicalVolume()->GetName()
          << " of density " << density/(g/cm3) << " g/cm3" ;
      theTrack->SetTrackStatus(fStopAndKill);
    }
  }
}

void SteppingAction::catchLowEnergyInVacuumNext

(

const G4Step *

aStep

)

[private]

Definition at line 72 of file SteppingAction.cc.

References g.

Referenced by UserSteppingAction().

                                                                    {
  G4Track * theTrack = aStep->GetTrack();
  double theKenergy = theTrack->GetKineticEnergy();
  if (theTrack->GetNextVolume()) {
    double density = theTrack->GetNextVolume()->GetLogicalVolume()->GetMaterial()->GetDensity();
    if (theKenergy <=  theCriticalEnergyForVacuum && theKenergy > 0.0 &&
        density <= theCriticalDensity && theTrack->GetDefinition()->GetPDGCharge() != 0 &&
        theTrack->GetTrackStatus() != fStopAndKill) {
      if (verbose>1)
        edm::LogInfo("SimG4CoreApplication") 
          << " SteppingAction: LoopCatchSteppingAction::catchLowEnergyInVacuumNext: "
          << " Track from " << theTrack->GetDefinition()->GetParticleName()
          << " of kinetic energy " << theKenergy/MeV << " MeV "
          << " stopped in " << theTrack->GetVolume()->GetLogicalVolume()->GetName()
          << " before going into "<< theTrack->GetNextVolume()->GetLogicalVolume()->GetName()
          << " of density " << density/(g/cm3) << " g/cm3" ;
      theTrack->SetTrackStatus(fStopButAlive);
    }
  }
}

bool SteppingAction::initPointer

(

)

[private]

Definition at line 93 of file SteppingAction.cc.

References calo, LogDebug, and tracker.

Referenced by UserSteppingAction().

                                 {

  const G4PhysicalVolumeStore * pvs = G4PhysicalVolumeStore::GetInstance();
  if (pvs) {
    std::vector<G4VPhysicalVolume*>::const_iterator pvcite;
    for (pvcite = pvs->begin(); pvcite != pvs->end(); pvcite++) {
      if ((*pvcite)->GetName() == "Tracker") tracker = (*pvcite);
      if ((*pvcite)->GetName() == "CALO")    calo    = (*pvcite);
      if (tracker && calo) break;
    }
    edm::LogInfo("SimG4CoreApplication") << "Pointer for Tracker " << tracker
                                         << " and for Calo " << calo;
    if (tracker) LogDebug("SimG4CoreApplication") << "Tracker vol name "
                                                  << tracker->GetName();
    if (calo)    LogDebug("SimG4CoreApplication") << "Calorimeter vol name "
                                                  << calo->GetName();
    return true;
  }
  return false;
}

bool SteppingAction::isThisVolume	(	const G4VTouchable *	touch,
		G4VPhysicalVolume *	pv
	)			[private]

Definition at line 114 of file SteppingAction.cc.

References level.

Referenced by UserSteppingAction().

                                                         {

  int level = ((touch->GetHistoryDepth())+1);
  if (level > 0 && level >= 3) {
    int ii = level - 3;
    return (touch->GetVolume(ii) == pv);
  }
  return false;
}

void SteppingAction::UserSteppingAction

(

const G4Step *

aStep

)

Definition at line 23 of file SteppingAction.cc.

References EventAction::addTkCaloStateInfo(), calo, catchLowEnergyInVacuumHere(), catchLowEnergyInVacuumNext(), eventAction_, initialized, initPointer(), isThisVolume(), killBeamPipe, m_g4StepSignal, p, tracker, x, y, and z.

                                                            {
  if (!initialized) initialized = initPointer();
  m_g4StepSignal(aStep);
  if (killBeamPipe){
    catchLowEnergyInVacuumHere(aStep);
    catchLowEnergyInVacuumNext(aStep);
  }
  if (aStep->GetPostStepPoint()->GetPhysicalVolume() != 0) {
    bool ok = (isThisVolume(aStep->GetPreStepPoint()->GetTouchable(),tracker)&&
               isThisVolume(aStep->GetPostStepPoint()->GetTouchable(),calo));
    if (ok) {

      math::XYZVectorD pos((aStep->GetPreStepPoint()->GetPosition()).x(),
                           (aStep->GetPreStepPoint()->GetPosition()).y(),
                           (aStep->GetPreStepPoint()->GetPosition()).z());
      
      math::XYZTLorentzVectorD mom((aStep->GetPreStepPoint()->GetMomentum()).x(),
                                   (aStep->GetPreStepPoint()->GetMomentum()).y(),
                                   (aStep->GetPreStepPoint()->GetMomentum()).z(),
                                   aStep->GetPreStepPoint()->GetTotalEnergy());
      
      uint32_t id = aStep->GetTrack()->GetTrackID();
      
      std::pair<math::XYZVectorD,math::XYZTLorentzVectorD> p(pos,mom);
      eventAction_->addTkCaloStateInfo(id,p);
    }
  }
}

---

Member Data Documentation

G4VPhysicalVolume * SteppingAction::calo [private]

Definition at line 33 of file SteppingAction.h.

Referenced by initPointer(), and UserSteppingAction().

EventAction* SteppingAction::eventAction_ [private]

Definition at line 27 of file SteppingAction.h.

Referenced by UserSteppingAction().

bool SteppingAction::initialized [private]

Definition at line 32 of file SteppingAction.h.

Referenced by UserSteppingAction().

bool SteppingAction::killBeamPipe [private]

Definition at line 28 of file SteppingAction.h.

Referenced by UserSteppingAction().

SimActivityRegistry::G4StepSignal SteppingAction::m_g4StepSignal

Definition at line 20 of file SteppingAction.h.

Referenced by RunManager::initializeUserActions(), and UserSteppingAction().

double SteppingAction::theCriticalDensity [private]

Definition at line 30 of file SteppingAction.h.

double SteppingAction::theCriticalEnergyForVacuum [private]

Definition at line 29 of file SteppingAction.h.

G4VPhysicalVolume* SteppingAction::tracker [private]

Definition at line 33 of file SteppingAction.h.

Referenced by initPointer(), and UserSteppingAction().

int SteppingAction::verbose [private]

Definition at line 31 of file SteppingAction.h.

---

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

• SimG4Core/Application/interface/SteppingAction.h
• SimG4Core/Application/src/SteppingAction.cc

---