LowEnergyFastSimModel Class Reference

#include <LowEnergyFastSimModel.h>

Inheritance diagram for LowEnergyFastSimModel:

Public Member Functions
void	DoIt (const G4FastTrack &fastTrack, G4FastStep &fastStep) override
G4bool	IsApplicable (const G4ParticleDefinition &particle) override
	LowEnergyFastSimModel (const G4String &name, G4Region *region, const edm::ParameterSet &parSet)
G4bool	ModelTrigger (const G4FastTrack &fastTrack) override

Private Attributes
G4double	fEmax
GFlashHitMaker	fHitMaker
const G4Envelope *	fRegion
const TrackingAction *	fTrackingAction
LowEnergyFastSimParam	param

Detailed Description

Definition at line 15 of file LowEnergyFastSimModel.h.

Constructor & Destructor Documentation

LowEnergyFastSimModel()

LowEnergyFastSimModel::LowEnergyFastSimModel	(	const G4String &	name,
		G4Region *	region,
		const edm::ParameterSet &	parSet
	)

Definition at line 16 of file LowEnergyFastSimModel.cc.

    : G4VFastSimulationModel(name, region), fRegion(region), fTrackingAction(nullptr) {
  fEmax = parSet.getParameter<double>("LowEnergyGflashEcalEmax") * CLHEP::GeV;
}

References fEmax, edm::ParameterSet::getParameter(), and GeV.

Member Function Documentation

DoIt()

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

override

Definition at line 37 of file LowEnergyFastSimModel.cc.

                                                                                   {
  fastStep.KillPrimaryTrack();
  fastStep.SetPrimaryTrackPathLength(0.0);
  double energy = fastTrack.GetPrimaryTrack()->GetKineticEnergy();
 
  const G4ThreeVector& pos = fastTrack.GetPrimaryTrack()->GetPosition();
 
  G4double inPointEnergy = param.GetInPointEnergyFraction(energy) * energy;
 
  // take into account positron annihilation (not included in in-point)
  if (-11 == fastTrack.GetPrimaryTrack()->GetDefinition()->GetPDGEncoding())
    energy += twomass;
 
  const G4ThreeVector& momDir = fastTrack.GetPrimaryTrack()->GetMomentumDirection();
  const G4ThreeVector& ortho = momDir.orthogonal();
  const G4ThreeVector& cross = momDir.cross(ortho);
 
  // in point energy deposition
  GFlashEnergySpot spot;
  spot.SetEnergy(inPointEnergy);
  spot.SetPosition(pos);
  fHitMaker.make(&spot, &fastTrack);
 
  // tail energy deposition
  G4double etail = energy - inPointEnergy;
  const G4int nspots = int(etail) + 1;
  const G4double tailEnergy = etail / (G4double)nspots;
  for (G4int i = 0; i < nspots; ++i) {
    const G4double radius = param.GetRadius(energy);
    const G4double z = param.GetZ();
 
    const G4double phi = CLHEP::twopi * G4UniformRand();
    const G4ThreeVector tailPos = pos + z * momDir + radius * std::cos(phi) * ortho + radius * std::sin(phi) * cross;
 
    spot.SetEnergy(tailEnergy);
    spot.SetPosition(tailPos);
    fHitMaker.make(&spot, &fastTrack);
  }
}

References funct::cos(), Basic3DVector< T >::cross(), cross(), HCALHighEnergyHPDFilter_cfi::energy, fHitMaker, LowEnergyFastSimParam::GetInPointEnergyFraction(), LowEnergyFastSimParam::GetRadius(), LowEnergyFastSimParam::GetZ(), mps_fire::i, createfilelist::int, param, phi, CosmicsPD_Skims::radius, funct::sin(), twomass, and z.

IsApplicable()

G4bool LowEnergyFastSimModel::IsApplicable ( const G4ParticleDefinition &  particle )

override

Definition at line 21 of file LowEnergyFastSimModel.cc.

                                                                               {
  return (11 == std::abs(particle.GetPDGEncoding()));
}

References funct::abs().

ModelTrigger()

G4bool LowEnergyFastSimModel::ModelTrigger ( const G4FastTrack &  fastTrack )

override

Definition at line 25 of file LowEnergyFastSimModel.cc.

                                                                       {
  const G4Track* track = fastTrack.GetPrimaryTrack();
  if (nullptr == fTrackingAction) {
    fTrackingAction = static_cast<const TrackingAction*>(G4EventManager::GetEventManager()->GetUserTrackingAction());
  }
  int pdgMother = std::abs(fTrackingAction->geant4Track()->GetDefinition()->GetPDGEncoding());
  if (pdgMother == 11 || pdgMother == 22)
    return false;
  G4double energy = track->GetKineticEnergy();
  return energy < fEmax && fRegion == fastTrack.GetEnvelope();
}

References funct::abs(), HCALHighEnergyHPDFilter_cfi::energy, fEmax, fRegion, fTrackingAction, TrackingAction::geant4Track(), and HLT_FULL_cff::track.

Member Data Documentation

fEmax

G4double LowEnergyFastSimModel::fEmax

private

Definition at line 24 of file LowEnergyFastSimModel.h.

Referenced by LowEnergyFastSimModel(), and ModelTrigger().

fHitMaker

GFlashHitMaker LowEnergyFastSimModel::fHitMaker

private

Definition at line 27 of file LowEnergyFastSimModel.h.

Referenced by DoIt().

fRegion

const G4Envelope* LowEnergyFastSimModel::fRegion

private

Definition at line 25 of file LowEnergyFastSimModel.h.

Referenced by ModelTrigger().

fTrackingAction

const TrackingAction* LowEnergyFastSimModel::fTrackingAction

private

Definition at line 26 of file LowEnergyFastSimModel.h.

Referenced by ModelTrigger().

param

LowEnergyFastSimParam LowEnergyFastSimModel::param

private

Definition at line 28 of file LowEnergyFastSimModel.h.

Referenced by DoIt().