ALFA_MDTracking.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ALFA_MDTRACKING_H
#define ALFA_MDTRACKING_H
 
#include <iostream>
#include <vector>
 
#include "AthenaBaseComps/AthMessaging.h"
 
#include "GaudiKernel/StatusCode.h"
 
#include "Rtypes.h"
#include "TMath.h"
#include "TH1.h"
#include "TH2.h"
 
#include "ALFA_LocRec/ALFA_UserObjects.h"
#include "ALFA_Geometry/ALFA_constants.h"
 
class ALFA_MDTracking : public AthMessaging
{
    public:
        ALFA_MDTracking();
        ~ALFA_MDTracking();
 
    private:
        Int_t m_iMultiplicityCut;
        Int_t m_iUVCut;
        Int_t m_iRPot;
        Int_t m_iNumU, m_iNumV;
        Float_t m_fOverlapCut;
        Float_t m_fRecXPos;
        Float_t m_fRecYPos;
        Float_t m_faMD[RPOTSCNT][ALFALAYERSCNT*ALFAPLATESCNT][ALFAFIBERSCNT];   //slope, offset and Z-pos for MD fibers [2*10][64]
        Float_t m_fbMD[RPOTSCNT][ALFALAYERSCNT*ALFAPLATESCNT][ALFAFIBERSCNT];
        double  m_fOverlapU, m_fOverlapV;
        Int_t m_iFibSel[ALFALAYERSCNT*ALFAPLATESCNT];
 
    private:
        std::map<int, FIBERS> m_MapLayers;
 
    public:
        StatusCode Initialize(Float_t faMD[RPOTSCNT][ALFALAYERSCNT*ALFAPLATESCNT][ALFAFIBERSCNT], Float_t fbMD[RPOTSCNT][ALFALAYERSCNT*ALFAPLATESCNT][ALFAFIBERSCNT], Int_t iMultiplicityCut, Int_t iUVCut, Float_t fOverlapCut);
        StatusCode Execute(Int_t iRPot, const std::list<MDHIT> &ListMDHits);
        StatusCode Finalize(Float_t &fRecXPos, Float_t &fRecYPos);
 
        void GetData(Int_t &NumU, Int_t &NumV, Float_t &OL_U, Float_t &OL_V, Int_t (&iFibSel)[ALFALAYERSCNT*ALFAPLATESCNT]);
 
    private:
        void HistFill(Float_t &b_p, Float_t &b_n, Float_t &Ov_p, Float_t &Ov_n, Int_t &NumU, Int_t &NumV, Int_t iFlag);
        void OverLap();
        void SetData(Int_t NumU, Int_t NumV, Float_t OL_U, Float_t OL_V);
 
};
#endif // ALFA_MDTRACKING_H