ZSEnergy_impl Namespace Reference

Functions
template<class DIGI >
bool	keepMe (const HcalDbService &db, const DIGI &inp, int threshold, int firstSample, int samples, bool twosided)

Function Documentation

template<class DIGI >

bool ZSEnergy_impl::keepMe	(	const HcalDbService &	db,
		const DIGI &	inp,
		int	threshold,
		int	firstSample,
		int	samples,
		bool	twosided
	)

Definition at line 13 of file HcalZSAlgoEnergy.cc.

References HcalQIECoder::charge(), CastorSimpleReconstructor_cfi::firstSample, HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalDbService::getPedestal(), HcalPedestal::getValue(), createfilelist::int, PFRecoTauDiscriminationByIsolation_cfi::offset, and slope.

                                                                                                                    {
    bool keepIt = false;
    const HcalQIECoder *channelCoder = db.getHcalCoder(inp.id());
    const HcalQIEShape *shape = db.getHcalShape(channelCoder);

    // determine average pedestal for channel
    float pedsum = 0, pedave = 0, offset = 0, slope = 0;
    const HcalPedestal *ped = db.getPedestal(inp.id());
    for (int j = 0; j < 4; j++) {
      pedave += ped->getValue(j) / 4.0;
      offset += channelCoder->charge(*shape, 0, j) / 4.0;
      slope += channelCoder->charge(*shape, 1, j) / 4.0;
    }
    slope -= offset;
    pedave -= offset;
    pedave /= slope;

    int sum = 0;

    for (int j = 0; j < samples && j + firstSample < (int)inp.size(); j++) {
      sum += inp[j + firstSample].adc();
      pedsum += pedave;
    }
    //    int presum=sum;
    sum -= (int)(pedsum);

    if (sum >= threshold)
      keepIt = true;
    else if (sum <= (-threshold) && twosided)
      keepIt = true;
    /*
    else
   std::cout << inp.id() << " " << sum << ":" << presum << " " << threshold
    << " " << pedsum << " " << pedave
    << " " << offset << " " << slope
    << std::endl;
  */
    return keepIt;
  }