#include <ZdcSimpleReconstructor.h>

Inheritance diagram for ZdcSimpleReconstructor:

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

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

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

	ZdcSimpleReconstructor (const edm::ParameterSet &ps)

	~ZdcSimpleReconstructor () override

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
DetId::Detector	det_

bool	dropZSmarkedPassed_

const HcalTimeSlew *	hcalTimeSlew_delay_

HcalLongRecoParams *	myobject

ZdcSimpleRecAlgo	reco_

int	subdet_

HcalOtherSubdetector	subdetOther_

edm::EDGetTokenT< ZDCDigiCollection >	tok_input_castor

edm::EDGetTokenT< ZDCDigiCollection >	tok_input_hcal

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer< T >
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

Detailed Description

Definition at line 23 of file ZdcSimpleReconstructor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file ZdcSimpleReconstructor.cc.

References DetId::Calo, gather_cfg::cout, det_, edm::ParameterSet::getParameter(), HcalCalibration, HcalOther, hcalTimeSlew_delay_, AlCaHLTBitMon_QueryRunRegistry::string, subdet_, HcalZDCDetId::SubdetectorId, subdetOther_, tok_input_castor, and tok_input_hcal.

     : reco_(conf.getParameter<bool>("correctForTimeslew"),
             conf.getParameter<bool>("correctForPhaseContainment"),
             conf.getParameter<double>("correctionPhaseNS"),
             conf.getParameter<int>("recoMethod"),
             conf.getParameter<int>("lowGainOffset"),
             conf.getParameter<double>("lowGainFrac")),
       det_(DetId::Hcal),
       dropZSmarkedPassed_(conf.getParameter<bool>("dropZSmarkedPassed")) {
   tok_input_castor = consumes<ZDCDigiCollection>(conf.getParameter<edm::InputTag>("digiLabelcastor"));
   tok_input_hcal = consumes<ZDCDigiCollection>(conf.getParameter<edm::InputTag>("digiLabelhcal"));
 
   std::string subd = conf.getParameter<std::string>("Subdetector");
   if (!strcasecmp(subd.c_str(), "ZDC")) {
     det_ = DetId::Calo;
     subdet_ = HcalZDCDetId::SubdetectorId;
     produces<ZDCRecHitCollection>();
   } else if (!strcasecmp(subd.c_str(), "CALIB")) {
     subdet_ = HcalOther;
     subdetOther_ = HcalCalibration;
     produces<HcalCalibRecHitCollection>();
   } else {
     std::cout << "ZdcSimpleReconstructor is not associated with a specific subdetector!" << std::endl;
   }
 
   hcalTimeSlew_delay_ = nullptr;
 }

ZdcSimpleReconstructor::~ZdcSimpleReconstructor ( )

override

Definition at line 41 of file ZdcSimpleReconstructor.cc.

41 {}

Member Function Documentation

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

final

Definition at line 42 of file ZdcSimpleReconstructor.cc.

References edm::EventSetup::get(), hcalTimeSlew_delay_, myobject, AlCaHLTBitMon_ParallelJobs::p, and edm::ESHandle< T >::product().

                                                                               {
   edm::ESHandle<HcalLongRecoParams> p;
   es.get<HcalLongRecoParamsRcd>().get(p);
   myobject = new HcalLongRecoParams(*p.product());
 
   edm::ESHandle<HcalTimeSlew> delay;
   es.get<HcalTimeSlewRecord>().get("HBHE", delay);
   hcalTimeSlew_delay_ = &*delay;
 }

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

final

Definition at line 52 of file ZdcSimpleReconstructor.cc.

References myobject.

                                                                             {
   delete myobject;
   myobject = nullptr;
 }

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

override

Definition at line 56 of file ZdcSimpleReconstructor.cc.

References edm::SortedCollection< T, SORT >::begin(), AlignmentProducer_cff::calibrations, DetId::Calo, det_, dropZSmarkedPassed_, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), edm::EventSetup::get(), edm::Event::getByToken(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalCondObjectContainer< Item >::getValues(), hcalTimeSlew_delay_, mps_fire::i, ZdcSimpleRecAlgo::initPulseCorr(), eostools::move(), myobject, HcalLongRecoParam::noiseTS(), edm::Event::put(), reco_, ZdcSimpleRecAlgo::reconstruct(), HcalLongRecoParam::signalTS(), edm::SortedCollection< T, SORT >::size(), subdet_, HcalZDCDetId::SubdetectorId, tok_input_castor, and tok_input_hcal.

                                                                                  {
   // get conditions
   edm::ESHandle<HcalDbService> conditions;
   eventSetup.get<HcalDbRecord>().get(conditions);
   // define vectors to pass noiseTS and signalTS
   std::vector<unsigned int> mySignalTS;
   std::vector<unsigned int> myNoiseTS;
 
   if (det_ == DetId::Calo && subdet_ == HcalZDCDetId::SubdetectorId) {
     edm::Handle<ZDCDigiCollection> digi;
     e.getByToken(tok_input_hcal, digi);
 
     if (digi->empty()) {
       e.getByToken(tok_input_castor, digi);
       if (digi->empty())
         edm::LogInfo("ZdcHitReconstructor") << "No ZDC info found in either castorDigis or hcalDigis." << std::endl;
     }
 
     // create empty output
     auto rec = std::make_unique<ZDCRecHitCollection>();
     rec->reserve(digi->size());
     // run the algorithm
     unsigned int toaddMem = 0;
 
     ZDCDigiCollection::const_iterator i;
     for (i = digi->begin(); i != digi->end(); i++) {
       HcalZDCDetId cell = i->id();
       DetId detcell = (DetId)cell;
       // rof 27.03.09: drop ZS marked and passed digis:
       if (dropZSmarkedPassed_)
         if (i->zsMarkAndPass())
           continue;
 
       // get db values for signalTSs and noiseTSs
       const HcalLongRecoParam* myParams = myobject->getValues(detcell);
       mySignalTS.clear();
       myNoiseTS.clear();
       mySignalTS = myParams->signalTS();
       myNoiseTS = myParams->noiseTS();
       // warning: the PulseCorrection is not used by ZDC. If it gets a non-contingious set of
       // signal TS, it may not work properly. Assume contiguous here....
       unsigned int toadd = mySignalTS.size();
       if (toaddMem != toadd) {
         reco_.initPulseCorr(toadd, hcalTimeSlew_delay_);
         toaddMem = toadd;
       }
       const HcalCalibrations& calibrations = conditions->getHcalCalibrations(cell);
       const HcalQIECoder* channelCoder = conditions->getHcalCoder(cell);
       const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
       HcalCoderDb coder(*channelCoder, *shape);
       rec->push_back(reco_.reconstruct(*i, myNoiseTS, mySignalTS, coder, calibrations));
     }
     // return result
     e.put(std::move(rec));
   }
 }