MTCCHLTrigger.cc

Go to the documentation of this file.

#include "EventFilter/SiStripChannelChargeFilter/interface/MTCCHLTrigger.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
#include "DataFormats/SiStripDigi/interface/SiStripDigi.h"
#include "DataFormats/Common/interface/DetSetVector.h"

namespace cms {
  MTCCHLTrigger::MTCCHLTrigger(const edm::ParameterSet& ps) {
    selOnDigiCharge = ps.getParameter<bool>("SelOnDigiCharge");
    ChargeThreshold = ps.getParameter<int>("ChargeThreshold");
    clusterProducer = ps.getParameter<std::string>("ClusterProducer");
    produces<int>();
    produces<unsigned int>();
  }

  bool MTCCHLTrigger::filter(edm::Event& e, edm::EventSetup const& c) {
    //get data
    //StripCluster
    edm::Handle<edm::DetSetVector<SiStripCluster> > h;
    e.getByLabel(clusterProducer, h);

    //StripDigi from RawToDigi and ZeroSuppressor
    std::vector<edm::Handle<edm::DetSetVector<SiStripDigi> > > di;
    e.getManyByType(di);

    if (selOnDigiCharge) {
      unsigned int digiadc = 0;
      for (std::vector<edm::Handle<edm::DetSetVector<SiStripDigi> > >::const_iterator mi = di.begin(); mi != di.end();
           mi++) {
        for (edm::DetSetVector<SiStripDigi>::const_iterator it = (*mi)->begin(); it != (*mi)->end(); it++) {
          for (std::vector<SiStripDigi>::const_iterator vit = (it->data).begin(); vit != (it->data).end(); vit++)
            digiadc += vit->adc();
        }
      }
      return (digiadc > ChargeThreshold) ? true : false;
    } else {
      unsigned int amplclus = 0;
      for (edm::DetSetVector<SiStripCluster>::const_iterator it = h->begin(); it != h->end(); it++) {
        for (std::vector<SiStripCluster>::const_iterator vit = (it->data).begin(); vit != (it->data).end(); vit++) {
          for (auto ia = vit->amplitudes().begin(); ia != vit->amplitudes().end(); ia++) {
            if ((*ia) > 0) {
              amplclus += (*ia);
            }
          }
        }
      }
      bool decision = (amplclus > ChargeThreshold) ? true : false;
      e.put(std::make_unique<unsigned int>(amplclus));
      e.put(std::make_unique<int>(decision));
      return decision;
    }
  }
}  // namespace cms