ISFTruthIncident.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_ISFEVENT_ISFTRUTHINCIDENT_H
#define ISF_ISFEVENT_ISFTRUTHINCIDENT_H
 
// stl includes
#include <vector>
 
// ISF includes
#include "ISF_Event/ISFParticleContainer.h"
#include "ISF_Event/ITruthIncident.h"
// DetectorDescription
#include "AtlasDetDescr/AtlasRegion.h" //enum
 
namespace ISF {
  class ISFParticle;

  enum KillPrimary {
    fPrimarySurvives = false,
    fKillsPrimary    = true
  };

 
  class ISFTruthIncident : public ITruthIncident {
  public:
    ISFTruthIncident( ISF::ISFParticle &parent,
                      const ISFParticleVector &children,
                      int process,
                      AtlasDetDescr::AtlasRegion geoID,
                      ISF::KillPrimary killsPrimary = ISF::fPrimarySurvives,
                      const HepMC::FourVector *position = 0);
    virtual ~ISFTruthIncident();
 
    ISFTruthIncident( const ISFTruthIncident& ) = delete;
    ISFTruthIncident& operator= ( const ISFTruthIncident& ) = delete;

    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;
    HepMC::GenParticlePtr      parentParticle() override final;
    int      parentStatus() override final;
    int  parentBarcode() override final; // TODO Remove this method
    int  parentUniqueID() override final;
    bool                      parentSurvivesIncident() const override final;
    HepMC::GenParticlePtr     parentParticleAfterIncident(int newBC) override final;

    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) const override final;
    HepMC::GenParticlePtr     childParticle(unsigned short index,
                                            int bc) override final;
    void updateParentAfterIncidentProperties();
    void updateChildParticleProperties();
 
  private:
    ISFTruthIncident();

    inline HepMC::GenParticlePtr getHepMCTruthParticle( ISF::ISFParticle& particle ) const;

    inline HepMC::GenParticlePtr updateHepMCTruthParticle( ISF::ISFParticle& particle, ISF::ISFParticle* parent=nullptr ) const;
 
    ISF::ISFParticle&                  m_parent;
    const ISFParticleVector&           m_children;
    const int  m_process;
    const ISF::KillPrimary             m_killsPrimary;
    const HepMC::FourVector*           m_position;
  };
 
} // end of namespace
 
#endif // ISF_ISFEVENT_ISFTRUTHINCIDENT_H