CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Namespaces | Functions | Variables
HcalSimpleRecAlgo.cc File Reference
#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 <algorithm>
#include <cmath>

Go to the source code of this file.

Namespaces

 HcalSimpleRecAlgoImpl
 

Functions

static float eCorr (int ieta, int iphi, double ampl, 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 PulseShapeFitOOTPileupCorrection *psFitOOTpuCorr)
 
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)
 
static float timeshift_ns_hbheho (float wpksamp)
 
static float timeshift_ns_hf (float wpksamp)
 Same as above, but for the HF PMTs. More...
 

Variables

static const float actual_ns_hbheho [num_bins_hbheho]
 
static const float actual_ns_hf [num_bins_hf]
 
int HPDShapev3DataNum = 105
 
int HPDShapev3MCNum = 105
 
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 eCorr ( int  ieta,
int  iphi,
double  ampl,
int  runnum 
)
static

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

Definition at line 503 of file HcalSimpleRecAlgo.cc.

References corr, relval_parameters_module::energy, create_public_lumi_plots::exp, and slope.

Referenced by PFClusterProducer::PFClusterProducer(), PFMultiDepthClusterProducer::PFMultiDepthClusterProducer(), and HcalSimpleRecAlgoImpl::reco().

503  {
504 // return energy correction factor for HBM channels
505 // iphi=6 ieta=(-1,-15) and iphi=32 ieta=(-1,-7)
506 // I.Vodopianov 28 Feb. 2011
507  static const float low32[7] = {0.741,0.721,0.730,0.698,0.708,0.751,0.861};
508  static const float high32[7] = {0.973,0.925,0.900,0.897,0.950,0.935,1};
509  static const float low6[15] = {0.635,0.623,0.670,0.633,0.644,0.648,0.600,
510  0.570,0.595,0.554,0.505,0.513,0.515,0.561,0.579};
511  static const float high6[15] = {0.875,0.937,0.942,0.900,0.922,0.925,0.901,
512  0.850,0.852,0.818,0.731,0.717,0.782,0.853,0.778};
513 
514 
515  double slope, mid, en;
516  double corr = 1.0;
517 
518  if (!(iphi==6 && ieta<0 && ieta>-16) && !(iphi==32 && ieta<0 && ieta>-8))
519  return corr;
520 
521  int jeta = -ieta-1;
522  double xeta = (double) ieta;
523  if (energy > 0.) en=energy;
524  else en = 0.;
525 
526  if (iphi == 32) {
527  slope = 0.2272;
528  mid = 17.14 + 0.7147*xeta;
529  if (en > 100.) corr = high32[jeta];
530  else corr = low32[jeta]+(high32[jeta]-low32[jeta])/(1.0+exp(-(en-mid)*slope));
531  }
532  else if (iphi == 6 && runnum < 216091 ) {
533  slope = 0.1956;
534  mid = 15.96 + 0.3075*xeta;
535  if (en > 100.0) corr = high6[jeta];
536  else corr = low6[jeta]+(high6[jeta]-low6[jeta])/(1.0+exp(-(en-mid)*slope));
537  }
538 
539  // std::cout << "HBHE cell: ieta, iphi = " << ieta << " " << iphi
540  // << " -> energy = " << en << " corr = " << corr << std::endl;
541 
542  return corr;
543 }
static const double slope[3]
JetCorrectorParameters corr
Definition: classes.h:5
tuple runnum
Definition: summaryLumi.py:210
float leakCorr ( double  energy)
static

Leak correction.

Definition at line 547 of file HcalSimpleRecAlgo.cc.

References corr.

Referenced by HcalSimpleRecAlgoImpl::reco().

547  {
548  double corr = 1.0;
549  return corr;
550 }
JetCorrectorParameters corr
Definition: classes.h:5
float timeshift_ns_hbheho ( float  wpksamp)
static

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 622 of file HcalSimpleRecAlgo.cc.

References actual_ns_hbheho, cmsHarvester::index, num_bins_hbheho, and wpksamp0_hbheho.

Referenced by HcalSimpleRecAlgoImpl::reco().

622  {
623  float flx = (num_bins_hbheho-1)*(wpksamp - wpksamp0_hbheho);
624  int index = (int)flx;
625  float yval;
626 
627  if (index < 0) return actual_ns_hbheho[0];
628  else if (index >= num_bins_hbheho-1) return actual_ns_hbheho[num_bins_hbheho-1];
629 
630  // else interpolate:
631  float y1 = actual_ns_hbheho[index];
632  float y2 = actual_ns_hbheho[index+1];
633 
634  yval = y1 + (y2-y1)*(flx-(float)index);
635 
636  return yval;
637 }
return((rh^lh)&mask)
static const float wpksamp0_hbheho
static const float actual_ns_hbheho[num_bins_hbheho]
static const int num_bins_hbheho
float timeshift_ns_hf ( float  wpksamp)
static

Same as above, but for the HF PMTs.

Definition at line 746 of file HcalSimpleRecAlgo.cc.

References actual_ns_hf, cmsHarvester::index, num_bins_hf, and wpksamp0_hf.

Referenced by HcalSimpleRecAlgoImpl::recoHFTime().

746  {
747  float flx = (num_bins_hf-1)*(wpksamp-wpksamp0_hf);
748  int index = (int)flx;
749  float yval;
750 
751  if (index < 0) return actual_ns_hf[0];
752  else if (index >= num_bins_hf-1) return actual_ns_hf[num_bins_hf-1];
753 
754  // else interpolate:
755  float y1 = actual_ns_hf[index];
756  float y2 = actual_ns_hf[index+1];
757 
758  // float delta_x = 1/(float)num_bins_hf;
759  // yval = y1 + (y2-y1)*(flx-(float)index)/delta_x;
760 
761  yval = y1 + (y2-y1)*(flx-(float)index);
762  return yval;
763 }
static const float actual_ns_hf[num_bins_hf]
return((rh^lh)&mask)
static const float wpksamp0_hf
static const int num_bins_hf

Variable Documentation

const float actual_ns_hbheho[num_bins_hbheho]
static

Definition at line 558 of file HcalSimpleRecAlgo.cc.

Referenced by timeshift_ns_hbheho().

const float actual_ns_hf[num_bins_hf]
static

Definition at line 642 of file HcalSimpleRecAlgo.cc.

Referenced by timeshift_ns_hf().

int HPDShapev3DataNum = 105

Definition at line 14 of file HcalSimpleRecAlgo.cc.

Referenced by HcalSimpleRecAlgo::setForData().

int HPDShapev3MCNum = 105

Definition at line 15 of file HcalSimpleRecAlgo.cc.

Referenced by HcalSimpleRecAlgo::setForData().

double MaximumFractionalError = 0.002

Definition at line 13 of file HcalSimpleRecAlgo.cc.

const int num_bins_hbheho = 61
static

Definition at line 556 of file HcalSimpleRecAlgo.cc.

Referenced by timeshift_ns_hbheho().

const int num_bins_hf = 101
static

Definition at line 639 of file HcalSimpleRecAlgo.cc.

Referenced by timeshift_ns_hf().

const float wpksamp0_hbheho = 0.5
static

Definition at line 555 of file HcalSimpleRecAlgo.cc.

Referenced by timeshift_ns_hbheho().

const float wpksamp0_hf = 0.5
static

Definition at line 640 of file HcalSimpleRecAlgo.cc.

Referenced by timeshift_ns_hf().