GenFilterEfficiencyProducer.cc

Go to the documentation of this file.

#include "GeneratorInterface/Core/interface/GenFilterEfficiencyProducer.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

GenFilterEfficiencyProducer::GenFilterEfficiencyProducer(const edm::ParameterSet& iConfig) :
  filterPath(iConfig.getParameter<std::string>("filterPath")),
  tns_(),
  thisProcess(),pathIndex(100000),
  numEventsPassPos_(0),
  numEventsPassNeg_(0),
  numEventsTotalPos_(0),
  numEventsTotalNeg_(0),
  sumpass_w_(0.),
  sumpass_w2_(0.),
  sumtotal_w_(0.),
  sumtotal_w2_(0.)
{
   //now do what ever initialization is needed
  if (edm::Service<edm::service::TriggerNamesService>().isAvailable()) {
    // get tns pointer
    tns_ = edm::Service<edm::service::TriggerNamesService>().operator->();
    if (tns_!=0) {
      thisProcess = tns_->getProcessName();
      std::vector<std::string> theNames = tns_->getTrigPaths();
      for ( unsigned int i = 0; i < theNames.size(); i++ ) {
        if ( theNames[i] == filterPath ) { pathIndex = i; continue; }
      }
    }
    else
      edm::LogError("ServiceNotAvailable") << "TriggerNamesServive not available, no filter information stored";
  }

  produces<GenFilterInfo, edm::InLumi>(); 

}


GenFilterEfficiencyProducer::~GenFilterEfficiencyProducer()
{

} 


//
// member functions
//

// ------------ method called to for each event  ------------
void
GenFilterEfficiencyProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{

  edm::InputTag theTrig("TriggerResults","",thisProcess);
  edm::Handle<edm::TriggerResults> trigR;
  iEvent.getByLabel(theTrig,trigR); 
  edm::Handle<GenEventInfoProduct>    genEventScale;
  if (!iEvent.getByLabel("generator", genEventScale))return;
  double weight = genEventScale->weight();
  
  unsigned int nSize = (*trigR).size();
  // std::cout << "Number of paths in TriggerResults = " << nSize  << std::endl;
  if ( nSize >= pathIndex ) {

    if (!trigR->wasrun(pathIndex))return;
    if ( trigR->accept(pathIndex) ) { 
      sumpass_w_ += weight;
      sumpass_w2_+= weight*weight;
      
      sumtotal_w_ += weight;
      sumtotal_w2_+= weight*weight;

      if(weight > 0)
        {
          numEventsPassPos_++;
          numEventsTotalPos_++;
        }
      else
        {
          numEventsPassNeg_++;
          numEventsTotalNeg_++;
        }

    }
    else // if fail the filter
      {
        sumtotal_w_ += weight;
        sumtotal_w2_+= weight*weight;

        if(weight > 0)
          numEventsTotalPos_++;
        else
          numEventsTotalNeg_++;
      }
    //    std::cout << "Total events = " << numEventsTotal << " passed = " << numEventsPassed << std::endl;
  }

}


void
GenFilterEfficiencyProducer::beginLuminosityBlock(edm::LuminosityBlock & iLumi, const edm::EventSetup&) {
 
  numEventsPassPos_=0;
  numEventsPassNeg_=0;
  numEventsTotalPos_=0;
  numEventsTotalNeg_=0;
  sumpass_w_=0;
  sumpass_w2_=0;
  sumtotal_w_=0;
  sumtotal_w2_=0;
  
}

void
GenFilterEfficiencyProducer::endLuminosityBlock(edm::LuminosityBlock & iLumi, const edm::EventSetup&) {

  std::auto_ptr<GenFilterInfo> thisProduct(new GenFilterInfo(
                                                             numEventsPassPos_,
                                                             numEventsPassNeg_,
                                                             numEventsTotalPos_,
                                                             numEventsTotalNeg_,
                                                             sumpass_w_,
                                                             sumpass_w2_,
                                                             sumtotal_w_,
                                                             sumtotal_w2_
                                                             ));
  iLumi.put(thisProduct);
  
}