#include <PATHemisphereProducer.h>

Inheritance diagram for PATHemisphereProducer:

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

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

	~PATHemisphereProducer ()

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	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 &	itemsToGetFromEvent () 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

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

virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef std::vector< float >	HemiAxis

Private Attributes
const int	_combinationMethod

const float	_maxElectronEta

const float	_maxJetEta

const float	_maxMuonEta

const float	_maxPhotonEta

const float	_maxTauEta

const float	_minElectronEt

const float	_minJetEt

const float	_minMuonEt

const float	_minPhotonEt

const float	_minTauEt

const edm::EDGetTokenT< reco::CandidateView >	_patElectronsToken

const edm::EDGetTokenT< reco::CandidateView >	_patJetsToken
	Input: All PAT objects that are to cross-clean or needed for that. More...

const edm::EDGetTokenT< reco::CandidateView >	_patMuonsToken

const edm::EDGetTokenT< reco::CandidateView >	_patPhotonsToken

const edm::EDGetTokenT< reco::CandidateView >	_patTausToken

const int	_seedMethod

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 40 of file PATHemisphereProducer.h.

Member Typedef Documentation

typedef std::vector<float> PATHemisphereProducer::HemiAxis

private

Definition at line 72 of file PATHemisphereProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 55 of file PATHemisphereProducer.cc.

                                                                            :
   _patJetsToken       ( consumes<reco::CandidateView> ( iConfig.getParameter<edm::InputTag>( "patJets" ) ) ),
   _patMuonsToken      ( consumes<reco::CandidateView> ( iConfig.getParameter<edm::InputTag>( "patMuons" ) ) ),
   _patElectronsToken  ( consumes<reco::CandidateView> ( iConfig.getParameter<edm::InputTag>( "patElectrons" ) ) ),
   _patPhotonsToken    ( consumes<reco::CandidateView> ( iConfig.getParameter<edm::InputTag>( "patPhotons" ) ) ),
   _patTausToken       ( consumes<reco::CandidateView> ( iConfig.getParameter<edm::InputTag>( "patTaus" ) ) ),
 
   _minJetEt       ( iConfig.getParameter<double>("minJetEt") ),
   _minMuonEt       ( iConfig.getParameter<double>("minMuonEt") ),
   _minElectronEt       ( iConfig.getParameter<double>("minElectronEt") ),
   _minTauEt       ( iConfig.getParameter<double>("minTauEt") ),
   _minPhotonEt       ( iConfig.getParameter<double>("minPhotonEt") ),
 
   _maxJetEta       ( iConfig.getParameter<double>("maxJetEta") ),
   _maxMuonEta       ( iConfig.getParameter<double>("maxMuonEta") ),
   _maxElectronEta       ( iConfig.getParameter<double>("maxElectronEta") ),
   _maxTauEta       ( iConfig.getParameter<double>("maxTauEta") ),
   _maxPhotonEta       ( iConfig.getParameter<double>("maxPhotonEta") ),
 
   _seedMethod    ( iConfig.getParameter<int>("seedMethod") ),
   _combinationMethod ( iConfig.getParameter<int>("combinationMethod") )
 
 {
 
 
   produces< std::vector<pat::Hemisphere> >();
 }

PATHemisphereProducer::~PATHemisphereProducer ( )

Definition at line 84 of file PATHemisphereProducer.cc.

 {
 
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

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

overridevirtual

Definition at line 99 of file PATHemisphereProducer.cc.

References _combinationMethod, _maxElectronEta, _maxJetEta, _maxMuonEta, _maxPhotonEta, _maxTauEta, _minElectronEt, _minJetEt, _minMuonEt, _minPhotonEt, _minTauEt, _patElectronsToken, _patJetsToken, _patMuonsToken, _patPhotonsToken, _patTausToken, _seedMethod, DEFINE_FWK_MODULE, stringResolutionProvider_cfi::eta, HemisphereAlgo::getAxis1(), HemisphereAlgo::getAxis2(), edm::Event::getByToken(), HemisphereAlgo::getGrouping(), i, createfilelist::int, eostools::move(), p1, p2, edm::View< T >::ptrAt(), and edm::Event::put().

                                                                                                {
   using namespace edm;
   using namespace std;
   
   std::vector<float> vPx, vPy, vPz, vE;
   std::vector<float> vA1, vA2;
   std::vector<int> vgroups;
   std::vector<reco::CandidatePtr> componentPtrs;
 
   //Jets
   Handle<reco::CandidateView> pJets;
   iEvent.getByToken(_patJetsToken,pJets);
 
   //Muons
   Handle<reco::CandidateView> pMuons;
   iEvent.getByToken(_patMuonsToken,pMuons);
   
   //Electrons
   Handle<reco::CandidateView> pElectrons;
   iEvent.getByToken(_patElectronsToken,pElectrons);
   
   //Photons
   Handle<reco::CandidateView> pPhotons;
   iEvent.getByToken(_patPhotonsToken,pPhotons);
   
   //Taus
   Handle<reco::CandidateView> pTaus;
   iEvent.getByToken(_patTausToken,pTaus);
   
   
   //fill e,p vector with information from all objects (hopefully cleaned before)
   for(int i = 0; i < (int) (*pJets).size() ; i++){
     if((*pJets)[i].pt() <  _minJetEt || fabs((*pJets)[i].eta()) >  _maxJetEta) continue;
     
     componentPtrs.push_back(pJets->ptrAt(i));
   }
   
   for(int i = 0; i < (int) (*pMuons).size() ; i++){
     if((*pMuons)[i].pt() <  _minMuonEt || fabs((*pMuons)[i].eta()) >  _maxMuonEta) continue;
     
     componentPtrs.push_back(pMuons->ptrAt(i));
   }
   
   for(int i = 0; i < (int) (*pElectrons).size() ; i++){
     if((*pElectrons)[i].pt() <  _minElectronEt || fabs((*pElectrons)[i].eta()) >  _maxElectronEta) continue;
     
     componentPtrs.push_back(pElectrons->ptrAt(i));
   }
   
   for(int i = 0; i < (int) (*pPhotons).size() ; i++){
     if((*pPhotons)[i].pt() <  _minPhotonEt || fabs((*pPhotons)[i].eta()) >  _maxPhotonEta) continue;
     
     componentPtrs.push_back(pPhotons->ptrAt(i));
   }
   
   //aren't taus included in jets?
   for(int i = 0; i < (int) (*pTaus).size() ; i++){
     if((*pTaus)[i].pt() <  _minTauEt || fabs((*pTaus)[i].eta()) >  _maxTauEta) continue;
     
     componentPtrs.push_back(pTaus->ptrAt(i));
   }
   
   // create product
   auto hemispheres = std::make_unique<std::vector<Hemisphere>>();
   hemispheres->reserve(2);
   
   //calls HemiAlgorithm for seed method 3 (transv. inv. Mass) and association method 3 (Lund algo)
   HemisphereAlgo myHemi(componentPtrs,_seedMethod,_combinationMethod);
   
   //get Hemisphere Axis
   vA1 = myHemi.getAxis1();
   vA2 = myHemi.getAxis2();
   
   reco::Particle::LorentzVector p1(vA1[0]*vA1[3],vA1[1]*vA1[3],vA1[2]*vA1[3],vA1[4]);
   hemispheres->push_back(Hemisphere(p1));
 
   reco::Particle::LorentzVector p2(vA2[0]*vA2[3],vA2[1]*vA2[3],vA2[2]*vA2[3],vA2[4]);
   hemispheres->push_back(Hemisphere(p2));
 
   //get information to which Hemisphere each object belongs
   vgroups = myHemi.getGrouping();
 
   for ( unsigned int i=0; i<vgroups.size(); ++i ) {
     if ( vgroups[i]==1 ) {
       (*hemispheres)[0].addDaughter(componentPtrs[i]);
     }
     else {
       (*hemispheres)[1].addDaughter(componentPtrs[i]);
     }
   }
   
   
   iEvent.put(std::move(hemispheres));  
 }

Member Data Documentation

const int PATHemisphereProducer::_combinationMethod

private

Definition at line 70 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_maxElectronEta

private

Definition at line 65 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_maxJetEta

private

Definition at line 63 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_maxMuonEta

private

Definition at line 64 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_maxPhotonEta

private

Definition at line 67 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_maxTauEta

private

Definition at line 66 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_minElectronEt

private

Definition at line 59 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_minJetEt

private

Definition at line 57 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_minMuonEt

private

Definition at line 58 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_minPhotonEt

private

Definition at line 61 of file PATHemisphereProducer.h.

Referenced by produce().

const float PATHemisphereProducer::_minTauEt

private

Definition at line 60 of file PATHemisphereProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::CandidateView> PATHemisphereProducer::_patElectronsToken

private

Definition at line 53 of file PATHemisphereProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::CandidateView> PATHemisphereProducer::_patJetsToken

private

Input: All PAT objects that are to cross-clean or needed for that.

Definition at line 50 of file PATHemisphereProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::CandidateView> PATHemisphereProducer::_patMuonsToken

private

Definition at line 52 of file PATHemisphereProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::CandidateView> PATHemisphereProducer::_patPhotonsToken

private

Definition at line 54 of file PATHemisphereProducer.h.

Referenced by produce().

const edm::EDGetTokenT<reco::CandidateView> PATHemisphereProducer::_patTausToken

private

Definition at line 55 of file PATHemisphereProducer.h.

Referenced by produce().

const int PATHemisphereProducer::_seedMethod

private

Definition at line 69 of file PATHemisphereProducer.h.

Referenced by produce().

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation