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 More...

	~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 More...

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] {}