HcalCorrectionFunctions.h File Reference

Go to the source code of this file.

Functions
float	hbminus_special_ecorr (int ieta, int iphi, double energy, int runnum)
	Special energy correction for some HB- cells. More...
float	timeshift_ns_hbheho (float wpksamp)

Function Documentation

hbminus_special_ecorr()

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

Special energy correction for some HB- cells.

Special energy correction for some HB- cells.

Definition at line 463 of file HcalSimpleRecAlgo.cc.

References alignCSCRings::corr, HCALHighEnergyHPDFilter_cfi::energy, JetChargeProducer_cfi::exp, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, METSignificanceParams_cfi::jeta, runnum, 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;
}

timeshift_ns_hbheho()

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

References actual_ns_hbheho, createfilelist::int, num_bins_hbheho, wpksamp0_hbheho, testProducerWithPsetDescEmpty_cfi::y1, testProducerWithPsetDescEmpty_cfi::y2, and makeMEIFBenchmarkPlots::yval.

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;
}