#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

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

Definition at line 46 of file ZdcHitReconstructor.h.

Constructor & Destructor Documentation

◆ ZdcHitReconstructor()

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

explicit

Definition at line 17 of file ZdcHitReconstructor.cc.

     : 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>();
 }

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.

◆ ~ZdcHitReconstructor()

ZdcHitReconstructor::~ZdcHitReconstructor ( )

override

Definition at line 63 of file ZdcHitReconstructor.cc.

63 { delete saturationFlagSetter_; }

References saturationFlagSetter_.

Member Function Documentation

◆ beginRun()

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

final

Definition at line 65 of file ZdcHitReconstructor.cc.

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

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

◆ endRun()

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

final

Definition at line 72 of file ZdcHitReconstructor.cc.

72 {}

◆ produce()

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

final

Definition at line 74 of file ZdcHitReconstructor.cc.

                                                                               {
   // 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(...)

References AuxTSvec_, edm::SortedCollection< T, SORT >::begin(), AlignmentProducer_cff::calibrations, DetId::Calo, submitPVValidationJobs::conditions, conditionsToken_, det_, HcalSeverityLevelComputer::dropChannel(), dropZSmarkedPassed_, MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), edm::EventSetup::getData(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), mps_fire::i, edm::HandleBase::isValid(), longRecoParams_, eostools::move(), HcalLongRecoParam::noiseTS(), qualToken_, DetId::rawId(), reco_, ZdcSimpleRecAlgo::reconstruct(), saturationFlagSetter_, HcalADCSaturationFlag::setSaturationFlag(), setSaturationFlags_, sevToken_, HcalLongRecoParam::signalTS(), edm::SortedCollection< T, SORT >::size(), subdet_, HcalZDCDetId::SubdetectorId, tok_input_castor, tok_input_hcal, and geometryCSVtoXML::xx.

Member Data Documentation

◆ AuxTSvec_

std::vector<int> ZdcHitReconstructor::AuxTSvec_

private

Definition at line 71 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ conditionsToken_

edm::ESGetToken<HcalDbService, HcalDbRecord> ZdcHitReconstructor::conditionsToken_

private

Definition at line 82 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ correctTiming_

bool ZdcHitReconstructor::correctTiming_

private

Definition at line 64 of file ZdcHitReconstructor.h.

◆ det_

DetId::Detector ZdcHitReconstructor::det_

private

Definition at line 58 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ dropZSmarkedPassed_

bool ZdcHitReconstructor::dropZSmarkedPassed_

private

Definition at line 70 of file ZdcHitReconstructor.h.

Referenced by produce().

◆ htopoToken_

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> ZdcHitReconstructor::htopoToken_

private

Definition at line 80 of file ZdcHitReconstructor.h.

Referenced by beginRun(), and ZdcHitReconstructor().

◆ longRecoParams_

std::unique_ptr<HcalLongRecoParams> ZdcHitReconstructor::longRecoParams_

private

Definition at line 77 of file ZdcHitReconstructor.h.

Referenced by beginRun(), and produce().

◆ lowGainFrac_

double ZdcHitReconstructor::lowGainFrac_

private

Definition at line 75 of file ZdcHitReconstructor.h.

◆ lowGainOffset_

int ZdcHitReconstructor::lowGainOffset_

private

Definition at line 74 of file ZdcHitReconstructor.h.

◆ paramsToken_

edm::ESGetToken<HcalLongRecoParams, HcalLongRecoParamsRcd> ZdcHitReconstructor::paramsToken_

private

Definition at line 81 of file ZdcHitReconstructor.h.

Referenced by beginRun(), and ZdcHitReconstructor().

◆ qualToken_

edm::ESGetToken<HcalChannelQuality, HcalChannelQualityRcd> ZdcHitReconstructor::qualToken_

private

Definition at line 83 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ reco_

ZdcSimpleRecAlgo ZdcHitReconstructor::reco_

private

Definition at line 55 of file ZdcHitReconstructor.h.

Referenced by produce().

◆ saturationFlagSetter_

HcalADCSaturationFlag* ZdcHitReconstructor::saturationFlagSetter_

private

Definition at line 56 of file ZdcHitReconstructor.h.

Referenced by produce(), ZdcHitReconstructor(), and ~ZdcHitReconstructor().

◆ setHSCPFlags_

bool ZdcHitReconstructor::setHSCPFlags_

private

Definition at line 66 of file ZdcHitReconstructor.h.

◆ setNoiseFlags_

bool ZdcHitReconstructor::setNoiseFlags_

private

Definition at line 65 of file ZdcHitReconstructor.h.

◆ setSaturationFlags_

bool ZdcHitReconstructor::setSaturationFlags_

private

Definition at line 67 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ setTimingTrustFlags_

bool ZdcHitReconstructor::setTimingTrustFlags_

private

Definition at line 68 of file ZdcHitReconstructor.h.

◆ sevToken_

edm::ESGetToken<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd> ZdcHitReconstructor::sevToken_

private

Definition at line 84 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ subdet_

int ZdcHitReconstructor::subdet_

private

Definition at line 59 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ subdetOther_

HcalOtherSubdetector ZdcHitReconstructor::subdetOther_

private

Definition at line 60 of file ZdcHitReconstructor.h.

Referenced by ZdcHitReconstructor().

◆ tok_input_castor

edm::EDGetTokenT<ZDCDigiCollection> ZdcHitReconstructor::tok_input_castor

private

Definition at line 62 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

◆ tok_input_hcal

edm::EDGetTokenT<ZDCDigiCollection> ZdcHitReconstructor::tok_input_hcal

private

Definition at line 61 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ZdcHitReconstructor()

◆ ~ZdcHitReconstructor()

Member Function Documentation

◆ beginRun()

◆ endRun()

◆ produce()

Member Data Documentation

◆ AuxTSvec_

◆ conditionsToken_

◆ correctTiming_

◆ det_

◆ dropZSmarkedPassed_

◆ htopoToken_

◆ longRecoParams_

◆ lowGainFrac_

◆ lowGainOffset_

◆ paramsToken_

◆ qualToken_

◆ reco_

◆ saturationFlagSetter_

◆ setHSCPFlags_

◆ setNoiseFlags_

◆ setSaturationFlags_

◆ setTimingTrustFlags_

◆ sevToken_

◆ subdet_

◆ subdetOther_

◆ tok_input_castor

◆ tok_input_hcal