ZdcHitReconstructor Class Reference

#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::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_, jetUpdater_cfi::sort, 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::~ZdcHitReconstructor ( )

override

Definition at line 63 of file ZdcHitReconstructor.cc.

References saturationFlagSetter_.

{ delete 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.

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

endRun()

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

final

Definition at line 72 of file ZdcHitReconstructor.cc.

{}

produce()

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

final

Definition at line 74 of file ZdcHitReconstructor.cc.

References AuxTSvec_, edm::SortedCollection< T, SORT >::begin(), AlignmentProducer_cff::calibrations, DetId::Calo, PixelBaryCentreAnalyzer_cfg::conditions, conditionsToken_, det_, HcalSeverityLevelComputer::dropChannel(), dropZSmarkedPassed_, MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), options_cfi::eventSetup, 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_, l1trig_cff::shape, HcalLongRecoParam::signalTS(), edm::SortedCollection< T, SORT >::size(), subdet_, HcalZDCDetId::SubdetectorId, tok_input_castor, tok_input_hcal, and geometryCSVtoXML::xx.

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

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