HcalTimingCorrector_impl Namespace Reference

Functions
template<class T , class V >
void	makeCorrection (T &rechit, const V &digi, int favorite_capid)

Function Documentation

template<class T , class V >

void HcalTimingCorrector_impl::makeCorrection	(	T &	rechit,
		const V &	digi,
		int	favorite_capid
	)

Definition at line 8 of file HcalTimingCorrector.cc.

References HcalCaloFlagLabels::TimingAddedBit, HcalCaloFlagLabels::TimingErrorBit, and HcalCaloFlagLabels::TimingSubtractedBit.

                                                                    {
    if (digi.fiberIdleOffset() == -1000) return; //reserved for bad stuff

    bool bad = false;
    //convention, -1 means digi is early, need to add
    //            +1 means digi is late, need to subtract
    int capid_shift = (favorite_capid - digi[0].capid());
    capid_shift %= 4;
    while (capid_shift < 0) capid_shift += 4;
    
    int idle_shift = digi.fiberIdleOffset();
    idle_shift %= 4;
    while (idle_shift < 0) idle_shift += 4;
    
    if (capid_shift != idle_shift) {
      bad = true;
    }
    else if (capid_shift == 0) {
      // do nothing, all is well
    }
    else if (capid_shift == 2) {
      bad = true;
    }
    else {
      /*
      std::cout << "capid_shift " << capid_shift 
        << " fiberOffset " << digi.fiberIdleOffset() << std::endl;
      */
      rechit = T(rechit.id(), rechit.energy(), rechit.time() +
         ( capid_shift == 1 ? -25 : +25 )
         );
      rechit.setFlagField(1,(capid_shift == 1 ? 
                 HcalCaloFlagLabels::TimingSubtractedBit : 
                 HcalCaloFlagLabels::TimingAddedBit ));
    }

    if (bad) {
      /*
    std::cout << "BAD capid_shift " << capid_shift 
    << " fiberOffset " << digi.fiberIdleOffset() << std::endl;
      */      
      rechit.setFlagField(1,HcalCaloFlagLabels::TimingErrorBit);
    }
  }