#include <ZdcHitReconstructor.h>

Inheritance diagram for ZdcHitReconstructor:

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) final

	ZdcHitReconstructor (const edm::ParameterSet &ps)

	~ZdcHitReconstructor () override

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
std::vector< int >	AuxTSvec_

edm::ESGetToken< HcalDbService, HcalDbRecord >	conditionsToken_

bool	correctTiming_

DetId::Detector	det_

bool	dropZSmarkedPassed_

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	htopoToken_

std::unique_ptr < HcalLongRecoParams >	longRecoParams_

double	lowGainFrac_

int	lowGainOffset_

edm::ESGetToken < HcalLongRecoParams, HcalLongRecoParamsRcd >	paramsToken_

edm::ESGetToken < HcalChannelQuality, HcalChannelQualityRcd >	qualToken_

ZdcSimpleRecAlgo	reco_

HcalADCSaturationFlag *	saturationFlagSetter_

bool	setHSCPFlags_

bool	setNoiseFlags_

bool	setSaturationFlags_

bool	setTimingTrustFlags_

edm::ESGetToken < HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd >	sevToken_

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<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

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

Definition at line 46 of file ZdcHitReconstructor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file ZdcHitReconstructor.cc.

References AuxTSvec_, DetId::Calo, conditionsToken_, gather_cfg::cout, det_, edm::ParameterSet::getParameter(), HcalCalibration, HcalOther, htopoToken_, paramsToken_, qualToken_, saturationFlagSetter_, setSaturationFlags_, sevToken_, 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")),
       saturationFlagSetter_(nullptr),
       det_(DetId::Hcal),
       correctTiming_(conf.getParameter<bool>("correctTiming")),
       setNoiseFlags_(conf.getParameter<bool>("setNoiseFlags")),
       setHSCPFlags_(conf.getParameter<bool>("setHSCPFlags")),
       setSaturationFlags_(conf.getParameter<bool>("setSaturationFlags")),
       setTimingTrustFlags_(conf.getParameter<bool>("setTimingTrustFlags")),
       dropZSmarkedPassed_(conf.getParameter<bool>("dropZSmarkedPassed")),
       AuxTSvec_(conf.getParameter<std::vector<int> >("AuxTSvec")) {
   tok_input_hcal = consumes<ZDCDigiCollection>(conf.getParameter<edm::InputTag>("digiLabelhcal"));
   tok_input_castor = consumes<ZDCDigiCollection>(conf.getParameter<edm::InputTag>("digiLabelcastor"));
 
   std::sort(AuxTSvec_.begin(), AuxTSvec_.end());  // sort vector in ascending TS order
   std::string subd = conf.getParameter<std::string>("Subdetector");
 
   if (setSaturationFlags_) {
     const edm::ParameterSet& pssat = conf.getParameter<edm::ParameterSet>("saturationParameters");
     saturationFlagSetter_ = new HcalADCSaturationFlag(pssat.getParameter<int>("maxADCvalue"));
   }
   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 << "ZdcHitReconstructor is not associated with a specific subdetector!" << std::endl;
   }
 
   // ES tokens
   htopoToken_ = esConsumes<HcalTopology, HcalRecNumberingRecord, edm::Transition::BeginRun>();
   paramsToken_ = esConsumes<HcalLongRecoParams, HcalLongRecoParamsRcd, edm::Transition::BeginRun>();
   conditionsToken_ = esConsumes<HcalDbService, HcalDbRecord>();
   qualToken_ = esConsumes<HcalChannelQuality, HcalChannelQualityRcd>(edm::ESInputTag("", "withTopo"));
   sevToken_ = esConsumes<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd>();
 }

ZdcHitReconstructor::~ZdcHitReconstructor ( )

override

Definition at line 63 of file ZdcHitReconstructor.cc.

References saturationFlagSetter_.

63 { delete saturationFlagSetter_; }

ZdcHitReconstructor::saturationFlagSetter_

HcalADCSaturationFlag * saturationFlagSetter_

Definition: ZdcHitReconstructor.h:56

Member Function Documentation

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

final

Definition at line 65 of file ZdcHitReconstructor.cc.

References edm::EventSetup::getData(), htopoToken_, longRecoParams_, AlCaHLTBitMon_ParallelJobs::p, and paramsToken_.

                                                                            {
   const HcalTopology& htopo = es.getData(htopoToken_);
   const HcalLongRecoParams& p = es.getData(paramsToken_);
   longRecoParams_ = std::make_unique<HcalLongRecoParams>(p);
   longRecoParams_->setTopo(&htopo);
 }

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

final

Definition at line 72 of file ZdcHitReconstructor.cc.

72 {}

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

final

Definition at line 74 of file ZdcHitReconstructor.cc.

References AuxTSvec_, DetId::Calo, submitPVValidationJobs::conditions, conditionsToken_, det_, HcalSeverityLevelComputer::dropChannel(), dropZSmarkedPassed_, edm::Event::getByToken(), edm::EventSetup::getData(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), mps_fire::i, edm::HandleBase::isValid(), longRecoParams_, eostools::move(), HcalLongRecoParam::noiseTS(), edm::Event::put(), qualToken_, DetId::rawId(), reco_, ZdcSimpleRecAlgo::reconstruct(), saturationFlagSetter_, HcalADCSaturationFlag::setSaturationFlag(), setSaturationFlags_, sevToken_, HcalLongRecoParam::signalTS(), subdet_, HcalZDCDetId::SubdetectorId, tok_input_castor, and tok_input_hcal.

                                                                               {
   // get conditions
   const HcalDbService* conditions = &eventSetup.getData(conditionsToken_);
   const HcalChannelQuality* myqual = &eventSetup.getData(qualToken_);
   const HcalSeverityLevelComputer* mySeverity = &eventSetup.getData(sevToken_);
 
   // 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()) {
       edm::Handle<ZDCDigiCollection> digi_castor;
       e.getByToken(tok_input_castor, digi_castor);
       if (!digi_castor.isValid() || digi_castor->empty())
         edm::LogInfo("ZdcHitReconstructor") << "No ZDC info found in either castorDigis or hcalDigis." << std::endl;
       if (digi_castor.isValid())
         e.getByToken(tok_input_castor, digi);
     }
 
     // create empty output
     auto rec = std::make_unique<ZDCRecHitCollection>();
     rec->reserve(digi->size());
     // run the algorithm
     ZDCDigiCollection::const_iterator i;
     for (i = digi->begin(); i != digi->end(); i++) {
       HcalZDCDetId cell = i->id();
       DetId detcell = (DetId)cell;
       // check on cells to be ignored and dropped: (rof,20.Feb.09)
       const HcalChannelStatus* mydigistatus = myqual->getValues(detcell.rawId());
       if (mySeverity->dropChannel(mydigistatus->getValue()))
         continue;
       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);
 
       // get db values for signalTSs and noiseTSs
       const HcalLongRecoParam* myParams = longRecoParams_->getValues(detcell);
       mySignalTS.clear();
       myNoiseTS.clear();
       mySignalTS = myParams->signalTS();
       myNoiseTS = myParams->noiseTS();
 
       rec->push_back(reco_.reconstruct(*i, myNoiseTS, mySignalTS, coder, calibrations));
       (rec->back()).setFlags(0);
       if (setSaturationFlags_)
         saturationFlagSetter_->setSaturationFlag(rec->back(), *i);
 
       // Set auxiliary flag with subset of digi information
       // ZDC aux flag can store non-contiguous set of values
       int auxflag = 0;
       for (unsigned int xx = 0; xx < AuxTSvec_.size() && xx < 4; ++xx) {
         if (AuxTSvec_[xx] < 0 || AuxTSvec_[xx] > 9)
           continue;  // don't allow
         auxflag += (i->sample(AuxTSvec_[xx]).adc())
                    << (7 * xx);  // store the time slices in the first 28 bits of aux, a set of 4 7-bit a dc values
       }
       // bits 28 and 29 are reserved for capid of the first time slice saved in aux
       if (!AuxTSvec_.empty())
         auxflag += ((i->sample(AuxTSvec_[0]).capid()) << 28);
       (rec->back()).setAux(auxflag);
     }
     // return result
     e.put(std::move(rec));
   }  // else if (det_==DetId::Calo...)
 
 }  // void HcalHitReconstructor::produce(...)