AFP_ProtonRecoAnalytical.h

Go to the documentation of this file.

/*
Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef AFP_GLOBRECO_AFP_PROTONRECOANALYTICAL_H
#define AFP_GLOBRECO_AFP_PROTONRECOANALYTICAL_H
 
#include "AFP_GlobReco/AFP_ProtonRecoBase.h"
#include "AFP_GlobReco/Parameterization.h"
 
// FrameWork includes
#include "Gaudi/Property.h"
#include "PathResolver/PathResolver.h"
 
// xAOD includes
#include "xAODForward/AFPProton.h"
#include "xAODForward/AFPProtonContainer.h"
#include "xAODForward/AFPProtonAuxContainer.h"
#include "xAODForward/AFPTrack.h"
#include "xAODForward/AFPTrackContainer.h"
#include "xAODForward/AFPTrackAuxContainer.h"
#include "xAODForward/AFPProtonRecoAlgID.h"
 
#include <string>
#include <vector>
#include <memory>
#include <algorithm>
 

class AFP_ProtonRecoAnalytical : public AFP_ProtonRecoBase {
 
  public:
    AFP_ProtonRecoAnalytical(const std::string& type, const std::string& name, const IInterface* parent);
    
    ~AFP_ProtonRecoAnalytical() = default;
    
    StatusCode initialize() override;
 
  private:
    
    StatusCode configInfo() const;
    
    virtual
    xAOD::AFPProton * reco(const xAOD::AFPTrack* trkNear, const xAOD::AFPTrack* trkFar, std::unique_ptr<xAOD::AFPProtonContainer>& outputContainer) const override;

    virtual
    xAOD::AFPProton * reco (const xAOD::AFPTrack* trkFar, std::unique_ptr<xAOD::AFPProtonContainer>& outputContainer) const override;

    virtual
    double chi2(double energy, double sx, double sy, const Measurement& my_measAFP) const override;

    double bisection(double (AFP_ProtonRecoAnalytical::*fun)(double,const Measurement&,std::vector<double>&, std::vector<double>&) const, const Measurement& my_measAFP, std::vector<double>& my_slopeCalculated, std::vector<double>& my_positionCalculated) const;

    double bothStations(double energy,const  Measurement& my_measAFP, std::vector<double>& my_slopeCalculated, std::vector<double>& my_positionCalculated) const;

    double singleStation(double energy, const Measurement& my_measAFP, std::vector<double>& my_slopeCalculated, std::vector<double>& my_positionCalculated) const;

    double calculateSlope(double energy, int XorY, std::vector<double>& my_slopeCalculated) const;

    double calculateXslope(double energy, std::vector<double>& my_slopeCalculated) const;

    double calculateYslope(double energy, std::vector<double>& my_slopeCalculated) const;

    std::unique_ptr<AFP::Parameterization> m_parametrization;

    double m_distanceBetweenStations = 0.0; // [m]

    Gaudi::Property<std::string> m_parametrizationFileName{this, "parametrizationFileName", {}, "Name of the file containing parameterization"};

    double m_parametrizationPosition = 0.0; // [m]

    double m_parametrizationEnergy = 0.0; // [GeV]
};
 
 
#endif