#include <DiJetVarProducer.h>

Inheritance diagram for DiJetVarProducer:

Public Member Functions
	DiJetVarProducer (const edm::ParameterSet &)

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

	~DiJetVarProducer () 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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::InputTag	inputJetTag_

edm::EDGetTokenT< reco::CaloJetCollection >	inputJetTagToken_

double	wideJetDeltaR_

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	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

Detailed Description

Definition at line 19 of file DiJetVarProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 24 of file DiJetVarProducer.cc.

References edm::InputTag::encode(), edm::ParameterSet::getParameter(), inputJetTag_, inputJetTagToken_, LogDebug, and wideJetDeltaR_.

     : inputJetTag_(iConfig.getParameter<edm::InputTag>("inputJetTag")),
       wideJetDeltaR_(iConfig.getParameter<double>("wideJetDeltaR")) {
   // register your products
   // produces<std::vector<double> >("dijetvariables");
   produces<std::vector<math::PtEtaPhiMLorentzVector>>("widejets");
 
   // set Token(-s)
   inputJetTagToken_ = consumes<reco::CaloJetCollection>(iConfig.getParameter<edm::InputTag>("inputJetTag"));
 
   LogDebug("") << "Input Jet Tag: " << inputJetTag_.encode() << " ";
   LogDebug("") << "Radius Parameter Wide Jet: " << wideJetDeltaR_ << ".";
 }

DiJetVarProducer::~DiJetVarProducer ( )

override

Definition at line 38 of file DiJetVarProducer.cc.

38 {}

Member Function Documentation

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

override

Definition at line 41 of file DiJetVarProducer.cc.

References DEFINE_FWK_MODULE, edm::Event::getByToken(), inputJetTagToken_, eostools::move(), edm::Event::put(), and wideJetDeltaR_.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                             {
   using namespace std;
   using namespace edm;
   using namespace reco;
 
   // ## The output collections
   // std::unique_ptr<std::vector<double> > dijetvariables(new
   // std::vector<double>);
   std::unique_ptr<std::vector<math::PtEtaPhiMLorentzVector>> widejets(new std::vector<math::PtEtaPhiMLorentzVector>);
 
   // ## Get jet collection
   edm::Handle<reco::CaloJetCollection> calojets_handle;
   iEvent.getByToken(inputJetTagToken_, calojets_handle);
   // cout << "size: " << calojets_handle->size() << endl;
 
   // ## Wide Jet Algorithm
   // At least two jets
   if (calojets_handle->size() >= 2) {
     TLorentzVector wj1_tmp;
     TLorentzVector wj2_tmp;
     TLorentzVector wj1;
     TLorentzVector wj2;
     TLorentzVector wdijet;
 
     //        // Loop over all the input jets
     //        for(reco::CaloJetCollection::const_iterator it =
     //        calojets_handle->begin(); it != calojets_handle->end(); ++it)
     //           {
     //     cout << "jet: " << it->pt() << " " << it->eta() << " " << it->phi()
     // << endl;
     //           }
 
     // Find two leading jets
     TLorentzVector jet1, jet2;
 
     reco::CaloJetCollection::const_iterator j1 = calojets_handle->begin();
     reco::CaloJetCollection::const_iterator j2 = j1;
     ++j2;
 
     jet1.SetPtEtaPhiM(j1->pt(), j1->eta(), j1->phi(), j1->mass());
     jet2.SetPtEtaPhiM(j2->pt(), j2->eta(), j2->phi(), j2->mass());
 
     // cout << "j1: " << jet1.Pt() << " " << jet1.Eta() << " " << jet1.Phi() <<
     // endl; cout << "j2: " << jet2.Pt() << " " << jet2.Eta() << " " <<
     // jet2.Phi() << endl;
 
     // Create wide jets (radiation recovery algorithm)
     for (reco::CaloJetCollection::const_iterator it = calojets_handle->begin(); it != calojets_handle->end(); ++it) {
       TLorentzVector currentJet;
       currentJet.SetPtEtaPhiM(it->pt(), it->eta(), it->phi(), it->mass());
 
       double DeltaR1 = currentJet.DeltaR(jet1);
       double DeltaR2 = currentJet.DeltaR(jet2);
 
       if (DeltaR1 < DeltaR2 && DeltaR1 < wideJetDeltaR_) {
         wj1_tmp += currentJet;
       } else if (DeltaR2 < wideJetDeltaR_) {
         wj2_tmp += currentJet;
       }
     }
 
     // Re-order the wide jets in pT
     if (wj1_tmp.Pt() > wj2_tmp.Pt()) {
       wj1 = wj1_tmp;
       wj2 = wj2_tmp;
     } else {
       wj1 = wj2_tmp;
       wj2 = wj1_tmp;
     }
 
     // Create dijet system
     wdijet = wj1 + wj2;
 
     //        cout << "j1 wide: " << wj1.Pt() << " " << wj1.Eta() << " " <<
     //        wj1.Phi() << " " << wj1.M() << endl; cout << "j2 wide: " <<
     //        wj2.Pt() << " " << wj2.Eta() << " " << wj2.Phi() << " " << wj2.M()
     //        << endl; cout << "MJJWide: " << wdijet.M() << endl; cout <<
     //        "DeltaEtaJJWide: " << fabs(wj1.Eta()-wj2.Eta()) << endl; cout <<
     //        "DeltaPhiJJWide: " << fabs(wj1.DeltaPhi(wj2)) << endl;
 
     //        // Put variables in the container
     //        dijetvariables->push_back( wdijet.M() );                 //0 =
     //        MJJWide dijetvariables->push_back( fabs(wj1.Eta()-wj2.Eta()) );
     //        //1 = DeltaEtaJJWide dijetvariables->push_back(
     //        fabs(wj1.DeltaPhi(wj2)) );    //2 = DeltaPhiJJWide
 
     // Put widejets in the container
     math::PtEtaPhiMLorentzVector wj1math(wj1.Pt(), wj1.Eta(), wj1.Phi(), wj1.M());
     math::PtEtaPhiMLorentzVector wj2math(wj2.Pt(), wj2.Eta(), wj2.Phi(), wj2.M());
     widejets->push_back(wj1math);
     widejets->push_back(wj2math);
   }
   //    else
   //      {
   //        // Put variables in the container
   //        dijetvariables->push_back( -1 );                //0 = MJJWide
   //        dijetvariables->push_back( -1 );                //1 = DeltaEtaJJWide
   //        dijetvariables->push_back( -1 );                //2 = DeltaPhiJJWide
   //      }
 
   // ## Put objects in the Event
   // iEvent.put(std::move(dijetvariables), "dijetvariables");
   iEvent.put(std::move(widejets), "widejets");
 }