EgammaHLTRechitInRegionsProducer Class Reference

#include <EgammaHLTRechitInRegionsProducer.h>

Inheritance diagram for EgammaHLTRechitInRegionsProducer:

Public Member Functions
	EgammaHLTRechitInRegionsProducer (const edm::ParameterSet &ps)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~EgammaHLTRechitInRegionsProducer ()
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
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
bool	counterExceeded () const

Private Attributes
bool	doIsolated_
std::string	hitcollection_
std::vector< edm::InputTag >	hitLabels
edm::InputTag	hitproducer_
std::vector< edm::EDGetTokenT < EcalRecHitCollection > >	hitTokens
double	l1LowerThr_
double	l1LowerThrIgnoreIsolation_
edm::InputTag	l1TagIsolated_
edm::InputTag	l1TagNonIsolated_
double	l1UpperThr_
int	nEvt_
int	nMaxPrintout_
double	regionEtaMargin_
double	regionPhiMargin_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
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 16 of file EgammaHLTRechitInRegionsProducer.h.

Constructor & Destructor Documentation

EgammaHLTRechitInRegionsProducer::EgammaHLTRechitInRegionsProducer

(

const edm::ParameterSet &

ps

)

Definition at line 52 of file EgammaHLTRechitInRegionsProducer.cc.

References doIsolated_, edm::ParameterSet::getParameter(), hitLabels, hitproducer_, hitTokens, i, l1LowerThr_, l1LowerThrIgnoreIsolation_, l1TagIsolated_, l1TagNonIsolated_, l1UpperThr_, nEvt_, regionEtaMargin_, and regionPhiMargin_.

                                                                                            {

  hitproducer_   = ps.getParameter<edm::InputTag>("ecalhitproducer");

  l1TagIsolated_ = ps.getParameter< edm::InputTag > ("l1TagIsolated");
  l1TagNonIsolated_ = ps.getParameter< edm::InputTag > ("l1TagNonIsolated");
  doIsolated_   = ps.getParameter<bool>("doIsolated");
  
  l1LowerThr_ = ps.getParameter<double> ("l1LowerThr");
  l1UpperThr_ = ps.getParameter<double> ("l1UpperThr");
  l1LowerThrIgnoreIsolation_ = ps.getParameter<double> ("l1LowerThrIgnoreIsolation");

  regionEtaMargin_   = ps.getParameter<double>("regionEtaMargin");
  regionPhiMargin_   = ps.getParameter<double>("regionPhiMargin");

  const std::vector<std::string> flagnames = ps.getParameter<std::vector<std::string> >("RecHitFlagToBeExcluded");
  const std::vector<int> flagsexcl = StringToEnumValue<EcalRecHit::Flags>(flagnames);
  
  const std::vector<std::string> severitynames = ps.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcluded");
  const std::vector<int> severitiesexcl = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynames);

  hitLabels = ps.getParameter<std::vector<edm::InputTag>>("ecalhitLabels");
  for (unsigned int i=0; i<hitLabels.size(); i++) 
    hitTokens.push_back(consumes<EcalRecHitCollection>(hitLabels[i]));
    
  produces<EcalRecHitCollection> ("EcalRegionalRecHitsEB");
  produces<EcalRecHitCollection> ("EcalRegionalRecHitsEE");
  produces<EcalRecHitCollection> ("EcalRegionalRecHitsES");

  nEvt_ = 0;
}

EgammaHLTRechitInRegionsProducer::~EgammaHLTRechitInRegionsProducer

(

)

Definition at line 85 of file EgammaHLTRechitInRegionsProducer.cc.

{}

Member Function Documentation

bool EgammaHLTRechitInRegionsProducer::counterExceeded ( ) const

inlineprivate

Definition at line 48 of file EgammaHLTRechitInRegionsProducer.h.

References nEvt_, and nMaxPrintout_.

{ return ((nEvt_ > nMaxPrintout_) || (nMaxPrintout_ < 0));}

void EgammaHLTRechitInRegionsProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 88 of file EgammaHLTRechitInRegionsProducer.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                                  {
  edm::ParameterSetDescription desc;
  std::vector<edm::InputTag> inputTags;
  inputTags.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEB"));
  inputTags.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEE"));
  inputTags.push_back(edm::InputTag("hltESRegionalEgammaRecHit:EcalRecHitsES"));
  desc.add<std::vector<edm::InputTag>>("ecalhitLabels", inputTags);
  desc.add<edm::InputTag>("ecalhitproducer", edm::InputTag("ecalRecHit"));
  desc.add<edm::InputTag>("l1TagIsolated", edm::InputTag("l1extraParticles","Isolated"));
  desc.add<edm::InputTag>("l1TagNonIsolated", edm::InputTag("l1extraParticles","NonIsolated"));
  desc.add<bool>("doIsolated", true);
  desc.add<double>("l1LowerThr", 5.0);
  desc.add<double>("l1UpperThr", 999.);
  desc.add<double>("l1LowerThrIgnoreIsolation", 0.0);
  desc.add<double>("regionEtaMargin", 0.14);
  desc.add<double>("regionPhiMargin", 0.4);
  desc.add<std::vector<std::string> >("RecHitFlagToBeExcluded", std::vector<std::string>());
  desc.add<std::vector<std::string> >("RecHitSeverityToBeExcluded", std::vector<std::string>());
  descriptions.add(("hltEgammaHLTRechitInRegionsProducer"), desc);  
}

void EgammaHLTRechitInRegionsProducer::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

virtual

Implements edm::EDProducer.

Definition at line 109 of file EgammaHLTRechitInRegionsProducer.cc.

References edm::SortedCollection< T, SORT >::begin(), doIsolated_, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, alcaDQMUpload::encode(), edm::SortedCollection< T, SORT >::end(), geometry, edm::EventSetup::get(), edm::Event::getByLabel(), edm::Event::getByToken(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), hitLabels, hitTokens, i, l1LowerThr_, l1LowerThrIgnoreIsolation_, l1TagIsolated_, l1TagNonIsolated_, l1UpperThr_, nEvt_, or, position, edm::Handle< T >::product(), edm::Event::put(), regionEtaMargin_, and regionPhiMargin_.

                                                                                     {

  std::auto_ptr<EcalRecHitCollection> hitsEB(new EcalRecHitCollection);
  std::auto_ptr<EcalRecHitCollection> hitsEE(new EcalRecHitCollection);
  std::auto_ptr<EcalRecHitCollection> hitsES(new EcalRecHitCollection);
  
  edm::Handle<EcalRecHitCollection> rhcH[3];

  for (unsigned int i=0; i<hitLabels.size(); i++) {
    evt.getByToken(hitTokens[i], rhcH[i]);  
    if (!(rhcH[i].isValid())) {
      edm::LogError("ProductNotFound")<< "could not get a handle on the EcalRecHitCollection! (" << hitLabels[i].encode() << ")" << std::endl;
      return;
    }
  }
  
  // get the collection geometry:
  edm::ESHandle<CaloGeometry> geoHandle;
  es.get<CaloGeometryRecord>().get(geoHandle);
  const CaloGeometry& geometry = *geoHandle;
  const CaloSubdetectorGeometry *geometry_p;
  std::auto_ptr<const CaloSubdetectorTopology> topology;

  //edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
  //es.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
  //const EcalSeverityLevelAlgo* sevLevel = sevlv.product();
  
  //Get the L1 EM Particle Collection
  //Get the L1 EM Particle Collection
  edm::Handle< l1extra::L1EmParticleCollection > emIsolColl ;
  if(doIsolated_) 
    evt.getByLabel(l1TagIsolated_, emIsolColl);

  //Get the L1 EM Particle Collection
  edm::Handle< l1extra::L1EmParticleCollection > emNonIsolColl ;
  evt.getByLabel(l1TagNonIsolated_, emNonIsolColl);
  
  // Get the CaloGeometry
  edm::ESHandle<L1CaloGeometry> l1CaloGeom ;
  es.get<L1CaloGeometryRecord>().get(l1CaloGeom) ;

  std::vector<EcalEtaPhiRegion> regions;

  if(doIsolated_) {
    for( l1extra::L1EmParticleCollection::const_iterator emItr = emIsolColl->begin(); emItr != emIsolColl->end() ;++emItr ) {
      if ((emItr->et() > l1LowerThr_) and (emItr->et() < l1UpperThr_)) {

    // Access the GCT hardware object corresponding to the L1Extra EM object.
    int etaIndex = emItr->gctEmCand()->etaIndex();
    int phiIndex = emItr->gctEmCand()->phiIndex();

    // Use the L1CaloGeometry to find the eta, phi bin boundaries.
    double etaLow  = l1CaloGeom->etaBinLowEdge(etaIndex);
    double etaHigh = l1CaloGeom->etaBinHighEdge(etaIndex);
    double phiLow  = l1CaloGeom->emJetPhiBinLowEdge( phiIndex ) ;
    double phiHigh = l1CaloGeom->emJetPhiBinHighEdge( phiIndex ) ;

    etaLow -= regionEtaMargin_;
    etaHigh += regionEtaMargin_;
    phiLow -= regionPhiMargin_;
    phiHigh += regionPhiMargin_;

    regions.push_back(EcalEtaPhiRegion(etaLow,etaHigh,phiLow,phiHigh));
      }
    }
  }

  if(!doIsolated_ or (l1LowerThrIgnoreIsolation_ < 64)) {
    for( l1extra::L1EmParticleCollection::const_iterator emItr = emNonIsolColl->begin(); emItr != emNonIsolColl->end() ;++emItr ) {
      
      if(doIsolated_ and (emItr->et() < l1LowerThrIgnoreIsolation_)) 
    continue;

      if ((emItr->et() > l1LowerThr_) and (emItr->et() < l1UpperThr_)) {
    
    // Access the GCT hardware object corresponding to the L1Extra EM object.
    int etaIndex = emItr->gctEmCand()->etaIndex();
    int phiIndex = emItr->gctEmCand()->phiIndex();
    
    // Use the L1CaloGeometry to find the eta, phi bin boundaries.
    double etaLow  = l1CaloGeom->etaBinLowEdge(etaIndex);
    double etaHigh = l1CaloGeom->etaBinHighEdge(etaIndex);
    double phiLow  = l1CaloGeom->emJetPhiBinLowEdge(phiIndex);
    double phiHigh = l1CaloGeom->emJetPhiBinHighEdge(phiIndex);

    etaLow -= regionEtaMargin_;
    etaHigh += regionEtaMargin_;
    phiLow -= regionPhiMargin_;
    phiHigh += regionPhiMargin_;
    
    regions.push_back(EcalEtaPhiRegion(etaLow,etaHigh,phiLow,phiHigh));
      }
    }
  }
  
  for (unsigned int i=0; i<hitLabels.size(); i++) {
    const EcalRecHitCollection *recHits = rhcH[i].product();
    if(hitLabels[i].encode() == "hltEcalRegionalEgammaRecHit:EcalRecHitsEB") {
      geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
      topology.reset(new EcalBarrelTopology(geoHandle));
    } else if(hitLabels[i].encode() == "hltEcalRegionalEgammaRecHit:EcalRecHitsEE") {
      geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
      topology.reset(new EcalEndcapTopology(geoHandle));
    } else if(hitLabels[i].encode() == "hltESRegionalEgammaRecHit:EcalRecHitsES") {
      geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
      topology.reset(new EcalPreshowerTopology (geoHandle));
    } else throw(std::runtime_error("\n\nProducer encountered invalied ecalhitcollection type.\n\n"));

    if(regions.size() != 0) {
      EcalRecHitCollection::const_iterator it;

      for (it = recHits->begin(); it != recHits->end(); it++){
    //Make the vector of seeds that we're going to use.
    //One of the few places position is used, needed for ET calculation.    
    const CaloCellGeometry *this_cell = (*geometry_p).getGeometry(it->id());
    GlobalPoint position = this_cell->getPosition();

    std::vector<EcalEtaPhiRegion>::const_iterator region;
    for (region=regions.begin(); region!=regions.end(); region++) {
      //std::cout << region->etaLow() << " " << region->etaHigh() <<  " " << region->phiLow() << " " << region->phiHigh() << std::endl;
      //std::cout << i << " -  " << position.eta() << " " << position.phi() << std::endl;
      if (region->inRegion(position)) {
        if (i == 0)
          hitsEB->push_back(*it);
        if (i == 1)
          hitsEE->push_back(*it);
        if (i == 2)
          hitsES->push_back(*it);
      }
    }
      }
    }
  }

  //std::cout << hitsEB->size() << " " << hitsEE->size() << " " << hitsES->size() << std::endl;
  evt.put(hitsEB, "EcalRegionalRecHitsEB");
  evt.put(hitsEE, "EcalRegionalRecHitsEE");
  evt.put(hitsES, "EcalRegionalRecHitsES");

  nEvt_++;
}

Member Data Documentation

bool EgammaHLTRechitInRegionsProducer::doIsolated_

private

Definition at line 31 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

std::string EgammaHLTRechitInRegionsProducer::hitcollection_

private

Definition at line 33 of file EgammaHLTRechitInRegionsProducer.h.

std::vector<edm::InputTag> EgammaHLTRechitInRegionsProducer::hitLabels

private

Definition at line 45 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

edm::InputTag EgammaHLTRechitInRegionsProducer::hitproducer_

private

Definition at line 32 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer().

std::vector<edm::EDGetTokenT<EcalRecHitCollection> > EgammaHLTRechitInRegionsProducer::hitTokens

private

Definition at line 46 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

double EgammaHLTRechitInRegionsProducer::l1LowerThr_

private

Definition at line 38 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

double EgammaHLTRechitInRegionsProducer::l1LowerThrIgnoreIsolation_

private

Definition at line 40 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

edm::InputTag EgammaHLTRechitInRegionsProducer::l1TagIsolated_

private

Definition at line 35 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

edm::InputTag EgammaHLTRechitInRegionsProducer::l1TagNonIsolated_

private

Definition at line 36 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

double EgammaHLTRechitInRegionsProducer::l1UpperThr_

private

Definition at line 39 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

int EgammaHLTRechitInRegionsProducer::nEvt_

private

Definition at line 29 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by counterExceeded(), EgammaHLTRechitInRegionsProducer(), and produce().

int EgammaHLTRechitInRegionsProducer::nMaxPrintout_

private

Definition at line 28 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by counterExceeded().

double EgammaHLTRechitInRegionsProducer::regionEtaMargin_

private

Definition at line 42 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().

double EgammaHLTRechitInRegionsProducer::regionPhiMargin_

private

Definition at line 43 of file EgammaHLTRechitInRegionsProducer.h.

Referenced by EgammaHLTRechitInRegionsProducer(), and produce().