#include <HcalSimpleReconstructor.h>

Inheritance diagram for HcalSimpleReconstructor:

Public Member Functions
virtual void	beginRun (edm::Run const &r, edm::EventSetup const &es) overridefinal

virtual void	endRun (edm::Run const &r, edm::EventSetup const &es) overridefinal

	HcalSimpleReconstructor (const edm::ParameterSet &ps)

virtual void	produce (edm::Event &e, const edm::EventSetup &c) overridefinal

virtual	~HcalSimpleReconstructor ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	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
std::vector< ConsumesInfo >	consumesInfo () const

	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
template<class DIGICOLL , class RECHITCOLL >
void	process (edm::Event &e, const edm::EventSetup &c, const edm::EDGetTokenT< DIGICOLL > &tok)

void	processUpgrade (edm::Event &e, const edm::EventSetup &c)

Private Attributes
DetId::Detector	det_

bool	dropZSmarkedPassed_

int	firstSample_

bool	HFQIE10_

edm::InputTag	inputLabel_

HcalRecoParams *	paramTS

HcalSimpleRecAlgo	reco_

int	samplesToAdd_

int	subdet_

HcalOtherSubdetector	subdetOther_

HcalTopology *	theTopology

edm::EDGetTokenT < HcalCalibDigiCollection >	tok_calib_

edm::EDGetTokenT < HBHEDigiCollection >	tok_hbhe_

edm::EDGetTokenT < HBHEUpgradeDigiCollection >	tok_hbheUp_

edm::EDGetTokenT < HFDigiCollection >	tok_hf_

edm::EDGetTokenT < QIE10DigiCollection >	tok_hfQIE10_

edm::EDGetTokenT < HFUpgradeDigiCollection >	tok_hfUp_

edm::EDGetTokenT < HODigiCollection >	tok_ho_

bool	tsFromDB_

bool	upgradeHBHE_

bool	upgradeHF_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Author: J. Mans - Minnesota; E. Garcia - CSU Based on HcalSimpleReconstructor.h by J. Mans

Definition at line 26 of file HcalSimpleReconstructor.h.

Constructor & Destructor Documentation

HcalSimpleReconstructor::HcalSimpleReconstructor ( const edm::ParameterSet & ps )

explicit

Definition at line 16 of file HcalSimpleReconstructor.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), HcalBarrel, HcalForward, HcalOuter, HFQIE10_, inputLabel_, reco_, HcalSimpleRecAlgo::setMeth3Params(), AlCaHLTBitMon_QueryRunRegistry::string, subdet_, tok_calib_, tok_hbhe_, tok_hbheUp_, tok_hf_, tok_hfQIE10_, tok_hfUp_, tok_ho_, upgradeHBHE_, and upgradeHF_.

                                                                            :
   reco_(conf.getParameter<bool>("correctForTimeslew"),
     conf.getParameter<bool>("correctForPhaseContainment"),conf.getParameter<double>("correctionPhaseNS")),
   det_(DetId::Hcal),
   inputLabel_(conf.getParameter<edm::InputTag>("digiLabel")),
   dropZSmarkedPassed_(conf.getParameter<bool>("dropZSmarkedPassed")),
   firstSample_(conf.getParameter<int>("firstSample")),
   samplesToAdd_(conf.getParameter<int>("samplesToAdd")),
   tsFromDB_(conf.getParameter<bool>("tsFromDB")),
   upgradeHBHE_(false),
   upgradeHF_(false),
   HFQIE10_(false),
   paramTS(0),
   theTopology(0)
 {
   // Intitialize "method 3"
   reco_.setMeth3Params(
             conf.getParameter<int>     ("pedestalSubtractionType"),
             conf.getParameter<double>  ("pedestalUpperLimit"),
             conf.getParameter<int>     ("timeSlewParsType"),
             conf.getParameter<std::vector<double> >("timeSlewPars"),
             conf.getParameter<double>  ("respCorrM3")
       );
 
   // register for data access
   tok_hbheUp_ = consumes<HBHEUpgradeDigiCollection>(inputLabel_);
   tok_hfUp_ = consumes<HFUpgradeDigiCollection>(inputLabel_);
   tok_hfQIE10_ = consumes<QIE10DigiCollection>(inputLabel_);
 
   tok_hbhe_ = consumes<HBHEDigiCollection>(inputLabel_);
   tok_hf_ = consumes<HFDigiCollection>(inputLabel_);
   tok_ho_ = consumes<HODigiCollection>(inputLabel_);
   tok_calib_ = consumes<HcalCalibDigiCollection>(inputLabel_);
 
   std::string subd=conf.getParameter<std::string>("Subdetector");
   if(!strcasecmp(subd.c_str(),"upgradeHBHE")) {
      upgradeHBHE_ = true;
      produces<HBHERecHitCollection>();
   }
   else if (!strcasecmp(subd.c_str(),"upgradeHF")) {
      upgradeHF_ = true;
      produces<HFRecHitCollection>();
   }
   else if (!strcasecmp(subd.c_str(),"HFQIE10")) {
      HFQIE10_ = true;
      produces<HFRecHitCollection>();
   }
   else if (!strcasecmp(subd.c_str(),"HO")) {
     subdet_=HcalOuter;
     produces<HORecHitCollection>();
   }  
   else if (!strcasecmp(subd.c_str(),"HBHE")) {
     if( !upgradeHBHE_) {
       subdet_=HcalBarrel;
       produces<HBHERecHitCollection>();
     }
   } 
   else if (!strcasecmp(subd.c_str(),"HF")) {
     if( !upgradeHF_) {
     subdet_=HcalForward;
     produces<HFRecHitCollection>();
     }
   } 
   else {
     std::cout << "HcalSimpleReconstructor is not associated with a specific subdetector!" << std::endl;
   }       
   
 }

HcalSimpleReconstructor::~HcalSimpleReconstructor ( )

virtual

Definition at line 96 of file HcalSimpleReconstructor.cc.

References paramTS, and theTopology.

                                                   { 
   delete paramTS;
   delete theTopology;
 }

Member Function Documentation

void HcalSimpleReconstructor::beginRun	(	edm::Run const &	r,
		edm::EventSetup const &	es
	)

finaloverridevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 101 of file HcalSimpleReconstructor.cc.

References HcalSimpleRecAlgo::beginRun(), edm::EventSetup::get(), AlCaHLTBitMon_ParallelJobs::p, paramTS, edm::ESHandle< class >::product(), reco_, HcalCondObjectContainerBase::setTopo(), theTopology, and tsFromDB_.

                                                                               {
   if(tsFromDB_) {
     edm::ESHandle<HcalRecoParams> p;
     es.get<HcalRecoParamsRcd>().get(p);
     paramTS = new HcalRecoParams(*p.product());
 
     edm::ESHandle<HcalTopology> htopo;
     es.get<HcalRecNumberingRecord>().get(htopo);
     theTopology=new HcalTopology(*htopo);
     paramTS->setTopo(theTopology);
 
   }
   reco_.beginRun(es);
 }

void HcalSimpleReconstructor::endRun	(	edm::Run const &	r,
		edm::EventSetup const &	es
	)

finaloverridevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 116 of file HcalSimpleReconstructor.cc.

References HcalSimpleRecAlgo::endRun(), paramTS, reco_, and tsFromDB_.

                                                                             {
   if(tsFromDB_ && paramTS) {
     delete paramTS;
     paramTS = 0;
     reco_.endRun();
   }
 }

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

static

Definition at line 85 of file HcalSimpleReconstructor.cc.

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

                                                                                          {
   edm::ParameterSetDescription desc;
   desc.setAllowAnything();
   desc.add<int>("pedestalSubtractionType", 1); 
   desc.add<double>("pedestalUpperLimit", 2.7); 
   desc.add<int>("timeSlewParsType",3);
   desc.add<std::vector<double>>("timeSlewPars", { 12.2999, -2.19142, 0, 12.2999, -2.19142, 0, 12.2999, -2.19142, 0 });
   desc.add<double>("respCorrM3", 0.95);
   descriptions.add("simpleHbhereco",desc);
 }

template<class DIGICOLL , class RECHITCOLL >

void HcalSimpleReconstructor::process	(	edm::Event &	e,
		const edm::EventSetup &	c,
		const edm::EDGetTokenT< DIGICOLL > &	tok
	)

private

Definition at line 126 of file HcalSimpleReconstructor.cc.

References dropZSmarkedPassed_, plotBeamSpotDB::first, HcalRecoParam::firstSample(), firstSample_, edm::EventSetup::get(), edm::Event::getByToken(), HcalCondObjectContainer< Item >::getValues(), i, paramTS, edm::Event::put(), DetId::rawId(), reco_, HcalSimpleRecAlgo::reconstruct(), HcalRecoParam::samplesToAdd(), samplesToAdd_, and tsFromDB_.

 {
   // get conditions
   edm::ESHandle<HcalDbService> conditions;
   eventSetup.get<HcalDbRecord>().get(conditions);
 
   edm::Handle<DIGICOLL> digi;
   e.getByToken(tok,digi);
 
   // create empty output
   std::auto_ptr<RECHITCOLL> rec(new RECHITCOLL);
   rec->reserve(digi->size());
   // run the algorithm
   int first = firstSample_;
   int toadd = samplesToAdd_;
   typename DIGICOLL::const_iterator i;
   for (i=digi->begin(); i!=digi->end(); i++) {
     HcalDetId cell = i->id();
     DetId detcell=(DetId)cell;
     // rof 27.03.09: drop ZS marked and passed digis:
     if (dropZSmarkedPassed_)
       if (i->zsMarkAndPass()) continue;
 
     const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
     const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
     const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); 
     HcalCoderDb coder (*channelCoder, *shape);
 
     //>>> firstSample & samplesToAdd
     if(tsFromDB_) {
       const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
       first = param_ts->firstSample();
       toadd = param_ts->samplesToAdd();
     }
     rec->push_back(reco_.reconstruct(*i,first,toadd,coder,calibrations));   
   }
   // return result
   e.put(rec);
 }

void HcalSimpleReconstructor::processUpgrade	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

private

Definition at line 167 of file HcalSimpleReconstructor.cc.

References dropZSmarkedPassed_, plotBeamSpotDB::first, HcalRecoParam::firstSample(), firstSample_, edm::EventSetup::get(), edm::Event::getByToken(), HcalCondObjectContainer< Item >::getValues(), HFQIE10_, i, paramTS, edm::Event::put(), DetId::rawId(), reco_, HcalSimpleRecAlgo::reconstructHBHEUpgrade(), HcalSimpleRecAlgo::reconstructHFUpgrade(), HcalSimpleRecAlgo::reconstructQIE10(), HcalRecoParam::samplesToAdd(), samplesToAdd_, tok_hbheUp_, tok_hfQIE10_, tok_hfUp_, tsFromDB_, upgradeHBHE_, upgradeHF_, and QIE10DataFrame::zsMarkAndPass().

Referenced by produce().

 {
   // get conditions
   edm::ESHandle<HcalDbService> conditions;
   eventSetup.get<HcalDbRecord>().get(conditions);
 
   if(upgradeHBHE_){
    
     edm::Handle<HBHEUpgradeDigiCollection> digi;
     e.getByToken(tok_hbheUp_, digi);
 
     // create empty output
     std::auto_ptr<HBHERecHitCollection> rec(new HBHERecHitCollection);
     rec->reserve(digi->size()); 
 
     // run the algorithm
     int first = firstSample_;
     int toadd = samplesToAdd_;
     HBHEUpgradeDigiCollection::const_iterator i;
     for (i=digi->begin(); i!=digi->end(); i++) {
       HcalDetId cell = i->id();
       DetId detcell=(DetId)cell;
       // rof 27.03.09: drop ZS marked and passed digis:
       if (dropZSmarkedPassed_)
       if (i->zsMarkAndPass()) continue;
       
       const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
       const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
       const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); 
       HcalCoderDb coder (*channelCoder, *shape);
 
       //>>> firstSample & samplesToAdd
       if(tsFromDB_) {
     const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
     first = param_ts->firstSample();
     toadd = param_ts->samplesToAdd();
       }
       rec->push_back(reco_.reconstructHBHEUpgrade(*i,first,toadd,coder,calibrations));
 
     }
 
     e.put(rec); // put results
   }// End of upgradeHBHE
 
   if(upgradeHF_){
 
     edm::Handle<HFUpgradeDigiCollection> digi;
     e.getByToken(tok_hfUp_, digi);
 
     // create empty output
     std::auto_ptr<HFRecHitCollection> rec(new HFRecHitCollection);
     rec->reserve(digi->size()); 
 
     // run the algorithm
     int first = firstSample_;
     int toadd = samplesToAdd_;
     HFUpgradeDigiCollection::const_iterator i;
     for (i=digi->begin(); i!=digi->end(); i++) {
       HcalDetId cell = i->id();
       DetId detcell=(DetId)cell;
       // rof 27.03.09: drop ZS marked and passed digis:
       if (dropZSmarkedPassed_)
       if (i->zsMarkAndPass()) continue;
       
       const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
       const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
       const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); 
       HcalCoderDb coder (*channelCoder, *shape);
 
       //>>> firstSample & samplesToAdd
       if(tsFromDB_) {
     const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
     first = param_ts->firstSample();
     toadd = param_ts->samplesToAdd();
       }
       rec->push_back(reco_.reconstructHFUpgrade(*i,first,toadd,coder,calibrations));
 
     }  
     e.put(rec); // put results
   }// End of upgradeHF
 
   if(HFQIE10_){
 
     edm::Handle<QIE10DigiCollection> digi;
     e.getByToken(tok_hfQIE10_, digi);
 
     // create empty output
     std::auto_ptr<HFRecHitCollection> rec(new HFRecHitCollection);
     rec->reserve(digi->size()); 
 
     // run the algorithm
     int first = firstSample_;
     int toadd = samplesToAdd_;
     for (QIE10DigiCollection::const_iterator i=digi->begin(); i!=digi->end(); i++) {
       HcalDetId cell = i->id();
       DetId detcell=(DetId)cell;
       
       //make dataframe
       QIE10DataFrame frame(*i);
       
       // rof 27.03.09: drop ZS marked and passed digis:
       if (dropZSmarkedPassed_){
         if (frame.zsMarkAndPass()) continue;
       }
 
       const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
       const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
       const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); 
       HcalCoderDb coder (*channelCoder, *shape);
 
       //>>> firstSample & samplesToAdd
       if(tsFromDB_) {
     const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
     first = param_ts->firstSample();
     toadd = param_ts->samplesToAdd();
       }
       rec->push_back(reco_.reconstructQIE10(frame,first,toadd,coder,calibrations));
 
     }  
     e.put(rec); // put results
   }// End of upgradeHF
   
 }

void HcalSimpleReconstructor::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

finaloverridevirtual

Implements edm::stream::EDProducerBase.

Definition at line 293 of file HcalSimpleReconstructor.cc.

References det_, alignCSCRings::e, DetId::Hcal, HcalBarrel, HcalCalibration, HcalEndcap, HcalForward, HcalOther, HcalOuter, HFQIE10_, edm::EventBase::isRealData(), processUpgrade(), reco_, edm::Event::run(), HcalSimpleRecAlgo::setForData(), subdet_, subdetOther_, tok_calib_, tok_hbhe_, tok_hf_, tok_ho_, upgradeHBHE_, and upgradeHF_.

 {
   // HACK related to HB- corrections
   if(e.isRealData()) reco_.setForData(e.run());
  
   // What to produce, better to avoid the same subdet Upgrade and regular 
   // rechits "clashes"
   if(upgradeHBHE_ || upgradeHF_ || HFQIE10_) {
       processUpgrade(e, eventSetup);
   } else if (det_==DetId::Hcal) {
     if ((subdet_==HcalBarrel || subdet_==HcalEndcap) && !upgradeHBHE_) {
       process<HBHEDigiCollection, HBHERecHitCollection>(e, eventSetup, tok_hbhe_);
     } else if (subdet_==HcalForward && !upgradeHF_) {
       process<HFDigiCollection, HFRecHitCollection>(e, eventSetup, tok_hf_);
     } else if (subdet_==HcalOuter) {
       process<HODigiCollection, HORecHitCollection>(e, eventSetup, tok_ho_);
     } else if (subdet_==HcalOther && subdetOther_==HcalCalibration) {
       process<HcalCalibDigiCollection, HcalCalibRecHitCollection>(e, eventSetup, tok_calib_);
     }
   } 
 }