37 std::vector<double> shortEnergyParams,
38 std::vector<double> shortETParams,
39 std::vector<double> longR,
40 std::vector<double> longEnergyParams,
41 std::vector<double> longETParams,
42 int HcalAcceptSeverityLevel,
43 std::vector<double> shortR29,
44 std::vector<double> longR29)
83 double ET = energy/fabs(cosh(fEta));
86 double ETthresh=0, Energythresh=0;
99 if (energy<Energythresh || ET < ETthresh)
124 HcalDetId neighbor(part->detid().rawId());
131 Ratio=(energy-part->energy())/(energy+part->energy());
134 double RatioThresh=0;
146 if (Ratio<=RatioThresh)
171 threshold+=params[
i]*
pow(abs_x, (
int)
i);
std::vector< double > long_ET_Thresh
void setFlagField(uint32_t value, int base, int width=1)
std::vector< T >::const_iterator const_iterator
std::vector< double > long_Energy_Thresh
std::vector< double > short_R
uint32_t rawId() const
get the raw id
void HFSetFlagFromPET(HFRecHit &hf, HFRecHitCollection &rec, HcalChannelQuality *myqual, const HcalSeverityLevelComputer *mySeverity)
int depth() const
get the tower depth
std::vector< double > long_R
int ieta() const
get the cell ieta
bool dropChannel(const uint32_t &mystatus) const
std::vector< double > long_R_29
const_iterator end() const
int iphi() const
get the cell iphi
static const double theHFEtaBounds[]
int HcalAcceptSeverityLevel_
std::vector< double > short_ET_Thresh
std::vector< double > short_Energy_Thresh
int getSeverityLevel(const DetId &myid, const uint32_t &myflag, const uint32_t &mystatus) const
iterator find(key_type k)
HcalDetId id() const
get the id
uint32_t getValue() const
const Item * getValues(DetId fId) const
Power< A, B >::type pow(const A &a, const B &b)
double CalcThreshold(double abs_energy, std::vector< double > params)
std::vector< double > short_R_29