Geant4TruthIncident.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_GEANT4TOOLS_Geant4TruthIncident_H
#define ISF_GEANT4TOOLS_Geant4TruthIncident_H
 
// ISF includes
#include "ISF_Event/ITruthIncident.h" //inheritance
 
//G4 includes
#include "G4ThreeVector.hh" //typedef
 
// HepMC includes
#include "AtlasHepMC/SimpleVector.h" //typedef FourVector
#include "AtlasHepMC/GenParticle_fwd.h" //typedef GenParticlePtr
 
#include "AtlasDetDescr/AtlasRegion.h" //enum
 
#include "CxxUtils/checker_macros.h" //ATLAS_THREAD_SAFE
// std
#include <vector>
 
 
// forward declarations
class G4Step;
class G4Track;
class AtlasG4EventUserInfo;
 
namespace ISF {
  class ISFParticle;
}
 
 
namespace iGeant4 {

 
  class Geant4TruthIncident : public ISF::ITruthIncident {
    public:
      Geant4TruthIncident( const G4Step*,
                           const ISF::ISFParticle& baseISP,
                           AtlasDetDescr::AtlasRegion geoID,
                           AtlasG4EventUserInfo* atlasG4EvtUserInfo);
      virtual ~Geant4TruthIncident() {};

      const HepMC::FourVector&  position() const override final;

      int                       physicsProcessCategory() const override final;
      int physicsProcessCode() const override final;

      double                    parentP2() const override final;
      double                    parentPt2() const override final;
      double                    parentEkin() const override final;
      int                       parentPdgCode() const override final;
      int  parentBarcode() override final; // TODO Remove this method
      int  parentUniqueID() override final;
      int  parentStatus() override final;
      bool                      parentSurvivesIncident() const override final;
      HepMC::GenParticlePtr     parentParticleAfterIncident(int newBC) override final;

      const G4ThreeVector       childP(unsigned short index) const;
      double                    childP2(unsigned short index) const override final;
      double                    childPt2(unsigned short index) const override final;
      double                    childEkin(unsigned short index) const override final;
      int                       childPdgCode(unsigned short index) const override final;
      int  childBarcode(unsigned short index) const override final; // TODO Remove - only used in one place in TruthSvc

      ISF::InteractionClass_t interactionClassification() const override final;
 
      // only called once accepted

      HepMC::GenParticlePtr      parentParticle() override final;
      HepMC::GenParticlePtr   childParticle(unsigned short index,
                                            int bc) override final;
    private:
      Geant4TruthIncident();
      inline void prepareChildren();

      inline bool particleAlive(const G4Track *track) const;
 
      HepMC::GenParticlePtr convert(const G4Track *particle, const int barcode, const bool secondary) const;
 
      bool                          m_positionSet;
      HepMC::FourVector             m_position;
      const G4Step*                 m_step{};
      const ISF::ISFParticle&       m_baseISP;
 
      AtlasG4EventUserInfo*         m_atlasG4EvtUserInfo ATLAS_THREAD_SAFE{};
      std::vector<const G4Track*>   m_children;
 
      HepMC::GenParticlePtr         m_parentParticleAfterIncident{};
   };
 
}
 
#endif // ISF_GEANT4TOOLS_Geant4TruthIncident_H