#include <CaloRecHitCandidateProducer.h>

Inheritance diagram for CaloRecHitCandidateProducer:

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

double	cellTresholdAndWeight (const CaloRecHit &, const HcalTopology &) const

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

	~CaloRecHitCandidateProducer ()

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 Attributes
bool	mAllowMissingInputs

double	mEBthreshold

double	mEBweight

std::vector< edm::InputTag >	mEcalLabels

double	mEEthreshold

double	mEEweight

edm::InputTag	mHBHELabel
	source collection tag More...

double	mHBthreshold

double	mHBweight

double	mHEDthreshold

double	mHEDweight

double	mHESthreshold

double	mHESweight

double	mHF1threshold

double	mHF1weight

double	mHF2threshold

double	mHF2weight

edm::InputTag	mHFLabel

edm::InputTag	mHOLabel

double	mHOthreshold

double	mHOweight

bool	mUseHO

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 13 of file CaloRecHitCandidateProducer.h.

Constructor & Destructor Documentation

CaloRecHitCandidateProducer::CaloRecHitCandidateProducer ( const edm::ParameterSet & fConfig )

Definition at line 67 of file CaloRecHitCandidateProducer.cc.

   :  mHBHELabel (fConfig.getParameter<edm::InputTag>("hbheInput")),
      mHOLabel (fConfig.getParameter<edm::InputTag>("hoInput")),
      mHFLabel (fConfig.getParameter<edm::InputTag>("hfInput")),
      mEcalLabels (fConfig.getParameter<std::vector<edm::InputTag> >("ecalInputs")),
      mAllowMissingInputs (fConfig.getUntrackedParameter<bool>("AllowMissingInputs",false)),
      mUseHO (fConfig.getParameter<bool>("UseHO")),
 
      mEBthreshold (fConfig.getParameter<double>("EBThreshold")),
      mEEthreshold  (fConfig.getParameter<double>("EEThreshold")),
      mHBthreshold  (fConfig.getParameter<double>("HBThreshold")),
      mHESthreshold  (fConfig.getParameter<double>("HESThreshold")),
      mHEDthreshold  (fConfig.getParameter<double>("HEDThreshold")),
      mHOthreshold (fConfig.getParameter<double>("HOThreshold")),
      mHF1threshold (fConfig.getParameter<double>("HF1Threshold")),
      mHF2threshold (fConfig.getParameter<double>("HF2Threshold")),
      mEBweight (fConfig.getParameter<double>("EBWeight")),
      mEEweight (fConfig.getParameter<double>("EEWeight")),
      mHBweight (fConfig.getParameter<double>("HBWeight")),
      mHESweight (fConfig.getParameter<double>("HESWeight")),
      mHEDweight (fConfig.getParameter<double>("HEDWeight")),
      mHOweight (fConfig.getParameter<double>("HOWeight")),
      mHF1weight (fConfig.getParameter<double>("HF1Weight")),
      mHF2weight (fConfig.getParameter<double>("HF2Weight"))
 {
   produces<CandidateCollection>();
 }

CaloRecHitCandidateProducer::~CaloRecHitCandidateProducer ( )

inline

Definition at line 16 of file CaloRecHitCandidateProducer.h.

16 { }

Member Function Documentation

double CaloRecHitCandidateProducer::cellTresholdAndWeight	(	const CaloRecHit &	fHit,
		const HcalTopology &	fTopology
	)		const

Definition at line 158 of file CaloRecHitCandidateProducer.cc.

References HcalDetId::depth(), DetId::det(), CaloRecHit::detid(), DetId::Ecal, EcalBarrel, EcalEndcap, CaloRecHit::energy(), HcalTopology::firstHEDoublePhiRing(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HcalDetId::ietaAbs(), mEBthreshold, mEBweight, mEEthreshold, mEEweight, mHBthreshold, mHBweight, mHEDthreshold, mHEDweight, mHESthreshold, mHESweight, mHF1threshold, mHF1weight, mHF2threshold, mHF2weight, mHOthreshold, mHOweight, HcalDetId::subdet(), DetId::subdetId(), and crabWrap::threshold.

                                                                                                                       {
   double weight = 0;
   double threshold = 0;
   DetId detId = fHit.detid ();
   DetId::Detector det = detId.det ();
   if(det == DetId::Ecal) {
     // may or may not be EB.  We'll find out.
     
     EcalSubdetector subdet = (EcalSubdetector)(detId.subdetId());
     if(subdet == EcalBarrel) {
       threshold = mEBthreshold;
       weight = mEBweight;
     }
     else if(subdet == EcalEndcap) {
       threshold = mEEthreshold;
       weight = mEEweight;
     }
   }
   else if(det == DetId::Hcal) {
     HcalDetId hcalDetId(detId);
     HcalSubdetector subdet = hcalDetId.subdet();
     
     if(subdet == HcalBarrel) {
       threshold = mHBthreshold;
       weight = mHBweight;
     }
     
     else if(subdet == HcalEndcap) {
       // check if it's single or double tower
       if(hcalDetId.ietaAbs() < fTopology.firstHEDoublePhiRing()) {
         threshold = mHESthreshold;
         weight = mHESweight;
       }
       else {
         threshold = mHEDthreshold;
         weight = mHEDweight;
       }
     } else if(subdet == HcalOuter) {
       threshold = mHOthreshold;
       weight = mHOweight;
     } else if(subdet == HcalForward) {
       if(hcalDetId.depth() == 1) {
         threshold = mHF1threshold;
         weight = mHF1weight;
       } else {
         threshold = mHF2threshold;
         weight = mHF2weight;
       }
     }
   }
   return fHit.energy () >= threshold ? weight : 0; 
 }

void CaloRecHitCandidateProducer::produce	(	edm::Event &	fEvent,
		const edm::EventSetup &	fSetup
	)

virtual

Implements edm::EDProducer.

Definition at line 95 of file CaloRecHitCandidateProducer.cc.

References geometry, edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), edm::Event::getByLabel(), i, edm::HandleBase::isValid(), mAllowMissingInputs, mEcalLabels, mHBHELabel, mHFLabel, mHOLabel, mUseHO, convertSQLitetoXML_cfg::output, edm::Event::put(), and record.

                                                                                           {
   // get geometry
   //  const IdealGeometryRecord& record = fSetup.template get<IdealGeometryRecord>();
   const CaloGeometryRecord& caloRecord = fSetup.get<CaloGeometryRecord>();
   ESHandle<CaloGeometry> geometry;
   caloRecord.get (geometry);
   const IdealGeometryRecord& record = fSetup.get<IdealGeometryRecord>();
   ESHandle<HcalTopology> topology;
   record.get (topology);
   // set Output
   auto_ptr<CandidateCollection> output ( new CandidateCollection );
   // get and process Inputs
   edm::Handle<HBHERecHitCollection> hbhe;
   fEvent.getByLabel(mHBHELabel,hbhe);
   if (!hbhe.isValid()) {
     // can't find it!
     if (!mAllowMissingInputs) {
       *hbhe;  // will throw the proper exception
     }
   } else {
     processHits (hbhe, *this, *geometry, *topology, &*output);
   }
 
   if (mUseHO) {
     edm::Handle<HORecHitCollection> ho;
     fEvent.getByLabel(mHOLabel,ho);
     if (!ho.isValid()) {
       // can't find it!
       if (!mAllowMissingInputs) {
     *ho;  // will throw the proper exception    
       }
     } else {
       processHits (ho, *this, *geometry, *topology, &*output);
     }
   }
 
   edm::Handle<HFRecHitCollection> hf;
   fEvent.getByLabel(mHFLabel,hf);
   if (!hf.isValid()) {
     // can't find it!
     if (!mAllowMissingInputs) {
       *hf;  // will throw the proper exception
     }
   } else {
     processHits (hf, *this, *geometry, *topology, &*output);
   }
 
   std::vector<edm::InputTag>::const_iterator i;
   for (i=mEcalLabels.begin(); i!=mEcalLabels.end(); i++) {
     edm::Handle<EcalRecHitCollection> ec;
     fEvent.getByLabel(*i,ec);
     if (!ec.isValid()) {
       // can't find it!
       if (!mAllowMissingInputs) {
     *ec;  // will throw the proper exception
       }
     } else {
       processHits (ec, *this, *geometry, *topology, &*output);
     }
   }
   fEvent.put(output);
 }