TtDilepEvtSolutionMaker.h

Go to the documentation of this file.

#include <string>
#include <vector>

#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "DataFormats/Candidate/interface/Candidate.h"
#include "TopQuarkAnalysis/TopKinFitter/interface/TtFullLepKinSolver.h"

class TtDilepLRSignalSelObservables;

class TtDilepEvtSolutionMaker : public edm::EDProducer {
public:
  explicit TtDilepEvtSolutionMaker(const edm::ParameterSet& iConfig);
  ~TtDilepEvtSolutionMaker() override;

  void beginJob() override;
  void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override;

private:
  // next methods are avoidable but they make the code legible
  inline bool PTComp(const reco::Candidate*, const reco::Candidate*) const;
  inline bool LepDiffCharge(const reco::Candidate*, const reco::Candidate*) const;
  inline bool HasPositiveCharge(const reco::Candidate*) const;

private:
  edm::EDGetTokenT<std::vector<pat::Electron> > electronSourceToken_;
  edm::EDGetTokenT<std::vector<pat::Muon> > muonSourceToken_;
  edm::EDGetTokenT<std::vector<pat::Tau> > tauSourceToken_;
  edm::EDGetTokenT<std::vector<pat::MET> > metSourceToken_;
  edm::EDGetTokenT<std::vector<pat::Jet> > jetSourceToken_;
  edm::EDGetTokenT<TtGenEvent> evtSourceToken_;
  int jetCorrScheme_;
  unsigned int nrCombJets_;
  bool matchToGenEvt_, calcTopMass_, useMCforBest_;
  bool eeChannel_, emuChannel_, mumuChannel_, etauChannel_, mutauChannel_, tautauChannel_;
  double tmassbegin_, tmassend_, tmassstep_;
  std::vector<double> nupars_;

  TtDilepLRSignalSelObservables* myLRSignalSelObservables;
  TtFullLepKinSolver* solver;
};

inline bool TtDilepEvtSolutionMaker::PTComp(const reco::Candidate* l1, const reco::Candidate* l2) const {
  return (l1->pt() > l2->pt());
}

inline bool TtDilepEvtSolutionMaker::LepDiffCharge(const reco::Candidate* l1, const reco::Candidate* l2) const {
  return (l1->charge() != l2->charge());
}

inline bool TtDilepEvtSolutionMaker::HasPositiveCharge(const reco::Candidate* l) const { return (l->charge() > 0); }