Shortcut to the ordinary tracking for tracking detectors. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par02FastSimModelTracker:

Public Member Functions
	Par02FastSimModelTracker (G4String aModelName, G4Region *aEnvelope, Par02DetectorParametrisation::Parametrisation aParamType)
	A constructor.

	Par02FastSimModelTracker (G4String aModelName, G4Region *aEnvelope)
	A constructor.

	Par02FastSimModelTracker (G4String aModelName)
	A constructor.

	~Par02FastSimModelTracker ()

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)
	Calculates the final position (at the outer boundary of the tracking detector) of a particle with the momentum at the entrance of the tracking detector.

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

Par02DetectorParametrisation::Parametrisation	fParametrisation
	A parametrisation type.

Detailed Description

Shortcut to the ordinary tracking for tracking detectors.

The fast simulation model describes what should be done instead of a normal tracking. Instead of the ordinary tracking, a particle momentum at the entrance of the tracking detector is smeared (by Par02Smearer::SmearMomentum()) and the particle is placed at the tracking detector exit, at the place it would reach without the change of its momentum. Based on G4 examples/extended/parametrisations/Par01/include/Par01EMShowerModel.hh .

Author: Anna Zaborowska

Definition at line 43 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

◆ Par02FastSimModelTracker() [1/3]

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

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

◆ Par02FastSimModelTracker() [2/3]

Par02FastSimModelTracker::Par02FastSimModelTracker	(	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 61 of file Par02FastSimModelTracker.cc.

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

◆ Par02FastSimModelTracker() [3/3]

Par02FastSimModelTracker::Par02FastSimModelTracker ( G4String aModelName )

A constructor.

Parameters

aModelName A name of the fast simulation model.

Definition at line 68 of file Par02FastSimModelTracker.cc.

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

◆ ~Par02FastSimModelTracker()

Par02FastSimModelTracker::~Par02FastSimModelTracker ( )

default

Member Function Documentation

◆ IsApplicable()

G4bool Par02FastSimModelTracker::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 78 of file Par02FastSimModelTracker.cc.

                                                                                   {
  return aParticleType.GetPDGCharge() != 0;  // Applicable for all charged particles
}

◆ ModelTrigger()

G4bool Par02FastSimModelTracker::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 85 of file Par02FastSimModelTracker.cc.

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

◆ DoIt()

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

virtual

Calculates the final position (at the outer boundary of the tracking detector) of a particle with the momentum at the entrance of the tracking detector.

Smears the particle momentum and saves it, together with the tracking detector resolution and efficiency to the Par02PrimaryParticleInformation.

Parameters

aFastTrack	A track.
aFastStep	A step.

Definition at line 91 of file Par02FastSimModelTracker.cc.

                                                             {
 
  G4cout << " ________Tracker model triggered _________" << G4endl;
 
  // Calculate the final position (at the outer boundary of the tracking detector)
  // of the particle with the momentum at the entrance of the tracking detector.
 
  G4Track track = * aFastTrack.GetPrimaryTrack();
  G4FieldTrack aFieldTrack( '0' );
  G4FieldTrackUpdator::Update( &aFieldTrack, &track );
 
  G4double retSafety = -1.0;
  ELimited retStepLimited;
  G4FieldTrack endTrack( 'a' );
  G4double currentMinimumStep = 10.0*m;  // Temporary: change that to sth connected
                                         // to particle momentum.
  G4PathFinder* fPathFinder = G4PathFinder::GetInstance();
  /*G4double lengthAlongCurve = */ 
  fPathFinder->ComputeStep( aFieldTrack,
                            currentMinimumStep,
                            0,
                            aFastTrack.GetPrimaryTrack()->GetCurrentStepNumber(),
                            retSafety,
                            retStepLimited,
                            endTrack,
                            aFastTrack.GetPrimaryTrack()->GetVolume() );
 
  // Place the particle at the tracking detector exit 
  // (at the place it would reach without the change of its momentum).
  aFastStep.ProposePrimaryTrackFinalPosition( endTrack.GetPosition() );
 
  // Consider only primary tracks (do nothing else for secondary charged particles)
  G4ThreeVector Porg = aFastTrack.GetPrimaryTrack()->GetMomentum();
  if ( ! aFastTrack.GetPrimaryTrack()->GetParentID() ) {
    auto  info = (Par02EventInformation*) 
                            G4EventManager::GetEventManager()->GetUserInformation();
    if ( info->GetDoSmearing() ) {
      // Smearing according to the tracking detector resolution
      G4double res = fCalculateParametrisation->
        GetResolution( Par02DetectorParametrisation::eTRACKER, 
                       fParametrisation, Porg.mag() );
      G4double eff = fCalculateParametrisation->
        GetEfficiency( Par02DetectorParametrisation::eTRACKER,
                       fParametrisation, Porg.mag() );
      G4ThreeVector Psm;
      Psm = Par02Smearer::Instance()->
                                 SmearMomentum( aFastTrack.GetPrimaryTrack(), res );
      Par02Output::Instance()->FillHistogram( 0, ((Psm.mag()/MeV) / (Porg.mag()/MeV)) );
      // Setting the values of Psm, res and eff
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetTrackerMomentum( Psm );
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetTrackerResolution( res );
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetTrackerEfficiency( eff );
    } else {
      // No smearing: simply setting the value of Porg
      ( (Par02PrimaryParticleInformation*) ( const_cast< G4PrimaryParticle* >
          ( aFastTrack.GetPrimaryTrack()->GetDynamicParticle()->GetPrimaryParticle() )->
            GetUserInformation() ) )->SetTrackerMomentum( Porg );
    }
  }
}

Member Data Documentation

◆ fCalculateParametrisation

Par02DetectorParametrisation* Par02FastSimModelTracker::fCalculateParametrisation

private

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

Definition at line 92 of file Par02FastSimModelTracker.hh.