TtFullHadSignalSelMVAComputer.cc

Go to the documentation of this file.

#include "PhysicsTools/JetMCUtils/interface/combination.h"
 
#include "TopQuarkAnalysis/TopEventSelection/plugins/TtFullHadSignalSelMVAComputer.h"
#include "TopQuarkAnalysis/TopEventSelection/interface/TtFullHadSignalSelEval.h"
 
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
 
#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Framework/interface/TriggerNamesService.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DataFormats/PatCandidates/interface/Flags.h"
 
TtFullHadSignalSelMVAComputer::TtFullHadSignalSelMVAComputer(const edm::ParameterSet& cfg)
    : jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))) {
  produces<double>("DiscSel");
}
 
TtFullHadSignalSelMVAComputer::~TtFullHadSignalSelMVAComputer() {}
 
void TtFullHadSignalSelMVAComputer::produce(edm::Event& evt, const edm::EventSetup& setup) {
  std::unique_ptr<double> pOutDisc(new double);
 
  mvaComputer.update<TtFullHadSignalSelMVARcd>(setup, "ttFullHadSignalSelMVA");
 
  // read name of the last processor in the MVA calibration
  // (to be used as meta information)
  edm::ESHandle<PhysicsTools::Calibration::MVAComputerContainer> calibContainer;
  setup.get<TtFullHadSignalSelMVARcd>().get(calibContainer);
  std::vector<PhysicsTools::Calibration::VarProcessor*> processors =
      (calibContainer->find("ttFullHadSignalSelMVA")).getProcessors();
 
  edm::Handle<std::vector<pat::Jet> > jets;
  evt.getByToken(jetsToken_, jets);
 
  //calculation of InputVariables
  //see TopQuarkAnalysis/TopTools/interface/TtFullHadSignalSel.h
  //                             /src/TtFullHadSignalSel.cc
  //all objects, jets, which are needed for the calculation
  //of the input-variables have to be passed to this class
  TtFullHadSignalSel selection(*jets);
 
  double discrim = evaluateTtFullHadSignalSel(mvaComputer, selection);
 
  *pOutDisc = discrim;
 
  evt.put(std::move(pOutDisc), "DiscSel");
 
  DiscSel = discrim;
}
 
void TtFullHadSignalSelMVAComputer::beginJob() {}
 
void TtFullHadSignalSelMVAComputer::endJob() {}
 
// implement the plugins for the computer container
// -> register TtFullHadSignalSelMVARcd
// -> define TtFullHadSignalSelMVAFileSource
MVA_COMPUTER_CONTAINER_IMPLEMENT(TtFullHadSignalSelMVA);