TBPlaneTrackingAlgo.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TBREC_TBPLANETRACKINGALGO_H
#define TBREC_TBPLANETRACKINGALGO_H
 
#include "AthenaBaseComps/AthAlgorithm.h"
 
#include "TBEvent/TBHitPlaneCont.h"
 
#include <utility>
#include <vector>

 
class TBPlaneTrackingAlgo : public AthAlgorithm {
 public:
  
  // Contructor & Destructor //
  TBPlaneTrackingAlgo(const std::string& name, ISvcLocator* pSvcLocator);
  virtual ~TBPlaneTrackingAlgo(){};
 
  // public methods //
  virtual StatusCode initialize() override;
  virtual StatusCode execute() override;
  virtual StatusCode finalize() override;
 
 private:
  void FillRandomHit();
  StatusCode getnewcalib();
 
  
  // Private function members //

  bool fitPlane(const TBHitPlaneCont * hitPlaneCont, double &a1, double &a2,
        double &chi2, std::vector<double> &residual);

  bool fitHits(const std::vector<double> & u, const std::vector<double> & w,
           const std::vector<double> & eu, // const std::vector<double> & ew,
           double &a1, double &a2);

  double getResidual(const double & u, const double & w,const double &a1,
             const double &a2);

  double getChi2(const std::vector<double> &v_u, 
         const std::vector<double> &v_w,
         const std::vector<double> &v_eu,
//       const std::vector<double> &v_ew,
         const double &a1, const double &a2);

  StatusCode buildHits();
    
  //Used to test this algo set to true if you want to debug
  bool m_testAlgo;
  TBHitPlaneCont m_hitPlaneCont_u;
  TBHitPlaneCont m_hitPlaneCont_v;
 
  // Positions of BPC needed to construct hits  
  std::vector<float> m_bpc_posZX,m_bpc_posZY,m_bpc_posX,m_bpc_posY;
  std::vector<float> m_bpc_errposZX,m_bpc_errposZY,m_bpc_errposX,m_bpc_errposY;
  std::vector<float> m_bpc_errmeasX,m_bpc_errmeasY;
 
  std::vector<std::string> m_bpc_names;
 
  // filename containing positions constants
  std::string m_calib_filename;
  unsigned int m_runnumber;
 
};
 
#endif