#include <StEvtSolutionMaker.h>

Inheritance diagram for StEvtSolutionMaker:

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

	StEvtSolutionMaker (const edm::ParameterSet &)

	~StEvtSolutionMaker () 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 Attributes
bool	addLRJetComb_

std::vector< int >	constraints_

bool	doKinFit_

edm::EDGetTokenT< std::vector < pat::Electron > >	electronSrcToken_

edm::EDGetTokenT< StGenEvent >	genEvtSrcToken_

int	jetCorrScheme_

int	jetParam_

edm::EDGetTokenT< std::vector < pat::Jet > >	jetSrcToken_

int	lepParam_

std::string	leptonFlavour_

bool	matchToGenEvt_

double	maxDeltaS_

double	maxF_

int	maxNrIter_

int	metParam_

edm::EDGetTokenT< std::vector < pat::MET > >	metSrcToken_

edm::EDGetTokenT< std::vector < pat::Muon > >	muonSrcToken_

StKinFitter *	myKinFitter

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

Definition at line 16 of file StEvtSolutionMaker.h.

Constructor & Destructor Documentation

StEvtSolutionMaker::StEvtSolutionMaker ( const edm::ParameterSet & iConfig )

explicit

Definition at line 8 of file StEvtSolutionMaker.cc.

References addLRJetComb_, constraints_, doKinFit_, electronSrcToken_, genEvtSrcToken_, edm::ParameterSet::getParameter(), HLT_FULL_cff::InputTag, jetCorrScheme_, jetParam_, jetSrcToken_, lepParam_, leptonFlavour_, matchToGenEvt_, maxDeltaS_, maxF_, maxNrIter_, metParam_, metSrcToken_, muonSrcToken_, myKinFitter, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                      {
   // configurables
   electronSrcToken_ = mayConsume<std::vector<pat::Electron> >(iConfig.getParameter<edm::InputTag>("electronSource"));
   muonSrcToken_ = mayConsume<std::vector<pat::Muon> >(iConfig.getParameter<edm::InputTag>("muonSource"));
   metSrcToken_ = consumes<std::vector<pat::MET> >(iConfig.getParameter<edm::InputTag>("metSource"));
   jetSrcToken_ = consumes<std::vector<pat::Jet> >(iConfig.getParameter<edm::InputTag>("jetSource"));
   genEvtSrcToken_ = mayConsume<StGenEvent>(edm::InputTag("genEvt"));
   leptonFlavour_ = iConfig.getParameter<std::string>("leptonFlavour");
   jetCorrScheme_ = iConfig.getParameter<int>("jetCorrectionScheme");
   //jetInput_        = iConfig.getParameter< std::string >    ("jetInput");
   doKinFit_ = iConfig.getParameter<bool>("doKinFit");
   addLRJetComb_ = iConfig.getParameter<bool>("addLRJetComb");
   maxNrIter_ = iConfig.getParameter<int>("maxNrIter");
   maxDeltaS_ = iConfig.getParameter<double>("maxDeltaS");
   maxF_ = iConfig.getParameter<double>("maxF");
   jetParam_ = iConfig.getParameter<int>("jetParametrisation");
   lepParam_ = iConfig.getParameter<int>("lepParametrisation");
   metParam_ = iConfig.getParameter<int>("metParametrisation");
   constraints_ = iConfig.getParameter<std::vector<int> >("constraints");
   matchToGenEvt_ = iConfig.getParameter<bool>("matchToGenEvt");
 
   // define kinfitter
   if (doKinFit_) {
     myKinFitter = new StKinFitter(jetParam_, lepParam_, metParam_, maxNrIter_, maxDeltaS_, maxF_, constraints_);
   }
   // define what will be produced
   produces<std::vector<StEvtSolution> >();
 }

StEvtSolutionMaker::~StEvtSolutionMaker ( )

override

Definition at line 37 of file StEvtSolutionMaker.cc.

References doKinFit_, and myKinFitter.

                                         {
   if (doKinFit_)
     delete myKinFitter;
 }

Member Function Documentation

void StEvtSolutionMaker::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 42 of file StEvtSolutionMaker.cc.

References StKinFitter::addKinFitInfo(), b, gather_cfg::cout, doKinFit_, HI_PhotonSkim_cff::electrons, electronSrcToken_, edm::Event::getByToken(), jetCorrScheme_, fwrapper::jets, jetSrcToken_, cmsLHEtoEOSManager::l, leptonFlavour_, matchToGenEvt_, metSrcToken_, eostools::move(), patZpeak::muons, muonSrcToken_, myKinFitter, edm::Event::put(), StEvtSolution::setBottom(), StEvtSolution::setElectron(), StEvtSolution::setJetCorrectionScheme(), StEvtSolution::setLight(), StEvtSolution::setMuon(), and StEvtSolution::setNeutrino().

                                                                               {
   //
   //  TopObject Selection
   //
 
   // select lepton (the TtLepton vectors are, for the moment, sorted on pT)
   bool leptonFound = false;
   edm::Handle<std::vector<pat::Muon> > muons;
   if (leptonFlavour_ == "muon") {
     iEvent.getByToken(muonSrcToken_, muons);
     if (!muons->empty())
       leptonFound = true;
   }
   edm::Handle<std::vector<pat::Electron> > electrons;
   if (leptonFlavour_ == "electron") {
     iEvent.getByToken(electronSrcToken_, electrons);
     if (!electrons->empty())
       leptonFound = true;
   }
 
   // select MET (TopMET vector is sorted on ET)
   bool metFound = false;
   edm::Handle<std::vector<pat::MET> > mets;
   iEvent.getByToken(metSrcToken_, mets);
   if (!mets->empty())
     metFound = true;
 
   // select Jets
   bool jetsFound = false;
   edm::Handle<std::vector<pat::Jet> > jets;
   iEvent.getByToken(jetSrcToken_, jets);
   unsigned int maxJets = 2;  //this has to become a custom-defined parameter (we may want 2 or 3 jets)
   if (jets->size() >= 2)
     jetsFound = true;
 
   std::vector<StEvtSolution>* evtsols = new std::vector<StEvtSolution>();
   if (leptonFound && metFound && jetsFound) {
     std::cout << "constructing solutions" << std::endl;
     for (unsigned int b = 0; b < maxJets; b++) {
       for (unsigned int l = 0; l < maxJets; l++) {
         if (b != l) {  // to avoid double counting
           StEvtSolution asol;
           asol.setJetCorrectionScheme(jetCorrScheme_);
           if (leptonFlavour_ == "muon")
             asol.setMuon(muons, 0);
           if (leptonFlavour_ == "electron")
             asol.setElectron(electrons, 0);
           asol.setNeutrino(mets, 0);
           asol.setBottom(jets, b);
           asol.setLight(jets, l);
 
           if (doKinFit_)
             asol = myKinFitter->addKinFitInfo(&asol);
 
           /* to be adapted to ST (Andrea)
          if(addLRJetComb_){
          asol.setPtrueCombExist(jetCombProbs[m].getPTrueCombExist(&afitsol));
          asol.setPtrueBJetSel(jetCombProbs[m].getPTrueBJetSel(&afitsol));
          asol.setPtrueBhadrSel(jetCombProbs[m].getPTrueBhadrSel(&afitsol));
          asol.setPtrueJetComb(afitsol.getPtrueCombExist()*afitsol.getPtrueBJetSel()*afitsol.getPtrueBhadrSel());
          }
       */
           evtsols->push_back(asol);
         }
       }
     }
 
     // if asked for, match the event solutions to the gen Event
     if (matchToGenEvt_) {
       /*
     edm::Handle<StGenEvent> genEvt;
     iEvent.getByToken (genEvtSrcToken_,genEvt);
     double bestSolDR = 9999.;
     int bestSol = 0;
     for(size_t s=0; s<evtsols->size(); s++) {
     (*evtsols)[s].setGenEvt(genEvt->particles());
     vector<double> bm = BestMatch((*evtsols)[s], false); //false to use DR, true SpaceAngles
     (*evtsols)[s].setSumDeltaRjp(bm[0]); // dRBB + dRLL
     (*evtsols)[s].setChangeBL((int) bm[1]); // has swapped or not
     (*evtsols)[s].setDeltaRB(bm[2]);
     (*evtsols)[s].setDeltaRL(bm[3]);
     if(bm[0]<bestSolDR) { bestSolDR =  bm[0]; bestSol = s; }
     }
     (*evtsols)[bestSol].setBestSol(true);
       */
     }
 
     //store the vector of solutions to the event
     std::unique_ptr<std::vector<StEvtSolution> > pOut(evtsols);
     iEvent.put(std::move(pOut));
   } else {
     std::cout << "@@@ No calibrated solutions built, because:  " << std::endl;
     ;
     if (jets->size() < maxJets)
       std::cout << "@ nr jets = " << jets->size() << " < " << maxJets << std::endl;
     if (leptonFlavour_ == "muon" && !leptonFound)
       std::cout << "@ no good muon candidate" << std::endl;
     if (leptonFlavour_ == "electron" && !leptonFound)
       std::cout << "@ no good electron candidate" << std::endl;
     if (mets->empty())
       std::cout << "@ no MET reconstruction" << std::endl;
 
     StEvtSolution asol;
     evtsols->push_back(asol);
     std::unique_ptr<std::vector<StEvtSolution> > pOut(evtsols);
     iEvent.put(std::move(pOut));
   }
 }