PATHemisphereProducer Class Reference

#include <PATHemisphereProducer.h>

Inheritance diagram for PATHemisphereProducer:

Public Member Functions
	PATHemisphereProducer (const edm::ParameterSet &)
	~PATHemisphereProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
typedef std::vector< float >	HemiAxis

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

Private Attributes
int	_combinationMethod
float	_maxElectronEta
float	_maxJetEta
float	_maxMuonEta
float	_maxPhotonEta
float	_maxTauEta
float	_minElectronEt
float	_minJetEt
float	_minMuonEt
float	_minPhotonEt
float	_minTauEt
edm::EDGetTokenT < reco::CandidateView >	_patElectronsToken
edm::EDGetTokenT < reco::CandidateView >	_patJetsToken
	Input: All PAT objects that are to cross-clean or needed for that. More...
edm::EDGetTokenT < reco::CandidateView >	_patMuonsToken
edm::EDGetTokenT < reco::CandidateView >	_patPhotonsToken
edm::EDGetTokenT < reco::CandidateView >	_patTausToken
int	_seedMethod
std::vector< reco::CandidatePtr >	componentPtrs_
HemisphereAlgo *	myHemi
std::vector< float >	vA1
std::vector< float >	vA2
std::vector< float >	vE
std::vector< int >	vgroups
std::vector< float >	vPx
std::vector< float >	vPy
std::vector< float >	vPz

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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 81 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::endJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 203 of file PATHemisphereProducer.cc.

                              {

}

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

overrideprivatevirtual

Implements edm::EDProducer.

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, componentPtrs_, eta(), HemisphereAlgo::getAxis1(), HemisphereAlgo::getAxis2(), edm::Event::getByToken(), HemisphereAlgo::getGrouping(), razorHemispheres_cff::hemispheres, i, myHemi, p1, p2, edm::Event::put(), vA1, vA2, vE, vgroups, vPx, vPy, and vPz.

{
   using namespace edm;
   using namespace std;

   //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
   std::auto_ptr< std::vector<Hemisphere> > hemispheres(new 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(hemispheres);

  //clean up

    vPx.clear();
    vPy.clear();
    vPz.clear();
    vE.clear();
    vgroups.clear();
    componentPtrs_.clear();
}

Member Data Documentation

int PATHemisphereProducer::_combinationMethod

private

Definition at line 71 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_maxElectronEta

private

Definition at line 66 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_maxJetEta

private

Definition at line 64 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_maxMuonEta

private

Definition at line 65 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_maxPhotonEta

private

Definition at line 68 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_maxTauEta

private

Definition at line 67 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_minElectronEt

private

Definition at line 60 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_minJetEt

private

Definition at line 58 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_minMuonEt

private

Definition at line 59 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_minPhotonEt

private

Definition at line 62 of file PATHemisphereProducer.h.

Referenced by produce().

float PATHemisphereProducer::_minTauEt

private

Definition at line 61 of file PATHemisphereProducer.h.

Referenced by produce().

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

private

Definition at line 54 of file PATHemisphereProducer.h.

Referenced by produce().

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

private

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

Definition at line 51 of file PATHemisphereProducer.h.

Referenced by produce().

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

private

Definition at line 53 of file PATHemisphereProducer.h.

Referenced by produce().

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

private

Definition at line 55 of file PATHemisphereProducer.h.

Referenced by produce().

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

private

Definition at line 56 of file PATHemisphereProducer.h.

Referenced by produce().

int PATHemisphereProducer::_seedMethod

private

Definition at line 70 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<reco::CandidatePtr> PATHemisphereProducer::componentPtrs_

private

Definition at line 78 of file PATHemisphereProducer.h.

Referenced by produce().

HemisphereAlgo* PATHemisphereProducer::myHemi

private

Definition at line 73 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<float> PATHemisphereProducer::vA1

private

Definition at line 76 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<float> PATHemisphereProducer::vA2

private

Definition at line 76 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<float> PATHemisphereProducer::vE

private

Definition at line 75 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<int> PATHemisphereProducer::vgroups

private

Definition at line 77 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<float> PATHemisphereProducer::vPx

private

Definition at line 75 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<float> PATHemisphereProducer::vPy

private

Definition at line 75 of file PATHemisphereProducer.h.

Referenced by produce().

std::vector<float> PATHemisphereProducer::vPz

private

Definition at line 75 of file PATHemisphereProducer.h.

Referenced by produce().