Loading [MathJax]/extensions/tex2jax.js
 |
ATLAS Offline Software
|
Go to the documentation of this file.
11 #ifndef JETCALIBTOOLS_JESCALIBSTEP_H
12 #define JETCALIBTOOLS_JESCALIBSTEP_H 1
45 double getLogPolN(
const double *factors,
double x)
const;
53 double getJES(
const double X,
const double Y=0,
const double Emax=-1)
const;
59 double getLowPtJES(
double E_uncorr,
double eta_det)
const ;
65 Gaudi::Property< std::string >
m_constantFileName {
this,
"CalibConstantFile",
"/afs/cern.ch/work/s/stapiaar/JetDev4/athena/JetToolHelpers/data/file_JES.config",
"text file containing constants" };
67 Gaudi::Property< std::string >
m_jetAlgo {
this,
"JetAlgo",
"AntiKt4EMPFlow",
"jet collection" };
68 Gaudi::Property< float >
m_minPt_JES = {
this,
"MinPtForETAJES",15,
"min pT"};
69 Gaudi::Property< bool >
m_freezeJESatHighE = {
this,
"FreezeJEScorrectionatHighE",
false,
" freeze at high e"};
70 Gaudi::Property< float >
m_lowPtExtrap = {
this,
"LowPtJESExtrapolationMethod", 0,
" low pt etrap"};
71 Gaudi::Property< float >
m_lowPtMinR = {
this,
"LowPtJESExtrapolationMinimumResponse", 0.25,
" low pt etrap min"};
74 Gaudi::Property< float >
m_maxE_EtaCorr = {
this,
"MaxEForEtaCorr" ,2500. ,
""};
76 Gaudi::Property< std::string >
m_histoFileName {
this,
"HistoFile",
"none",
"root file containing histos for spline calib" };
77 Gaudi::Property< bool >
m_useSpline = {
this,
"UseSpline",
false,
" use spline"};
80 ToolHandle<JetHelper::IVarTool>
m_vartoolE{
this,
"VarToolE",
"VarTool",
"InputVariable instance E (or pT?)" };
82 ToolHandle<JetHelper::IVarTool>
m_vartoolEta{
this,
"VarToolEta",
"VarTool",
"InputVariable instance eta (or rapididty?)" };
virtual StatusCode initialize() override
Dummy implementation of the initialisation function.
int getEtaBin(double eta_det) const
double m_JES_MinPt_R[s_nEtaBins]
double m_etaCorrFactors[s_nEtaBins][s_nParMax]
EtaMassJESCalibStep(const std::string &name="EtaMassJESCalibStep")
double m_energyFreezeJES[s_nEtaBins]
double getLowPtJES(double E_uncorr, double eta_det) const
deal with low pt jets
double getEtaCorr(double X, double Y=0) const
return Eta correction
static const unsigned int s_nEtaBins
Gaudi::Property< std::string > m_jetAlgo
jet collection to be calibrated
Gaudi::Property< float > m_lowPtExtrap
double m_JESFactors[s_nEtaBins][s_nParMax]
void loadSplineHists(const std::string &fileName, const std::string &etajes_name)
virtual StatusCode calibrate(xAOD::JetContainer &) const override
Apply calibration to a jet container.
double getSplineCorr(const int etaBin, double E) const
double m_JES_MinPt_E[s_nEtaBins]
setSAddress setEtaMS setDirPhiMS setDirZMS setBarrelRadius setEndcapAlpha setEndcapRadius setInterceptInner setEtaMap etaBin
double getJES(const double X, const double Y=0, const double Emax=-1) const
return MCJES calibration factor
double getLogPolNSlope(const double *factors, double x) const
::StatusCode StatusCode
StatusCode definition for legacy code.
double m_JES_MinPt_Slopes[s_nEtaBins]
Gaudi::Property< float > m_minPt_JES
ToolHandle< JetHelper::IVarTool > m_vartoolEta
interface for xAOD::jet variable to be defined by user, this must correspond to jet Eta in currect ve...
Gaudi::Property< float > m_maxE_EtaCorr
bool readMCJESFromHists()
ToolHandle< JetHelper::IVarTool > m_vartoolE
double getLogPolN(const double *factors, double x) const
Gaudi::Property< float > m_minPt_EtaCorr
Gaudi::Property< std::string > m_histoFileName
static const unsigned int s_nParMax
Gaudi::Property< float > m_lowPtMinR
Gaudi::Property< bool > m_useSpline
double getEmaxJES(const double Y) const
return Emax
double getSplineSlope(const int ieta, const double minE) const
Gaudi::Property< std::string > m_constantFileName
name of the text file
std::vector< std::unique_ptr< TH1 > > m_etajesFactors
Gaudi::Property< bool > m_freezeJESatHighE