#include <L2TauNarrowConeIsolationProducer.h>

Inheritance diagram for L2TauNarrowConeIsolationProducer:

Classes
struct	CrystalPtComparator

Public Types
typedef reco::CaloJet	CaloJet

typedef reco::CaloJetCollection	CaloJetCollection

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

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

	~L2TauNarrowConeIsolationProducer ()

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

math::PtEtaPhiELorentzVectorCollection	getECALHits (const CaloJet &, const edm::Event &, const edm::EventSetup &iSetup)

math::PtEtaPhiELorentzVectorCollection	getHCALHits (const CaloJet &, const edm::Event &)

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

Private Attributes
double	associationRadius_

edm::EDGetTokenT < CaloTowerCollection >	CaloTowers_

CrystalPtComparator	comparePt

double	crystalThresholdB_

double	crystalThresholdE_

edm::EDGetTokenT < EBRecHitCollection >	EBRecHits_

double	ECALClustering_clusterRadius_

bool	ECALClustering_run_

double	ECALIsolation_innerCone_

double	ECALIsolation_outerCone_

bool	ECALIsolation_run_

edm::EDGetTokenT < EERecHitCollection >	EERecHits_

edm::EDGetTokenT < reco::CaloJetCollection >	l2CaloJets_

double	TowerIsolation_innerCone_

double	TowerIsolation_outerCone_

bool	TowerIsolation_run_

double	towerThreshold_

Additional Inherited Members
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 52 of file L2TauNarrowConeIsolationProducer.h.

Member Typedef Documentation

typedef reco::CaloJet L2TauNarrowConeIsolationProducer::CaloJet

Definition at line 54 of file L2TauNarrowConeIsolationProducer.h.

typedef reco::CaloJetCollection L2TauNarrowConeIsolationProducer::CaloJetCollection

Definition at line 55 of file L2TauNarrowConeIsolationProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 8 of file L2TauNarrowConeIsolationProducer.cc.

References ECALClustering_clusterRadius_, ECALClustering_run_, ECALIsolation_innerCone_, ECALIsolation_outerCone_, ECALIsolation_run_, edm::ParameterSet::getParameter(), TowerIsolation_innerCone_, TowerIsolation_outerCone_, and TowerIsolation_run_.

                                                                                                 :
   l2CaloJets_(consumes<CaloJetCollection>(iConfig.getParameter<edm::InputTag>("L2TauJetCollection"))),
   EBRecHits_(consumes<EBRecHitCollection>(iConfig.getParameter<edm::InputTag>("EBRecHits"))),
   EERecHits_(consumes<EERecHitCollection>(iConfig.getParameter<edm::InputTag>("EERecHits"))),
   CaloTowers_(consumes<CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("CaloTowers"))),
   associationRadius_(iConfig.getParameter<double>("associationRadius")),
   crystalThresholdE_(iConfig.getParameter<double>("crystalThresholdEE")),
   crystalThresholdB_(iConfig.getParameter<double>("crystalThresholdEB")),
   towerThreshold_(iConfig.getParameter<double>("towerThreshold"))
  {
         
   //ECAL Isolation
   edm::ParameterSet ECALIsolParams = iConfig.getParameter<edm::ParameterSet>("ECALIsolation") ;
   ECALIsolation_innerCone_ =  ECALIsolParams.getParameter<double>( "innerCone" );
   ECALIsolation_outerCone_ =  ECALIsolParams.getParameter<double>( "outerCone" );
   ECALIsolation_run_       =  ECALIsolParams.getParameter<bool>( "runAlgorithm" );
 
   //ECAL Clustering
   edm::ParameterSet ECALClusterParams = iConfig.getParameter<edm::ParameterSet>("ECALClustering") ;
   ECALClustering_run_                 =  ECALClusterParams.getParameter<bool>( "runAlgorithm" );
   ECALClustering_clusterRadius_       =  ECALClusterParams.getParameter<double>( "clusterRadius" );
     
   //Tower Isolation
   edm::ParameterSet TowerIsolParams = iConfig.getParameter<edm::ParameterSet>("TowerIsolation") ;
   TowerIsolation_innerCone_         =  TowerIsolParams.getParameter<double>( "innerCone" );
   TowerIsolation_outerCone_         =  TowerIsolParams.getParameter<double>( "outerCone" );
   TowerIsolation_run_               =  TowerIsolParams.getParameter<bool>( "runAlgorithm" );
  
 
   //Add the products
   produces<L2TauInfoAssociation>();
 
 }

L2TauNarrowConeIsolationProducer::~L2TauNarrowConeIsolationProducer ( )

Definition at line 43 of file L2TauNarrowConeIsolationProducer.cc.

 {
   //Destruction
 
 }

Member Function Documentation

void L2TauNarrowConeIsolationProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 119 of file L2TauNarrowConeIsolationProducer.cc.

120 {

121 }

void L2TauNarrowConeIsolationProducer::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 125 of file L2TauNarrowConeIsolationProducer.cc.

125 {

126 }

math::PtEtaPhiELorentzVectorCollection L2TauNarrowConeIsolationProducer::getECALHits	(	const CaloJet &	jet,
		const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 156 of file L2TauNarrowConeIsolationProducer.cc.

References associationRadius_, comparePt, crystalThresholdB_, crystalThresholdE_, egHLT::errCodes::EBRecHits, EBRecHits_, DetId::Ecal, EcalBarrel, EcalEndcap, HWWFunctions::EE, egHLT::errCodes::EERecHits, EERecHits_, relval_parameters_module::energy, PV3DBase< T, PVType, FrameType >::eta(), eta(), geometry, edm::EventSetup::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), M_PI, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::p4(), phi, PV3DBase< T, PVType, FrameType >::phi(), funct::sin(), python.multivaluedict::sort(), PV3DBase< T, PVType, FrameType >::theta(), and theta().

Referenced by produce().

 {
   //Init Geometry
   edm::ESHandle<CaloGeometry> geometry;
   iSetup.get<CaloGeometryRecord>().get(geometry);
 
   //Create ECAL Geometry
   const CaloSubdetectorGeometry* EB = geometry->getSubdetectorGeometry(DetId::Ecal,EcalBarrel);
   const CaloSubdetectorGeometry* EE = geometry->getSubdetectorGeometry(DetId::Ecal,EcalEndcap);
 
   //Handle To the ECAL
   edm::Handle<EBRecHitCollection> EBRecHits;
   edm::Handle<EERecHitCollection> EERecHits;
 
   //Create a container for the hits
   math::PtEtaPhiELorentzVectorCollection jetRecHits;
 
   //Loop on the barrel hits
   if(iEvent.getByToken( EBRecHits_, EBRecHits))
      for(EBRecHitCollection::const_iterator hit = EBRecHits->begin();hit!=EBRecHits->end();++hit)
        {
      //get Detector Geometry
      const CaloCellGeometry* this_cell = EB->getGeometry(hit->detid());
      GlobalPoint posi = this_cell->getPosition();
      double energy = hit->energy();
      double eta = posi.eta();
      double phi = posi.phi();
      double theta = posi.theta();
      if(theta > M_PI) theta = 2 * M_PI- theta;
      double et = energy * sin(theta);
      math::PtEtaPhiELorentzVector p(et, eta, phi, energy);
      if(ROOT::Math::VectorUtil::DeltaR(p,jet.p4()) <associationRadius_)
        if(p.energy()>crystalThresholdB_)
          jetRecHits.push_back(p);
        }
 
  if(iEvent.getByToken( EERecHits_, EERecHits))
      for(EERecHitCollection::const_iterator hit = EERecHits->begin();hit!=EERecHits->end();++hit)
        {
      //get Detector Geometry
      const CaloCellGeometry* this_cell = EE->getGeometry(hit->detid());
      GlobalPoint posi = this_cell->getPosition();
      double energy = hit->energy();
      double eta = posi.eta();
      double phi = posi.phi();
      double theta = posi.theta();
      if(theta > M_PI) theta = 2 * M_PI- theta;
      double et = energy * sin(theta);
      math::PtEtaPhiELorentzVector p(et, eta, phi, energy);
      if(ROOT::Math::VectorUtil::DeltaR(p,jet.p4()) < associationRadius_)
        if(p.energy()>crystalThresholdE_)
          jetRecHits.push_back(p);
        }
 
 
  std::sort(jetRecHits.begin(),jetRecHits.end(),comparePt);
 
  return jetRecHits;
 }

math::PtEtaPhiELorentzVectorCollection L2TauNarrowConeIsolationProducer::getHCALHits	(	const CaloJet &	jet,
		const edm::Event &	iEvent
	)

private

Definition at line 131 of file L2TauNarrowConeIsolationProducer.cc.

References associationRadius_, CaloTowers_, comparePt, relval_parameters_module::energy, eta(), edm::Event::getByToken(), i, reco::LeafCandidate::p4(), phi, python.multivaluedict::sort(), and towerThreshold_.

Referenced by produce().

 {
   edm::Handle<CaloTowerCollection> towers;
   math::PtEtaPhiELorentzVectorCollection towers2;
 
   if(iEvent.getByToken(CaloTowers_,towers))
     for(size_t i=0;i<towers->size();++i)
       {
     math::PtEtaPhiELorentzVector tower((*towers)[i].et(),(*towers)[i].eta(),(*towers)[i].phi(),(*towers)[i].energy());
     if(ROOT::Math::VectorUtil::DeltaR(tower,jet.p4()) <associationRadius_)
       {
         if(tower.energy()>towerThreshold_)
           towers2.push_back(tower);
 
       }   
 
       }
 
  std::sort(towers2.begin(),towers2.end(),comparePt);
 
   return towers2;
 }

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 52 of file L2TauNarrowConeIsolationProducer.cc.

References L2TauSimpleClustering::clusterize(), L2TauIsolationAlgs::clusterShape(), ECALClustering_clusterRadius_, ECALClustering_run_, ECALIsolation_innerCone_, ECALIsolation_outerCone_, ECALIsolation_run_, edm::Event::getByToken(), getECALHits(), getHCALHits(), spr::hitsHCAL(), info(), L2TauIsolationAlgs::isolatedEt(), l2CaloJets_, edm::Event::put(), reco::L2TauIsolationInfo::setEcalClusterShape(), reco::L2TauIsolationInfo::setEcalIsolEt(), reco::L2TauIsolationInfo::setHcalIsolEt(), reco::L2TauIsolationInfo::setNEcalHits(), reco::L2TauIsolationInfo::setSeedEcalHitEt(), reco::L2TauIsolationInfo::setSeedHcalHitEt(), TowerIsolation_innerCone_, TowerIsolation_outerCone_, and TowerIsolation_run_.

 {
 
 
    edm::Handle<CaloJetCollection> l2CaloJets; //Handle to the input (L2TauCaloJets);
    iEvent.getByToken(l2CaloJets_ ,l2CaloJets);//get the handle
 
    //Create the Association
    std::auto_ptr<L2TauInfoAssociation> l2InfoAssoc( new L2TauInfoAssociation);
 
    //If the JetCrystalsAssociation exists -> RUN The Producer
    if(l2CaloJets->size()>0)
      {
       CaloJetCollection::const_iterator jcStart = l2CaloJets->begin();
        //Loop on Jets
        for(CaloJetCollection::const_iterator jc = jcStart ;jc!=l2CaloJets->end();++jc)
      {
        //Create Algorithm Object
        L2TauIsolationAlgs alg;
        
        //Create Info Object
        L2TauIsolationInfo info; 
 
        //get Hits
        math::PtEtaPhiELorentzVectorCollection hitsECAL = getECALHits(*jc,iEvent,iSetup);
        math::PtEtaPhiELorentzVectorCollection hitsHCAL = getHCALHits(*jc,iEvent);
 
 
 
        //Run ECALIsolation 
        if(ECALIsolation_run_)
          {
            info.setEcalIsolEt( alg.isolatedEt(hitsECAL , jc->p4().Vect(), ECALIsolation_innerCone_,ECALIsolation_outerCone_) );
            if(hitsECAL.size()>0)
          info.setSeedEcalHitEt(hitsECAL[0].pt());
          }
 
        //Run ECALClustering 
        if(ECALClustering_run_)
          {
            //load simple clustering algorithm
            L2TauSimpleClustering clustering(ECALClustering_clusterRadius_);
            math::PtEtaPhiELorentzVectorCollection clusters = clustering.clusterize(hitsECAL);
            info.setEcalClusterShape(alg.clusterShape(clusters,jc->p4().Vect(),0,0.5) );
            info.setNEcalHits(clusters.size());
          }
 
        //Run CaloTower Isolation
            if(TowerIsolation_run_)
          {
            info.setHcalIsolEt( alg.isolatedEt(hitsHCAL , jc->p4().Vect(), TowerIsolation_innerCone_,TowerIsolation_outerCone_) );
            if(hitsHCAL.size()>0)
          info.setSeedHcalHitEt(hitsHCAL[0].pt());
          }
 
        //Store the info Class
        edm::Ref<CaloJetCollection> jcRef(l2CaloJets,jc-jcStart);
        l2InfoAssoc->insert(jcRef, info);
      }
 
      } //end of if(*jetCrystalsObj)
 
     iEvent.put(l2InfoAssoc);
 }