#include <L2TauModularIsolationSelector.h>

Inheritance diagram for L2TauModularIsolationSelector:

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

	~L2TauModularIsolationSelector ()

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

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

Private Attributes
edm::EDGetTokenT < reco::L2TauInfoAssociation >	associationInput_

std::vector< double >	ecalClusterDrRMS_

std::vector< double >	ecalClusterEtaRMS_

std::vector< double >	ecalClusterPhiRMS_

std::vector< double >	ecalIsolEt_

double	et_

std::vector< double >	hcalClusterDrRMS_

std::vector< double >	hcalClusterEtaRMS_

std::vector< double >	hcalClusterPhiRMS_

std::vector< double >	hcalIsolEt_

std::vector< double >	nEcalClusters_

std::vector< double >	nHcalClusters_

double	seedTowerEt_

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 22 of file L2TauModularIsolationSelector.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file L2TauModularIsolationSelector.cc.

                                                                                           :
   associationInput_(consumes<L2TauInfoAssociation>(iConfig.getParameter<edm::InputTag>("L2InfoAssociation"))),
 
   ecalIsolEt_(iConfig.getParameter<std::vector<double> >("EcalIsolationEt")),
   nEcalClusters_(iConfig.getParameter<std::vector<double> >("NumberOfECALClusters")),
   ecalClusterPhiRMS_(iConfig.getParameter<std::vector<double> >("ECALClusterPhiRMS")),
   ecalClusterEtaRMS_(iConfig.getParameter<std::vector<double> >("ECALClusterEtaRMS")),
   ecalClusterDrRMS_(iConfig.getParameter<std::vector<double> >("ECALClusterDRRMS")),
   hcalIsolEt_(iConfig.getParameter<std::vector<double> >("HcalIsolationEt")),
   nHcalClusters_(iConfig.getParameter<std::vector<double> >("NumberOfHCALClusters")),
   hcalClusterPhiRMS_(iConfig.getParameter<std::vector<double> >("HCALClusterPhiRMS")),
   hcalClusterEtaRMS_(iConfig.getParameter<std::vector<double> >("HCALClusterEtaRMS")),
   hcalClusterDrRMS_(iConfig.getParameter<std::vector<double> >("HCALClusterDRRMS")),
   et_(iConfig.getParameter<double>("MinJetEt")),
   seedTowerEt_(iConfig.getParameter<double>("SeedTowerEt"))
 {
   produces<CaloJetCollection>("Isolated");
 }

L2TauModularIsolationSelector::~L2TauModularIsolationSelector ( )

Definition at line 26 of file L2TauModularIsolationSelector.cc.

 {
 
 }

Member Function Documentation

void L2TauModularIsolationSelector::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 77 of file L2TauModularIsolationSelector.cc.

78 {

79 }

void L2TauModularIsolationSelector::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 83 of file L2TauModularIsolationSelector.cc.

83 {

84 }

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 37 of file L2TauModularIsolationSelector.cc.

References associationInput_, ecalClusterDrRMS_, ecalClusterEtaRMS_, ecalClusterPhiRMS_, reco::L2TauIsolationInfo::ecalClusterShape(), reco::L2TauIsolationInfo::ecalIsolEt(), ecalIsolEt_, reco::LeafCandidate::et(), et_, edm::Event::getByToken(), hcalClusterDrRMS_, hcalClusterEtaRMS_, hcalClusterPhiRMS_, reco::L2TauIsolationInfo::hcalClusterShape(), reco::L2TauIsolationInfo::hcalIsolEt(), hcalIsolEt_, metsig::jet, nEcalClusters_, reco::L2TauIsolationInfo::nEcalHits(), nHcalClusters_, reco::L2TauIsolationInfo::nHcalHits(), AlCaHLTBitMon_ParallelJobs::p, edm::Event::put(), reco::L2TauIsolationInfo::seedHcalHitEt(), and seedTowerEt_.

 {
    using namespace edm;
    edm::Handle<L2TauInfoAssociation> Imap;
 
    std::auto_ptr<CaloJetCollection> l2IsolCaloJets( new CaloJetCollection );
    iEvent.getByToken(associationInput_ ,Imap); 
 
    if(Imap->size()>0)
      for(L2TauInfoAssociation::const_iterator p = Imap->begin();p!=Imap->end();++p)
        {
          //Retrieve The L2TauIsolationInfo Class from the AssociationMap
          const L2TauIsolationInfo l2info = p->val;
          //Retrieve the Jet
          const CaloJet jet =*(p->key);
          
          //If The Cuts are Satisfied
        if(jet.et()>et_) 
          if(l2info.seedHcalHitEt()>seedTowerEt_)
            if(l2info.ecalIsolEt()< ecalIsolEt_[0]+ecalIsolEt_[1]*jet.et()+ecalIsolEt_[2]*jet.et()*jet.et())
          if(l2info.nEcalHits() <(int)(nEcalClusters_[0]+nEcalClusters_[1]*jet.et()+nEcalClusters_[2]*jet.et()*jet.et()))
            if(l2info.ecalClusterShape()[0] <ecalClusterEtaRMS_[0]+ecalClusterEtaRMS_[1]*jet.et()+ecalClusterEtaRMS_[2]*jet.et()*jet.et())
              if(l2info.ecalClusterShape()[1] <ecalClusterPhiRMS_[0]+ecalClusterPhiRMS_[1]*jet.et()+ecalClusterPhiRMS_[2]*jet.et()*jet.et())
                if(l2info.ecalClusterShape()[2] <ecalClusterDrRMS_[0]+ecalClusterDrRMS_[1]*jet.et()+ecalClusterDrRMS_[2]*jet.et()*jet.et())
              if(l2info.hcalIsolEt()<hcalIsolEt_[0]+hcalIsolEt_[1]*jet.et()+hcalIsolEt_[2]*jet.et()*jet.et())
                if(l2info.nHcalHits() <(int)(nHcalClusters_[0]+nHcalClusters_[1]*jet.et()+nHcalClusters_[2]*jet.et()*jet.et()))
                  if(l2info.hcalClusterShape()[0] <hcalClusterEtaRMS_[0]+hcalClusterEtaRMS_[1]*jet.et()+hcalClusterEtaRMS_[2]*jet.et()*jet.et())
                    if(l2info.hcalClusterShape()[1] <hcalClusterPhiRMS_[0]+hcalClusterPhiRMS_[1]*jet.et()+hcalClusterPhiRMS_[2]*jet.et()*jet.et())
                  if(l2info.hcalClusterShape()[2] <hcalClusterDrRMS_[0]+hcalClusterDrRMS_[1]*jet.et()+hcalClusterDrRMS_[2]*jet.et()*jet.et())
                    {
                      l2IsolCaloJets->push_back(jet);
                    }
 
        }
  
         iEvent.put(l2IsolCaloJets, "Isolated");
 }