TtSemiLepHitFitProducer< LeptonCollection > Class Template Reference

#include <TtSemiLepHitFitProducer.h>

Inheritance diagram for TtSemiLepHitFitProducer< LeptonCollection >:

Classes
struct	FitResult

Public Member Functions
	TtSemiLepHitFitProducer (const edm::ParameterSet &)
	~TtSemiLepHitFitProducer () override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef hitfit::RunHitFit < pat::Electron, pat::Muon, pat::Jet, pat::MET >	PatHitFit

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	bTagAlgo_
	input tag for b-tagging algorithm More...
hitfit::LeptonTranslatorBase < pat::Electron >	electronTranslator_
PatHitFit *	HitFit
edm::FileInPath	hitfitBJetResolution_
edm::FileInPath	hitfitDefault_
edm::FileInPath	hitfitElectronResolution_
edm::FileInPath	hitfitMETResolution_
edm::FileInPath	hitfitMuonResolution_
edm::FileInPath	hitfitUdscJetResolution_
double	jes_
	jet energy scale More...
double	jesB_
std::string	jetCorrectionLevel_
	jet correction level More...
edm::EDGetTokenT< std::vector < pat::Jet > >	jetsToken_
hitfit::JetTranslatorBase < pat::Jet >	jetTranslator_
edm::EDGetTokenT < LeptonCollection >	lepsToken_
double	maxBTagValueNonBJet_
	max value of bTag for a non-b-jet More...
double	maxEtaEle_
	maximum eta value for electrons, needed to limited range in which resolutions are provided More...
double	maxEtaJet_
	maximum eta value for jets, needed to limited range in which resolutions are provided More...
double	maxEtaMu_
	maximum eta value for muons, needed to limited range in which resolutions are provided More...
int	maxNComb_
	maximal number of combinations to be written to the event More...
int	maxNJets_
	maximal number of jets (-1 possible to indicate 'all') More...
edm::EDGetTokenT< std::vector < pat::MET > >	metsToken_
hitfit::METTranslatorBase < pat::MET >	metTranslator_
double	minBTagValueBJet_
	min value of bTag for a b-jet More...
double	mTop_
hitfit::LeptonTranslatorBase < pat::Muon >	muonTranslator_
double	mW_
	constraints More...
bool	useBTag_
	switch to tell whether to use b-tagging or not More...

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

template<typename LeptonCollection>
class TtSemiLepHitFitProducer< LeptonCollection >

Definition at line 20 of file TtSemiLepHitFitProducer.h.

Member Typedef Documentation

template<typename LeptonCollection >

typedef hitfit::RunHitFit<pat::Electron, pat::Muon, pat::Jet, pat::MET> TtSemiLepHitFitProducer< LeptonCollection >::PatHitFit

private

Definition at line 81 of file TtSemiLepHitFitProducer.h.

Constructor & Destructor Documentation

template<typename LeptonCollection >

TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer ( const edm::ParameterSet &  cfg )

explicit

Definition at line 99 of file TtSemiLepHitFitProducer.h.

References TtSemiLepHitFitProducer< LeptonCollection >::electronTranslator_, edm::FileInPath::fullPath(), TtSemiLepHitFitProducer< LeptonCollection >::HitFit, TtSemiLepHitFitProducer< LeptonCollection >::hitfitDefault_, TtSemiLepHitFitProducer< LeptonCollection >::jetTranslator_, TtSemiLepHitFitProducer< LeptonCollection >::maxEtaEle_, TtSemiLepHitFitProducer< LeptonCollection >::maxEtaJet_, TtSemiLepHitFitProducer< LeptonCollection >::maxEtaMu_, TtSemiLepHitFitProducer< LeptonCollection >::metTranslator_, TtSemiLepHitFitProducer< LeptonCollection >::mTop_, TtSemiLepHitFitProducer< LeptonCollection >::muonTranslator_, and TtSemiLepHitFitProducer< LeptonCollection >::mW_.

    : jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
      lepsToken_(consumes<LeptonCollection>(cfg.getParameter<edm::InputTag>("leps"))),
      metsToken_(consumes<std::vector<pat::MET> >(cfg.getParameter<edm::InputTag>("mets"))),
      maxNJets_(cfg.getParameter<int>("maxNJets")),
      maxNComb_(cfg.getParameter<int>("maxNComb")),
      bTagAlgo_(cfg.getParameter<std::string>("bTagAlgo")),
      minBTagValueBJet_(cfg.getParameter<double>("minBDiscBJets")),
      maxBTagValueNonBJet_(cfg.getParameter<double>("maxBDiscLightJets")),
      useBTag_(cfg.getParameter<bool>("useBTagging")),
      mW_(cfg.getParameter<double>("mW")),
      mTop_(cfg.getParameter<double>("mTop")),
      jetCorrectionLevel_(cfg.getParameter<std::string>("jetCorrectionLevel")),
      jes_(cfg.getParameter<double>("jes")),
      jesB_(cfg.getParameter<double>("jesB")),

      // The following five initializers read the config parameters for the
      // ASCII text files which contains the physics object resolutions.
      hitfitDefault_(cfg.getUntrackedParameter<edm::FileInPath>(
          std::string("hitfitDefault"),
          edm::FileInPath(std::string("TopQuarkAnalysis/TopHitFit/data/setting/RunHitFitConfiguration.txt")))),
      hitfitElectronResolution_(cfg.getUntrackedParameter<edm::FileInPath>(
          std::string("hitfitElectronResolution"),
          edm::FileInPath(std::string("TopQuarkAnalysis/TopHitFit/data/resolution/tqafElectronResolution.txt")))),
      hitfitMuonResolution_(cfg.getUntrackedParameter<edm::FileInPath>(
          std::string("hitfitMuonResolution"),
          edm::FileInPath(std::string("TopQuarkAnalysis/TopHitFit/data/resolution/tqafMuonResolution.txt")))),
      hitfitUdscJetResolution_(cfg.getUntrackedParameter<edm::FileInPath>(
          std::string("hitfitUdscJetResolution"),
          edm::FileInPath(std::string("TopQuarkAnalysis/TopHitFit/data/resolution/tqafUdscJetResolution.txt")))),
      hitfitBJetResolution_(cfg.getUntrackedParameter<edm::FileInPath>(
          std::string("hitfitBJetResolution"),
          edm::FileInPath(std::string("TopQuarkAnalysis/TopHitFit/data/resolution/tqafBJetResolution.txt")))),
      hitfitMETResolution_(cfg.getUntrackedParameter<edm::FileInPath>(
          std::string("hitfitMETResolution"),
          edm::FileInPath(std::string("TopQuarkAnalysis/TopHitFit/data/resolution/tqafKtResolution.txt")))),

      // The following four initializers instantiate the translator between PAT objects
      // and HitFit objects using the ASCII text files which contains the resolutions.
      electronTranslator_(hitfitElectronResolution_.fullPath()),
      muonTranslator_(hitfitMuonResolution_.fullPath()),
      jetTranslator_(
          hitfitUdscJetResolution_.fullPath(), hitfitBJetResolution_.fullPath(), jetCorrectionLevel_, jes_, jesB_),
      metTranslator_(hitfitMETResolution_.fullPath())

{
  // Create an instance of RunHitFit and initialize it.
  HitFit = new PatHitFit(
      electronTranslator_, muonTranslator_, jetTranslator_, metTranslator_, hitfitDefault_.fullPath(), mW_, mW_, mTop_);

  maxEtaMu_ = 2.4;
  maxEtaEle_ = 2.5;
  maxEtaJet_ = 2.5;

  edm::LogVerbatim("TopHitFit") << "\n"
                                << "+++++++++++ TtSemiLepHitFitProducer ++++++++++++ \n"
                                << " Due to the eta ranges for which resolutions     \n"
                                << " are provided in                                 \n"
                                << " TopQuarkAnalysis/TopHitFit/data/resolution/     \n"
                                << " so far, the following cuts are currently        \n"
                                << " implemented in the TtSemiLepHitFitProducer:     \n"
                                << " |eta(muons    )| <= " << maxEtaMu_ << " \n"
                                << " |eta(electrons)| <= " << maxEtaEle_ << " \n"
                                << " |eta(jets     )| <= " << maxEtaJet_ << " \n"
                                << "++++++++++++++++++++++++++++++++++++++++++++++++ \n";

  produces<std::vector<pat::Particle> >("PartonsHadP");
  produces<std::vector<pat::Particle> >("PartonsHadQ");
  produces<std::vector<pat::Particle> >("PartonsHadB");
  produces<std::vector<pat::Particle> >("PartonsLepB");
  produces<std::vector<pat::Particle> >("Leptons");
  produces<std::vector<pat::Particle> >("Neutrinos");

  produces<std::vector<std::vector<int> > >();
  produces<std::vector<double> >("Chi2");
  produces<std::vector<double> >("Prob");
  produces<std::vector<double> >("MT");
  produces<std::vector<double> >("SigMT");
  produces<std::vector<int> >("Status");
  produces<int>("NumberOfConsideredJets");
}

template<typename LeptonCollection >

TtSemiLepHitFitProducer< LeptonCollection >::~TtSemiLepHitFitProducer ( )

override

Definition at line 182 of file TtSemiLepHitFitProducer.h.

                                                                    {
  delete HitFit;
}

Member Function Documentation

template<typename LeptonCollection >

void TtSemiLepHitFitProducer< LeptonCollection >::produce ( edm::Event &  evt, const edm::EventSetup &  setup  )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 187 of file TtSemiLepHitFitProducer.h.

References funct::abs(), TtSemiLepHitFitProducer< LeptonCollection >::FitResult::Chi2, PVValHelper::eta, funct::exp(), gpuVertexFinder::fit, edm::Event::getByToken(), TtSemiLepEvtPartons::HadB, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::HadB, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::HadP, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::HadQ, mps_fire::i, metsig::jet, hitfit::Lepjets_Event::jet(), TtSemiLepHitFitProducer< LeptonCollection >::FitResult::JetCombi, fwrapper::jets, hitfit::Lepjets_Event::lep(), TtSemiLepEvtPartons::LepB, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::LepB, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::LepL, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::LepN, TtSemiLepEvtPartons::LightQ, TtSemiLepEvtPartons::LightQBar, hitfit::Lepjets_Event::met(), eostools::move(), TtSemiLepHitFitProducer< LeptonCollection >::FitResult::MT, nPartons, hitfit::Lepjets_Event_Lep::p(), TtSemiLepHitFitProducer< LeptonCollection >::FitResult::Prob, edm::Event::put(), mps_fire::result, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::SigMT, TtSemiLepHitFitProducer< LeptonCollection >::FitResult::Status, hitfit::Lepjets_Event_Lep::type(), and trackerHitRTTI::vector.

                                                                                                 {
  std::unique_ptr<std::vector<pat::Particle> > pPartonsHadP(new std::vector<pat::Particle>);
  std::unique_ptr<std::vector<pat::Particle> > pPartonsHadQ(new std::vector<pat::Particle>);
  std::unique_ptr<std::vector<pat::Particle> > pPartonsHadB(new std::vector<pat::Particle>);
  std::unique_ptr<std::vector<pat::Particle> > pPartonsLepB(new std::vector<pat::Particle>);
  std::unique_ptr<std::vector<pat::Particle> > pLeptons(new std::vector<pat::Particle>);
  std::unique_ptr<std::vector<pat::Particle> > pNeutrinos(new std::vector<pat::Particle>);

  std::unique_ptr<std::vector<std::vector<int> > > pCombi(new std::vector<std::vector<int> >);
  std::unique_ptr<std::vector<double> > pChi2(new std::vector<double>);
  std::unique_ptr<std::vector<double> > pProb(new std::vector<double>);
  std::unique_ptr<std::vector<double> > pMT(new std::vector<double>);
  std::unique_ptr<std::vector<double> > pSigMT(new std::vector<double>);
  std::unique_ptr<std::vector<int> > pStatus(new std::vector<int>);
  std::unique_ptr<int> pJetsConsidered(new int);

  edm::Handle<std::vector<pat::Jet> > jets;
  evt.getByToken(jetsToken_, jets);

  edm::Handle<std::vector<pat::MET> > mets;
  evt.getByToken(metsToken_, mets);

  edm::Handle<LeptonCollection> leps;
  evt.getByToken(lepsToken_, leps);

  // -----------------------------------------------------
  // skip events with no appropriate lepton candidate in
  // or empty MET or less jets than partons
  // -----------------------------------------------------

  const unsigned int nPartons = 4;

  // Clear the internal state
  HitFit->clear();

  // Add lepton into HitFit
  bool foundLepton = false;
  if (!leps->empty()) {
    double maxEtaLep = maxEtaMu_;
    if (!dynamic_cast<const reco::Muon*>(&((*leps)[0])))  // assume electron if it is not a muon
      maxEtaLep = maxEtaEle_;
    for (unsigned iLep = 0; iLep < (*leps).size() && !foundLepton; ++iLep) {
      if (std::abs((*leps)[iLep].eta()) <= maxEtaLep) {
        HitFit->AddLepton((*leps)[iLep]);
        foundLepton = true;
      }
    }
  }

  // Add jets into HitFit
  unsigned int nJetsFound = 0;
  for (unsigned iJet = 0; iJet < (*jets).size() && (int)nJetsFound != maxNJets_; ++iJet) {
    double jet_eta = (*jets)[iJet].eta();
    if ((*jets)[iJet].isCaloJet()) {
      jet_eta = ((reco::CaloJet*)((*jets)[iJet]).originalObject())->detectorP4().eta();
    }
    if (std::abs(jet_eta) <= maxEtaJet_) {
      HitFit->AddJet((*jets)[iJet]);
      nJetsFound++;
    }
  }
  *pJetsConsidered = nJetsFound;

  // Add missing transverse energy into HitFit
  if (!mets->empty())
    HitFit->SetMet((*mets)[0]);

  if (!foundLepton || mets->empty() || nJetsFound < nPartons) {
    // the kinFit getters return empty objects here
    pPartonsHadP->push_back(pat::Particle());
    pPartonsHadQ->push_back(pat::Particle());
    pPartonsHadB->push_back(pat::Particle());
    pPartonsLepB->push_back(pat::Particle());
    pLeptons->push_back(pat::Particle());
    pNeutrinos->push_back(pat::Particle());
    // indices referring to the jet combination
    std::vector<int> invalidCombi;
    for (unsigned int i = 0; i < nPartons; ++i)
      invalidCombi.push_back(-1);
    pCombi->push_back(invalidCombi);
    // chi2
    pChi2->push_back(-1.);
    // chi2 probability
    pProb->push_back(-1.);
    // fitted top mass
    pMT->push_back(-1.);
    pSigMT->push_back(-1.);
    // status of the fitter
    pStatus->push_back(-1);
    // feed out all products
    evt.put(std::move(pCombi));
    evt.put(std::move(pPartonsHadP), "PartonsHadP");
    evt.put(std::move(pPartonsHadQ), "PartonsHadQ");
    evt.put(std::move(pPartonsHadB), "PartonsHadB");
    evt.put(std::move(pPartonsLepB), "PartonsLepB");
    evt.put(std::move(pLeptons), "Leptons");
    evt.put(std::move(pNeutrinos), "Neutrinos");
    evt.put(std::move(pChi2), "Chi2");
    evt.put(std::move(pProb), "Prob");
    evt.put(std::move(pMT), "MT");
    evt.put(std::move(pSigMT), "SigMT");
    evt.put(std::move(pStatus), "Status");
    evt.put(std::move(pJetsConsidered), "NumberOfConsideredJets");
    return;
  }

  std::list<FitResult> FitResultList;

  //
  // BEGIN DECLARATION OF VARIABLES FROM KINEMATIC FIT
  //

  // In this part are variables from the
  // kinematic fit procedure

  // Number of all permutations of the event
  size_t nHitFit = 0;

  // Number of jets in the event
  size_t nHitFitJet = 0;

  // Results of the fit for all jet permutations of the event
  std::vector<hitfit::Fit_Result> hitFitResult;

  //
  // R U N   H I T F I T
  //
  // Run the kinematic fit and get how many permutations are possible
  // in the fit

  nHitFit = HitFit->FitAllPermutation();

  //
  // BEGIN PART WHICH EXTRACTS INFORMATION FROM HITFIT
  //

  // Get the number of jets
  nHitFitJet = HitFit->GetUnfittedEvent()[0].njets();

  // Get the fit results for all permutations
  hitFitResult = HitFit->GetFitAllPermutation();

  // Loop over all permutations and extract the information
  for (size_t fit = 0; fit != nHitFit; ++fit) {
    // Get the event after the fit
    hitfit::Lepjets_Event fittedEvent = hitFitResult[fit].ev();

    /*
        Get jet permutation according to TQAF convention
        11 : leptonic b
        12 : hadronic b
        13 : hadronic W
        14 : hadronic W
      */
    std::vector<int> hitCombi(4);
    for (size_t jet = 0; jet != nHitFitJet; ++jet) {
      int jet_type = fittedEvent.jet(jet).type();

      switch (jet_type) {
        case 11:
          hitCombi[TtSemiLepEvtPartons::LepB] = jet;
          break;
        case 12:
          hitCombi[TtSemiLepEvtPartons::HadB] = jet;
          break;
        case 13:
          hitCombi[TtSemiLepEvtPartons::LightQ] = jet;
          break;
        case 14:
          hitCombi[TtSemiLepEvtPartons::LightQBar] = jet;
          break;
      }
    }

    // Store the kinematic quantities in the corresponding containers.

    hitfit::Lepjets_Event_Jet hadP_ = fittedEvent.jet(hitCombi[TtSemiLepEvtPartons::LightQ]);
    hitfit::Lepjets_Event_Jet hadQ_ = fittedEvent.jet(hitCombi[TtSemiLepEvtPartons::LightQBar]);
    hitfit::Lepjets_Event_Jet hadB_ = fittedEvent.jet(hitCombi[TtSemiLepEvtPartons::HadB]);
    hitfit::Lepjets_Event_Jet lepB_ = fittedEvent.jet(hitCombi[TtSemiLepEvtPartons::LepB]);
    hitfit::Lepjets_Event_Lep lepL_ = fittedEvent.lep(0);

    /*
  std::string bTagAlgo_;
  double minBTagValueBJet_;
  double maxBTagValueNonBJet_;
  bool useBTag_;
      */

    if (hitFitResult[fit].chisq() > 0  // only take into account converged fits
        && (!useBTag_ ||
            (useBTag_  // use btag information if chosen
             && jets->at(hitCombi[TtSemiLepEvtPartons::LightQ]).bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
             jets->at(hitCombi[TtSemiLepEvtPartons::LightQBar]).bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
             jets->at(hitCombi[TtSemiLepEvtPartons::HadB]).bDiscriminator(bTagAlgo_) > minBTagValueBJet_ &&
             jets->at(hitCombi[TtSemiLepEvtPartons::LepB]).bDiscriminator(bTagAlgo_) > minBTagValueBJet_))) {
      FitResult result;
      result.Status = 0;
      result.Chi2 = hitFitResult[fit].chisq();
      result.Prob = exp(-1.0 * (hitFitResult[fit].chisq()) / 2.0);
      result.MT = hitFitResult[fit].mt();
      result.SigMT = hitFitResult[fit].sigmt();
      result.HadB = pat::Particle(reco::LeafCandidate(
          0, math::XYZTLorentzVector(hadB_.p().x(), hadB_.p().y(), hadB_.p().z(), hadB_.p().t()), math::XYZPoint()));
      result.HadP = pat::Particle(reco::LeafCandidate(
          0, math::XYZTLorentzVector(hadP_.p().x(), hadP_.p().y(), hadP_.p().z(), hadP_.p().t()), math::XYZPoint()));
      result.HadQ = pat::Particle(reco::LeafCandidate(
          0, math::XYZTLorentzVector(hadQ_.p().x(), hadQ_.p().y(), hadQ_.p().z(), hadQ_.p().t()), math::XYZPoint()));
      result.LepB = pat::Particle(reco::LeafCandidate(
          0, math::XYZTLorentzVector(lepB_.p().x(), lepB_.p().y(), lepB_.p().z(), lepB_.p().t()), math::XYZPoint()));
      result.LepL = pat::Particle(reco::LeafCandidate(
          0, math::XYZTLorentzVector(lepL_.p().x(), lepL_.p().y(), lepL_.p().z(), lepL_.p().t()), math::XYZPoint()));
      result.LepN = pat::Particle(reco::LeafCandidate(
          0,
          math::XYZTLorentzVector(
              fittedEvent.met().x(), fittedEvent.met().y(), fittedEvent.met().z(), fittedEvent.met().t()),
          math::XYZPoint()));
      result.JetCombi = hitCombi;

      FitResultList.push_back(result);
    }
  }

  // sort results w.r.t. chi2 values
  FitResultList.sort();

  // -----------------------------------------------------
  // feed out result
  // starting with the JetComb having the smallest chi2
  // -----------------------------------------------------

  if (((unsigned)FitResultList.size()) < 1) {  // in case no fit results were stored in the list (all fits aborted)
    pPartonsHadP->push_back(pat::Particle());
    pPartonsHadQ->push_back(pat::Particle());
    pPartonsHadB->push_back(pat::Particle());
    pPartonsLepB->push_back(pat::Particle());
    pLeptons->push_back(pat::Particle());
    pNeutrinos->push_back(pat::Particle());
    // indices referring to the jet combination
    std::vector<int> invalidCombi;
    for (unsigned int i = 0; i < nPartons; ++i)
      invalidCombi.push_back(-1);
    pCombi->push_back(invalidCombi);
    // chi2
    pChi2->push_back(-1.);
    // chi2 probability
    pProb->push_back(-1.);
    // fitted top mass
    pMT->push_back(-1.);
    pSigMT->push_back(-1.);
    // status of the fitter
    pStatus->push_back(-1);
  } else {
    unsigned int iComb = 0;
    for (typename std::list<FitResult>::const_iterator result = FitResultList.begin(); result != FitResultList.end();
         ++result) {
      if (maxNComb_ >= 1 && iComb == (unsigned int)maxNComb_)
        break;
      iComb++;
      // partons
      pPartonsHadP->push_back(result->HadP);
      pPartonsHadQ->push_back(result->HadQ);
      pPartonsHadB->push_back(result->HadB);
      pPartonsLepB->push_back(result->LepB);
      // lepton
      pLeptons->push_back(result->LepL);
      // neutrino
      pNeutrinos->push_back(result->LepN);
      // indices referring to the jet combination
      pCombi->push_back(result->JetCombi);
      // chi2
      pChi2->push_back(result->Chi2);
      // chi2 probability
      pProb->push_back(result->Prob);
      // fitted top mass
      pMT->push_back(result->MT);
      pSigMT->push_back(result->SigMT);
      // status of the fitter
      pStatus->push_back(result->Status);
    }
  }
  evt.put(std::move(pCombi));
  evt.put(std::move(pPartonsHadP), "PartonsHadP");
  evt.put(std::move(pPartonsHadQ), "PartonsHadQ");
  evt.put(std::move(pPartonsHadB), "PartonsHadB");
  evt.put(std::move(pPartonsLepB), "PartonsLepB");
  evt.put(std::move(pLeptons), "Leptons");
  evt.put(std::move(pNeutrinos), "Neutrinos");
  evt.put(std::move(pChi2), "Chi2");
  evt.put(std::move(pProb), "Prob");
  evt.put(std::move(pMT), "MT");
  evt.put(std::move(pSigMT), "SigMT");
  evt.put(std::move(pStatus), "Status");
  evt.put(std::move(pJetsConsidered), "NumberOfConsideredJets");
}

Member Data Documentation

template<typename LeptonCollection >

std::string TtSemiLepHitFitProducer< LeptonCollection >::bTagAlgo_

private

input tag for b-tagging algorithm

Definition at line 46 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

hitfit::LeptonTranslatorBase<pat::Electron> TtSemiLepHitFitProducer< LeptonCollection >::electronTranslator_

private

Definition at line 90 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

PatHitFit* TtSemiLepHitFitProducer< LeptonCollection >::HitFit

private

Definition at line 95 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

edm::FileInPath TtSemiLepHitFitProducer< LeptonCollection >::hitfitBJetResolution_

private

Definition at line 87 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

edm::FileInPath TtSemiLepHitFitProducer< LeptonCollection >::hitfitDefault_

private

Definition at line 83 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

edm::FileInPath TtSemiLepHitFitProducer< LeptonCollection >::hitfitElectronResolution_

private

Definition at line 84 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

edm::FileInPath TtSemiLepHitFitProducer< LeptonCollection >::hitfitMETResolution_

private

Definition at line 88 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

edm::FileInPath TtSemiLepHitFitProducer< LeptonCollection >::hitfitMuonResolution_

private

Definition at line 85 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

edm::FileInPath TtSemiLepHitFitProducer< LeptonCollection >::hitfitUdscJetResolution_

private

Definition at line 86 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::jes_

private

jet energy scale

Definition at line 62 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::jesB_

private

Definition at line 63 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

std::string TtSemiLepHitFitProducer< LeptonCollection >::jetCorrectionLevel_

private

jet correction level

Definition at line 59 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

edm::EDGetTokenT<std::vector<pat::Jet> > TtSemiLepHitFitProducer< LeptonCollection >::jetsToken_

private

Definition at line 29 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

hitfit::JetTranslatorBase<pat::Jet> TtSemiLepHitFitProducer< LeptonCollection >::jetTranslator_

private

Definition at line 92 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

edm::EDGetTokenT<LeptonCollection> TtSemiLepHitFitProducer< LeptonCollection >::lepsToken_

private

Definition at line 30 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::maxBTagValueNonBJet_

private

max value of bTag for a non-b-jet

Definition at line 50 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::maxEtaEle_

private

maximum eta value for electrons, needed to limited range in which resolutions are provided

Definition at line 41 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::maxEtaJet_

private

maximum eta value for jets, needed to limited range in which resolutions are provided

Definition at line 43 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::maxEtaMu_

private

maximum eta value for muons, needed to limited range in which resolutions are provided

Definition at line 39 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

int TtSemiLepHitFitProducer< LeptonCollection >::maxNComb_

private

maximal number of combinations to be written to the event

Definition at line 36 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

int TtSemiLepHitFitProducer< LeptonCollection >::maxNJets_

private

maximal number of jets (-1 possible to indicate 'all')

Definition at line 34 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

edm::EDGetTokenT<std::vector<pat::MET> > TtSemiLepHitFitProducer< LeptonCollection >::metsToken_

private

Definition at line 31 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

hitfit::METTranslatorBase<pat::MET> TtSemiLepHitFitProducer< LeptonCollection >::metTranslator_

private

Definition at line 93 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::minBTagValueBJet_

private

min value of bTag for a b-jet

Definition at line 48 of file TtSemiLepHitFitProducer.h.

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::mTop_

private

Definition at line 56 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

hitfit::LeptonTranslatorBase<pat::Muon> TtSemiLepHitFitProducer< LeptonCollection >::muonTranslator_

private

Definition at line 91 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

double TtSemiLepHitFitProducer< LeptonCollection >::mW_

private

constraints

Definition at line 55 of file TtSemiLepHitFitProducer.h.

Referenced by TtSemiLepHitFitProducer< LeptonCollection >::TtSemiLepHitFitProducer().

template<typename LeptonCollection >

bool TtSemiLepHitFitProducer< LeptonCollection >::useBTag_

private

switch to tell whether to use b-tagging or not

Definition at line 52 of file TtSemiLepHitFitProducer.h.