get resolution systematics from material effect (mc12c setup, MVA calibration) More...

#include <get_MaterialResolutionEffect.h>

Inheritance diagram for get_MaterialResolutionEffect:

Collaboration diagram for get_MaterialResolutionEffect:

Public Member Functions
	get_MaterialResolutionEffect ()
	constructor (initialization done there reading root files with resolution fit parameters More...

	~get_MaterialResolutionEffect ()

double	getDelta (int particle_type, double energy, double eta, int response_type, int isyst) const
	get material effect on resolution from distorted geometry as difference to 40 GeV Et electrons smearing More...

double	interpolateTH1 (TH1 *hist, double x, bool abs_bins) const

void	store_IBL_PP0_YProjections ()

void	setInterpolate (bool interpolate)

void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
std::array< std::array< std::array< std::unique_ptr< TH1 >, 3 >, s_nEtaBins >, s_nSys >	m_hSystPeak

std::array< std::array< std::array< std::unique_ptr< TH1 >, 3 >, s_nEtaBins >, s_nSys >	m_hSystResol

std::array< std::unique_ptr< TH2 >, 3 >	m_hsyst_IBL_PP0

std::array< std::vector< std::unique_ptr< TH1 > >, 3 >	m_hsyst_IBL_PP0_ProjectionY

bool	m_interpolate = false

const TArrayD *	m_etBins

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 int	s_nEtaBins = 8

static const int	s_nSys = 4

Detailed Description

get resolution systematics from material effect (mc12c setup, MVA calibration)

Different parameterizations (gaussian core, sigma eff 90% and 80% from crystal ball fits) are available

Definition at line 28 of file get_MaterialResolutionEffect.h.

Constructor & Destructor Documentation

◆ get_MaterialResolutionEffect()

get_MaterialResolutionEffect::get_MaterialResolutionEffect ( )

constructor (initialization done there reading root files with resolution fit parameters

Definition at line 36 of file get_MaterialResolutionEffect.cxx.

     : asg::AsgMessaging("get_MaterialResolutionEffect") {
  
   const std::string pType[3] = {"Elec", "Unconv", "Conv"};
   const std::string sType[s_nSys] = {"A", "CD", "EL", "FMX"};
  
   const std::string filename = PathResolverFindCalibFile(
       "ElectronPhotonFourMomentumCorrection/v8/"
       "histos-systematics-material.root");
  
   std::unique_ptr<TFile> file0(TFile::Open(filename.c_str(), "READ"));
  
   for (Int_t isys = 0; isys < s_nSys; isys++) {  // 0=A, 1=CD, 2=EL, 3=FMX
     for (Int_t ieta = 0; ieta < s_nEtaBins; ieta++) {
       for (Int_t iconv = 0; iconv < 3;
            iconv++) {  // 0=electron, 1=unconverted, 2=converted
         char name[31];
         sprintf(name, "syst%s_%s_etaBin_%d", pType[iconv].c_str(),
                 sType[isys].c_str(), ieta);
         m_hSystPeak.at(isys).at(ieta).at(iconv).reset(
             checked_cast<TH1*>(file0->Get(name)));
         m_hSystPeak.at(isys).at(ieta).at(iconv)->SetDirectory(nullptr);
  
         sprintf(name, "syst%s_sigmaG_%s_etaBin_%d", pType[iconv].c_str(),
                 sType[isys].c_str(), ieta);
         m_hSystResol.at(isys).at(ieta).at(iconv).reset(
             checked_cast<TH1*>(file0->Get(name)));
         m_hSystResol.at(isys).at(ieta).at(iconv)->SetDirectory(nullptr);
       }
     }
   }
  
   // IBL+PP0 material systematics stored in 2D histograms
   m_hsyst_IBL_PP0.at(0).reset(
       checked_cast<TH2*>(file0->Get("systElec_IBLPP0")));
   m_hsyst_IBL_PP0.at(0)->SetDirectory(nullptr);
  
   m_hsyst_IBL_PP0.at(1).reset(
       checked_cast<TH2*>(file0->Get("systUnconv_IBLPP0")));
   m_hsyst_IBL_PP0.at(1)->SetDirectory(nullptr);
  
   m_hsyst_IBL_PP0.at(2).reset(
       checked_cast<TH2*>(file0->Get("systConv_IBLPP0")));
   m_hsyst_IBL_PP0.at(2)->SetDirectory(nullptr);
  
   m_etBins = m_hSystResol.at(0).at(0).at(1)->GetXaxis()->GetXbins();
 }

◆ ~get_MaterialResolutionEffect()

get_MaterialResolutionEffect::~get_MaterialResolutionEffect ( )

inline

Definition at line 34 of file get_MaterialResolutionEffect.h.

34 {};

Member Function Documentation

◆ getDelta()

double get_MaterialResolutionEffect::getDelta	(	int	particle_type,
		double	energy,
		double	eta,
		int	response_type,
		int	isyst
	)		const

get material effect on resolution from distorted geometry as difference to 40 GeV Et electrons smearing

particle_type: 0=electron, 1=reco unconverted photon, 2=reco converted photon

energy : energy in MeV

eta : eta

response_type: 0=gaussian core, 1=sigma eff 80%, 2=sigma eff 90%

isyst : 0 configA, 1 configCD, 2 configE, 3 configFMX

returned value is sigmaE/E change in quadrature to resolution

Definition at line 94 of file get_MaterialResolutionEffect.cxx.

                                                                {
  
   // cout << " in getDelta " << endl;
   if (particle_type < 0 || particle_type > 2)
     return -999;
   if (response_type < 0 || response_type > 1)
     return -999;
  
   float aeta = std::fabs(eta);
   double energyGeV = energy * 0.001;
   double et = energyGeV / cosh(eta);
  
   // IBL+PP0
   if (isyst == 5) {
     double et2 = et;
     if (et < 5.) {
       et2 = 5.1;
     }
     if (et > 2000) {
       et2 = 1999.;
     }
     if (aeta>=2.5) aeta=2.49; 
  
     int ieta = m_hsyst_IBL_PP0.at(particle_type)->GetXaxis()->FindBin(aeta);
     int iet = m_hsyst_IBL_PP0.at(particle_type)->GetYaxis()->FindBin(et2);
     if (m_interpolate) {
       return 0.01*interpolateTH1(m_hsyst_IBL_PP0_ProjectionY[particle_type][ieta - 1].get(), et2, true);
     }
     else {
       return 0.01*m_hsyst_IBL_PP0.at(particle_type)->GetBinContent(ieta, iet);
     }
   }
  
   int ieta = 0;
   if (aeta < 0.4) {
     ieta = 0;
   } else if (aeta < 0.8) {
     ieta = 1;
   } else if (aeta < 1.1) {
     ieta = 2;
   } else if (aeta < 1.37) {
     ieta = 3;
   } else if (aeta < 1.52) {
     ieta = 4;
   } else if (aeta < 1.80) {
     ieta = 5;
   } else if (aeta < 2.10) {
     ieta = 6;
   } else {
     ieta = 7;
   }
  
   int ibinEt = m_etBins->GetSize() - 2;
   for (int i = 1; i < m_etBins->GetSize(); i++) {
     if (et < m_etBins->GetAt(i)) {
       ibinEt = i - 1;
       break;
     }
   }
  
   auto& hist = response_type==0 ? m_hSystPeak : m_hSystResol;
   if (m_interpolate) {
     return 0.01*interpolateTH1(hist.at(isyst).at(ieta).at(particle_type).get(), et, true);
   }
   else {
     return 0.01*hist.at(isyst).at(ieta).at(particle_type)->GetBinContent(ibinEt+1);
   }
  
 }

◆ 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();
   m_lvl = m_imsg ?
     static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
     MSG::INFO;
 }

◆ interpolateTH1()

double get_MaterialResolutionEffect::interpolateTH1	(	TH1 *	hist,
		double	x,
		bool	abs_bins
	)		const

Definition at line 183 of file get_MaterialResolutionEffect.cxx.

 {
   if (!hist) return 0.;
   int nbins = hist->GetNbinsX();
  
   if(x<=hist->GetBinCenter(1)) {
     return abs_bins? std::abs(hist->GetBinContent(1)) : hist->GetBinContent(1);
   } 
   else if(x>=hist->GetBinCenter(nbins)) {
     return abs_bins? std::abs(hist->GetBinContent(nbins)) : hist->GetBinContent(nbins);
   } 
   else {
     int xbin = hist->FindBin(x);
     int xbin0, xbin1;
     if(x<=hist->GetBinCenter(xbin)) {
       xbin0 = xbin - 1;
       xbin1 = xbin;
     } 
     else {
       xbin0 = xbin;
       xbin1 = xbin + 1;
     }
     double y0 = abs_bins? std::abs(hist->GetBinContent(xbin0)) : hist->GetBinContent(xbin0);
     double x0 = hist->GetBinCenter(xbin0);
     double y1 = abs_bins? std::abs(hist->GetBinContent(xbin1)) : hist->GetBinContent(xbin1);
     double x1 = hist->GetBinCenter(xbin1);
     return y0 + (x-x0)*((y1-y0)/(x1-x0));
   }
 }

◆ 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
    }

◆ setInterpolate()

void get_MaterialResolutionEffect::setInterpolate ( bool interpolate )

inline

Definition at line 50 of file get_MaterialResolutionEffect.h.

                                         {
     m_interpolate = interpolate;
     store_IBL_PP0_YProjections();
   }

◆ 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;
 }

◆ store_IBL_PP0_YProjections()

void get_MaterialResolutionEffect::store_IBL_PP0_YProjections ( )

Definition at line 166 of file get_MaterialResolutionEffect.cxx.

 {
   for (size_t i = 0; i < m_hsyst_IBL_PP0.size(); i++) {
     if (m_hsyst_IBL_PP0[i]) {
       int nEtaBins = m_hsyst_IBL_PP0[i]->GetNbinsX();
       // Resize the inner vector to match the size of TH2 x-axis bins
       m_hsyst_IBL_PP0_ProjectionY[i].resize(nEtaBins);
       for (int ieta : std::views::iota(1, nEtaBins + 1)) {
         std::string histName = std::format("h1d_IBL_PP0_ptype{}_EtaBin_{}", i, ieta);
         m_hsyst_IBL_PP0_ProjectionY[i][ieta - 1].reset( dynamic_cast<TH1*>(m_hsyst_IBL_PP0[i]->ProjectionY(histName.c_str(), ieta, ieta)) );
         m_hsyst_IBL_PP0_ProjectionY[i][ieta - 1]->SetDirectory( nullptr );
       }
     }      
   }
 }

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_etBins

const TArrayD* get_MaterialResolutionEffect::m_etBins

private

Definition at line 66 of file get_MaterialResolutionEffect.h.

◆ m_hsyst_IBL_PP0

std::array<std::unique_ptr<TH2>, 3> get_MaterialResolutionEffect::m_hsyst_IBL_PP0

private

Definition at line 61 of file get_MaterialResolutionEffect.h.

◆ m_hsyst_IBL_PP0_ProjectionY

std::array<std::vector<std::unique_ptr<TH1> >, 3> get_MaterialResolutionEffect::m_hsyst_IBL_PP0_ProjectionY

private

Definition at line 64 of file get_MaterialResolutionEffect.h.

◆ m_hSystPeak

std::array<std::array<std::array<std::unique_ptr<TH1>, 3>, s_nEtaBins>, s_nSys> get_MaterialResolutionEffect::m_hSystPeak

private

Definition at line 59 of file get_MaterialResolutionEffect.h.

◆ m_hSystResol

std::array<std::array<std::array<std::unique_ptr<TH1>, 3>, s_nEtaBins>, s_nSys> get_MaterialResolutionEffect::m_hSystResol

private

Definition at line 60 of file get_MaterialResolutionEffect.h.

◆ m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_interpolate

bool get_MaterialResolutionEffect::m_interpolate = false

private

Definition at line 65 of file get_MaterialResolutionEffect.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_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_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ s_nEtaBins

const int get_MaterialResolutionEffect::s_nEtaBins = 8

staticprivate

Definition at line 57 of file get_MaterialResolutionEffect.h.

◆ s_nSys

const int get_MaterialResolutionEffect::s_nSys = 4

staticprivate

Definition at line 58 of file get_MaterialResolutionEffect.h.

The documentation for this class was generated from the following files:

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ get_MaterialResolutionEffect()

◆ ~get_MaterialResolutionEffect()

Member Function Documentation

◆ getDelta()

◆ initMessaging()

◆ interpolateTH1()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ setInterpolate()

◆ setLevel()

◆ store_IBL_PP0_YProjections()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_etBins

◆ m_hsyst_IBL_PP0

◆ m_hsyst_IBL_PP0_ProjectionY

◆ m_hSystPeak

◆ m_hSystResol

◆ m_imsg

◆ m_interpolate

◆ m_lvl

◆ m_msg_tls

◆ m_nm

◆ s_nEtaBins

◆ s_nSys