Par02FastSimModelEMCal Class Reference

Shortcut to the ordinary tracking for electromagnetic calorimeters. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par02FastSimModelEMCal:

Public Member Functions
	Par02FastSimModelEMCal (G4String aModelName, G4Region *aEnvelope, Par02DetectorParametrisation::Parametrisation aParamType)
	A constructor.
	Par02FastSimModelEMCal (G4String aModelName, G4Region *aEnvelope)
	A constructor.
	Par02FastSimModelEMCal (G4String aModelName)
	A constructor.
	~Par02FastSimModelEMCal ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &aParticle)
	Checks if this model should be applied to this particle type.
virtual G4bool	ModelTrigger (const G4FastTrack &aFastTrack)
	Checks if the model should be applied, taking into account the kinematics of a track.
virtual void	DoIt (const G4FastTrack &aFastTrack, G4FastStep &aFastStep)
	Smears the energy deposit and saves it, together with the position of the deposit, the electromagnetic calorimeter resolution and efficiency to the Par02PrimaryParticleInformation.

Private Attributes
Par02DetectorParametrisation *	fCalculateParametrisation
	A pointer to Par02DetectorParametrisation used to get the efficiency and resolution of the detector for a given particle and parametrisation type.
Par02DetectorParametrisation::Parametrisation	fParametrisation
	A parametrisation type.

Detailed Description

Shortcut to the ordinary tracking for electromagnetic calorimeters.

The fast simulation model describes what should be done instead of a normal tracking. Instead of the ordinary tracking, a particle deposits its energy at the entrance to the electromagnetic calorimeter and its value is smeared (by Par02Smearer::SmearMomentum()). Based on G4 examples/extended/parametrisations/Par01/include/Par01EMShowerModel.hh .

Author

Anna Zaborowska

Definition at line 41 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

Par02FastSimModelEMCal() [1/3]

Par02FastSimModelEMCal::Par02FastSimModelEMCal	(	G4String	aModelName,
		G4Region *	aEnvelope,
		Par02DetectorParametrisation::Parametrisation	aParamType
	)

A constructor.

Parameters

aModelName	A name of the fast simulation model.
aEnvelope	A region where the model can take over the ordinary tracking.
aParamType	A parametrisation type.

Definition at line 50 of file Par02FastSimModelEMCal.cc.

                                                                           :
  G4VFastSimulationModel( aModelName, aEnvelope ), fCalculateParametrisation(),
  fParametrisation( aType ) {}

Par02FastSimModelEMCal() [2/3]

Par02FastSimModelEMCal::Par02FastSimModelEMCal	(	G4String	aModelName,
		G4Region *	aEnvelope
	)

A constructor.

Parameters

aModelName	A name of the fast simulation model.
aEnvelope	A region where the model can take over the ordinary tracking.

Definition at line 57 of file Par02FastSimModelEMCal.cc.

                                                                      : 
  G4VFastSimulationModel( aModelName, aEnvelope ), fCalculateParametrisation(),
  fParametrisation( Par02DetectorParametrisation::eCMS ) {}

Par02FastSimModelEMCal() [3/3]

Par02FastSimModelEMCal::Par02FastSimModelEMCal

(

G4String

aModelName

)

A constructor.

Parameters

aModelName

A name of the fast simulation model.

Definition at line 64 of file Par02FastSimModelEMCal.cc.

                                                                    :
  G4VFastSimulationModel( aModelName ), fCalculateParametrisation(), 
  fParametrisation( Par02DetectorParametrisation::eCMS ) {}

~Par02FastSimModelEMCal()

Par02FastSimModelEMCal::~Par02FastSimModelEMCal ( )

default

Member Function Documentation

IsApplicable()

G4bool Par02FastSimModelEMCal::IsApplicable ( const G4ParticleDefinition &  aParticle )

virtual

Checks if this model should be applied to this particle type.

Parameters

aParticle

A particle definition (type).

Definition at line 74 of file Par02FastSimModelEMCal.cc.

                                              {
  // Applicable for electrons, positrons, and gammas
  return &aParticleType == G4Electron::Definition()  ||
         &aParticleType == G4Positron::Definition()  ||
         &aParticleType == G4Gamma::Definition();
}

ModelTrigger()

G4bool Par02FastSimModelEMCal::ModelTrigger ( const G4FastTrack &  aFastTrack )

virtual

Checks if the model should be applied, taking into account the kinematics of a track.

Parameters

aFastTrack

A track.

Definition at line 84 of file Par02FastSimModelEMCal.cc.

                                                                  {
  return true;  // No kinematical restrictions to apply the parametrisation
}

DoIt()

void Par02FastSimModelEMCal::DoIt ( const G4FastTrack &  aFastTrack, G4FastStep &  aFastStep  )

virtual

Smears the energy deposit and saves it, together with the position of the deposit, the electromagnetic calorimeter resolution and efficiency to the Par02PrimaryParticleInformation.

Parameters

aFastTrack	A track.
aFastStep	A step.

Definition at line 90 of file Par02FastSimModelEMCal.cc.

                                                           {
  //G4cout << " ________EMCal model triggered _________" << G4endl;
 
  // Kill the parameterised particle at the entrance of the electromagnetic calorimeter
  aFastStep.KillPrimaryTrack();
  aFastStep.ProposePrimaryTrackPathLength( 0.0 );
  G4double Edep = aFastTrack.GetPrimaryTrack()->GetKineticEnergy();
 
  // Consider only primary tracks (do nothing else for secondary e-, e+, gammas)
  G4ThreeVector Pos = aFastTrack.GetPrimaryTrack()->GetPosition();
  if ( ! aFastTrack.GetPrimaryTrack()->GetParentID() ) {
    auto  info = (Par02EventInformation*) 
                            G4EventManager::GetEventManager()->GetUserInformation();
    if ( info->GetDoSmearing() ) {
      // Smearing according to the electromagnetic calorimeter resolution
      G4ThreeVector Porg = aFastTrack.GetPrimaryTrack()->GetMomentum();
      G4double res = fCalculateParametrisation->GetResolution( 
        Par02DetectorParametrisation::eEMCAL, fParametrisation, Porg.mag() );
      G4double eff = fCalculateParametrisation->GetEfficiency( 
        Par02DetectorParametrisation::eEMCAL, fParametrisation, Porg.mag() );
      G4double Esm;
      Esm = std::abs( Par02Smearer::Instance()->
                        SmearEnergy( aFastTrack.GetPrimaryTrack(), res ) );
      Par02Output::Instance()->FillHistogram( 1, (Esm/MeV) / (Edep/MeV) );
      // Setting the values of Pos, Esm, res and eff
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetEMCalPosition( Pos );
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetEMCalEnergy( Esm );
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetEMCalResolution( res );
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetEMCalEfficiency( eff );
      // The (smeared) energy of the particle is deposited in the step
      // (which corresponds to the entrance of the electromagnetic calorimeter)
      aFastStep.ProposeTotalEnergyDeposited( Esm );
    } else {
      // No smearing: simply setting the value of Edep
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetEMCalEnergy( Edep );
      // The (initial) energy of the particle is deposited in the step
      // (which corresponds to the entrance of the electromagnetic calorimeter)
      aFastStep.ProposeTotalEnergyDeposited( Edep );
    }
  }
}

Member Data Documentation

fCalculateParametrisation

Par02DetectorParametrisation* Par02FastSimModelEMCal::fCalculateParametrisation

private

A pointer to Par02DetectorParametrisation used to get the efficiency and resolution of the detector for a given particle and parametrisation type.

Definition at line 87 of file Par02FastSimModelEMCal.hh.

fParametrisation

Par02DetectorParametrisation::Parametrisation Par02FastSimModelEMCal::fParametrisation

private

A parametrisation type.

Definition at line 90 of file Par02FastSimModelEMCal.hh.

---

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

• .doxygen/Doxymodules_parameterisations.h
• extended/parameterisations/Par02/include/Par02FastSimModelEMCal.hh
• extended/parameterisations/Par02/src/Par02FastSimModelEMCal.cc