Par02FastSimModelEMCal Class Reference

#include <Par02FastSimModelEMCal.hh>

Inheritance diagram for Par02FastSimModelEMCal:

Collaboration diagram for Par02FastSimModelEMCal:

Public Member Functions
	Par02FastSimModelEMCal (G4String aModelName, G4Region *aEnvelope, Par02DetectorParametrisation::Parametrisation aParamType)
	Par02FastSimModelEMCal (G4String aModelName, G4Region *aEnvelope)
	Par02FastSimModelEMCal (G4String aModelName)
	~Par02FastSimModelEMCal ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &aParticle)
virtual G4bool	ModelTrigger (const G4FastTrack &aFastTrack)
virtual void	DoIt (const G4FastTrack &aFastTrack, G4FastStep &aFastStep)
Public Member Functions inherited from G4VFastSimulationModel
	G4VFastSimulationModel (const G4String &aName)
	G4VFastSimulationModel (const G4String &aName, G4Envelope *, G4bool IsUnique=FALSE)
virtual	~G4VFastSimulationModel ()
virtual G4bool	AtRestModelTrigger (const G4FastTrack &)
virtual void	AtRestDoIt (const G4FastTrack &, G4FastStep &)
const G4String	GetName () const
G4bool	operator== (const G4VFastSimulationModel &) const

Private Attributes
Par02DetectorParametrisation *	fCalculateParametrisation
Par02DetectorParametrisation::Parametrisation	fParametrisation
	A parametrisation type. More...

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 47 of file Par02FastSimModelEMCal.hh.

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 51 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 58 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 65 of file Par02FastSimModelEMCal.cc.

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

~Par02FastSimModelEMCal()

Par02FastSimModelEMCal::~Par02FastSimModelEMCal

(

)

Definition at line 71 of file Par02FastSimModelEMCal.cc.

{}

Member Function Documentation

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.

Implements G4VFastSimulationModel.

Definition at line 91 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() ) {
    Par02EventInformation* 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 );
    }
  }
}

Here is the call graph for this function:

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).

Implements G4VFastSimulationModel.

Definition at line 75 of file Par02FastSimModelEMCal.cc.

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

Here is the call graph for this function:

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.

Implements G4VFastSimulationModel.

Definition at line 85 of file Par02FastSimModelEMCal.cc.

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

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 88 of file Par02FastSimModelEMCal.hh.

fParametrisation

Par02DetectorParametrisation::Parametrisation Par02FastSimModelEMCal::fParametrisation

private

A parametrisation type.

Definition at line 91 of file Par02FastSimModelEMCal.hh.

---

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

• Par02FastSimModelEMCal.hh
• Par02FastSimModelEMCal.cc