L2TauModularIsolationProducer Class Reference

#include <L2TauModularIsolationProducer.h>

Inheritance diagram for L2TauModularIsolationProducer:

Classes
struct	RecHitPtComparator

Public Member Functions
	L2TauModularIsolationProducer (const edm::ParameterSet &)
	~L2TauModularIsolationProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Member Functions
virtual void	beginJob ()
virtual void	endJob ()
math::PtEtaPhiELorentzVectorCollection	getECALHits (const reco::CaloJet &, const edm::Event &, const edm::EventSetup &iSetup)
math::PtEtaPhiELorentzVectorCollection	getHCALHits (const reco::CaloJet &, const edm::Event &)
math::PtEtaPhiELorentzVectorCollection	getPFClusters (const reco::CaloJet &, const edm::Event &, const edm::InputTag &)
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
double	associationRadius_
edm::InputTag	caloTowers_
RecHitPtComparator	comparePt
double	crystalThresholdB_
double	crystalThresholdE_
edm::InputTag	EBRecHits_
std::string	ecalClusteringAlg_
std::string	ecalIsolationAlg_
edm::InputTag	EERecHits_
std::string	hcalClusteringAlg_
std::string	hcalIsolationAlg_
double	innerConeECAL_
double	innerConeHCAL_
edm::InputTag	l2CaloJets_
double	outerConeECAL_
double	outerConeHCAL_
edm::InputTag	pfClustersECAL_
edm::InputTag	pfClustersHCAL_
double	simpleClusterRadiusECAL_
double	simpleClusterRadiusHCAL_
double	towerThreshold_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
typedef WorkerT< EDProducer >	WorkerType
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)
Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 53 of file L2TauModularIsolationProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 10 of file L2TauModularIsolationProducer.cc.

                                                                                          :
  l2CaloJets_(iConfig.getParameter<edm::InputTag>("L2TauJetCollection")),
  EBRecHits_(iConfig.getParameter<edm::InputTag>("EBRecHits")),
  EERecHits_(iConfig.getParameter<edm::InputTag>("EERecHits")),
  caloTowers_(iConfig.getParameter<edm::InputTag>("CaloTowers")),
  pfClustersECAL_(iConfig.getParameter<edm::InputTag>("pfClustersECAL")),
  pfClustersHCAL_(iConfig.getParameter<edm::InputTag>("pfClustersHCAL")),
  ecalIsolationAlg_(iConfig.getParameter<std::string>("ecalIsolationAlgorithm")),
  hcalIsolationAlg_(iConfig.getParameter<std::string>("hcalIsolationAlgorithm")),
  ecalClusteringAlg_(iConfig.getParameter<std::string>("ecalClusteringAlgorithm")),
  hcalClusteringAlg_(iConfig.getParameter<std::string>("hcalClusteringAlgorithm")),
  associationRadius_(iConfig.getParameter<double>("associationRadius")),
  simpleClusterRadiusECAL_(iConfig.getParameter<double>("simpleClusterRadiusEcal")),
  simpleClusterRadiusHCAL_(iConfig.getParameter<double>("simpleClusterRadiusHcal")),
  innerConeECAL_(iConfig.getParameter<double>("innerConeECAL")),
  outerConeECAL_(iConfig.getParameter<double>("outerConeECAL")),
  innerConeHCAL_(iConfig.getParameter<double>("innerConeHCAL")),
  outerConeHCAL_(iConfig.getParameter<double>("outerConeHCAL")),
  crystalThresholdE_(iConfig.getParameter<double>("crystalThresholdEE")),
  crystalThresholdB_(iConfig.getParameter<double>("crystalThresholdEB")),
  towerThreshold_(iConfig.getParameter<double>("towerThreshold"))

 {
       
  //Add the products
  produces<L2TauInfoAssociation>();

}

L2TauModularIsolationProducer::~L2TauModularIsolationProducer

(

)

Definition at line 40 of file L2TauModularIsolationProducer.cc.

{
  //Destruction

}

Member Function Documentation

void L2TauModularIsolationProducer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 206 of file L2TauModularIsolationProducer.cc.

{
}

void L2TauModularIsolationProducer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 212 of file L2TauModularIsolationProducer.cc.

                                      {
}

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

private

Definition at line 244 of file L2TauModularIsolationProducer.cc.

References associationRadius_, comparePt, crystalThresholdB_, crystalThresholdE_, egHLT::errCodes::EBRecHits, EBRecHits_, DetId::Ecal, EcalBarrel, EcalEndcap, egHLT::errCodes::EERecHits, EERecHits_, relval_parameters_module::energy, ExpressReco_HICollisions_FallBack::et, PV3DBase< T, PVType, FrameType >::eta(), eta(), geometry, edm::EventSetup::get(), edm::Event::getByLabel(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), M_PI, L1TEmulatorMonitor_cff::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
  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.getByLabel( 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.pt()>crystalThresholdB_)
         jetRecHits.push_back(p);
       }

 if(iEvent.getByLabel( 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.pt()>crystalThresholdE_)
         jetRecHits.push_back(p);
       }

 if(jetRecHits.size()>0)
   std::sort(jetRecHits.begin(),jetRecHits.end(),comparePt);

 return jetRecHits;
}

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

private

Definition at line 219 of file L2TauModularIsolationProducer.cc.

References associationRadius_, caloTowers_, comparePt, relval_parameters_module::energy, ExpressReco_HICollisions_FallBack::et, eta(), edm::Event::getByLabel(), i, reco::LeafCandidate::p4(), phi, python.multivaluedict::sort(), and towerThreshold_.

Referenced by produce().

{
  edm::Handle<CaloTowerCollection> towers;

  math::PtEtaPhiELorentzVectorCollection towers2;

  if(iEvent.getByLabel(caloTowers_,towers))
    if(towers->size()>0)
    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.pt()>towerThreshold_)
          towers2.push_back(tower);
      }   
      }

  if(towers2.size()>0)
    std::sort(towers2.begin(),towers2.end(),comparePt);
  return towers2;
}

math::PtEtaPhiELorentzVectorCollection L2TauModularIsolationProducer::getPFClusters ( const reco::CaloJet &  jet, const edm::Event &  iEvent, const edm::InputTag &  input  )

private

Definition at line 307 of file L2TauModularIsolationProducer.cc.

References associationRadius_, trackerHits::c, comparePt, relval_parameters_module::energy, ExpressReco_HICollisions_FallBack::et, eta(), edm::Event::getByLabel(), M_PI, L1TEmulatorMonitor_cff::p, reco::LeafCandidate::p4(), phi, funct::sin(), python.multivaluedict::sort(), and theta().

Referenced by produce().

{
  edm::Handle<PFClusterCollection> clusters;
  math::PtEtaPhiELorentzVectorCollection clusters2;

  //get Clusters near the jet
  if(iEvent.getByLabel(input,clusters))
    if(clusters->size()>0)
    for(PFClusterCollection::const_iterator c = clusters->begin();c!=clusters->end();++c)
    {
      double energy = c->energy();
      double eta = c->eta();
      double phi = c->phi();
      double theta = c->position().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_)
    clusters2.push_back(p);
    }

  if(clusters2.size()>0)
    std::sort(clusters2.begin(),clusters2.end(),comparePt);
  return clusters2;
}

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

privatevirtual

Implements edm::EDProducer.

Definition at line 49 of file L2TauModularIsolationProducer.cc.

References L2TauSimpleClustering::clusterize(), L2TauIsolationAlgs::clusterShape(), ecalClusteringAlg_, ecalIsolationAlg_, edm::Event::getByLabel(), getECALHits(), getHCALHits(), getPFClusters(), hcalClusteringAlg_, hcalIsolationAlg_, spr::hitsHCAL(), info, innerConeECAL_, innerConeHCAL_, L2TauIsolationAlgs::isolatedEt(), l2CaloJets_, L2TauIsolationAlgs::nClustersAnnulus(), outerConeECAL_, outerConeHCAL_, pfClustersECAL_, pfClustersHCAL_, edm::Event::put(), reco::L2TauIsolationInfo::setEcalClusterShape(), reco::L2TauIsolationInfo::setEcalIsolEt(), reco::L2TauIsolationInfo::setHcalClusterShape(), reco::L2TauIsolationInfo::setHcalIsolEt(), reco::L2TauIsolationInfo::setNEcalHits(), reco::L2TauIsolationInfo::setNHcalHits(), reco::L2TauIsolationInfo::setSeedEcalHitEt(), reco::L2TauIsolationInfo::setSeedHcalHitEt(), simpleClusterRadiusECAL_, and simpleClusterRadiusHCAL_.

{
   edm::Handle<CaloJetCollection> l2CaloJets; //Handle to the input (L2TauCaloJets);
   iEvent.getByLabel(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; 

       //Define objects to be loaded from file:
       math::PtEtaPhiELorentzVectorCollection hitsECAL;
       math::PtEtaPhiELorentzVectorCollection hitsHCAL;
       math::PtEtaPhiELorentzVectorCollection pfClustersECAL;
       math::PtEtaPhiELorentzVectorCollection pfClustersHCAL;

       if(ecalIsolationAlg_=="recHits" || ecalClusteringAlg_=="recHits" || ecalIsolationAlg_ =="simpleClusters" ||ecalClusteringAlg_ =="simpleClusters")
         hitsECAL =getECALHits(*jc,iEvent,iSetup); 

       if(hcalIsolationAlg_=="recHits" || hcalClusteringAlg_=="recHits" || hcalIsolationAlg_ =="simpleClusters" ||hcalClusteringAlg_ =="simpleClusters")
         hitsHCAL =getHCALHits(*jc,iEvent); 

       if(ecalIsolationAlg_=="particleFlow" || ecalClusteringAlg_=="particleFlow")
         pfClustersECAL =getPFClusters(*jc,iEvent,pfClustersECAL_);

       if(hcalIsolationAlg_=="particleFlow" || hcalClusteringAlg_=="particleFlow")
         pfClustersHCAL =getPFClusters(*jc,iEvent,pfClustersHCAL_);



       //Do ECAL Isolation

       if(ecalIsolationAlg_ == "recHits")
         {
           //Use Rechits
               info.setEcalIsolEt( alg.isolatedEt(hitsECAL , jc->p4().Vect(), innerConeECAL_,outerConeECAL_) );
         }
       else if(ecalIsolationAlg_ == "simpleClusters")
         {
           //create the simple clusters
           L2TauSimpleClustering clustering(simpleClusterRadiusECAL_);
           math::PtEtaPhiELorentzVectorCollection clusters = clustering.clusterize(hitsECAL);
               info.setEcalIsolEt( alg.isolatedEt(clusters , jc->p4().Vect(), innerConeECAL_,outerConeECAL_) );

         }
       else if(ecalIsolationAlg_ == "particleFlow")
         {
           //Use ParticleFlow
               info.setEcalIsolEt( alg.isolatedEt(pfClustersECAL , jc->p4().Vect(), innerConeECAL_,outerConeECAL_) );
         }

       //Do ECAL Clustering

       if(ecalClusteringAlg_ == "recHits")
         {
           //Use Rechits
               info.setEcalClusterShape(alg.clusterShape(hitsECAL,jc->p4().Vect(),0.,outerConeECAL_) );
           info.setNEcalHits( alg.nClustersAnnulus(hitsECAL , jc->p4().Vect(), innerConeECAL_,outerConeECAL_));
           if(hitsECAL.size()>0)
         info.setSeedEcalHitEt(hitsECAL[0].pt());

         }
       else if(ecalClusteringAlg_ == "simpleClusters")
         {
           //create the simple clusters
           L2TauSimpleClustering clustering(simpleClusterRadiusECAL_);
           math::PtEtaPhiELorentzVectorCollection clusters = clustering.clusterize(hitsECAL);
               info.setEcalClusterShape(alg.clusterShape(clusters,jc->p4().Vect(),0.,outerConeECAL_) );
           info.setNEcalHits( alg.nClustersAnnulus(clusters, jc->p4().Vect(), innerConeECAL_,outerConeECAL_));

           if(clusters.size()>0)
         info.setSeedEcalHitEt(clusters[0].pt());
         }
       else if(ecalClusteringAlg_ == "particleFlow")
         {
           //Use ParticleFlow
               info.setEcalClusterShape(alg.clusterShape(pfClustersECAL,jc->p4().Vect(),0.,outerConeECAL_) );
           info.setNEcalHits( alg.nClustersAnnulus(pfClustersECAL, jc->p4().Vect(), innerConeECAL_,outerConeECAL_));
           if(pfClustersECAL.size()>0)
         info.setSeedEcalHitEt(pfClustersECAL[0].pt());
         }


       //Do HCAL Isolation

       if(hcalIsolationAlg_ == "recHits")
         {
           //Use Rechits
               info.setHcalIsolEt( alg.isolatedEt(hitsHCAL , jc->p4().Vect(), innerConeHCAL_,outerConeHCAL_) );
         }
       else if(hcalIsolationAlg_ == "simpleClusters")
         {
           //create the simple clusters
           L2TauSimpleClustering clustering(simpleClusterRadiusHCAL_);
           math::PtEtaPhiELorentzVectorCollection clusters = clustering.clusterize(hitsHCAL);
               info.setHcalIsolEt( alg.isolatedEt(clusters , jc->p4().Vect(), innerConeHCAL_,outerConeHCAL_) );

         }
       else if(hcalIsolationAlg_ == "particleFlow")
         {
           //Use ParticleFlow
               info.setHcalIsolEt( alg.isolatedEt(pfClustersHCAL , jc->p4().Vect(), innerConeHCAL_,outerConeHCAL_) );
         }



       //Do HCAL Clustering

       if(hcalClusteringAlg_ == "recHits")
         {
           //Use Rechits
               info.setHcalClusterShape(alg.clusterShape(hitsHCAL,jc->p4().Vect(),0.,outerConeHCAL_) );
           info.setNHcalHits( alg.nClustersAnnulus(hitsHCAL, jc->p4().Vect(), innerConeHCAL_,outerConeHCAL_));
           if(hitsHCAL.size()>0)
         info.setSeedHcalHitEt(hitsHCAL[0].pt());
         }
       else if(hcalClusteringAlg_ == "simpleClusters")
         {
           //create the simple clusters
           L2TauSimpleClustering clustering(simpleClusterRadiusHCAL_);
           math::PtEtaPhiELorentzVectorCollection clusters = clustering.clusterize(hitsHCAL);
               info.setHcalClusterShape(alg.clusterShape(clusters,jc->p4().Vect(),0.,outerConeHCAL_) );
           info.setNHcalHits( alg.nClustersAnnulus(clusters, jc->p4().Vect(), innerConeHCAL_,outerConeHCAL_));
           if(clusters.size()>0)
         info.setSeedHcalHitEt(clusters[0].pt());
         }
       else if(hcalClusteringAlg_ == "particleFlow")
         {
           //Use ParticleFlow
               info.setHcalClusterShape(alg.clusterShape(pfClustersHCAL,jc->p4().Vect(),0.,outerConeHCAL_) );
           info.setNHcalHits( alg.nClustersAnnulus(pfClustersHCAL, jc->p4().Vect(), innerConeHCAL_,outerConeHCAL_));
           if(pfClustersHCAL.size()>0)
         info.setSeedHcalHitEt(pfClustersHCAL[0].pt());
         }

       //Store the info Class
       edm::Ref<CaloJetCollection> jcRef(l2CaloJets,jc-jcStart);
       l2InfoAssoc->insert(jcRef, info);
     }

     } //end of if(*jetCrystalsObj)

    iEvent.put(l2InfoAssoc);
}

Member Data Documentation

double L2TauModularIsolationProducer::associationRadius_

private

Definition at line 81 of file L2TauModularIsolationProducer.h.

Referenced by getECALHits(), getHCALHits(), and getPFClusters().

edm::InputTag L2TauModularIsolationProducer::caloTowers_

private

Definition at line 71 of file L2TauModularIsolationProducer.h.

Referenced by getHCALHits().

RecHitPtComparator L2TauModularIsolationProducer::comparePt

private

Definition at line 104 of file L2TauModularIsolationProducer.h.

Referenced by getECALHits(), getHCALHits(), and getPFClusters().

double L2TauModularIsolationProducer::crystalThresholdB_

private

Definition at line 92 of file L2TauModularIsolationProducer.h.

Referenced by getECALHits().

double L2TauModularIsolationProducer::crystalThresholdE_

private

Definition at line 91 of file L2TauModularIsolationProducer.h.

Referenced by getECALHits().

edm::InputTag L2TauModularIsolationProducer::EBRecHits_

private

Definition at line 69 of file L2TauModularIsolationProducer.h.

Referenced by getECALHits().

std::string L2TauModularIsolationProducer::ecalClusteringAlg_

private

Definition at line 78 of file L2TauModularIsolationProducer.h.

Referenced by produce().

std::string L2TauModularIsolationProducer::ecalIsolationAlg_

private

Definition at line 76 of file L2TauModularIsolationProducer.h.

Referenced by produce().

edm::InputTag L2TauModularIsolationProducer::EERecHits_

private

Definition at line 70 of file L2TauModularIsolationProducer.h.

Referenced by getECALHits().

std::string L2TauModularIsolationProducer::hcalClusteringAlg_

private

Definition at line 79 of file L2TauModularIsolationProducer.h.

Referenced by produce().

std::string L2TauModularIsolationProducer::hcalIsolationAlg_

private

Definition at line 77 of file L2TauModularIsolationProducer.h.

Referenced by produce().

double L2TauModularIsolationProducer::innerConeECAL_

private

Definition at line 85 of file L2TauModularIsolationProducer.h.

Referenced by produce().

double L2TauModularIsolationProducer::innerConeHCAL_

private

Definition at line 87 of file L2TauModularIsolationProducer.h.

Referenced by produce().

edm::InputTag L2TauModularIsolationProducer::l2CaloJets_

private

Definition at line 68 of file L2TauModularIsolationProducer.h.

Referenced by produce().

double L2TauModularIsolationProducer::outerConeECAL_

private

Definition at line 86 of file L2TauModularIsolationProducer.h.

Referenced by produce().

double L2TauModularIsolationProducer::outerConeHCAL_

private

Definition at line 88 of file L2TauModularIsolationProducer.h.

Referenced by produce().

edm::InputTag L2TauModularIsolationProducer::pfClustersECAL_

private

Definition at line 72 of file L2TauModularIsolationProducer.h.

Referenced by produce().

edm::InputTag L2TauModularIsolationProducer::pfClustersHCAL_

private

Definition at line 73 of file L2TauModularIsolationProducer.h.

Referenced by produce().

double L2TauModularIsolationProducer::simpleClusterRadiusECAL_

private

Definition at line 83 of file L2TauModularIsolationProducer.h.

Referenced by produce().

double L2TauModularIsolationProducer::simpleClusterRadiusHCAL_

private

Definition at line 84 of file L2TauModularIsolationProducer.h.

Referenced by produce().

double L2TauModularIsolationProducer::towerThreshold_

private

Definition at line 93 of file L2TauModularIsolationProducer.h.

Referenced by getHCALHits().