Correction for wrong HV EMEC presampler correction for 20.7 processed data (2015,2016) More...

#include <corr_HV_EMECPS.h>

Collaboration diagram for corr_HV_EMECPS:

Public Member Functions
	corr_HV_EMECPS ()
	constructor (initialization done there reading a root file for the HV maps per period
	~corr_HV_EMECPS ()
float	getCorr (int run, float eta, float phi) const
	get correction factor to apply to raw EMECPS energy : corrected raw EMECPS energy = correction factor * raw uncorrected EMECPS energy

Private Attributes
TH2D *	m_HV [2] {}
TFile *	m_file

Detailed Description

Correction for wrong HV EMEC presampler correction for 20.7 processed data (2015,2016)

Perform correction of E0_RAW (EMEC PS) as function of run number and phi

Definition at line 14 of file corr_HV_EMECPS.h.

Constructor & Destructor Documentation

◆ corr_HV_EMECPS()

corr_HV_EMECPS::corr_HV_EMECPS ( )

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

Definition at line 16 of file corr_HV_EMECPS.cxx.

                              {
  std::string filename = PathResolverFindCalibFile("egammaLayerRecalibTool/v2/emecps_hv_corr_2015-2016_207X.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 {
    m_HV[0]=(TH2D*) (m_file->Get("hc"));
        if (!m_HV[0]) {
            std::cerr << " cannot read histogram hc " << std::endl;
        }
    m_HV[1]=(TH2D*) (m_file->Get("ha"));
        if (!m_HV[1]) {
            std::cerr << " cannot read histogram hc " << std::endl;
        }
  }
}

◆ ~corr_HV_EMECPS()

corr_HV_EMECPS::~corr_HV_EMECPS ( )

Definition at line 35 of file corr_HV_EMECPS.cxx.

{
    m_file->Close();
}

Member Function Documentation

◆ getCorr()

float corr_HV_EMECPS::getCorr	(	int	run,
		float	eta,
		float	phi ) const

get correction factor to apply to raw EMECPS energy : corrected raw EMECPS energy = correction factor * raw uncorrected EMECPS energy

inputs: run number, eta and phi

Definition at line 41 of file corr_HV_EMECPS.cxx.

{
   float newCorr=1.;
   if (eta>-1.5 && eta<1.5) return newCorr;
   if (eta<-1.85 || eta > 1.85) return newCorr;
 
   int iside=0;
   if (eta>0.) iside=1;
 
   if (!m_HV[iside]) {
      std::cerr << " cannot find histogram to apply correction " << std::endl;
   } else {
     int ibin = std::as_const(*m_HV[iside]).GetXaxis()->FindFixBin(((Double_t)(run+0.1)));
      if (ibin<1 || ibin>7) return newCorr;
      if (phi>3.14159) phi=phi-6.283185;
      int jbin = std::as_const(*m_HV[iside]).GetYaxis()->FindFixBin(phi);
      if (jbin<1 || jbin>2048) return newCorr;
      newCorr =  m_HV[iside]->GetBinContent(ibin,jbin); 
   }
   //std::cout << " in  corr_HV_EMECPS run,eta,phi,corr " << run << " " << eta << " " << phi << " " << newCorr << std::endl;
   return newCorr;
 
}

Member Data Documentation

◆ m_file

TFile* corr_HV_EMECPS::m_file

private

Definition at line 29 of file corr_HV_EMECPS.h.

◆ m_HV

TH2D* corr_HV_EMECPS::m_HV[2] {}

private

Definition at line 28 of file corr_HV_EMECPS.h.

28{};

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

Public Member Functions

Private Attributes