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.
Functions providing the same interface as AthMessaging
bool	msgLvl (const MSG::Level lvl) const
	Test the output level of the object.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.

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.
int	getEtaBin (double eta_det) const
void	initMessaging () const
	Initialize our message level and MessageSvc.

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 {}
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.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels).
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging).

Static Private Attributes
static const int	s_nEtaBins =90
static const int	s_nParMin =7
static const 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 13 of file EtaJESCorrection.cxx.

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

EtaJESCorrection() [2/2]

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

Definition at line 21 of file EtaJESCorrection.cxx.

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

~EtaJESCorrection()

EtaJESCorrection::~EtaJESCorrection ( )

virtual

Definition at line 29 of file EtaJESCorrection.cxx.

                                    {
 
  if (m_etaBinAxis) delete m_etaBinAxis;
 
}

Member Function Documentation

calibrate()

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

overridevirtual

Implements JetCalibrationStep.

Definition at line 224 of file EtaJESCorrection.cxx.

                                                                        {
 
  xAOD::JetFourMom_t jetStartP4;
  ATH_CHECK( setStartP4(jet) );
  jetStartP4 = jet.jetP4();
 
  //Apply the JES calibration scale factor
  //Takes the uncorrected jet eta (in case the origin and/or 4vector jet area corrections were applied
  float detectorEta = jet.getAttribute<float>("DetectorEta");
 
  xAOD::JetFourMom_t calibP4 = jetStartP4*getJES( jetStartP4.e(), detectorEta );
 
  const double etaCorr = calibP4.eta() + getEtaCorr( calibP4.e(), detectorEta );
  double massCorr;
  if(!m_applyMassCorrection) massCorr = calibP4.mass();
  else{ massCorr = jetStartP4.mass()*getMassCorr(calibP4.e(), detectorEta); }
  TLorentzVector TLVjet;
  TLVjet.SetPtEtaPhiM( calibP4.P()/cosh(etaCorr),etaCorr,calibP4.phi(),massCorr );
  calibP4.SetPxPyPzE( TLVjet.Px(), TLVjet.Py(), TLVjet.Pz(), TLVjet.E() );
 
  if(m_dev){
    float JESFactor = calibP4.e()/jetStartP4.e();
    jet.setAttribute<float>("JetJESCalibFactor",JESFactor);
  }
  
  //Transfer calibrated jet properties to the Jet object
  jet.setAttribute<xAOD::JetFourMom_t>("JetEtaJESScaleMomentum",calibP4);
  jet.setJetP4( calibP4 );
 
  return StatusCode::SUCCESS; 
}

getEtaBin()

int EtaJESCorrection::getEtaBin ( double eta_det ) const

private

Definition at line 391 of file EtaJESCorrection.cxx.

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

getEtaCorr()

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

private

Definition at line 344 of file EtaJESCorrection.cxx.

                                                                       {
  int ieta = getEtaBin(eta_det);
  const double *factors = m_etaCorrFactors[ieta];
  
  double E = E_corr/m_GeV;
  if ( E < m_minPt_EtaCorr*std::cosh(eta_det) ) 
    E = m_minPt_EtaCorr*std::cosh(eta_det);
  if ( E>m_maxE_EtaCorr ) E=m_maxE_EtaCorr;
  
  double etaCorr = getLogPolN(factors,E);
  
  // This is ( reco_eta - truth_eta )
  // to make it an additive correction return the negative value
  return -etaCorr;
}

getJES()

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

private

Definition at line 258 of file EtaJESCorrection.cxx.

                                                                     {
 
  double E = E_uncorr/m_GeV; // E in GeV
  //Check if the Pt goes below the minimum value, if so use the special GetLowPtJES method
  const double pt = E / std::cosh(eta_det);
 
  const double minPt =
    m_useSecondaryminPt_JES && std::abs(eta_det) >= m_etaSecondaryminPt_JES
     ? m_secondaryminPt_JES : m_minPt_JES;
 
  if (pt < minPt) {
    double R = getLowPtJES(E, eta_det);
    if (R == 0.)[[unlikely]]{
      ATH_MSG_ERROR("EtaJESCorrection::getJES: Divisor 'R' is zero.");
      return 0.;
    }
    return 1.0 / R;
  }
  // Get the factors
  int ieta = getEtaBin(eta_det);
 
  // Freeze correction
  if(m_freezeJESatHighE && E>m_energyFreezeJES[ieta] && m_energyFreezeJES[ieta]!=-1) E = m_energyFreezeJES[ieta];
 
 
  // The low pT extrapolation doesn't work for the spline yet, so putting this code here.
  if(m_isSpline){
    double R = getSplineCorr(ieta, E);
    return 1.0/R;
  }
  const double *factors = m_JESFactors[ieta];
 
  
  // Calculate the jet response and then the JES as 1/R
  double R = getLogPolN(factors,E);
  if (R == 0.)[[unlikely]]{
      ATH_MSG_ERROR("EtaJESCorrection::getJES: Divisor 'R' is zero.");
      return 0.;
  }
  return 1.0/R;
}

getLogPolN()

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

private

Definition at line 376 of file EtaJESCorrection.cxx.

                                                                         {
  double y=0;
  for ( int i=0; i<m_nPar; ++i )
    y += factors[i]*std::pow(std::log(x),i);
  return y;
}

getLogPolNSlope()

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

private

Definition at line 383 of file EtaJESCorrection.cxx.

                                                                              {
  double y=0;
  const double inv_x = 1. / x;
  for ( int i=0; i<m_nPar; ++i )
    y += i*factors[i]*std::pow(std::log(x),i-1)*inv_x;
  return y;
}

getLowPtJES()

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

private

Definition at line 300 of file EtaJESCorrection.cxx.

                                                                          {
  int ieta = getEtaBin(eta_det);
  if (m_lowPtExtrap == 0) {
    const double *factors = m_JESFactors[ieta];
    double E = m_minPt_JES*cosh(eta_det);
    double R = getLogPolN(factors,E);
    return R;
  }
  else if (m_lowPtExtrap == 1) {
    double Ecutoff = m_JES_MinPt_E[ieta];
    double Rcutoff = m_JES_MinPt_R[ieta];
    double slope = m_JES_MinPt_Slopes[ieta];
    double R = slope*(E_uncorr-Ecutoff)+Rcutoff;
    return R;
  }
  else if(m_lowPtExtrap == 2) {
    double minR = m_lowPtMinR;
    double R = minR + m_JES_MinPt_Param1[ieta]*E_uncorr + m_JES_MinPt_Param2[ieta]*E_uncorr*E_uncorr;
    return R;
  }
  else ATH_MSG_WARNING("Incorrect specification of low Pt JES extrapolation, please check the value of the LowPtJESExtrapolationMethod config flag.");
  return 1;
}

getMassCorr()

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

private

Definition at line 360 of file EtaJESCorrection.cxx.

                                                                        {
 
  if (!m_applyMassCorrection) { ATH_MSG_FATAL( "You can't apply the mass correction unless you specify ApplyMassCorrection: true in the configuration file!" ); return 0; }
 
  int ieta = getEtaBin(eta_det);
  const double *factors = m_JMSFactors[ieta];
  double E = ( E_corr/std::cosh(eta_det)<5.0*m_GeV ? 5.0*std::cosh(eta_det) : E_corr/m_GeV ); // E in GeV
 
  double massR = getLogPolN(factors,E);
  if (massR == 0.) [[unlikely]]{
    ATH_MSG_ERROR( "EtaJESCorrection::getMassCorr: Divisor 'massR' is zero");
    return 0.;
  }
  return 1.0/massR;
}

getNominalResolutionData()

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

virtualinherited

Reimplemented in JetSmearingCorrection.

Definition at line 33 of file JetCalibrationStep.cxx.

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

getNominalResolutionMC()

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

virtualinherited

Reimplemented in JetSmearingCorrection.

Definition at line 38 of file JetCalibrationStep.cxx.

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

getSplineCorr()

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

private

Definition at line 337 of file EtaJESCorrection.cxx.

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

getSplineSlope()

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

private

Definition at line 324 of file EtaJESCorrection.cxx.

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

initialize()

StatusCode EtaJESCorrection::initialize ( )

overridevirtual

Implements JetCalibrationStep.

Definition at line 35 of file EtaJESCorrection.cxx.

                                        {
 
  ATH_MSG_DEBUG("Initializing JES correction.");
 
  if(!m_config){
    ATH_MSG_ERROR("EtaJES tool received a null config pointer.");
    return StatusCode::FAILURE;
  }
 
  m_jetStartScale = m_config->GetValue("EtaJESStartingScale","JetPileupScaleMomentum");
 
  TString absoluteJESCalibFile = m_config->GetValue("AbsoluteJES.CalibFile","");
  if(m_dev){
    absoluteJESCalibFile.Remove(0,33);
    absoluteJESCalibFile.Insert(0,"JetCalibTools/");
  }
  else{absoluteJESCalibFile.Insert(14,m_calibAreaTag);}
  TString calibFile = PathResolverFindCalibFile(absoluteJESCalibFile.Data());
  m_config->ReadFile(calibFile, kEnvLocal);
  ATH_MSG_DEBUG("Reading absolute calibration factors from: " << calibFile);
  m_jesDesc = m_config->GetValue("AbsoluteJES.Description","");
  ATH_MSG_DEBUG("Description: " << m_jesDesc);
 
  // minPt_JES (always in GeV) determines at which point we stop using the correction curve and switch to an extrapolated value
  m_minPt_JES = m_config->GetValue(m_jetAlgo+".MinPtForETAJES",10);
  //Which extrapolation method to use at low Et (Et < _minPt_JES)
  m_lowPtExtrap = m_config->GetValue("LowPtJESExtrapolationMethod",0);
  //For order 2 extrapolation only, set the minimum value of the response for Et = 0
  m_lowPtMinR = m_config->GetValue("LowPtJESExtrapolationMinimumResponse",0.25);
  //Allowing to use different minPt_JES depending on eta (only for extrapolation methon 1)
  m_useSecondaryminPt_JES = m_config->GetValue(m_jetAlgo+".UseSecondaryMinPtForETAJES", false);
  //Starting eta for secondary minPt_JES (Default |eta|>=1.9) (Used only if UseSecondaryMinPtForETAJES is true)
  m_etaSecondaryminPt_JES = m_config->GetValue(m_jetAlgo+".EtaSecondaryMinPtForETAJES", 1.9);
  //SecondaryminPt_JES (Default 7 GeV) (Used only if UseSecondaryMinPtForETAJES is true)
  m_secondaryminPt_JES = m_config->GetValue(m_jetAlgo+".SecondaryMinPtForETAJES",7);
  // Freeze JES correction at maximum values of energy for each eta bin
  m_freezeJESatHighE = m_config->GetValue(m_jetAlgo+".FreezeJEScorrectionatHighE", false);
  m_isSpline = m_config->GetValue(m_jetAlgo+".isSpline", false);
 
  // From mswiatlo, help from dag: variable eta binning
  std::vector<double> etaBins = JetCalibUtils::VectorizeD(m_config->GetValue("JES.EtaBins",""));
  if (etaBins.empty()){ // default binning
    for (int i=0;i<=90; i++) 
      etaBins.push_back(0.1*i-4.5);
  }
  else if (etaBins.empty()) { ATH_MSG_FATAL("JES.EtaBins incorrectly specified"); return StatusCode::FAILURE; }
  else if (etaBins.size()>s_nEtaBins+1) {
    ATH_MSG_FATAL( "JES.EtaBins has " << etaBins.size()-1 << " bins, can be maximally 90!" );
    return StatusCode::FAILURE;
  }
  m_etaBinAxis = new TAxis(etaBins.size()-1,&etaBins[0]);
 
  m_applyMassCorrection = m_config->GetValue("ApplyMassCorrection",false);
 
  if(m_mass){ // Only for the calibration sequence: EtaMassJES
    if(m_applyMassCorrection) ATH_MSG_INFO("Jet mass correction will be applied.");
    else { ATH_MSG_FATAL( "You can't apply the mass correction unless you specify ApplyMassCorrection: true in the configuration file!"); return StatusCode::FAILURE; }
  }
 
  for (uint ieta=0;ieta<etaBins.size()-1;++ieta) {
    if(!m_isSpline){
      // Read in absolute JES calibration factors
      TString key=Form("JES.%s_Bin%d",m_jetAlgo.Data(),ieta);
      std::vector<double> params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
      m_nPar = std::ssize(params);
      if (m_nPar<s_nParMin || m_nPar>s_nParMax) { ATH_MSG_FATAL( "Cannot read JES calib constants " << key ); return StatusCode::FAILURE; }
      for (int ipar=0;ipar<m_nPar;++ipar) m_JESFactors[ieta][ipar] = params[ipar];
 
        //Protections for high order extrapolation methods at low Et (Et < _minPt_JES)
        if(m_lowPtExtrap > 0) {
          //Calculate the slope of the response curve at the minPt for each eta bin
          //Used in the GetLowPtJES method when Pt < minPt
          const double *factors = m_JESFactors[ieta];
          double Ecutoff;
          if(!m_useSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
          else {
            if(fabs(etaBins[ieta]) < m_etaSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
            else{ Ecutoff = m_secondaryminPt_JES*cosh(etaBins[ieta]);}
          }
          const double Rcutoff = getLogPolN(factors,Ecutoff);
          const double Slope = getLogPolNSlope(factors,Ecutoff);
          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] );
 
          m_JES_MinPt_E[ieta] = Ecutoff;
          m_JES_MinPt_R[ieta] = Rcutoff;
          m_JES_MinPt_Slopes[ieta] = Slope;
      
          //Calculate the parameters for a 2nd order polynomial extension to the calibration curve below minimum ET
          //Used in the GetLowPtJES method when Pt < minPt
          if(m_lowPtExtrap == 2) {
            const double h = m_lowPtMinR;
            const double Param1 = (2/Ecutoff)*(Rcutoff-h)-Slope;
            const double Param2 = (0.5/Ecutoff)*(Slope-Param1);
            //Slope of the calibration curve should always be positive
            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] );
            m_JES_MinPt_Param1[ieta] = Param1;
            m_JES_MinPt_Param2[ieta] = Param2;
          }
        }
      }
 
    // Read in jet eta calibration factors
    TString key=Form("EtaCorr.%s_Bin%d",m_jetAlgo.Data(),ieta);
    std::vector<double> params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
    m_nPar = std::ssize(params);
 
    for (int ipar=0;ipar<m_nPar;++ipar) m_etaCorrFactors[ieta][ipar] = params[ipar];
 
    if(m_freezeJESatHighE){ // Read starting energy values to freeze JES correction
      key=Form("EmaxJES.%s_Bin%d",m_jetAlgo.Data(),ieta);
      params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
      if (params.size()!=1) { ATH_MSG_FATAL( "Cannot read starting energy for the freezing of JES correction " << key ); return StatusCode::FAILURE; }
      for (uint ipar=0;ipar<1;++ipar) m_energyFreezeJES[ieta] = params[ipar];
    }
 
    if(m_applyMassCorrection) {
        // Read in absolute JMS calibration factors
        key=Form("MassCorr.%s_Bin%d",m_jetAlgo.Data(),ieta);
        params = JetCalibUtils::VectorizeD(m_config->GetValue(key,""));
        if (std::ssize(params)!=m_nPar) {ATH_MSG_FATAL( "Cannot read JMS calib constants " << key ); return StatusCode::FAILURE;}
        for (int ipar=0;ipar<m_nPar;++ipar) m_JMSFactors[ieta][ipar] = params[ipar];
    }
 
  }
  if(m_isSpline){
    TString absoluteJESCalibHists = m_config->GetValue("AbsoluteJES.CalibHists","");
    if(m_dev){
      absoluteJESCalibHists.Remove(0,33);
      absoluteJESCalibHists.Insert(0,"JetCalibTools/");
    }
    else{
      absoluteJESCalibHists.Insert(14,m_calibAreaTag);
    }
    TString calibHistFile = PathResolverFindCalibFile(absoluteJESCalibHists.Data());
    loadSplineHists(calibHistFile, "etaJes");
 
    //Protections for high order extrapolation methods at low Et (Et < _minPt_JES)
    if(m_lowPtExtrap != 1) {
      ATH_MSG_ERROR("Only linear extrapolations are supported for p-splines currently. Please change the config file to reflect this");
      return StatusCode::FAILURE;
    }
 
    if(m_lowPtExtrap > 0) {
      for (uint ieta=0;ieta<etaBins.size()-1;++ieta) {
        //Calculate the slope of the response curve at the minPt for each eta bin
        //Used in the GetLowPtJES method when Pt < minPt
        double Ecutoff;
        if(!m_useSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
        else {
          if(std::abs(etaBins[ieta]) < m_etaSecondaryminPt_JES) Ecutoff = m_minPt_JES*cosh(etaBins[ieta]);
          else{ Ecutoff = m_secondaryminPt_JES*cosh(etaBins[ieta]);}
        }
        const double Rcutoff = getSplineCorr(ieta, Ecutoff);
        const double Slope = getSplineSlope(ieta, Ecutoff);
        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] );
  
        m_JES_MinPt_E[ieta] = Ecutoff;
        m_JES_MinPt_R[ieta] = Rcutoff;
        m_JES_MinPt_Slopes[ieta] = Slope;
      }
    }
  }
 
 
  return StatusCode::SUCCESS;
}

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.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

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.

{
  std::unique_ptr<TFile> tmpF(TFile::Open( fileName ));
  TList *etajes_l = static_cast<TList*>( tmpF->Get(etajes_name.c_str()));
 
  m_etajesFactors.resize( etajes_l->GetSize() );
  if(etajes_l->GetSize() != m_etaBinAxis->GetNbins()+1){
    ATH_MSG_WARNING("Do not have the correct number of eta bins for " << fileName << "\t" << etajes_name << "\t" << etajes_l->GetSize() );
  }
 
  for(int i=0 ; i<m_etaBinAxis->GetNbins(); i++){
    auto *pTH1 = dynamic_cast<TH1*>(etajes_l->At(i));
    if (not pTH1) continue;
    m_etajesFactors[i].reset(pTH1);
    m_etajesFactors[i]->SetDirectory(nullptr);
  }
  tmpF->Close();
}

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.

                                      {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg();
#else // not XAOD_STANDALONE
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msg() [2/2]

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

inherited

The standard message stream.

Parameters

lvl	The 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.

                                                          {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg( lvl );
#else // not XAOD_STANDALONE
      m_msg << lvl;
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msgLvl()

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

inherited

Test the output level of the object.

Parameters

lvl	The 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.

                                                       {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msgLvl( lvl );
#else // not XAOD_STANDALONE
      return m_msg.msgLevel( lvl );
#endif // not XAOD_STANDALONE
   }

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.

{
  m_lvl = lvl;
}

setStartP4()

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

protectedvirtualinherited

Definition at line 21 of file JetCalibrationStep.cxx.

                                                            {
  xAOD::JetFourMom_t tmp;
  if ( m_jetStartScale.compare("Default") == 0 ) {
    return StatusCode::SUCCESS;
  } else if ( jet.getAttribute<xAOD::JetFourMom_t>(m_jetStartScale.c_str(),tmp) ) {
    jet.setJetP4(tmp);
    return StatusCode::SUCCESS;
  }
  ATH_MSG_WARNING("Jet does not have the requested momentum state: " << m_jetStartScale);
  return StatusCode::FAILURE;
}

setUnitsGeV()

virtual void JetCalibrationStep::setUnitsGeV ( bool useGeV )

inlinevirtualinherited

Definition at line 28 of file JetCalibrationStep.h.

{ 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 84 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 82 of file EtaJESCorrection.h.

m_etajesFactors

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

private

Definition at line 88 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 38 of file JetCalibrationStep.h.

m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

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 80 of file EtaJESCorrection.h.

m_JES_MinPt_Param2

double EtaJESCorrection::m_JES_MinPt_Param2[s_nEtaBins] {}

private

Definition at line 81 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 39 of file JetCalibrationStep.h.

m_JMSFactors

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

private

Definition at line 83 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.

{ MSG::NIL };

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{}, m_maxE_EtaCorr{};

m_minPt_EtaCorr

double EtaJESCorrection::m_minPt_EtaCorr {}

private

Definition at line 61 of file EtaJESCorrection.h.

{}, m_minPt_EtaCorr{}, m_maxE_EtaCorr{};

m_minPt_JES

double EtaJESCorrection::m_minPt_JES {}

private

Definition at line 61 of file EtaJESCorrection.h.

{}, m_minPt_EtaCorr{}, m_maxE_EtaCorr{};

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 40 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

int EtaJESCorrection::m_nPar {}

private

Definition at line 75 of file EtaJESCorrection.h.

{}; // number of parameters in config file

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 int EtaJESCorrection::s_nEtaBins =90

staticprivate

Definition at line 72 of file EtaJESCorrection.h.

s_nParMax

const int EtaJESCorrection::s_nParMax =9

staticprivate

Definition at line 74 of file EtaJESCorrection.h.

s_nParMin

const 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.h
• EtaJESCorrection.cxx