ElectronPATIdMVAProducer.cc

Go to the documentation of this file.

// system include files
#include <memory>
 
// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"
 
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "EgammaAnalysis/ElectronTools/interface/EGammaMvaEleEstimator.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
#include "DataFormats/PatCandidates/interface/Electron.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "TrackingTools/Records/interface/TransientTrackRecord.h"
#include "TrackingTools/IPTools/interface/IPTools.h"
//
// class declaration
//
 
class ElectronPATIdMVAProducer : public edm::EDProducer {
public:
  explicit ElectronPATIdMVAProducer(const edm::ParameterSet&);
  ~ElectronPATIdMVAProducer() override;
 
private:
  void produce(edm::Event&, const edm::EventSetup&) override;
 
  // ----------member data ---------------------------
  bool verbose_;
  edm::EDGetTokenT<pat::ElectronCollection> electronToken_;
  double _Rho;
  edm::EDGetTokenT<double> eventrhoToken_;
  std::string method_;
  std::vector<std::string> mvaWeightFiles_;
 
  EGammaMvaEleEstimator* mvaID_;
};
 
//
// constants, enums and typedefs
//
 
//
// static data member definitions
//
 
//
// constructors and destructor
//
ElectronPATIdMVAProducer::ElectronPATIdMVAProducer(const edm::ParameterSet& iConfig) {
  verbose_ = iConfig.getUntrackedParameter<bool>("verbose", false);
  electronToken_ = consumes<pat::ElectronCollection>(iConfig.getParameter<edm::InputTag>("electronTag"));
  method_ = iConfig.getParameter<std::string>("method");
  std::vector<std::string> fpMvaWeightFiles = iConfig.getParameter<std::vector<std::string> >("mvaWeightFile");
  eventrhoToken_ = consumes<double>(iConfig.getParameter<edm::InputTag>("Rho"));
 
  produces<edm::ValueMap<float> >();
 
  mvaID_ = new EGammaMvaEleEstimator();
 
  EGammaMvaEleEstimator::MVAType type_;
 
  type_ = EGammaMvaEleEstimator::kTrigNoIP;
 
  bool manualCat_ = true;
 
  std::string path_mvaWeightFileEleID;
  for (unsigned ifile = 0; ifile < fpMvaWeightFiles.size(); ++ifile) {
    path_mvaWeightFileEleID = edm::FileInPath(fpMvaWeightFiles[ifile].c_str()).fullPath();
    mvaWeightFiles_.push_back(path_mvaWeightFileEleID);
  }
 
  mvaID_->initialize(method_, type_, manualCat_, mvaWeightFiles_);
}
 
ElectronPATIdMVAProducer::~ElectronPATIdMVAProducer() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}
 
//
// member functions
//
 
// ------------ method called on each new Event  ------------
void ElectronPATIdMVAProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  std::unique_ptr<edm::ValueMap<float> > out(new edm::ValueMap<float>());
 
  edm::Handle<pat::ElectronCollection> egCollection;
  iEvent.getByToken(electronToken_, egCollection);
  const pat::ElectronCollection egCandidates = (*egCollection.product());
 
  _Rho = 0;
  edm::Handle<double> rhoPtr;
  //const edm::InputTag eventrhoToken_("kt6PFJets", "rho");
  iEvent.getByToken(eventrhoToken_, rhoPtr);
  _Rho = *rhoPtr;
 
  std::vector<float> values;
  values.reserve(egCollection->size());
 
  for (pat::ElectronCollection::const_iterator egIter = egCandidates.begin(); egIter != egCandidates.end(); ++egIter) {
    double mvaVal = -999999;
 
    mvaVal = mvaID_->mvaValue(*egIter, _Rho, verbose_);
 
    values.push_back(mvaVal);
  }
 
  edm::ValueMap<float>::Filler filler(*out);
  filler.insert(egCollection, values.begin(), values.end());
  filler.fill();
 
  iEvent.put(std::move(out));
}
 
//define this as a plug-in
DEFINE_FWK_MODULE(ElectronPATIdMVAProducer);