GFlashEMShowerModel Class Reference

#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
bool	theWatcherOn

Detailed Description

Definition at line 29 of file GFlashEMShowerModel.h.

Constructor & Destructor Documentation

GFlashEMShowerModel()

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

Definition at line 25 of file GFlashEMShowerModel.cc.

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

References edm::ParameterSet::getParameter(), theGflashNavigator, theGflashStep, theGflashTouchableHandle, theProfile, theRegion, and theWatcherOn.

~GFlashEMShowerModel()

GFlashEMShowerModel::~GFlashEMShowerModel ( )

override

Definition at line 41 of file GFlashEMShowerModel.cc.

                                          {
  delete theProfile;
  delete theGflashStep;
}

References theGflashStep, and theProfile.

Member Function Documentation

DoIt()

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

override

Definition at line 88 of file GFlashEMShowerModel.cc.

                                                                                 {
  // 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);
}

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

excludeDetectorRegion()

G4bool GFlashEMShowerModel::excludeDetectorRegion ( const G4FastTrack &  fastTrack )

private

Definition at line 172 of file GFlashEMShowerModel.cc.

                                                                              {
  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;
}

References PVValHelper::eta.

Referenced by ModelTrigger().

IsApplicable()

G4bool GFlashEMShowerModel::IsApplicable ( const G4ParticleDefinition &  particleType )

override

Definition at line 46 of file GFlashEMShowerModel.cc.

                                                                                 {
  return (&particleType == G4Electron::Electron() || &particleType == G4Positron::Positron());
}

References nanoDQM_cff::Electron, and PbPb_ZMuSkimMuonDPG_cff::particleType.

makeHits()

void GFlashEMShowerModel::makeHits ( const G4FastTrack &  fastTrack )

private

Definition at line 112 of file GFlashEMShowerModel.cc.

                                                               {
  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());
 
    // If there is a watcher defined in a job and the flag is turned on
    if (theWatcherOn) {
      SteppingAction* userSteppingAction = (SteppingAction*)G4EventManager::GetEventManager()->GetUserSteppingAction();
      userSteppingAction->m_g4StepSignal(theGflashStep);
    }
 
    // 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);
  }
}

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

Referenced by DoIt().

ModelTrigger()

G4bool GFlashEMShowerModel::ModelTrigger ( const G4FastTrack &  fastTrack )

override

Definition at line 51 of file GFlashEMShowerModel.cc.

                                                                     {
  // Mininum energy cutoff to parameterize
  if (fastTrack.GetPrimaryTrack()->GetKineticEnergy() < Gflash::energyCutOff) {
    return false;
  }
  if (excludeDetectorRegion(fastTrack)) {
    return false;
  }
 
  // This will be changed accordingly when the way
  // dealing with CaloRegion changes later.
  const G4TouchableHistory* touch = static_cast<const G4TouchableHistory*>(fastTrack.GetPrimaryTrack()->GetTouchable());
  const G4VPhysicalVolume* pCurrentVolume = touch->GetVolume();
  if (pCurrentVolume == nullptr) {
    return false;
  }
 
  const G4LogicalVolume* lv = pCurrentVolume->GetLogicalVolume();
  if (lv->GetRegion() != theRegion) {
    return false;
  }
  //std::cout << "GFlashEMShowerModel::ModelTrigger: LV "
  //        << lv->GetRegion()->GetName() << std::endl;
 
  // The parameterization starts inside crystals
  //std::size_t pos1 = lv->GetName().find("EBRY");
  //std::size_t pos2 = lv->GetName().find("EFRY");
 
  //std::size_t pos3 = lv->GetName().find("HVQ");
  //std::size_t pos4 = lv->GetName().find("HF");
  //if(pos1 == std::string::npos && pos2 == std::string::npos &&
  //   pos3 == std::string::npos && pos4 == std::string::npos) return false;
 
  return true;
}

References Gflash::energyCutOff, excludeDetectorRegion(), and theRegion.

updateGflashStep()

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

private

Definition at line 164 of file GFlashEMShowerModel.cc.

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

References theGflashStep, and theGflashTouchableHandle.

Referenced by makeHits().

Member Data Documentation

theGflashNavigator

G4Navigator* GFlashEMShowerModel::theGflashNavigator

private

Definition at line 52 of file GFlashEMShowerModel.h.

Referenced by GFlashEMShowerModel(), and makeHits().

theGflashStep

G4Step* GFlashEMShowerModel::theGflashStep

private

Definition at line 51 of file GFlashEMShowerModel.h.

Referenced by GFlashEMShowerModel(), makeHits(), updateGflashStep(), and ~GFlashEMShowerModel().

theGflashTouchableHandle

G4TouchableHandle GFlashEMShowerModel::theGflashTouchableHandle

private

Definition at line 53 of file GFlashEMShowerModel.h.

Referenced by GFlashEMShowerModel(), makeHits(), and updateGflashStep().

theParSet

edm::ParameterSet GFlashEMShowerModel::theParSet

private

Definition at line 44 of file GFlashEMShowerModel.h.

theProfile

GflashEMShowerProfile* GFlashEMShowerModel::theProfile

private

Definition at line 47 of file GFlashEMShowerModel.h.

Referenced by DoIt(), GFlashEMShowerModel(), makeHits(), and ~GFlashEMShowerModel().

theRegion

const G4Region* GFlashEMShowerModel::theRegion

private

Definition at line 49 of file GFlashEMShowerModel.h.

Referenced by GFlashEMShowerModel(), makeHits(), and ModelTrigger().

theWatcherOn

bool GFlashEMShowerModel::theWatcherOn

private

Definition at line 45 of file GFlashEMShowerModel.h.

Referenced by GFlashEMShowerModel(), and makeHits().