RiddersAlgorithm.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// RiddersAlgorithm.h, (c) ATLAS Detector software
 
#ifndef TRKRIDDERSALGORITHM_H
#define TRKRIDDERSALGORITHM_H
 
// Gaudi includes
#include "AthenaBaseComps/AthAlgorithm.h"
#include "GaudiKernel/ToolHandle.h"
#include "GaudiKernel/IRndmGenSvc.h"
#include "GaudiKernel/RndmGenerators.h"
#include "GeoPrimitives/GeoPrimitives.h"
#include <string>
 
class AtlasDetectorID;
class Identifier;
class TTree;
 
#ifndef RIDDLERSSTEPS
#define RIDDLERSSTEPS 10
#endif
 
namespace Trk 
{
 
  class IPropagator;
  class MagneticFieldProperties;
 
  class RiddersAlgorithm : public AthAlgorithm
    {
    public:
 
       RiddersAlgorithm(const std::string& name, ISvcLocator* pSvcLocator);
       ~RiddersAlgorithm();
 
       StatusCode          initialize();
       StatusCode          execute();
       StatusCode          finalize();
 
    private:
      static Amg::Transform3D createTransform(double x, double y, double z, double phi=0., double theta=0., double alphaZ=0.);      
     
      static double parabolicInterpolation(double y0, double y1, double y2,
                                           double x0, double x1, double x2);
 
      ToolHandle<IPropagator>           m_propagator;       
      bool                              m_useCustomField;
      bool                              m_useAlignedSurfaces;
      double                            m_fieldValue;
      MagneticFieldProperties*          m_magFieldProperties;
 
      double                    m_sigmaLoc;       
      double                    m_sigmaR;         
      double                    m_minPhi;         
      double                    m_maxPhi;         
      double                    m_minEta;         
      double                    m_maxEta;         
      double                    m_minP;           
      double                    m_maxP;           
 
      double                            m_minimumR;
      double                            m_maximumR;
 
      std::vector<double>               m_localVariations;
      std::vector<double>               m_angularVariations;
      std::vector<double>               m_qOpVariations;
 
      TTree*                    m_validationTree;            
 
      std::string               m_validationTreeName;        
      std::string               m_validationTreeDescription; 
      std::string               m_validationTreeFolder;      
 
      int               m_steps;
      float             m_loc1loc1[RIDDLERSSTEPS];
      float             m_loc1loc2[RIDDLERSSTEPS];
      float             m_loc1phi[RIDDLERSSTEPS];
      float             m_loc1theta[RIDDLERSSTEPS];
      float             m_loc1qop[RIDDLERSSTEPS];
      float             m_loc1steps[RIDDLERSSTEPS];
 
      float             m_loc2loc1[RIDDLERSSTEPS];
      float             m_loc2loc2[RIDDLERSSTEPS];
      float             m_loc2phi[RIDDLERSSTEPS];
      float             m_loc2theta[RIDDLERSSTEPS];
      float             m_loc2qop[RIDDLERSSTEPS];
      float             m_loc2steps[RIDDLERSSTEPS];
 
      float             m_philoc1[RIDDLERSSTEPS];
      float             m_philoc2[RIDDLERSSTEPS];
      float             m_phiphi[RIDDLERSSTEPS];
      float             m_phitheta[RIDDLERSSTEPS];
      float             m_phiqop[RIDDLERSSTEPS];
      float             m_phisteps[RIDDLERSSTEPS];
 
      float             m_thetaloc1[RIDDLERSSTEPS];
      float             m_thetaloc2[RIDDLERSSTEPS];
      float             m_thetaphi[RIDDLERSSTEPS];
      float             m_thetatheta[RIDDLERSSTEPS];
      float             m_thetaqop[RIDDLERSSTEPS];
      float             m_thetasteps[RIDDLERSSTEPS];
 
      float             m_qoploc1[RIDDLERSSTEPS];
      float             m_qoploc2[RIDDLERSSTEPS];
      float             m_qopphi[RIDDLERSSTEPS];
      float             m_qoptheta[RIDDLERSSTEPS];
      float             m_qopqop[RIDDLERSSTEPS];
      float             m_qopsteps[RIDDLERSSTEPS];
 
 
      Rndm::Numbers*                    m_gaussDist;
      Rndm::Numbers*                    m_flatDist;
 
    }; 
} // end of namespace
 
#endif