InputConverter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_INPUTCONVERTER_H
#define ISF_INPUTCONVERTER_H 1
 
// STL includes
#include <string>
// ISF include
#include "ISF_Interfaces/IInputConverter.h"
// Geant4 includes
#include <G4Event.hh>
// FrameWork includes
#include "AthenaBaseComps/AthService.h"
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/ToolHandle.h"
 
#include "AtlasHepMC/GenEvent_fwd.h"
#include "AtlasHepMC/GenParticle_fwd.h"
#include "TruthUtils/MagicNumbers.h"
#include "GeneratorModules/GenData.h"
 
class McEventCollection;
namespace ISFTesting {
  class InputConverter_test;
}
namespace ISF {
  class ISFParticle;
  class IGenParticleFilter;
}
 
class G4ParticleDefinition;
class G4PrimaryParticle;
class G4VSolid;
 
namespace ISF {

  class InputConverter final: public extends<AthService, IInputConverter> {
 
    // allow test to access private data
    friend ISFTesting::InputConverter_test;
 
  public:
    InputConverter(const std::string& name, ISvcLocator* svc);
    virtual ~InputConverter();

    virtual StatusCode  initialize() override final;
    virtual StatusCode  finalize() override final;

    virtual StatusCode convert(McEventCollection& inputGenEvents,
                               ISF::ISFParticleContainer& simParticles) const override final;

    virtual StatusCode convertHepMCToG4Event(
        McEventCollection& inputGenEvents, G4Event& outputG4Event,
        McEventCollection& shadowGenEvents) const override final;
 
    virtual StatusCode convertHepMCToG4EventLegacy(
        McEventCollection& inputGenEvents,
        G4Event& outputG4Event) const override final;

    void ISF_to_G4Event(G4Event& event,
                        const std::vector<ISF::ISFParticle*>& isp,
                        HepMC::GenEvent* genEvent,
                        HepMC::GenEvent* shadowGenEvent = nullptr,
                        bool useHepMC = false) const override final;
 
   private:
    const G4ParticleDefinition* getG4ParticleDefinition(int pdgcode) const;
 
    G4PrimaryParticle* getDaughterG4PrimaryParticle(const HepMC::ConstGenParticlePtr& gp) const;
    G4PrimaryParticle* getDaughterG4PrimaryParticle(const HepMC::GenParticlePtr& gp, bool makeLinkToTruth=true) const;
 
    G4PrimaryParticle* getG4PrimaryParticle(ISF::ISFParticle& isp, bool useHepMC, HepMC::GenEvent *shadowGenEvent) const;
 
    void addG4PrimaryVertex(G4Event& g4evt, ISF::ISFParticle& isp,
                            bool useHepMC,
                            HepMC::GenEvent* shadowGenEvent) const;
 
    void processPredefinedDecays(const HepMC::ConstGenParticlePtr& genpart, ISF::ISFParticle& isp, G4PrimaryParticle* g4particle) const;
    void processPredefinedDecays(const HepMC::GenParticlePtr& genpart, ISF::ISFParticle& isp, G4PrimaryParticle* g4particle, bool makeLinkToTruth=true) const;
 
    bool matchedGenParticles(const HepMC::ConstGenParticlePtr& p1,
                             const HepMC::ConstGenParticlePtr& p2) const;
 
    HepMC::GenParticlePtr findShadowParticle(const HepMC::ConstGenParticlePtr& genParticle, HepMC::GenEvent *shadowGenEvent) const;

    bool isInsideG4WorldVolume(const ISF::ISFParticle& isp, const G4VSolid* worldSolid) const;

    double getParticleMass(const HepMC::ConstGenParticlePtr& p) const;

    std::vector<HepMC::GenParticlePtr > getSelectedParticles(HepMC::GenEvent& evnt, bool legacyOrdering=false) const;

    bool passesFilters(const HepMC::ConstGenParticlePtr& p) const;

    ISF::ISFParticle* convertParticle(const HepMC::GenParticlePtr& genPartPtr) const;

    std::shared_ptr<GenData> m_gendata{nullptr};
 
    bool                                  m_useGeneratedParticleMass; 
 
    ToolHandleArray<IGenParticleFilter>   m_genParticleFilters;       
 
    bool                                  m_quasiStableParticlesIncluded; //<! will quasi-stable particles be included in the simulation
    BooleanProperty m_useShadowEvent{this, "UseShadowEvent", false, "New approach to selecting particles for simulation" };
 
  };
 
}
 
 
#endif //> !ISF_INPUTCONVERTER_H