#include "RecoLocalCalo/HcalRecAlgos/interface/HcalSimpleRecAlgo.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CalibCalorimetry/HcalAlgos/interface/HcalTimeSlew.h"
#include "RecoLocalCalo/HcalRecAlgos/src/HcalTDCReco.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/rawEnergy.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HcalCorrectionFunctions.h"
#include "DataFormats/METReco/interface/HcalCaloFlagLabels.h"
#include "CondFormats/DataRecord/interface/HcalTimeSlewRecord.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include <algorithm>
#include <cmath>

Namespaces
	HcalSimpleRecAlgoImpl

Functions
float	hbminus_special_ecorr (int ieta, int iphi, double energy, int runnum)
	Ugly hack to apply energy corrections to some HB- cells. More...

static float	leakCorr (double energy)
	Leak correction. More...

template<class Digi , class RecHit >
RecHit	HcalSimpleRecAlgoImpl::reco (const Digi &digi, const HcalCoder &coder, const HcalCalibrations &calibs, const int ifirst, const int n, const bool slewCorrect, const bool pulseCorrect, const HcalPulseContainmentCorrection corr, const HcalTimeSlew::BiasSetting slewFlavor, const int runnum, const bool useLeak, const AbsOOTPileupCorrection pileupCorrection, const BunchXParameter bxInfo, const unsigned lenInfo, const int puCorrMethod, const HcalTimeSlew hcalTimeSlew_delay_)

template<class Digi >
float	HcalSimpleRecAlgoImpl::recoHFTime (const Digi &digi, const int maxI, const double amp_fC, const bool slewCorrect, double maxA, float t0, float t2)

template<class Digi >
void	HcalSimpleRecAlgoImpl::removePileup (const Digi &digi, const HcalCoder &coder, const HcalCalibrations &calibs, const int ifirst, const int n, const bool pulseCorrect, const HcalPulseContainmentCorrection corr, const AbsOOTPileupCorrection pileupCorrection, const BunchXParameter bxInfo, const unsigned lenInfo, double p_maxA, double p_ampl, double p_uncorr_ampl, double p_fc_ampl, int p_nRead, int p_maxI, bool leakCorrApplied, float p_t0, float p_t2)

float	timeshift_ns_hbheho (float wpksamp)

static float	timeshift_ns_hf (float wpksamp)
	Timeshift correction for the HF PMTs. More...

Variables
static const float	actual_ns_hbheho [num_bins_hbheho]

static const float	actual_ns_hf [num_bins_hf]

double	MaximumFractionalError = 0.002

static const int	num_bins_hbheho = 61

static const int	num_bins_hf = 101

static const float	wpksamp0_hbheho = 0.5

static const float	wpksamp0_hf = 0.5

Function Documentation

float hbminus_special_ecorr	(	int	ieta,
		int	iphi,
		double	energy,
		int	runnum
	)

Ugly hack to apply energy corrections to some HB- cells.

Special energy correction for some HB- cells.

Definition at line 412 of file HcalSimpleRecAlgo.cc.

References corr, JetChargeProducer_cfi::exp, METSignificanceParams_cfi::jeta, and slope.

Referenced by SimpleHBHEPhase1Algo::hbminusCorrectionFactor().

                                                                            {
 // return energy correction factor for HBM channels 
 // iphi=6 ieta=(-1,-15) and iphi=32 ieta=(-1,-7)
 // I.Vodopianov 28 Feb. 2011
   static const float low32[7]  = {0.741,0.721,0.730,0.698,0.708,0.751,0.861};
   static const float high32[7] = {0.973,0.925,0.900,0.897,0.950,0.935,1};
   static const float low6[15]  = {0.635,0.623,0.670,0.633,0.644,0.648,0.600,
                   0.570,0.595,0.554,0.505,0.513,0.515,0.561,0.579};
   static const float high6[15] = {0.875,0.937,0.942,0.900,0.922,0.925,0.901,
                   0.850,0.852,0.818,0.731,0.717,0.782,0.853,0.778};
 
   
   double slope, mid, en;
   double corr = 1.0;
 
   if (!(iphi==6 && ieta<0 && ieta>-16) && !(iphi==32 && ieta<0 && ieta>-8)) 
     return corr;
 
   int jeta = -ieta-1;
   double xeta = (double) ieta;
   if (energy > 0.) en=energy;
   else en = 0.;
 
   if (iphi == 32) {
     slope = 0.2272;
     mid = 17.14 + 0.7147*xeta;
     if (en > 100.) corr = high32[jeta];
     else corr = low32[jeta]+(high32[jeta]-low32[jeta])/(1.0+exp(-(en-mid)*slope));
   }
   else if (iphi == 6 && runnum < 216091 ) {
     slope = 0.1956;
     mid = 15.96 + 0.3075*xeta;
     if (en > 100.0) corr = high6[jeta];
     else corr = low6[jeta]+(high6[jeta]-low6[jeta])/(1.0+exp(-(en-mid)*slope));
   }
 
   //  std::cout << "HBHE cell:  ieta, iphi = " << ieta << "  " << iphi 
   //        << "  ->  energy = " << en << "   corr = " << corr << std::endl;
 
   return corr;
 }

float leakCorr ( double energy )

static

Leak correction.

Definition at line 456 of file HcalSimpleRecAlgo.cc.

References corr.

Referenced by HcalSimpleRecAlgoImpl::reco(), and HcalSimpleRecAlgo::setBXInfo().

                               {
   double corr = 1.0;
   return corr;
 }

float timeshift_ns_hbheho ( float wpksamp )

Timeshift correction for HPDs based on the position of the peak ADC measurement. Allows for an accurate determination of the relative phase of the pulse shape from the HPD. Calculated based on a weighted sum of the -1,0,+1 samples relative to the peak as follows: wpksamp = (0*sample[0] + 1*sample[1] + 2*sample[2]) / (sample[0] + sample[1] + sample[2]) where sample[1] is the maximum ADC sample value.

Definition at line 531 of file HcalSimpleRecAlgo.cc.

References actual_ns_hbheho, createfilelist::int, num_bins_hbheho, and wpksamp0_hbheho.

Referenced by HcalSimpleRecAlgoImpl::reco().

                                          {
   float flx = (num_bins_hbheho-1)*(wpksamp - wpksamp0_hbheho);
   int index = (int)flx;
   float yval;
 
   if      (index <    0)               return actual_ns_hbheho[0];
   else if (index >= num_bins_hbheho-1) return actual_ns_hbheho[num_bins_hbheho-1];
 
   // else interpolate:
   float y1 = actual_ns_hbheho[index];
   float y2 = actual_ns_hbheho[index+1];
 
   yval = y1 + (y2-y1)*(flx-(float)index);
 
   return yval;
 }

float timeshift_ns_hf ( float wpksamp )

static

Timeshift correction for the HF PMTs.

Definition at line 655 of file HcalSimpleRecAlgo.cc.

References actual_ns_hf, createfilelist::int, num_bins_hf, and wpksamp0_hf.

Referenced by HcalSimpleRecAlgoImpl::recoHFTime(), and HcalSimpleRecAlgo::setBXInfo().

                                      {
   float flx = (num_bins_hf-1)*(wpksamp-wpksamp0_hf);
   int index = (int)flx;
   float yval;
   
   if      (index <  0)             return actual_ns_hf[0];
   else if (index >= num_bins_hf-1) return actual_ns_hf[num_bins_hf-1];
 
   // else interpolate:
   float y1       = actual_ns_hf[index];
   float y2       = actual_ns_hf[index+1];
 
   // float delta_x  = 1/(float)num_bins_hf;
   // yval = y1 + (y2-y1)*(flx-(float)index)/delta_x;
 
   yval = y1 + (y2-y1)*(flx-(float)index);
   return yval;
 }