#include <GflashEMShowerModel.h>

Inheritance diagram for GflashEMShowerModel:

Public Member Functions
void	DoIt (const G4FastTrack &, G4FastStep &) override

	GflashEMShowerModel (const G4String &name, G4Envelope *env, const edm::ParameterSet &parSet)

G4bool	IsApplicable (const G4ParticleDefinition &) override

G4bool	ModelTrigger (const G4FastTrack &) override

	~GflashEMShowerModel () override

Private Member Functions
G4bool	excludeDetectorRegion (const G4FastTrack &fastTrack)

void	makeHits (const G4FastTrack &fastTrack)

void	updateGflashStep (const G4ThreeVector &position, G4double time)

Private Attributes
G4Navigator *	theGflashNavigator

G4Step *	theGflashStep

G4TouchableHandle	theGflashTouchableHandle

edm::ParameterSet	theParSet

GflashEMShowerProfile *	theProfile

const G4Region *	theRegion

Detailed Description

Definition at line 28 of file GflashEMShowerModel.h.

Constructor & Destructor Documentation

GflashEMShowerModel::GflashEMShowerModel	(	const G4String &	name,
		G4Envelope *	env,
		const edm::ParameterSet &	parSet
	)

Definition at line 24 of file GflashEMShowerModel.cc.

References theGflashNavigator, theGflashStep, theGflashTouchableHandle, theProfile, and theRegion.

     : G4VFastSimulationModel(modelName, envelope), theParSet(parSet) {
   theProfile = new GflashEMShowerProfile(parSet);
   theRegion = const_cast<const G4Region *>(envelope);
 
   theGflashStep = new G4Step();
   theGflashTouchableHandle = new G4TouchableHistory();
   theGflashNavigator = new G4Navigator();
 }

GflashEMShowerModel::~GflashEMShowerModel ( )

override

Definition at line 38 of file GflashEMShowerModel.cc.

References theGflashStep, and theProfile.

                                           {
   delete theProfile;
   delete theGflashStep;
 }

Member Function Documentation

void GflashEMShowerModel::DoIt	(	const G4FastTrack &	fastTrack,
		G4FastStep &	fastStep
	)

override

Definition at line 72 of file GflashEMShowerModel.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, randomXiThetaGunProducer_cfi::energy, Gflash::findShowerType(), GeV, GflashEMShowerProfile::initialize(), makeHits(), GflashEMShowerProfile::parameterization(), position, and theProfile.

                                                                                  {
   // Kill the parameterised particle:
   fastStep.KillPrimaryTrack();
   fastStep.ProposePrimaryTrackPathLength(0.0);
 
   // input variables for GflashEMShowerProfile with showerType = 1,5 (shower
   // starts inside crystals)
   G4double energy = fastTrack.GetPrimaryTrack()->GetKineticEnergy() / GeV;
   G4double globalTime = fastTrack.GetPrimaryTrack()->GetStep()->GetPostStepPoint()->GetGlobalTime();
   G4double charge = fastTrack.GetPrimaryTrack()->GetStep()->GetPreStepPoint()->GetCharge();
   G4ThreeVector position = fastTrack.GetPrimaryTrack()->GetPosition() / cm;
   G4ThreeVector momentum = fastTrack.GetPrimaryTrack()->GetMomentum() / GeV;
   G4int showerType = Gflash::findShowerType(position);
 
   // Do actual parameterization. The result of parameterization is gflashHitList
   theProfile->initialize(showerType, energy, globalTime, charge, position, momentum);
   theProfile->parameterization();
 
   // make hits
   makeHits(fastTrack);
 }

G4bool GflashEMShowerModel::excludeDetectorRegion ( const G4FastTrack & fastTrack )

private

Definition at line 144 of file GflashEMShowerModel.cc.

References PVValHelper::eta.

Referenced by ModelTrigger().

                                                                               {
   G4bool isExcluded = false;
 
   // exclude regions where geometry are complicated
   //+- one supermodule around the EB/EE boundary: 1.479 +- 0.0174*5
   G4double eta = fastTrack.GetPrimaryTrack()->GetPosition().pseudoRapidity();
   if (std::fabs(eta) > 1.392 && std::fabs(eta) < 1.566) {
     return true;
   }
 
   return isExcluded;
 }

G4bool GflashEMShowerModel::IsApplicable ( const G4ParticleDefinition & particleType )

override

Definition at line 43 of file GflashEMShowerModel.cc.

                                                                                  {
   return (&particleType == G4Electron::ElectronDefinition() || &particleType == G4Positron::PositronDefinition());
 }

void GflashEMShowerModel::makeHits ( const G4FastTrack & fastTrack )

private

Definition at line 94 of file GflashEMShowerModel.cc.

References GflashEMShowerProfile::getGflashHitList(), theGflashNavigator, theGflashStep, theGflashTouchableHandle, theProfile, theRegion, and updateGflashStep().

Referenced by DoIt().

                                                                {
   std::vector<GflashHit> &gflashHitList = theProfile->getGflashHitList();
 
   theGflashStep->SetTrack(const_cast<G4Track *>(fastTrack.GetPrimaryTrack()));
 
   theGflashStep->GetPostStepPoint()->SetProcessDefinedStep(
       const_cast<G4VProcess *>(fastTrack.GetPrimaryTrack()->GetStep()->GetPostStepPoint()->GetProcessDefinedStep()));
   theGflashNavigator->SetWorldVolume(
       G4TransportationManager::GetTransportationManager()->GetNavigatorForTracking()->GetWorldVolume());
 
   std::vector<GflashHit>::const_iterator spotIter = gflashHitList.begin();
   std::vector<GflashHit>::const_iterator spotIterEnd = gflashHitList.end();
 
   for (; spotIter != spotIterEnd; spotIter++) {
     // put touchable for each hit so that touchable history keeps track of each
     // step.
     theGflashNavigator->LocateGlobalPointAndUpdateTouchableHandle(
         spotIter->getPosition(), G4ThreeVector(0, 0, 0), theGflashTouchableHandle, false);
     updateGflashStep(spotIter->getPosition(), spotIter->getTime());
 
     // Send G4Step information to Hit/Digi if the volume is sensitive
     // Copied from G4SteppingManager.cc
 
     G4VPhysicalVolume *aCurrentVolume = theGflashStep->GetPreStepPoint()->GetPhysicalVolume();
     if (aCurrentVolume == nullptr)
       continue;
 
     G4LogicalVolume *lv = aCurrentVolume->GetLogicalVolume();
     if (lv->GetRegion() != theRegion)
       continue;
 
     theGflashStep->GetPreStepPoint()->SetSensitiveDetector(aCurrentVolume->GetLogicalVolume()->GetSensitiveDetector());
     G4VSensitiveDetector *aSensitive = theGflashStep->GetPreStepPoint()->GetSensitiveDetector();
 
     if (aSensitive == nullptr)
       continue;
 
     theGflashStep->SetTotalEnergyDeposit(spotIter->getEnergy());
     aSensitive->Hit(theGflashStep);
   }
 }

G4bool GflashEMShowerModel::ModelTrigger ( const G4FastTrack & fastTrack )

override

Definition at line 48 of file GflashEMShowerModel.cc.

References excludeDetectorRegion(), GetVolume(), GeV, and theRegion.

                                                                      {
   // Mininum energy cutoff to parameterize
   if (fastTrack.GetPrimaryTrack()->GetKineticEnergy() < GeV) {
     return false;
   }
   if (excludeDetectorRegion(fastTrack)) {
     return false;
   }
 
   // This will be changed accordingly when the way
   // dealing with CaloRegion changes later.
   G4VPhysicalVolume *pCurrentVolume = (fastTrack.GetPrimaryTrack()->GetTouchable())->GetVolume();
   if (pCurrentVolume == nullptr) {
     return false;
   }
 
   G4LogicalVolume *lv = pCurrentVolume->GetLogicalVolume();
   if (lv->GetRegion() != theRegion) {
     return false;
   }
   return true;
 }

void GflashEMShowerModel::updateGflashStep	(	const G4ThreeVector &	position,
		G4double	time
	)

private

Definition at line 136 of file GflashEMShowerModel.cc.

References theGflashStep, and theGflashTouchableHandle.

Referenced by makeHits().

                                                                                                  {
   theGflashStep->GetPostStepPoint()->SetGlobalTime(timeGlobal);
   theGflashStep->GetPreStepPoint()->SetPosition(spotPosition);
   theGflashStep->GetPostStepPoint()->SetPosition(spotPosition);
   theGflashStep->GetPreStepPoint()->SetTouchableHandle(theGflashTouchableHandle);
 }