Par02FastSimModelHCal Class Reference

#include <Par02FastSimModelHCal.hh>

Inheritance diagram for Par02FastSimModelHCal:

Collaboration diagram for Par02FastSimModelHCal:

Public Member Functions
	Par02FastSimModelHCal (G4String aModelName, G4Region *aEnvelope, Par02DetectorParametrisation::Parametrisation aParamType)
	Par02FastSimModelHCal (G4String aModelName, G4Region *aEnvelope)
	Par02FastSimModelHCal (G4String aModelName)
	~Par02FastSimModelHCal ()
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

Detailed Description

Shortcut to the ordinary tracking for hadronic calorimeters.

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 hadronic 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 Par02FastSimModelHCal.hh.

Constructor & Destructor Documentation

Par02FastSimModelHCal::Par02FastSimModelHCal	(	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 47 of file Par02FastSimModelHCal.cc.

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

Par02FastSimModelHCal::Par02FastSimModelHCal	(	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 54 of file Par02FastSimModelHCal.cc.

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

Par02FastSimModelHCal::Par02FastSimModelHCal

(

G4String

aModelName

)

A constructor.

Parameters

aModelName

A name of the fast simulation model.

Definition at line 61 of file Par02FastSimModelHCal.cc.

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

Par02FastSimModelHCal::~Par02FastSimModelHCal

(

)

Definition at line 67 of file Par02FastSimModelHCal.cc.

{}

Member Function Documentation

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

virtual

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

Parameters

aFastTrack	A track.
aFastStep	A step.

Implements G4VFastSimulationModel.

Definition at line 93 of file Par02FastSimModelHCal.cc.

                                                          {
  //G4cout << " ________HCal model triggered _________" << G4endl;

  // Kill the parameterised particle at the entrance of the hadronic calorimeter
  aFastStep.KillPrimaryTrack();
  aFastStep.ProposePrimaryTrackPathLength( 0.0 );
  G4double Edep = aFastTrack.GetPrimaryTrack()->GetKineticEnergy();

  // Consider only primary tracks (do nothing else for secondary hadrons)
  G4ThreeVector Pos = aFastTrack.GetPrimaryTrack()->GetPosition();
  if ( ! aFastTrack.GetPrimaryTrack()->GetParentID() ) {
    Par02EventInformation* info = (Par02EventInformation*) 
                            G4EventManager::GetEventManager()->GetUserInformation();
    if ( info->GetDoSmearing() ) {
      // Smearing according to the hadronic calorimeter resolution
      G4ThreeVector Porg = aFastTrack.GetPrimaryTrack()->GetMomentum();
      G4double res = fCalculateParametrisation->
        GetResolution( Par02DetectorParametrisation::eHCAL, 
                       fParametrisation, Porg.mag() );
      G4double eff = fCalculateParametrisation->
        GetEfficiency( Par02DetectorParametrisation::eHCAL, 
                       fParametrisation, Porg.mag() );
      G4double Esm;
      Esm = std::abs( Par02Smearer::Instance()->
                        SmearEnergy( aFastTrack.GetPrimaryTrack(), res ) );
      Par02Output::Instance()->FillHistogram( 2, (Esm/MeV) / (Edep/MeV) );
      // Setting the values of Pos, Esm, res and eff
      Par02PrimaryParticleInformation* primaryInfo= 
         static_cast<Par02PrimaryParticleInformation*>(
           ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) ;
      primaryInfo->SetHCalPosition( Pos );
      primaryInfo->SetHCalEnergy( Esm );
      primaryInfo->SetHCalResolution( res );
      primaryInfo->SetHCalEfficiency( eff );
      // The (smeared) energy of the particle is deposited in the step
      // (which corresponds to the entrance of the hadronic calorimeter)
      aFastStep.ProposeTotalEnergyDeposited( Esm );
    } else {
      // No smearing: simply setting the value of Edep
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetHCalEnergy( Edep );
      // The (initial) energy of the particle is deposited in the step
      // (which corresponds to the entrance of the hadronic calorimeter)
      aFastStep.ProposeTotalEnergyDeposited( Edep );
    }
  }
}

Here is the call graph for this function:

G4bool Par02FastSimModelHCal::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 71 of file Par02FastSimModelHCal.cc.

                                                                                      {
  G4bool isOk = false;
  // Applicable to all hadrons, i.e. any particle made of quarks
  if ( aParticleType.GetQuarkContent(1) +
       aParticleType.GetQuarkContent(2) +
       aParticleType.GetQuarkContent(3) +
       aParticleType.GetQuarkContent(4) +
       aParticleType.GetQuarkContent(5) +
       aParticleType.GetQuarkContent(6) != 0 ) {
    isOk = true;
  }
   return isOk;
}

Here is the call graph for this function:

G4bool Par02FastSimModelHCal::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 87 of file Par02FastSimModelHCal.cc.

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

---

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

• source/geant4.10.03.p02/examples/extended/parameterisations/Par02/include/Par02FastSimModelHCal.hh
• source/geant4.10.03.p02/examples/extended/parameterisations/Par02/src/Par02FastSimModelHCal.cc