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);
}