#include <HcalSimpleReconstructor.h>

Inheritance diagram for HcalSimpleReconstructor:

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

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

	HcalSimpleReconstructor (const edm::ParameterSet &ps)

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

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)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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)

Private Attributes
DetId::Detector	det_

bool	dropZSmarkedPassed_

int	firstSample_

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< HFDigiCollection >	tok_hf_

edm::EDGetTokenT< HODigiCollection >	tok_ho_

bool	tsFromDB_

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, inputLabel_, reco_, HcalSimpleRecAlgo::setMeth3Params(), AlCaHLTBitMon_QueryRunRegistry::string, subdet_, tok_calib_, tok_hbhe_, tok_hf_, and tok_ho_.

                                                                            :
   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")),
   paramTS(0),
   theTopology(0)
 {
   // Intitialize "method 3"
   reco_.setMeth3Params(
         conf.getParameter<bool>    ("applyTimeSlewM3"),
             conf.getParameter<double>  ("pedestalUpperLimit"),
             conf.getParameter<int>     ("timeSlewParsType"),
             conf.getParameter<std::vector<double> >("timeSlewPars"),
             conf.getParameter<double>  ("respCorrM3")
       );
 
   // register for data access
   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(),"HO")) {
     subdet_=HcalOuter;
     produces<HORecHitCollection>();
   }  
   else if (!strcasecmp(subd.c_str(),"HBHE")) {
     subdet_=HcalBarrel;
     produces<HBHERecHitCollection>();
   } 
   else if (!strcasecmp(subd.c_str(),"HF")) {
     subdet_=HcalForward;
     produces<HFRecHitCollection>();
   } 
   else {
     std::cout << "HcalSimpleReconstructor is not associated with a specific subdetector!" << std::endl;
   }       
   
 }

HcalSimpleReconstructor::~HcalSimpleReconstructor ( )

virtual

Definition at line 73 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 78 of file HcalSimpleReconstructor.cc.

References HcalSimpleRecAlgo::beginRun(), edm::EventSetup::get(), AlCaHLTBitMon_ParallelJobs::p, paramTS, edm::ESHandle< T >::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 93 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 62 of file HcalSimpleReconstructor.cc.

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

                                                                                          {
   edm::ParameterSetDescription desc;
   desc.setAllowAnything();
   desc.add<bool>("applyTimeSlewM3", true);
   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", 1.0);
   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 103 of file HcalSimpleReconstructor.cc.

References AlignmentProducer_cff::calibrations, dropZSmarkedPassed_, plotBeamSpotDB::first, HcalRecoParam::firstSample(), firstSample_, edm::EventSetup::get(), edm::Event::getByToken(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalCondObjectContainer< Item >::getValues(), mps_fire::i, eostools::move(), 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
   auto rec = std::make_unique<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(std::move(rec));
 }

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

finaloverridevirtual

Definition at line 143 of file HcalSimpleReconstructor.cc.

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

 {
   // HACK related to HB- corrections
   if(e.isRealData()) reco_.setForData(e.run());
  
   if (det_==DetId::Hcal) {
     if ((subdet_==HcalBarrel || subdet_==HcalEndcap)) {
       process<HBHEDigiCollection, HBHERecHitCollection>(e, eventSetup, tok_hbhe_);
     } else if (subdet_==HcalForward) {
       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_);
     }
   } 
 }