da/d43/CETmaterial_8h_source.html

/*

  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration

*/


// CETmaterial.h, (c) ATLAS Detector software


#ifndef TRKEXALGS_CETMATERIAL_H

#define TRKEXALGS_CETMATERIAL_H


// Gaudi includes

#include "AthenaBaseComps/AthAlgorithm.h"

#include "GaudiKernel/IRndmGenSvc.h"

#include "GaudiKernel/RndmGenerators.h"

#include "GaudiKernel/ToolHandle.h"

#include "GaudiKernel/IChronoStatSvc.h"

#include "TrkParameters/TrackParameters.h"

#include <string>


namespace Trk

{


  class IExtrapolator;

  class Surface;

  class TrackingVolume;

  class TrackingGeometry;


  class CETmaterial : public AthAlgorithm

    {

    public:


       CETmaterial(const std::string& name, ISvcLocator* pSvcLocator);

       ~CETmaterial();


       StatusCode          initialize();

       StatusCode          execute();

       StatusCode          finalize();


    private:

       void printMat(double th, double ph, double mat,double dtheta=0.,double dphi=0.) const;

       void printMatScan(double theta, double phi, double r, double z, double mat, const std::string& name) const;

       void printMatPrec(double theta, double phi, const Trk::TrackParameters*, const Trk::TrackParameters*,

             double mat, int id, const std::string& name);

       void printMatComp(double theta, double phi, const Trk::TrackParameters* currPar, const std::string& name, double mat, double matApp,double dx, double dy) const;


      ToolHandle<IExtrapolator>         m_extrapolator;

      ToolHandle<IExtrapolator>         m_extraprec;


      Rndm::Numbers*                    m_gaussDist;

      Rndm::Numbers*                    m_flatDist;


      double                   m_sigmaD0;

      double                   m_minZ0;

      double                   m_maxZ0;

      double                   m_minEta;

      double                   m_maxEta;

      double                   m_minTheta;

      double                   m_maxTheta;

      double                   m_minPhi;

      double                   m_maxPhi;

      double                   m_minP;

      double                   m_maxP;

      double                   m_charge;

      unsigned int             m_numScan;

      bool                     m_checkStepWise;

      bool                     m_printMaterial;

      bool                     m_printActive;

      const char*              m_matTotFile;

      const char*              m_matScanFile;

      const char*              m_matActiveFile;

      const char*              m_matCompFile;

      bool                     m_backward;

      bool                     m_domsentry;

      bool                     m_doprecision;

      double                   m_th;

      double                   m_ph;

      int                      m_id;

      double                   m_matSaved;

      Trk::TrackParameters*    m_next;

      Amg::MatrixX*            m_err;


      const Trk::Surface*                m_outerBoundary;

      const Trk::TrackingGeometry*       m_trackingGeometry;

      const Trk::TrackingVolume*         m_msentry;


      int                      m_particleType;


      typedef ServiceHandle<IChronoStatSvc> IChronoStatSvc_t;

      IChronoStatSvc_t              m_chronoStatSvc;


    };

} // end of namespace


#endif