SimulationEnergies.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef CaloG4_SimulationEnergies_H
#define CaloG4_SimulationEnergies_H
 
#include "G4TrackStatus.hh"
#include "G4ParticleDefinition.hh"
#include "G4DynamicParticle.hh"
#include "G4ThreeVector.hh"
#include "globals.hh"
 
#include <map>
#include <memory>
#include <vector>
 
// Forward declarations.
class G4Step;
class G4ParticleDefinition;
class G4Track;
 
namespace CaloG4
{

  class SimulationEnergies
  {
  public:
 
    SimulationEnergies();
    virtual ~SimulationEnergies();

    void Energies( const G4Step* , std::vector<G4double> & ) const;

    enum eEnergyCategory
      { kEm = 0,
        kNonEm,
        kInvisible0,
        /* 12-Aug-2009 we compared three different techniques
        kInvisible1, //  to determine the energy lost
        kInvisible2, // due to "invisible" processes in the detector
             and chose one of them for production runs
        */
        kEscaped,
        kNumberOfEnergyCategories };

    typedef struct {
      std::map<eEnergyCategory, G4double> energy;
    } ClassifyResult_t;
 
    /* 12-Aug-2009 We keep only info about energy in ClassifyResult_t
                   for production runs (see above).
    typedef struct {
      G4double time;
      std::map<eEnergyCategory, G4double> energy;
      G4TrackStatus trackStatus;
      const G4ParticleDefinition* particle;
      const G4DynamicParticle* dynamicParticle;
      G4String processName;
    } ClassifyResult_t;
    */

    ClassifyResult_t Classify( const G4Step*, const G4bool processEscaped = true ) const;
 
  private:

 
    G4double measurableEnergy(const G4ParticleDefinition *particleDef,
                              G4int PDGEncoding,
                              G4double totalEnergy,
                              G4double kineticEnergy) const;
 
    G4double measurableEnergyV2(const G4ParticleDefinition *particleDef,
                                G4int PDGEncoding,
                                G4double totalEnergy,
                                G4double kineticEnergy) const;
 
    // Escaped energy requires special processing.
    G4bool ProcessEscapedEnergy( G4Step* fakeStep ) const;
 
    std::unique_ptr<G4Step> CreateFakeStep(G4Track* a_track, G4double a_energy) const;
 
  }; // class SimulationEnergies
 
} // namespace CaloG4
 
#endif // CaloG4_SimulationEnergies_H