PhotonConversionTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_FATRASTOOLS_PHOTONCONVERSIONTOOL_H
#define ISF_FATRASTOOLS_PHOTONCONVERSIONTOOL_H
 
// Gaudi & Athena
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ServiceHandle.h"
#include "GaudiKernel/ToolHandle.h"
#include "AthenaKernel/IAtRndmGenSvc.h"
#include "CxxUtils/checker_macros.h"
// Trk
#include "TrkEventPrimitives/PropDirection.h"
#include "TrkExUtils/MaterialUpdateMode.h"
#include "TrkDetDescrUtils/GeometrySignature.h"
// Fatras
#include "ISF_FatrasInterfaces/IPhotonConversionTool.h"
 
class TTree;
 
class StoreGateSvc;
 
namespace Trk {
  class Layer;
  class CylinderVolumeBounds;
  class PdgToParticleHypothesis;
  class TrackingGeometry;
}
 
namespace ISF {
  class IParticleBroker;
  class ITruthSvc;
}
 
namespace iFatras {
 
  class IPhysicsValidationTool;
  
  class ATLAS_NOT_THREAD_SAFE PhotonConversionTool : public extends<AthAlgTool, iFatras::IPhotonConversionTool> {  // deprecated: ATLASSIM-6020
  public:      
      PhotonConversionTool(const std::string&,const std::string&,const IInterface*);
      virtual ~PhotonConversionTool();
       
      StatusCode initialize();
      StatusCode finalize();
 
    
      bool pairProduction(const Trk::MaterialProperties& mprop,
              double pathCorrection,
              double p) const;
      
      bool doConversion(double time, const Trk::NeutralParameters& parm,
            const Trk::ExtendedMaterialProperties* extMatProp=0) const;
      
      ISF::ISFParticleVector doConversionOnLayer(const ISF::ISFParticle* parent, 
                         double time, const Trk::NeutralParameters& parm,
                         const Trk::ExtendedMaterialProperties *ematprop=0) const;
 
 
    
   private:
      void recordChilds(double time, 
            const Amg::Vector3D& vertex,
            const Amg::Vector3D& photonDirection,
            double childEnergy, double photonEnergy,
            const Amg::Vector3D& childDirection,
            Trk::ParticleHypothesis childType=Trk::electron) const;
      
      ISF::ISFParticleVector getChilds(const ISF::ISFParticle* parent,
                       double time,
                       const Amg::Vector3D& vertex,
                       const Amg::Vector3D& photonDirection,
                       double childEnergy, double photonEnergy,
                       const Amg::Vector3D& childDirection,
                       Trk::ParticleHypothesis childType) const;
    
      double childEnergyFraction(const Trk::MaterialProperties& mprop, double gammaMom) const;
      double childEnergyFraction(double gammaMom) const;
      
      Amg::Vector3D childDirection(const Amg::Vector3D& gammaMom,
                      double childE) const;
      
      double phi1(double delta) const;
      
      double phi2(double delta) const;
 
      ServiceHandle<ISF::IParticleBroker>          m_particleBroker;
      ServiceHandle<ISF::ITruthSvc>                m_truthRecordSvc;
      
      int                  m_processCode;
      
      bool                                         m_referenceMaterial;
      
      double                                       m_minChildEnergy;
      double                                       m_childEnergyScaleFactor;
      double                                       m_conversionProbScaleFactor;
      
      
      ServiceHandle<IAtRndmGenSvc>                 m_rndGenSvc;
      CLHEP::HepRandomEngine*                             m_randomEngine;
      std::string                                  m_randomEngineName;         
      
      mutable int                                  m_recordedConversions;
      mutable int                                  m_droppedConversions;
      
      // ------------------------ Validation section ------------------------------------
      bool                                         m_validationMode;
      std::string                                  m_validationTreeName;        
      std::string                                  m_validationTreeDescription; 
      std::string                                  m_validationTreeFolder;      
      
      TTree*                                       m_validationTree;            
      ToolHandle<IPhysicsValidationTool>           m_validationTool;
      
      // ------ ntuple branches
      mutable float                                m_conversionPointX;
      mutable float                                m_conversionPointY;
      mutable float                                m_conversionPointR;
      mutable float                                m_conversionPointZ;
      mutable float                                m_conversionPhotonEnergy;
      mutable float                                m_conversionChildEnergy;
      mutable float                                m_conversionChildAngle;
      // --------------------------------------------------------------------------------
      
      static const Trk::PdgToParticleHypothesis          s_pdgToHypo;
   };
  
 
inline double PhotonConversionTool::phi1(double delta) const {
  if (delta <= 1.)
     return 20.867 - 3.242 * delta  + 0.625*delta*delta;
  else
    return 21.12 - 4.184*log(delta+0.952);
}
 
inline double PhotonConversionTool::phi2(double delta) const {
  if (delta <= 1.)
     return 20.209 - 1.930 * delta  + 0.086*delta*delta;
   return 21.12 - 4.184*log(delta+0.952);
}
 
}
 
#endif