e1hg_systematics Class Reference

#include <e1hg_systematics.h>

Collaboration diagram for e1hg_systematics:

Public Member Functions
	e1hg_systematics (const std::string &filename)
	constructor (initialization done there reading root files) More...
	~e1hg_systematics ()
double	getAlpha (int particle_type, double energy, double eta, bool interpolate=false) const
	particle_type = 0 (electrons), 1 (unconv photons), 2 (conv photons) More...

Private Attributes
TH1D *	m_helec [8]
TH1D *	m_hphot [8]
TH1D *	m_hphot2 [8]
TFile *	m_file0
const TArrayD *	m_etBins

Detailed Description

Definition at line 18 of file e1hg_systematics.h.

Constructor & Destructor Documentation

e1hg_systematics()

e1hg_systematics::e1hg_systematics

(

const std::string &

filename

)

constructor (initialization done there reading root files)

Definition at line 14 of file e1hg_systematics.cxx.

                                                            {
 
  m_file0 = TFile::Open(filename.c_str());
 
  for (Int_t ieta = 0; ieta < 8; ieta++) {
    char name[60];
    sprintf(name, "elec_%d", ieta);
    m_helec[ieta] = (TH1D*)m_file0->Get(name);
    sprintf(name, "unconv_%d", ieta);
    m_hphot[ieta] = (TH1D*)m_file0->Get(name);
    sprintf(name, "conv_%d", ieta);
    m_hphot2[ieta] = (TH1D*)m_file0->Get(name);
  }
 
  TAxis* aa = m_helec[0]->GetXaxis();
  m_etBins = aa->GetXbins();
}

~e1hg_systematics()

e1hg_systematics::~e1hg_systematics

(

)

Definition at line 33 of file e1hg_systematics.cxx.

                                    {
  m_file0->Close();
}

Member Function Documentation

getAlpha()

double e1hg_systematics::getAlpha	(	int	particle_type,
		double	energy,
		double	eta,
		bool	interpolate = false
	)		const

particle_type = 0 (electrons), 1 (unconv photons), 2 (conv photons)

energy = energy in MeV

eta

return value alpha Ebiased = E*(1 + alpha)

Definition at line 44 of file e1hg_systematics.cxx.

                                                          {
 
  // cout << " in getDelta " << endl;
  if (particle_type < 0 || particle_type > 2)
    return 0.;
 
  float aeta = fabs(eta);
  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;
 
  double energyGeV = energy * 0.001;
  double et = energyGeV / std::cosh(eta);
 
  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;
    }
  }
 
  //   cout << " energy, et , ibinEt " << energyGeV << " " << et << " " <<
  //   ibinEt << endl;
 
  Double_t scale = 0.;
  // HACK: some ES model dependency needs to be introduced
  /*Default up to es2017_summer
    if (aeta<1.80) scale=0.;
    else if (aeta<2.3) scale = 0.050*(aeta-1.8)/0.4 /0.05;
    else scale=0.025/0.05;
  */
 
  if (aeta < 1.80)
    scale = 0.;
  else if (aeta < 2.3)
    scale = (0.050 * (aeta - 1.8) / 0.4 + 0.025) / 0.05;
  else
    scale = 0.05 / 0.05;
 
  if (!interpolate) {
    if (particle_type == 0)
      return scale * m_helec[ieta]->GetBinContent(ibinEt + 1);
    else if (particle_type == 1)
      return scale * m_hphot[ieta]->GetBinContent(ibinEt + 1);
    else
      return scale * m_hphot2[ieta]->GetBinContent(
                         ibinEt + 1);  // This is 2, since if particle_type is
                                       // not 0,1,2 we have returned 0 above
  } else {
    if (particle_type == 0)
      return scale * m_helec[ieta]->Interpolate(et);
    else if (particle_type == 1)
      return scale * m_hphot[ieta]->Interpolate(et);
    else
      return scale * m_hphot2[ieta]->Interpolate(
                         et);  // This is 2, since if particle_type is not 0,1,2
                               // we have returned 0
  }
}

Member Data Documentation

m_etBins

const TArrayD* e1hg_systematics::m_etBins

private

Definition at line 40 of file e1hg_systematics.h.

m_file0

TFile* e1hg_systematics::m_file0

private

Definition at line 39 of file e1hg_systematics.h.

m_helec

TH1D* e1hg_systematics::m_helec[8]

private

Definition at line 36 of file e1hg_systematics.h.

m_hphot

TH1D* e1hg_systematics::m_hphot[8]

private

Definition at line 37 of file e1hg_systematics.h.

m_hphot2

TH1D* e1hg_systematics::m_hphot2[8]

private

Definition at line 38 of file e1hg_systematics.h.

---

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

• e1hg_systematics.h
• e1hg_systematics.cxx