d3/da4/ExtrapolationValidation_8h_source.html

/*

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

*/


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


#ifndef TRKEXALGS_EXTRAPOLATIONVALIDATION_H

#define TRKEXALGS_EXTRAPOLATIONVALIDATION_H


// Gaudi includes

#include "AthenaBaseComps/AthAlgorithm.h"

#include "GaudiKernel/MsgStream.h"

#include "GaudiKernel/ToolHandle.h"

#include "GaudiKernel/IRndmGenSvc.h"

#include "GaudiKernel/RndmGenerators.h"

#include "EventPrimitives/EventPrimitives.h"

#include "GeoPrimitives/GeoPrimitivesHelpers.h"

#include <string>


class TTree;


#ifndef TRKEXALGS_MAXPARAMETERS

#define TRKEXALGS_MAXPARAMETERS 10

#endif


namespace Trk

{


  class IExtrapolator;

  class Surface;

  class TrackingGeometry;

  class TrackingVolume;


  class ExtrapolationValidation : public AthAlgorithm

    {

    public:


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

       ~ExtrapolationValidation();


       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.);


      const TrackingVolume*     m_highestVolume;


      ToolHandle<IExtrapolator> m_extrapolator;


      Rndm::Numbers*            m_gaussDist;

      Rndm::Numbers*            m_flatDist;


      bool                      m_materialCollectionValidation;


      bool                      m_direct; //extrapolate directly


      TTree*                    m_validationTree;


      std::string               m_validationTreeName;

      std::string               m_validationTreeDescription;

      std::string               m_validationTreeFolder;


      double                    m_maximumR;

      double                    m_maximumZ;


      double                    m_sigmaLoc;

      double                    m_sigmaR;

      double                    m_sigmaZ;

      double                    m_minEta;

      double                    m_maxEta;

      double                    m_minP;

      double                    m_maxP;


      int                       m_particleType;


      int               m_parameters;

      float             m_parameterLoc1[TRKEXALGS_MAXPARAMETERS];

      float             m_parameterLoc2[TRKEXALGS_MAXPARAMETERS];

      float             m_parameterPhi[TRKEXALGS_MAXPARAMETERS];

      float             m_parameterTheta[TRKEXALGS_MAXPARAMETERS];

      float             m_parameterEta[TRKEXALGS_MAXPARAMETERS];

      float             m_parameterQoverP[TRKEXALGS_MAXPARAMETERS];


      float             m_covarianceLoc1[TRKEXALGS_MAXPARAMETERS];

      float             m_covarianceLoc2[TRKEXALGS_MAXPARAMETERS];

      float             m_covariancePhi[TRKEXALGS_MAXPARAMETERS];

      float             m_covarianceTheta[TRKEXALGS_MAXPARAMETERS];

      float             m_covarianceQoverP[TRKEXALGS_MAXPARAMETERS];

      float             m_covarianceDeterminant[TRKEXALGS_MAXPARAMETERS];


      int               m_destinationSurfaceType;

      float             m_startX;

      float             m_startY;

      float             m_startR;

      float             m_startZ;

      float             m_startP;


      float             m_estimationX;

      float             m_estimationY;

      float             m_estimationR;

      float             m_estimationZ;


      float             m_destinationX;

      float             m_destinationY;

      float             m_destinationR;

      float             m_destinationZ;


      // ---- output statistics

      unsigned int     m_triesFront;

      unsigned int     m_breaksFront;

      unsigned int     m_triesBack;

      unsigned int     m_breaksBack;


      unsigned int     m_collectedLayerFront;

      unsigned int     m_collectedLayerBack;


      //AmgSymMatrix(5) * m_covariance;


    };

} // end of namespace


#endif