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
G4bool	fCheck
G4double	fEmax
GFlashHitMaker	fHitMaker
const G4Material *	fMaterial
LowEnergyFastSimParam	fParam
const G4ParticleDefinition *	fPositron
const G4Envelope *	fRegion
G4ThreeVector	fTailPos
const TrackingAction *	fTrackingAction

Detailed Description

Definition at line 18 of file LowEnergyFastSimModel.h.

Constructor & Destructor Documentation

LowEnergyFastSimModel()

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

Definition at line 22 of file LowEnergyFastSimModel.cc.

References fEmax, fMaterial, fPositron, edm::ParameterSet::getParameter(), dqmiodumpmetadata::n, DD4hep2DDDName::noNameSpace(), and TableParser::table.

    : G4VFastSimulationModel(name, region),
      fRegion(region),
      fTrackingAction(nullptr),
      fCheck(true),
      fTailPos(0., 0., 0.) {
  fEmax = parSet.getParameter<double>("LowEnergyGflashEcalEmax") * CLHEP::GeV;
  fPositron = G4Positron::Positron();
  fMaterial = nullptr;
  auto table = G4Material::GetMaterialTable();
  for (auto const& mat : *table) {
    G4String nam = (G4String)(DD4hep2DDDName::noNameSpace(mat->GetName()));
    size_t n = nam.size();
    if (n > 4) {
      G4String sn = nam.substr(n - 5, 5);
      if (sn == "PbWO4") {
        fMaterial = mat;
        break;
      }
    }
  }
  G4String nm = (nullptr == fMaterial) ? "not found!" : (G4String)(DD4hep2DDDName::noNameSpace(fMaterial->GetName()));
  edm::LogVerbatim("LowEnergyFastSimModel") << "LowEGFlash material: <" << nm << ">";
}

Member Function Documentation

DoIt()

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

override

Definition at line 74 of file LowEnergyFastSimModel.cc.

References funct::cos(), hcalRecHitTable_cff::energy, fHitMaker, fParam, fPositron, fTailPos, LowEnergyFastSimParam::GetInPointEnergyFraction(), LowEnergyFastSimParam::GetRadius(), LowEnergyFastSimParam::GetZ(), mps_fire::i, phi, alignCSCRings::r, scaleFactor, funct::sin(), HLT_2023v12_cff::track, twomass, and z.

                                                                                   {
  fastStep.KillPrimaryTrack();
  fastStep.ProposePrimaryTrackPathLength(0.0);
  auto track = fastTrack.GetPrimaryTrack();
  G4double energy = track->GetKineticEnergy() * scaleFactor;

  const G4ThreeVector& pos = track->GetPosition();

  G4double inPointEnergy = fParam.GetInPointEnergyFraction(energy) * energy;

  // take into account positron annihilation (not included in in-point)
  if (fPositron == track->GetDefinition())
    energy += twomass;

  const G4ThreeVector& momDir = track->GetMomentumDirection();

  // in point energy deposition
  GFlashEnergySpot spot;
  spot.SetEnergy(inPointEnergy);
  spot.SetPosition(pos);
  fHitMaker.make(&spot, &fastTrack);

  // Russian roulette
  double wt2 = track->GetWeight();
  if (wt2 <= 0.0) {
    wt2 = 1.0;
  }

  // tail energy deposition
  const G4double etail = energy - inPointEnergy;
  const G4int nspots = etail;
  const G4double tailEnergy = etail * wt2 / (nspots + 1);
  /*  
  edm::LogVerbatim("LowEnergyFastSimModel") << track->GetDefinition()->GetParticleName()
                        << " Ekin(MeV)=" << energy << " material: <"
                                            << track->GetMaterial()->GetName() 
                        << "> Elocal=" << inPointEnergy
                        << " Etail=" << tailEnergy
                        << " Nspots=" << nspots+1;
  */
  for (G4int i = 0; i <= nspots; ++i) {
    const G4double r = fParam.GetRadius(energy);
    const G4double z = fParam.GetZ();

    const G4double phi = CLHEP::twopi * G4UniformRand();
    fTailPos.set(r * std::cos(phi), r * std::sin(phi), z);
    fTailPos.rotateUz(momDir);
    fTailPos += pos;

    spot.SetEnergy(tailEnergy);
    spot.SetPosition(fTailPos);
    fHitMaker.make(&spot, &fastTrack);
  }
}

IsApplicable()

G4bool LowEnergyFastSimModel::IsApplicable ( const G4ParticleDefinition &  particle )

override

Definition at line 47 of file LowEnergyFastSimModel.cc.

References funct::abs().

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

ModelTrigger()

G4bool LowEnergyFastSimModel::ModelTrigger ( const G4FastTrack &  fastTrack )

override

Definition at line 51 of file LowEnergyFastSimModel.cc.

References funct::abs(), hcalRecHitTable_cff::energy, fCheck, fEmax, fMaterial, fTrackingAction, TrackInformation::genParticlePID(), and HLT_2023v12_cff::track.

                                                                       {
  const G4Track* track = fastTrack.GetPrimaryTrack();
  G4double energy = track->GetKineticEnergy();
  if (fMaterial != track->GetMaterial() || energy >= fEmax)
    return false;

  /*
  edm::LogVerbatim("LowEnergyFastSimModel") << track->GetDefinition()->GetParticleName()
                        << " Ekin(MeV)=" << energy << " material: <"
                                            << track->GetMaterial()->GetName() << ">";
  */
  if (fCheck) {
    if (nullptr == fTrackingAction) {
      fTrackingAction = static_cast<const TrackingAction*>(G4EventManager::GetEventManager()->GetUserTrackingAction());
    }
    const TrackInformation* ptrti = static_cast<TrackInformation*>(track->GetUserInformation());
    int pdg = ptrti->genParticlePID();
    if (std::abs(pdg) == 11 || pdg == 22)
      return false;
  }
  return true;
}

Member Data Documentation

fCheck

G4bool LowEnergyFastSimModel::fCheck

private

Definition at line 32 of file LowEnergyFastSimModel.h.

Referenced by ModelTrigger().

fEmax

G4double LowEnergyFastSimModel::fEmax

private

Definition at line 27 of file LowEnergyFastSimModel.h.

Referenced by LowEnergyFastSimModel(), and ModelTrigger().

fHitMaker

GFlashHitMaker LowEnergyFastSimModel::fHitMaker

private

Definition at line 34 of file LowEnergyFastSimModel.h.

Referenced by DoIt().

fMaterial

const G4Material* LowEnergyFastSimModel::fMaterial

private

Definition at line 31 of file LowEnergyFastSimModel.h.

Referenced by LowEnergyFastSimModel(), and ModelTrigger().

fParam

LowEnergyFastSimParam LowEnergyFastSimModel::fParam

private

Definition at line 35 of file LowEnergyFastSimModel.h.

Referenced by DoIt().

fPositron

const G4ParticleDefinition* LowEnergyFastSimModel::fPositron

private

Definition at line 30 of file LowEnergyFastSimModel.h.

Referenced by DoIt(), and LowEnergyFastSimModel().

fRegion

const G4Envelope* LowEnergyFastSimModel::fRegion

private

Definition at line 28 of file LowEnergyFastSimModel.h.

fTailPos

G4ThreeVector LowEnergyFastSimModel::fTailPos

private

Definition at line 33 of file LowEnergyFastSimModel.h.

Referenced by DoIt().

fTrackingAction

const TrackingAction* LowEnergyFastSimModel::fTrackingAction

private

Definition at line 29 of file LowEnergyFastSimModel.h.

Referenced by ModelTrigger().