#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

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)

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 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 39 of file DiJetVarProducer.cc.

40 {

41 }

Member Function Documentation

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

override

Definition at line 45 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");
 
 }