/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/EventFilter/CastorRawToDigi/src/CastorDataFrameFilter.cc

Go to the documentation of this file.

#include "EventFilter/CastorRawToDigi/interface/CastorDataFrameFilter.h"

namespace CastorDataFrameFilter_impl {

  template <class DataFrame> 
  bool check(const DataFrame& df, bool capcheck, bool dvercheck) {
    if (capcheck || dvercheck) {
      int lastcapid=0, capid=0;
      for (int i=0; i<df.size(); i++) {
        capid=df[i].capid();
        if (capcheck && i!=0 && ((lastcapid+1)%4)!=capid) 
          return false;
        if (dvercheck && ( df[i].er() || !df[i].dv() ))
          return false;
        lastcapid=capid;      
      }
    }
    return true;
  }

  template <class DataFrame> 
  double energySum(const DataFrame& df, int fs, int ls) {
    double es=0;
    for (int i=fs; i<=ls && i<=df.size(); i++) 
      es+=df[i].nominal_fC();
    return es;
  }

}


CastorDataFrameFilter::CastorDataFrameFilter(bool requireCapid, bool requireDVER, bool energyFilter, int firstSample, int lastSample, double minAmpl) :
  requireCapid_(requireCapid), requireDVER_(requireDVER), energyFilter_(energyFilter),
  firstSample_(firstSample), lastSample_(lastSample), minimumAmplitude_(minAmpl) {
}

CastorDigiCollection CastorDataFrameFilter::filter(const CastorDigiCollection& incol, HcalUnpackerReport& r) {
  CastorDigiCollection output;
  for (CastorDigiCollection::const_iterator i=incol.begin(); i!=incol.end(); i++) {
    if (!CastorDataFrameFilter_impl::check(*i,requireCapid_,requireDVER_)) 
      r.countBadQualityDigi();
    else if (!energyFilter_ || minimumAmplitude_<CastorDataFrameFilter_impl::energySum(*i,firstSample_,lastSample_))
      output.push_back(*i);
  }
  return output;
}

bool CastorDataFrameFilter::active() const {
  return requireCapid_|requireDVER_|energyFilter_;
}