#include "RecoLocalCalo/HcalRecAlgos/interface/ZdcSimpleRecAlgo.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CalibCalorimetry/HcalAlgos/interface/HcalTimeSlew.h"
#include <algorithm>
#include <iostream>
#include <math.h>

Namespaces
	ZdcSimpleRecAlgoImpl

Functions
template<class Digi , class RecHit >
RecHit	ZdcSimpleRecAlgoImpl::reco1 (const Digi &digi, const HcalCoder &coder, const HcalCalibrations &calibs, const std::vector< unsigned int > &myNoiseTS, const std::vector< unsigned int > &mySignalTS, bool slewCorrect, const HcalPulseContainmentCorrection *corr, HcalTimeSlew::BiasSetting slewFlavor)

template<class Digi , class RecHit >
RecHit	ZdcSimpleRecAlgoImpl::reco2 (const Digi &digi, const HcalCoder &coder, const HcalCalibrations &calibs, const std::vector< unsigned int > &myNoiseTS, const std::vector< unsigned int > &mySignalTS, bool slewCorrect, const HcalPulseContainmentCorrection *corr, HcalTimeSlew::BiasSetting slewFlavor)

static float	timeshift_ns_zdc (float wpksamp)

Variables
static const float	actual_ns_zdc [num_bins_zdc]

static double	MaximumFractionalError = 0.0005

static const int	num_bins_zdc = 100

static const float	scale_zdc = 0.999683

static const float	wpksamp0_zdc = 0.500053

Function Documentation

float timeshift_ns_zdc ( float wpksamp )

static

Definition at line 305 of file ZdcSimpleRecAlgo.cc.

References actual_ns_zdc, getHLTprescales::index, num_bins_zdc, scale_zdc, and wpksamp0_zdc.

Referenced by ZdcSimpleRecAlgoImpl::reco1().

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