corr_pileupShift Class Reference

Correction for pileup induced energy shit as function of mu per layer for 2016 data. More...

#include <corr_pileupShift.h>

Collaboration diagram for corr_pileupShift:

Public Member Functions
	corr_pileupShift ()
	constructor (initialization done there reading a root file for the HV maps per period More...
	~corr_pileupShift ()
float	getCorr (int layer, int run, float mu, float eta) const
	get shift to subtract to raw layer energy : corrected layer energy = raw energy - shift More...

Private Attributes
int	m_runMin
int	m_runMax
TH1D *	m_haverage [4] {}
TFile *	m_file

Detailed Description

Correction for pileup induced energy shit as function of mu per layer for 2016 data.

Author

: Guillaume Unal

Definition at line 13 of file corr_pileupShift.h.

Constructor & Destructor Documentation

corr_pileupShift()

corr_pileupShift::corr_pileupShift

(

)

constructor (initialization done there reading a root file for the HV maps per period

Definition at line 9 of file corr_pileupShift.cxx.

                                  {
 
  std::string filename = PathResolverFindCalibFile("egammaLayerRecalibTool/v2/layer_average.root");
    
  m_file = TFile::Open(filename.c_str());
  if (not m_file or m_file->IsZombie()) {
    std::cerr << "FATAL: cannot open " << filename << std::endl;
  }
  else {
    for (int layer=0;layer<4;layer++) {
      char name[12];
      snprintf(name,sizeof(name),"average_%d",layer);
      m_haverage[layer]=(TH1D*) (m_file->Get(name));
      if (not m_haverage[layer]) {
      std::cerr << "FATAL: cannot find " << name << std::endl;
    }
    }
  }
  m_runMin=296939;
  m_runMax=311563;
}

~corr_pileupShift()

corr_pileupShift::~corr_pileupShift

(

)

Definition at line 31 of file corr_pileupShift.cxx.

{
    m_file->Close();
}

Member Function Documentation

getCorr()

float corr_pileupShift::getCorr	(	int	layer,
		int	run,
		float	mu,
		float	eta
	)		const

get shift to subtract to raw layer energy : corrected layer energy = raw energy - shift

inputs: layer (0 to 3), run (run number), mu (average number of interactions per crossing) and eta

Definition at line 37 of file corr_pileupShift.cxx.

{
   
   if (layer<0 || layer>3) return 0.;
   if (run<m_runMin || run>m_runMax) return 0.;
   double aeta=std::fabs(eta);
   if (aeta> 2.5) return 0.;
   if (aeta>2.37) aeta=2.37;
 
   int ibin = std::as_const(*m_haverage[layer]).GetXaxis()->FindFixBin(aeta);
   float corr = mu * m_haverage[layer]->GetBinContent(ibin);
 
   return corr;
 
}

Member Data Documentation

m_file

TFile* corr_pileupShift::m_file

private

Definition at line 31 of file corr_pileupShift.h.

m_haverage

TH1D* corr_pileupShift::m_haverage[4] {}

private

Definition at line 30 of file corr_pileupShift.h.

m_runMax

int corr_pileupShift::m_runMax

private

Definition at line 27 of file corr_pileupShift.h.

m_runMin

int corr_pileupShift::m_runMin

private

Definition at line 27 of file corr_pileupShift.h.

---

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

• corr_pileupShift.h
• corr_pileupShift.cxx