ATLAS Offline Software
Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
EtaJESCorrection Class Reference

#include <EtaJESCorrection.h>

Inheritance diagram for EtaJESCorrection:
Collaboration diagram for EtaJESCorrection:

Public Member Functions

 EtaJESCorrection ()
 
 EtaJESCorrection (const std::string &name, TEnv *config, TString jetAlgo, TString calibAreaTag, bool mass, bool dev)
 
virtual ~EtaJESCorrection ()
 
virtual StatusCode initialize () override
 
virtual StatusCode calibrate (xAOD::Jet &jet, JetEventInfo &) const override
 
virtual void setUnitsGeV (bool useGeV)
 
virtual StatusCode getNominalResolutionData (const xAOD::Jet &, double &) const
 
virtual StatusCode getNominalResolutionMC (const xAOD::Jet &, double &) const
 
void setLevel (MSG::Level lvl)
 Change the current logging level. More...
 

Protected Member Functions

virtual StatusCode setStartP4 (xAOD::Jet &jet) const
 

Protected Attributes

double m_GeV
 
std::string m_jetStartScale
 
std::string m_name
 

Private Member Functions

double getJES (double E_uncorr, double eta_det) const
 
double getLowPtJES (double E_uncorr, double eta_det) const
 
double getEtaCorr (double E_corr, double eta_det) const
 
double getMassCorr (double E_corr, double eta_det) const
 
double getLogPolN (const double *factors, double x) const
 
double getLogPolNSlope (const double *factors, double x) const
 
double getSplineCorr (const int etaBin, double E) const
 
double getSplineSlope (const int ieta, const double minE) const
 
void loadSplineHists (const TString &fileName, const std::string &etajes_name="etaJes")
 Loads the calib constants from histograms in TFile named fileName. More...
 
int getEtaBin (double eta_det) const
 
void initMessaging () const
 Initialize our message level and MessageSvc. More...
 

Private Attributes

TEnv * m_config {}
 
TString m_jetAlgo
 
TString m_calibAreaTag
 
bool m_mass {}
 
bool m_dev {}
 
bool m_freezeJESatHighE {}
 
bool m_isSpline {}
 
TString m_jesDesc
 
double m_minPt_JES {}
 
double m_minPt_EtaCorr {}
 
double m_maxE_EtaCorr {}
 
unsigned int m_lowPtExtrap {}
 
double m_lowPtMinR {}
 
bool m_applyMassCorrection {}
 
bool m_useSecondaryminPt_JES {}
 
double m_etaSecondaryminPt_JES {}
 
double m_secondaryminPt_JES {}
 
TAxis * m_etaBinAxis {}
 
unsigned int m_nPar {}
 
double m_JESFactors [s_nEtaBins][s_nParMax] ={}
 
double m_JES_MinPt_Slopes [s_nEtaBins] ={}
 
double m_JES_MinPt_E [s_nEtaBins] ={}
 
double m_JES_MinPt_R [s_nEtaBins] ={}
 
double m_JES_MinPt_Param1 [s_nEtaBins] ={}
 
double m_JES_MinPt_Param2 [s_nEtaBins]
 
double m_etaCorrFactors [s_nEtaBins][s_nParMax] ={}
 
double m_JMSFactors [s_nEtaBins][s_nParMax] ={}
 
double m_energyFreezeJES [s_nEtaBins] ={}
 
std::vector< std::unique_ptr< TH1 > > m_etajesFactors
 
std::string m_nm
 Message source name. More...
 
boost::thread_specific_ptr< MsgStream > m_msg_tls
 MsgStream instance (a std::cout like with print-out levels) More...
 
std::atomic< IMessageSvc * > m_imsg { nullptr }
 MessageSvc pointer. More...
 
std::atomic< MSG::Level > m_lvl { MSG::NIL }
 Current logging level. More...
 
std::atomic_flag m_initialized ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
 Messaging initialized (initMessaging) More...
 

Static Private Attributes

static const unsigned int s_nEtaBins =90
 
static const unsigned int s_nParMin =7
 
static const unsigned int s_nParMax =9
 

Detailed Description

Definition at line 26 of file EtaJESCorrection.h.

Constructor & Destructor Documentation

◆ EtaJESCorrection() [1/2]

EtaJESCorrection::EtaJESCorrection ( )

Definition at line 12 of file EtaJESCorrection.cxx.

13  : JetCalibrationStep::JetCalibrationStep(),
14  m_config(nullptr), m_jetAlgo(""), m_calibAreaTag(""), m_mass(false), m_dev(false),
15  m_minPt_JES(10), m_minPt_EtaCorr(8), m_maxE_EtaCorr(2500),
16  m_lowPtExtrap(0), m_lowPtMinR(0.25),
17  m_etaBinAxis(nullptr)
18 { }

◆ EtaJESCorrection() [2/2]

EtaJESCorrection::EtaJESCorrection ( const std::string &  name,
TEnv *  config,
TString  jetAlgo,
TString  calibAreaTag,
bool  mass,
bool  dev 
)

Definition at line 20 of file EtaJESCorrection.cxx.

21  : JetCalibrationStep::JetCalibrationStep(name.c_str()),
22  m_config(config), m_jetAlgo(std::move(jetAlgo)), m_calibAreaTag(std::move(calibAreaTag)), m_mass(mass), m_dev(dev),
23  m_minPt_JES(10), m_minPt_EtaCorr(8), m_maxE_EtaCorr(2500),
24  m_lowPtExtrap(0), m_lowPtMinR(0.25),
25  m_etaBinAxis(nullptr)
26 { }

◆ ~EtaJESCorrection()

EtaJESCorrection::~EtaJESCorrection ( )
virtual

Definition at line 28 of file EtaJESCorrection.cxx.

28  {
29 
30  if (m_etaBinAxis) delete m_etaBinAxis;
31 
32 }

Member Function Documentation

◆ calibrate()

StatusCode EtaJESCorrection::calibrate ( xAOD::Jet &  jet,
JetEventInfo &   
) const
overridevirtual

Implements JetCalibrationStep.

Definition at line 224 of file EtaJESCorrection.cxx.

224  {
225 
226  xAOD::JetFourMom_t jetStartP4;
227  ATH_CHECK( setStartP4(jet) );
228  jetStartP4 = jet.jetP4();
229 
230  //Apply the JES calibration scale factor
231  //Takes the uncorrected jet eta (in case the origin and/or 4vector jet area corrections were applied
232  float detectorEta = jet.getAttribute<float>("DetectorEta");
233 
234  xAOD::JetFourMom_t calibP4 = jetStartP4*getJES( jetStartP4.e(), detectorEta );
235 
236  const double etaCorr = calibP4.eta() + getEtaCorr( calibP4.e(), detectorEta );
237  double massCorr;
238  if(!m_applyMassCorrection) massCorr = calibP4.mass();
239  else{ massCorr = jetStartP4.mass()*getMassCorr(calibP4.e(), detectorEta); }
240  TLorentzVector TLVjet;
241  TLVjet.SetPtEtaPhiM( calibP4.P()/cosh(etaCorr),etaCorr,calibP4.phi(),massCorr );
242  calibP4.SetPxPyPzE( TLVjet.Px(), TLVjet.Py(), TLVjet.Pz(), TLVjet.E() );
243 
244  if(m_dev){
245  float JESFactor = calibP4.e()/jetStartP4.e();
246  jet.setAttribute<float>("JetJESCalibFactor",JESFactor);
247  }
248 
249  //Transfer calibrated jet properties to the Jet object
250  jet.setAttribute<xAOD::JetFourMom_t>("JetEtaJESScaleMomentum",calibP4);
251  jet.setJetP4( calibP4 );
252 
253  return StatusCode::SUCCESS;
254 }

◆ getEtaBin()

int EtaJESCorrection::getEtaBin ( double  eta_det) const
private

Definition at line 385 of file EtaJESCorrection.cxx.

385  {
386  int bin = std::as_const(m_etaBinAxis)->FindBin(eta_det);
387  if (bin<=0) return 0;
388  if (bin>m_etaBinAxis->GetNbins()) return bin-2; // overflow
389  return bin-1;
390 }

◆ getEtaCorr()

double EtaJESCorrection::getEtaCorr ( double  E_corr,
double  eta_det 
) const
private

Definition at line 342 of file EtaJESCorrection.cxx.

342  {
343  int ieta = getEtaBin(eta_det);
344  const double *factors = m_etaCorrFactors[ieta];
345 
346  double E = E_corr/m_GeV;
347  if ( E < m_minPt_EtaCorr*cosh(eta_det) )
348  E = m_minPt_EtaCorr*cosh(eta_det);
349  if ( E>m_maxE_EtaCorr ) E=m_maxE_EtaCorr;
350 
351  double etaCorr = getLogPolN(factors,E);
352 
353  // This is ( reco_eta - truth_eta )
354  // to make it an additive correction return the negative value
355  return -etaCorr;
356 }

◆ getJES()

double EtaJESCorrection::getJES ( double  E_uncorr,
double  eta_det 
) const
private

Definition at line 258 of file EtaJESCorrection.cxx.

258  {
259 
260  double E = E_uncorr/m_GeV; // E in GeV
261  //Check if the Pt goes below the minimum value, if so use the special GetLowPtJES method
262  if(m_useSecondaryminPt_JES){
263  if(fabs(eta_det) < m_etaSecondaryminPt_JES && E/cosh(eta_det) < m_minPt_JES){
264  double R = getLowPtJES(E,eta_det);
265  return 1.0/R;
266  }
267  if(fabs(eta_det) >= m_etaSecondaryminPt_JES && E/cosh(eta_det) < m_secondaryminPt_JES){
268  double R = getLowPtJES(E,eta_det);
269  return 1.0/R;
270  }
271  }else{
272  if ( E/cosh(eta_det) < m_minPt_JES ) {
273  double R = getLowPtJES(E,eta_det);
274  return 1.0/R;
275  }
276  }
277 
278  // Get the factors
279  int ieta = getEtaBin(eta_det);
280 
281  // Freeze correction
282  if(m_freezeJESatHighE && E>m_energyFreezeJES[ieta] && m_energyFreezeJES[ieta]!=-1) E = m_energyFreezeJES[ieta];
283 
284 
285  // The low pT extrapolation doesn't work for the spline yet, so putting this code here.
286  if(m_isSpline){
287  double R = getSplineCorr(ieta, E);
288  return 1.0/R;
289  }
290  const double *factors = m_JESFactors[ieta];
291 
292 
293  // Calculate the jet response and then the JES as 1/R
294  double R = getLogPolN(factors,E);
295  return 1.0/R;
296 }

◆ getLogPolN()

double EtaJESCorrection::getLogPolN ( const double *  factors,
double  x 
) const
private

Definition at line 370 of file EtaJESCorrection.cxx.

370  {
371  double y=0;
372  for ( uint i=0; i<m_nPar; ++i )
373  y += factors[i]*TMath::Power(log(x),Int_t(i));
374  return y;
375 }

◆ getLogPolNSlope()

double EtaJESCorrection::getLogPolNSlope ( const double *  factors,
double  x 
) const
private

Definition at line 377 of file EtaJESCorrection.cxx.

377  {
378  double y=0;
379  const double inv_x = 1. / x;
380  for ( uint i=0; i<m_nPar; ++i )
381  y += i*factors[i]*TMath::Power(log(x),Int_t(i-1))*inv_x;
382  return y;
383 }

◆ getLowPtJES()

double EtaJESCorrection::getLowPtJES ( double  E_uncorr,
double  eta_det 
) const
private

Definition at line 298 of file EtaJESCorrection.cxx.

298  {
299  int ieta = getEtaBin(eta_det);
300  if (m_lowPtExtrap == 0) {
301  const double *factors = m_JESFactors[ieta];
302  double E = m_minPt_JES*cosh(eta_det);
303  double R = getLogPolN(factors,E);
304  return R;
305  }
306  else if (m_lowPtExtrap == 1) {
307  double Ecutoff = m_JES_MinPt_E[ieta];
308  double Rcutoff = m_JES_MinPt_R[ieta];
309  double slope = m_JES_MinPt_Slopes[ieta];
310  double R = slope*(E_uncorr-Ecutoff)+Rcutoff;
311  return R;
312  }
313  else if(m_lowPtExtrap == 2) {
314  double minR = m_lowPtMinR;
315  double R = minR + m_JES_MinPt_Param1[ieta]*E_uncorr + m_JES_MinPt_Param2[ieta]*E_uncorr*E_uncorr;
316  return R;
317  }
318  else ATH_MSG_WARNING("Incorrect specification of low Pt JES extrapolation, please check the value of the LowPtJESExtrapolationMethod config flag.");
319  return 1;
320 }

◆ getMassCorr()

double EtaJESCorrection::getMassCorr ( double  E_corr,
double  eta_det 
) const
private

Definition at line 358 of file EtaJESCorrection.cxx.

358  {
359 
360  if (!m_applyMassCorrection) { ATH_MSG_FATAL( "You can't apply the mass correction unless you specify ApplyMassCorrection: true in the configuration file!" ); return 0; }
361 
362  int ieta = getEtaBin(eta_det);
363  const double *factors = m_JMSFactors[ieta];
364  double E = ( E_corr/cosh(eta_det)<5.0*m_GeV ? 5.0*cosh(eta_det) : E_corr/m_GeV ); // E in GeV
365 
366  double massR = getLogPolN(factors,E);
367  return 1.0/massR;
368 }

◆ getNominalResolutionData()

StatusCode JetCalibrationStep::getNominalResolutionData ( const xAOD::Jet &  ,
double &   
) const
virtualinherited

Reimplemented in JetSmearingCorrection.

Definition at line 33 of file JetCalibrationStep.cxx.

33  {
34  ATH_MSG_ERROR("Nominal data resolution requested from a jet calibration step that doesn't define it.");
35  return StatusCode::FAILURE;
36 }

◆ getNominalResolutionMC()

StatusCode JetCalibrationStep::getNominalResolutionMC ( const xAOD::Jet &  ,
double &   
) const
virtualinherited

Reimplemented in JetSmearingCorrection.

Definition at line 38 of file JetCalibrationStep.cxx.

38  {
39  ATH_MSG_ERROR("Nominal MC resolution requested from a jet calibration step that doesn't define it.");
40  return StatusCode::FAILURE;
41 }

◆ getSplineCorr()

double EtaJESCorrection::getSplineCorr ( const int  etaBin,
double  E 
) const
private

Definition at line 335 of file EtaJESCorrection.cxx.

335  {
336  double R = m_etajesFactors[ etaBin ]->Interpolate(E);
337  return R;
338 }

◆ getSplineSlope()

double EtaJESCorrection::getSplineSlope ( const int  ieta,
const double  minE 
) const
private

Definition at line 322 of file EtaJESCorrection.cxx.

322  {
323  // Don't want to use interpolation here, so instead just use the values at the bin centers near the cutoff
324  int minBin = m_etajesFactors[ieta]->FindBin(minE);
325 
326  double rFirst = m_etajesFactors[ ieta ]->GetBinContent(minBin);
327  double rSecond = m_etajesFactors[ ieta ]->GetBinContent(minBin+1);
328  double binWidth = m_etajesFactors[ ieta ]->GetBinCenter(minBin+1) - m_etajesFactors[ ieta ]->GetBinCenter(minBin);
329  double slope = (rSecond - rFirst) / binWidth;
330 
331  return slope;
332 }

◆ initialize()

StatusCode EtaJESCorrection::initialize ( )
overridevirtual

Implements JetCalibrationStep.

Definition at line 34 of file EtaJESCorrection.cxx.

34  {
35 
36  ATH_MSG_INFO("Initializing JES correction.");
37 
38  if(!m_config){
39  ATH_MSG_ERROR("EtaJES tool received a null config pointer.");
40  return StatusCode::FAILURE;
41  }
42 
43  m_jetStartScale = m_config->GetValue("EtaJESStartingScale","JetPileupScaleMomentum");
44 
45  TString absoluteJESCalibFile = m_config->GetValue("AbsoluteJES.CalibFile","");
46  if(m_dev){
47  absoluteJESCalibFile.Remove(0,33);
48  absoluteJESCalibFile.Insert(0,"JetCalibTools/");
49  }
50  else{absoluteJESCalibFile.Insert(14,m_calibAreaTag);}
51  TString calibFile = PathResolverFindCalibFile(absoluteJESCalibFile.Data());
52  m_config->ReadFile(calibFile, kEnvLocal);
53  ATH_MSG_INFO("Reading absolute calibration factors from: " << calibFile);
54  m_jesDesc = m_config->GetValue("AbsoluteJES.Description","");
55  ATH_MSG_INFO("Description: " << m_jesDesc);
56 
57  // minPt_JES (always in GeV) determines at which point we stop using the correction curve and switch to an extrapolated value
58  m_minPt_JES = m_config->GetValue(m_jetAlgo+".MinPtForETAJES",10);
59  //Which extrapolation method to use at low Et (Et < _minPt_JES)
60  m_lowPtExtrap = m_config->GetValue("LowPtJESExtrapolationMethod",0);
61  //For order 2 extrapolation only, set the minimum value of the response for Et = 0
62  m_lowPtMinR = m_config->GetValue("LowPtJESExtrapolationMinimumResponse",0.25);
63  //Allowing to use different minPt_JES depending on eta (only for extrapolation methon 1)
64  m_useSecondaryminPt_JES = m_config->GetValue(m_jetAlgo+".UseSecondaryMinPtForETAJES", false);
65  //Starting eta for secondary minPt_JES (Default |eta|>=1.9) (Used only if UseSecondaryMinPtForETAJES is true)
66  m_etaSecondaryminPt_JES = m_config->GetValue(m_jetAlgo+".EtaSecondaryMinPtForETAJES", 1.9);
67  //SecondaryminPt_JES (Default 7 GeV) (Used only if UseSecondaryMinPtForETAJES is true)
68  m_secondaryminPt_JES = m_config->GetValue(m_jetAlgo+".SecondaryMinPtForETAJES",7);
69  // Freeze JES correction at maximum values of energy for each eta bin
70  m_freezeJESatHighE = m_config->GetValue(m_jetAlgo+".FreezeJEScorrectionatHighE", false);
71  m_isSpline = m_config->GetValue(m_jetAlgo+".isSpline", false);
72 
73  // From mswiatlo, help from dag: variable eta binning
74  std::vector<double> etaBins = JetCalibUtils::VectorizeD(m_config->GetValue("JES.EtaBins",""));
75  if (etaBins.empty()){ // default binning
76  for (int i=0;i<=90; i++)
77  etaBins.push_back(0.1*i-4.5);
78  }
79  else if (etaBins.empty()) { ATH_MSG_FATAL("JES.EtaBins incorrectly specified"); return StatusCode::FAILURE; }
80  else if (etaBins.size()>s_nEtaBins+1) {
81  ATH_MSG_FATAL( "JES.EtaBins has " << etaBins.size()-1 << " bins, can be maximally 90!" );
82  return StatusCode::FAILURE;
83  }
84  m_etaBinAxis = new TAxis(etaBins.size()-1,&etaBins[0]);
85 
86  m_applyMassCorrection = m_config->GetValue("ApplyMassCorrection",false);
87 
88  if(m_mass){ // Only for the calibration sequence: EtaMassJES
89  if(m_applyMassCorrection) ATH_MSG_INFO("Jet mass correction will be applied.");
90  else { ATH_MSG_FATAL( "You can't apply the mass correction unless you specify ApplyMassCorrection: true in the configuration file!"); return StatusCode::FAILURE; }
91  }
92 
93  for (uint ieta=0;ieta<etaBins.size()-1;++ieta) {
94  if(!m_isSpline){
95  // Read in absolute JES calibration factors
96  TString key=Form("JES.%s_Bin%d",m_jetAlgo.Data(),ieta);
97  std::vector<double> params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
98  m_nPar = params.size();
99  if (m_nPar<s_nParMin || m_nPar>s_nParMax) { ATH_MSG_FATAL( "Cannot read JES calib constants " << key ); return StatusCode::FAILURE; }
100  for (uint ipar=0;ipar<m_nPar;++ipar) m_JESFactors[ieta][ipar] = params[ipar];
101 
102  //Protections for high order extrapolation methods at low Et (Et < _minPt_JES)
103  if(m_lowPtExtrap > 0) {
104  //Calculate the slope of the response curve at the minPt for each eta bin
105  //Used in the GetLowPtJES method when Pt < minPt
106  const double *factors = m_JESFactors[ieta];
107  double Ecutoff;
108  if(!m_useSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
109  else {
110  if(fabs(etaBins[ieta]) < m_etaSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
111  else{ Ecutoff = m_secondaryminPt_JES*cosh(etaBins[ieta]);}
112  }
113  const double Rcutoff = getLogPolN(factors,Ecutoff);
114  const double Slope = getLogPolNSlope(factors,Ecutoff);
115  if(Slope > Rcutoff/Ecutoff) ATH_MSG_FATAL("Slope of calibration curve at minimum ET is too steep for the JES factors of etabin " << ieta << ", eta = " << etaBins[ieta] );
116 
117  m_JES_MinPt_E[ieta] = Ecutoff;
118  m_JES_MinPt_R[ieta] = Rcutoff;
119  m_JES_MinPt_Slopes[ieta] = Slope;
120 
121  //Calculate the parameters for a 2nd order polynomial extension to the calibration curve below minimum ET
122  //Used in the GetLowPtJES method when Pt < minPt
123  if(m_lowPtExtrap == 2) {
124  const double h = m_lowPtMinR;
125  const double Param1 = (2/Ecutoff)*(Rcutoff-h)-Slope;
126  const double Param2 = (0.5/Ecutoff)*(Slope-Param1);
127  //Slope of the calibration curve should always be positive
128  if( Param1 < 0 || Param1 + 2*Param2*Ecutoff < 0) ATH_MSG_FATAL("Polynomial extension to calibration curve below minimum ET is not monotonically increasing for etabin " << ieta << ", eta = " << etaBins[ieta] );
129  m_JES_MinPt_Param1[ieta] = Param1;
130  m_JES_MinPt_Param2[ieta] = Param2;
131  }
132  }
133  }
134 
135  // Read in jet eta calibration factors
136  TString key=Form("EtaCorr.%s_Bin%d",m_jetAlgo.Data(),ieta);
137  std::vector<double> params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
138  m_nPar = params.size();
139 
140  if (params.size()!=m_nPar) { ATH_MSG_FATAL( "Cannot read jet eta calib constants " << key ); return StatusCode::FAILURE; }
141  for (uint ipar=0;ipar<m_nPar;++ipar) m_etaCorrFactors[ieta][ipar] = params[ipar];
142 
143  if(m_freezeJESatHighE){ // Read starting energy values to freeze JES correction
144  key=Form("EmaxJES.%s_Bin%d",m_jetAlgo.Data(),ieta);
145  params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
146  if (params.size()!=1) { ATH_MSG_FATAL( "Cannot read starting energy for the freezing of JES correction " << key ); return StatusCode::FAILURE; }
147  for (uint ipar=0;ipar<1;++ipar) m_energyFreezeJES[ieta] = params[ipar];
148  }
149 
150  if(m_applyMassCorrection) {
151  // Read in absolute JMS calibration factors
152  key=Form("MassCorr.%s_Bin%d",m_jetAlgo.Data(),ieta);
153  params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
154  if (params.size()!=m_nPar) {ATH_MSG_FATAL( "Cannot read JMS calib constants " << key ); return StatusCode::FAILURE;}
155  for (uint ipar=0;ipar<m_nPar;++ipar) m_JMSFactors[ieta][ipar] = params[ipar];
156  }
157 
158  }
159  if(m_isSpline){
160  TString absoluteJESCalibHists = m_config->GetValue("AbsoluteJES.CalibHists","");
161  if(m_dev){
162  absoluteJESCalibHists.Remove(0,33);
163  absoluteJESCalibHists.Insert(0,"JetCalibTools/");
164  }
165  else{
166  absoluteJESCalibHists.Insert(14,m_calibAreaTag);
167  }
168  TString calibHistFile = PathResolverFindCalibFile(absoluteJESCalibHists.Data());
169  loadSplineHists(calibHistFile, "etaJes");
170 
171  //Protections for high order extrapolation methods at low Et (Et < _minPt_JES)
172  if(m_lowPtExtrap != 1) {
173  ATH_MSG_ERROR("Only linear extrapolations are supported for p-splines currently. Please change the config file to reflect this");
174  return StatusCode::FAILURE;
175  }
176 
177  if(m_lowPtExtrap > 0) {
178  for (uint ieta=0;ieta<etaBins.size()-1;++ieta) {
179  //Calculate the slope of the response curve at the minPt for each eta bin
180  //Used in the GetLowPtJES method when Pt < minPt
181  double Ecutoff;
182  if(!m_useSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
183  else {
184  if(std::abs(etaBins[ieta]) < m_etaSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
185  else{ Ecutoff = m_secondaryminPt_JES*cosh(etaBins[ieta]);}
186  }
187  const double Rcutoff = getSplineCorr(ieta, Ecutoff);
188  const double Slope = getSplineSlope(ieta, Ecutoff);
189  if(Slope > Rcutoff/Ecutoff) ATH_MSG_WARNING("Slope of calibration curve at minimum ET is too steep for the JES factors of etabin " << ieta << ", eta = " << etaBins[ieta] );
190 
191  m_JES_MinPt_E[ieta] = Ecutoff;
192  m_JES_MinPt_R[ieta] = Rcutoff;
193  m_JES_MinPt_Slopes[ieta] = Slope;
194  }
195  }
196  }
197 
198 
199  return StatusCode::SUCCESS;
200 }

◆ initMessaging()

void AthMessaging::initMessaging ( ) const
privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

40 {
41  m_imsg = Athena::getMessageSvc();
42  m_lvl = m_imsg ?
43  static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
44  MSG::INFO;
45 }

◆ loadSplineHists()

void EtaJESCorrection::loadSplineHists ( const TString &  fileName,
const std::string &  etajes_name = "etaJes" 
)
private

Loads the calib constants from histograms in TFile named fileName.

Definition at line 204 of file EtaJESCorrection.cxx.

205 {
206  std::unique_ptr<TFile> tmpF(TFile::Open( fileName ));
207  TList *etajes_l = static_cast<TList*>( tmpF->Get(etajes_name.c_str()));
208 
209  m_etajesFactors.resize( etajes_l->GetSize() );
210  if(etajes_l->GetSize() != m_etaBinAxis->GetNbins()+1){
211  ATH_MSG_WARNING("Do not have the correct number of eta bins for " << fileName << "\t" << etajes_name << "\t" << etajes_l->GetSize() );
212  }
213 
214  for(int i=0 ; i<m_etaBinAxis->GetNbins(); i++){
215  auto *pTH1 = dynamic_cast<TH1*>(etajes_l->At(i));
216  if (not pTH1) continue;
217  m_etajesFactors[i].reset(pTH1);
218  m_etajesFactors[i]->SetDirectory(nullptr);
219  }
220  tmpF->Close();
221 }

◆ msg() [1/2]

MsgStream & asg::AsgMessaging::msg ( ) const
inherited

The standard message stream.

Returns
A reference to the default message stream of this object.

Definition at line 49 of file AsgMessaging.cxx.

49  {
50 #ifndef XAOD_STANDALONE
51  return ::AthMessaging::msg();
52 #else // not XAOD_STANDALONE
53  return m_msg;
54 #endif // not XAOD_STANDALONE
55  }

◆ msg() [2/2]

MsgStream & asg::AsgMessaging::msg ( const MSG::Level  lvl) const
inherited

The standard message stream.

Parameters
lvlThe message level to set the stream to
Returns
A reference to the default message stream, set to level "lvl"

Definition at line 57 of file AsgMessaging.cxx.

57  {
58 #ifndef XAOD_STANDALONE
59  return ::AthMessaging::msg( lvl );
60 #else // not XAOD_STANDALONE
61  m_msg << lvl;
62  return m_msg;
63 #endif // not XAOD_STANDALONE
64  }

◆ msgLvl()

bool asg::AsgMessaging::msgLvl ( const MSG::Level  lvl) const
inherited

Test the output level of the object.

Parameters
lvlThe message level to test against
Returns
boolean Indicting if messages at given level will be printed
true If messages at level "lvl" will be printed

Definition at line 41 of file AsgMessaging.cxx.

41  {
42 #ifndef XAOD_STANDALONE
43  return ::AthMessaging::msgLvl( lvl );
44 #else // not XAOD_STANDALONE
45  return m_msg.msgLevel( lvl );
46 #endif // not XAOD_STANDALONE
47  }

◆ setLevel()

void AthMessaging::setLevel ( MSG::Level  lvl)
inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

29 {
30  m_lvl = lvl;
31 }

◆ setStartP4()

StatusCode JetCalibrationStep::setStartP4 ( xAOD::Jet &  jet) const
protectedvirtualinherited

Definition at line 21 of file JetCalibrationStep.cxx.

21  {
22  xAOD::JetFourMom_t tmp;
23  if ( m_jetStartScale.compare("Default") == 0 ) {
24  return StatusCode::SUCCESS;
25  } else if ( jet.getAttribute<xAOD::JetFourMom_t>(m_jetStartScale.c_str(),tmp) ) {
26  jet.setJetP4(tmp);
27  return StatusCode::SUCCESS;
28  }
29  ATH_MSG_WARNING("Jet does not have the requested momentum state: " << m_jetStartScale);
30  return StatusCode::FAILURE;
31 }

◆ setUnitsGeV()

virtual void JetCalibrationStep::setUnitsGeV ( bool  useGeV)
inlinevirtualinherited

Definition at line 30 of file JetCalibrationStep.h.

30 { if (useGeV) m_GeV=1; else m_GeV=1000; }

Member Data Documentation

◆ ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

◆ m_applyMassCorrection

bool EtaJESCorrection::m_applyMassCorrection {}
private

Definition at line 64 of file EtaJESCorrection.h.

◆ m_calibAreaTag

TString EtaJESCorrection::m_calibAreaTag
private

Definition at line 54 of file EtaJESCorrection.h.

◆ m_config

TEnv* EtaJESCorrection::m_config {}
private

Definition at line 52 of file EtaJESCorrection.h.

◆ m_dev

bool EtaJESCorrection::m_dev {}
private

Definition at line 56 of file EtaJESCorrection.h.

◆ m_energyFreezeJES

double EtaJESCorrection::m_energyFreezeJES[s_nEtaBins] ={}
private

Definition at line 85 of file EtaJESCorrection.h.

◆ m_etaBinAxis

TAxis* EtaJESCorrection::m_etaBinAxis {}
private

Definition at line 69 of file EtaJESCorrection.h.

◆ m_etaCorrFactors

double EtaJESCorrection::m_etaCorrFactors[s_nEtaBins][s_nParMax] ={}
private

Definition at line 83 of file EtaJESCorrection.h.

◆ m_etajesFactors

std::vector<std::unique_ptr<TH1> > EtaJESCorrection::m_etajesFactors
private

Definition at line 89 of file EtaJESCorrection.h.

◆ m_etaSecondaryminPt_JES

double EtaJESCorrection::m_etaSecondaryminPt_JES {}
private

Definition at line 66 of file EtaJESCorrection.h.

◆ m_freezeJESatHighE

bool EtaJESCorrection::m_freezeJESatHighE {}
private

Definition at line 57 of file EtaJESCorrection.h.

◆ m_GeV

double JetCalibrationStep::m_GeV
protectedinherited

Definition at line 40 of file JetCalibrationStep.h.

◆ m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }
mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_isSpline

bool EtaJESCorrection::m_isSpline {}
private

Definition at line 58 of file EtaJESCorrection.h.

◆ m_JES_MinPt_E

double EtaJESCorrection::m_JES_MinPt_E[s_nEtaBins] ={}
private

Definition at line 78 of file EtaJESCorrection.h.

◆ m_JES_MinPt_Param1

double EtaJESCorrection::m_JES_MinPt_Param1[s_nEtaBins] ={}
private

Definition at line 81 of file EtaJESCorrection.h.

◆ m_JES_MinPt_Param2

double EtaJESCorrection::m_JES_MinPt_Param2[s_nEtaBins]
private

Definition at line 82 of file EtaJESCorrection.h.

◆ m_JES_MinPt_R

double EtaJESCorrection::m_JES_MinPt_R[s_nEtaBins] ={}
private

Definition at line 79 of file EtaJESCorrection.h.

◆ m_JES_MinPt_Slopes

double EtaJESCorrection::m_JES_MinPt_Slopes[s_nEtaBins] ={}
private

Definition at line 77 of file EtaJESCorrection.h.

◆ m_jesDesc

TString EtaJESCorrection::m_jesDesc
private

Definition at line 60 of file EtaJESCorrection.h.

◆ m_JESFactors

double EtaJESCorrection::m_JESFactors[s_nEtaBins][s_nParMax] ={}
private

Definition at line 76 of file EtaJESCorrection.h.

◆ m_jetAlgo

TString EtaJESCorrection::m_jetAlgo
private

Definition at line 53 of file EtaJESCorrection.h.

◆ m_jetStartScale

std::string JetCalibrationStep::m_jetStartScale
protectedinherited

Definition at line 41 of file JetCalibrationStep.h.

◆ m_JMSFactors

double EtaJESCorrection::m_JMSFactors[s_nEtaBins][s_nParMax] ={}
private

Definition at line 84 of file EtaJESCorrection.h.

◆ m_lowPtExtrap

unsigned int EtaJESCorrection::m_lowPtExtrap {}
private

Definition at line 62 of file EtaJESCorrection.h.

◆ m_lowPtMinR

double EtaJESCorrection::m_lowPtMinR {}
private

Definition at line 63 of file EtaJESCorrection.h.

◆ m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }
mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

◆ m_mass

bool EtaJESCorrection::m_mass {}
private

Definition at line 55 of file EtaJESCorrection.h.

◆ m_maxE_EtaCorr

double EtaJESCorrection::m_maxE_EtaCorr {}
private

Definition at line 61 of file EtaJESCorrection.h.

◆ m_minPt_EtaCorr

double EtaJESCorrection::m_minPt_EtaCorr {}
private

Definition at line 61 of file EtaJESCorrection.h.

◆ m_minPt_JES

double EtaJESCorrection::m_minPt_JES {}
private

Definition at line 61 of file EtaJESCorrection.h.

◆ m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls
mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

◆ m_name

std::string JetCalibrationStep::m_name
protectedinherited

Definition at line 42 of file JetCalibrationStep.h.

◆ m_nm

std::string AthMessaging::m_nm
privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ m_nPar

unsigned int EtaJESCorrection::m_nPar {}
private

Definition at line 75 of file EtaJESCorrection.h.

◆ m_secondaryminPt_JES

double EtaJESCorrection::m_secondaryminPt_JES {}
private

Definition at line 67 of file EtaJESCorrection.h.

◆ m_useSecondaryminPt_JES

bool EtaJESCorrection::m_useSecondaryminPt_JES {}
private

Definition at line 65 of file EtaJESCorrection.h.

◆ s_nEtaBins

const unsigned int EtaJESCorrection::s_nEtaBins =90
staticprivate

Definition at line 72 of file EtaJESCorrection.h.

◆ s_nParMax

const unsigned int EtaJESCorrection::s_nParMax =9
staticprivate

Definition at line 74 of file EtaJESCorrection.h.

◆ s_nParMin

const unsigned int EtaJESCorrection::s_nParMin =7
staticprivate

Definition at line 73 of file EtaJESCorrection.h.

The documentation for this class was generated from the following files:
EtaJESCorrection::loadSplineHists
void loadSplineHists(const TString &fileName, const std::string &etajes_name="etaJes")
Loads the calib constants from histograms in TFile named fileName.
Definition: EtaJESCorrection.cxx:204
AthMessaging::m_lvl
std::atomic< MSG::Level > m_lvl
Current logging level.
Definition: AthMessaging.h:138
EtaJESCorrection::m_jetAlgo
TString m_jetAlgo
Definition: EtaJESCorrection.h:53
EtaJESCorrection::m_freezeJESatHighE
bool m_freezeJESatHighE
Definition: EtaJESCorrection.h:57
EtaJESCorrection::m_mass
bool m_mass
Definition: EtaJESCorrection.h:55
EtaJESCorrection::getLogPolN
double getLogPolN(const double *factors, double x) const
Definition: EtaJESCorrection.cxx:370
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
EtaJESCorrection::m_maxE_EtaCorr
double m_maxE_EtaCorr
Definition: EtaJESCorrection.h:61
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
IDTPM::R
float R(const U &p)
Definition: TrackParametersHelper.h:101
JetCalibrationStep::setStartP4
virtual StatusCode setStartP4(xAOD::Jet &jet) const
Definition: JetCalibrationStep.cxx:21
EtaJESCorrection::m_minPt_EtaCorr
double m_minPt_EtaCorr
Definition: EtaJESCorrection.h:61
EtaJESCorrection::m_lowPtMinR
double m_lowPtMinR
Definition: EtaJESCorrection.h:63
EtaJESCorrection::getEtaCorr
double getEtaCorr(double E_corr, double eta_det) const
Definition: EtaJESCorrection.cxx:342
ConvertOldUJHistosToNewHistos.etaBins
list etaBins
Definition: ConvertOldUJHistosToNewHistos.py:145
EtaJESCorrection::getLogPolNSlope
double getLogPolNSlope(const double *factors, double x) const
Definition: EtaJESCorrection.cxx:377
EtaJESCorrection::m_minPt_JES
double m_minPt_JES
Definition: EtaJESCorrection.h:61
bin
Definition: BinsDiffFromStripMedian.h:43
AthMessaging::m_imsg
std::atomic< IMessageSvc * > m_imsg
MessageSvc pointer.
Definition: AthMessaging.h:135
x
#define x
Athena::getMessageSvc
IMessageSvc * getMessageSvc(bool quiet=false)
Definition: getMessageSvc.cxx:20
dqt_zlumi_pandas.mass
mass
Definition: dqt_zlumi_pandas.py:170
binWidth
void binWidth(TH1 *h)
Definition: listroot.cxx:80
EtaJESCorrection::m_etajesFactors
std::vector< std::unique_ptr< TH1 > > m_etajesFactors
Definition: EtaJESCorrection.h:89
config
Definition: PhysicsAnalysis/AnalysisCommon/AssociationUtils/python/config.py:1
uint
unsigned int uint
Definition: LArOFPhaseFill.cxx:20
JetCalibrationStep::JetCalibrationStep
JetCalibrationStep(const char *name="JetCalibrationStep")
Definition: JetCalibrationStep.cxx:12
xAOD::etaBin
setSAddress setEtaMS setDirPhiMS setDirZMS setBarrelRadius setEndcapAlpha setEndcapRadius setInterceptInner setEtaMap etaBin
Definition: L2StandAloneMuon_v1.cxx:148
EtaJESCorrection::m_JES_MinPt_Param1
double m_JES_MinPt_Param1[s_nEtaBins]
Definition: EtaJESCorrection.h:81
FortranAlgorithmOptions.fileName
fileName
Definition: FortranAlgorithmOptions.py:13
jet
Definition: JetCalibTools_PlotJESFactors.cxx:23
TrigConf::MSGTC::Level
Level
Definition: Trigger/TrigConfiguration/TrigConfBase/TrigConfBase/MsgStream.h:21
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
lumiFormat.i
int i
Definition: lumiFormat.py:92
EtaJESCorrection::getSplineCorr
double getSplineCorr(const int etaBin, double E) const
Definition: EtaJESCorrection.cxx:335
EtaJESCorrection::m_dev
bool m_dev
Definition: EtaJESCorrection.h:56
EtaJESCorrection::m_JES_MinPt_Param2
double m_JES_MinPt_Param2[s_nEtaBins]
Definition: EtaJESCorrection.h:82
EtaJESCorrection::m_isSpline
bool m_isSpline
Definition: EtaJESCorrection.h:58
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
DeMoUpdate.tmp
string tmp
Definition: DeMoUpdate.py:1167
EtaJESCorrection::m_energyFreezeJES
double m_energyFreezeJES[s_nEtaBins]
Definition: EtaJESCorrection.h:85
EtaJESCorrection::getSplineSlope
double getSplineSlope(const int ieta, const double minE) const
Definition: EtaJESCorrection.cxx:322
EtaJESCorrection::s_nEtaBins
static const unsigned int s_nEtaBins
Definition: EtaJESCorrection.h:72
EtaJESCorrection::m_jesDesc
TString m_jesDesc
Definition: EtaJESCorrection.h:60
xAOD::JetFourMom_t
ROOT::Math::LorentzVector< ROOT::Math::PtEtaPhiM4D< double > > JetFourMom_t
Base 4 Momentum type for Jet.
Definition: JetTypes.h:17
EtaJESCorrection::m_etaBinAxis
TAxis * m_etaBinAxis
Definition: EtaJESCorrection.h:69
EtaJESCorrection::s_nParMax
static const unsigned int s_nParMax
Definition: EtaJESCorrection.h:74
JetCalibrationStep::m_jetStartScale
std::string m_jetStartScale
Definition: JetCalibrationStep.h:41
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:195
VP1PartSpect::E
@ E
Definition: VP1PartSpectFlags.h:21
EtaJESCorrection::m_config
TEnv * m_config
Definition: EtaJESCorrection.h:52
EtaJESCorrection::m_JES_MinPt_E
double m_JES_MinPt_E[s_nEtaBins]
Definition: EtaJESCorrection.h:78
EtaJESCorrection::m_secondaryminPt_JES
double m_secondaryminPt_JES
Definition: EtaJESCorrection.h:67
EtaJESCorrection::m_etaCorrFactors
double m_etaCorrFactors[s_nEtaBins][s_nParMax]
Definition: EtaJESCorrection.h:83
EtaJESCorrection::m_JES_MinPt_Slopes
double m_JES_MinPt_Slopes[s_nEtaBins]
Definition: EtaJESCorrection.h:77
EtaJESCorrection::m_JMSFactors
double m_JMSFactors[s_nEtaBins][s_nParMax]
Definition: EtaJESCorrection.h:84
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:431
EtaJESCorrection::getLowPtJES
double getLowPtJES(double E_uncorr, double eta_det) const
Definition: EtaJESCorrection.cxx:298
EtaJESCorrection::getMassCorr
double getMassCorr(double E_corr, double eta_det) const
Definition: EtaJESCorrection.cxx:358
y
#define y
h
EtaJESCorrection::m_nPar
unsigned int m_nPar
Definition: EtaJESCorrection.h:75
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
EtaJESCorrection::m_calibAreaTag
TString m_calibAreaTag
Definition: EtaJESCorrection.h:54
TH1
Definition: rootspy.cxx:268
AthMessaging::m_nm
std::string m_nm
Message source name.
Definition: AthMessaging.h:129
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
JetCalibUtils::VectorizeD
VecD VectorizeD(const TString &str, const TString &sep=" ")
Definition: JetCalibUtils.cxx:25
PowhegControl_ttFCNC_NLO.params
params
Definition: PowhegControl_ttFCNC_NLO.py:226
EtaJESCorrection::m_etaSecondaryminPt_JES
double m_etaSecondaryminPt_JES
Definition: EtaJESCorrection.h:66
python.AutoConfigFlags.msg
msg
Definition: AutoConfigFlags.py:7
EtaJESCorrection::m_applyMassCorrection
bool m_applyMassCorrection
Definition: EtaJESCorrection.h:64
EtaJESCorrection::m_useSecondaryminPt_JES
bool m_useSecondaryminPt_JES
Definition: EtaJESCorrection.h:65
EtaJESCorrection::getJES
double getJES(double E_uncorr, double eta_det) const
Definition: EtaJESCorrection.cxx:258
EtaJESCorrection::getEtaBin
int getEtaBin(double eta_det) const
Definition: EtaJESCorrection.cxx:385
EtaJESCorrection::m_JES_MinPt_R
double m_JES_MinPt_R[s_nEtaBins]
Definition: EtaJESCorrection.h:79
EtaJESCorrection::m_lowPtExtrap
unsigned int m_lowPtExtrap
Definition: EtaJESCorrection.h:62
JetCalibrationStep::m_GeV
double m_GeV
Definition: JetCalibrationStep.h:40
EtaJESCorrection::m_JESFactors
double m_JESFactors[s_nEtaBins][s_nParMax]
Definition: EtaJESCorrection.h:76
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
Generated on Tue Jul 2 2024 21:44:51 for ATLAS Offline Software by   doxygen 1.8.18