HadIntProcessorParametric.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_FATRASTOOLS_HadIntProcessorParametric_H
#define ISF_FATRASTOOLS_HadIntProcessorParametric_H
 
// GaudiKernel & Athena
#include "AthenaBaseComps/AthAlgTool.h"
#include "AthenaKernel/IAtRndmGenSvc.h"
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/ToolHandle.h"
#include "CxxUtils/checker_macros.h"
 
// Fatras
#include "ISF_FatrasInterfaces/IHadronicInteractionProcessor.h"
 
// Trk
#include "TrkEventPrimitives/PropDirection.h"
 
// ISF
#include "ISF_Event/ISFParticleContainer.h"
 
#ifndef MAXHADINTCHILDREN
#define MAXHADINTCHILDREN 20
#endif
 
class StoreGateSvc;
 
class TTree;
  
namespace Trk {
  class Material;
  class MaterialProperties;
}
 
namespace ISF {
  class IParticleBroker;
  class ITruthSvc;    
}
namespace iFatras {
 
  class IPhysicsValidationTool;
      
  class ATLAS_NOT_THREAD_SAFE HadIntProcessorParametric : public extends<AthAlgTool, iFatras::IHadronicInteractionProcessor> {  // deprecated: ATLASSIM-6020
    public:      
      HadIntProcessorParametric(const std::string&,const std::string&,const IInterface*);
      virtual ~HadIntProcessorParametric();
       
      StatusCode initialize();
      StatusCode finalize();
 
      bool hadronicInteraction(const Amg::Vector3D& position, const Amg::Vector3D& momentum, 
                   double p, double E, double charge,
                               const Trk::MaterialProperties& mprop, double pathCorrection,
                               Trk::ParticleHypothesis particle=Trk::pion) const;
 
      bool recordHadState(double time, double p,
              const Amg::Vector3D& vertex,
              const Amg::Vector3D& particleDir,
              Trk::ParticleHypothesis particle ) const;
 
      bool doHadronicInteraction(double time, const Amg::Vector3D& position, const Amg::Vector3D& momentum, 
                 const Trk::Material* emat,
                 Trk::ParticleHypothesis particle,
                 bool processSecondaries) const;
 
      ISF::ISFParticleVector doHadIntOnLayer(const ISF::ISFParticle* parent, double time,
                         const Amg::Vector3D& position, const Amg::Vector3D& momentum, 
                         const Trk::Material* emat,
                         Trk::ParticleHypothesis particle) const;
 
   private:
      static double absorptionLength(const Trk::MaterialProperties* mat,
            double p, double q, Trk::ParticleHypothesis particle=Trk::pion) ;
 
      ISF::ISFParticleVector getHadState(const ISF::ISFParticle* parent,
                     double time, double p,
                     const Amg::Vector3D& vertex,
                     const Amg::Vector3D& particleDir,
                     Trk::ParticleHypothesis particle ) const;
                           
       double                       m_minimumHadOutEnergy;
       double                       m_childMomParam;
       bool                         m_cutChain;
       bool                         m_hadIntFromX0;
       double                       m_hadIntProbScale;
       double                       m_minimumHadInitialEnergy;
 
 
       ServiceHandle<ISF::IParticleBroker>  m_particleBroker;
       ServiceHandle<ISF::ITruthSvc>        m_truthRecordSvc;
                  
       ServiceHandle<IAtRndmGenSvc>         m_rndGenSvc;
       int          m_processCode;
 
       CLHEP::HepRandomEngine*                     m_randomEngine;
       std::string                          m_randomEngineName; 
 
                    
       bool                          m_validationMode;
       ToolHandle<IPhysicsValidationTool>  m_validationTool;
 
       bool                          m_hadIntValidation;
       std::string                   m_hadIntValidationTreeName;        
       std::string                   m_hadIntValidationTreeDescription; 
       std::string                   m_hadIntValidationTreeFolder;      
                                       
       TTree*                        m_hadIntValidationTree;            
                                       
       // ------ ntuple branches       
       mutable float                 m_hadIntPointX;               
       mutable float                 m_hadIntPointY;               
       mutable float                 m_hadIntPointR;               
       mutable float                 m_hadIntPointZ;               
       mutable int                   m_hadIntMotherPdg;            
       mutable int                   m_hadIntMotherBarcode;         
       mutable float                 m_hadIntMotherP;         
       mutable float                 m_hadIntMotherPt;         
       mutable float                 m_hadIntMotherPhi;            
       mutable float                 m_hadIntMotherEta;            
 
       mutable int                   m_hadIntChildren;                          
       mutable float                 m_hadIntChildE;                            
       mutable float                 m_hadIntChildP[MAXHADINTCHILDREN]{};     
       mutable float                 m_hadIntChildPcms[MAXHADINTCHILDREN]{};     
       mutable int                   m_hadIntChildPdg[MAXHADINTCHILDREN]{};        
       mutable float                 m_hadIntChildPhi[MAXHADINTCHILDREN]{};        
       mutable float                 m_hadIntChildEta[MAXHADINTCHILDREN]{};        
       mutable float                 m_hadIntChildTh[MAXHADINTCHILDREN]{};        
       mutable float                 m_hadIntChildThc[MAXHADINTCHILDREN]{};        
       mutable float                 m_hadIntChildDeltaPhi[MAXHADINTCHILDREN]{};   
       mutable float                 m_hadIntChildDeltaEta[MAXHADINTCHILDREN]{};   
 
   };
}
 
#endif // ISF_FATRASTOOLS_HadIntProcessorParametric_H