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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
std::vector< int >	AuxTSvec_
bool	correctTiming_
DetId::Detector	det_
bool	dropZSmarkedPassed_
HBHETimeProfileStatusBitSetter *	hbheHSCPFlagSetter_
HBHETimingShapedFlagSetter *	hbheTimingShapedFlagSetter_
HcalHFStatusBitFromDigis *	hfdigibit_
HcalHFStatusBitFromRecHits *	hfrechitbit_
HFTimingTrustFlag *	HFTimingTrustFlagSetter_
double	lowGainFrac_
int	lowGainOffset_
HcalLongRecoParams *	myobject
ZdcSimpleRecAlgo	reco_
HcalADCSaturationFlag *	saturationFlagSetter_
bool	setHSCPFlags_
bool	setNoiseFlags_
bool	setSaturationFlags_
bool	setTimingTrustFlags_
int	subdet_
HcalOtherSubdetector	subdetOther_
HcalTopology *	theTopology
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 35 of file ZdcHitReconstructor.h.

Constructor & Destructor Documentation

ZdcHitReconstructor()

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

explicit

Definition at line 18 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),
      HFTimingTrustFlagSetter_(nullptr),
      hbheHSCPFlagSetter_(nullptr),
      hbheTimingShapedFlagSetter_(nullptr),
      hfrechitbit_(nullptr),
      hfdigibit_(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")),
      myobject(nullptr),
      theTopology(nullptr)
 
{
  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;
  }
}

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

~ZdcHitReconstructor()

ZdcHitReconstructor::~ZdcHitReconstructor ( )

override

Definition at line 66 of file ZdcHitReconstructor.cc.

{ delete saturationFlagSetter_; }

References saturationFlagSetter_.

Member Function Documentation

beginRun()

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

final

Definition at line 68 of file ZdcHitReconstructor.cc.

                                                                           {
  edm::ESHandle<HcalLongRecoParams> p;
  es.get<HcalLongRecoParamsRcd>().get(p);
  myobject = new HcalLongRecoParams(*p.product());
 
  edm::ESHandle<HcalTopology> htopo;
  es.get<HcalRecNumberingRecord>().get(htopo);
  theTopology = new HcalTopology(*htopo);
  myobject->setTopo(theTopology);
}

References edm::EventSetup::get(), get, myobject, AlCaHLTBitMon_ParallelJobs::p, HcalCondObjectContainerBase::setTopo(), and theTopology.

endRun()

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

final

Definition at line 79 of file ZdcHitReconstructor.cc.

                                                                         {
  delete myobject;
  myobject = nullptr;
  delete theTopology;
  theTopology = nullptr;
}

References myobject, and theTopology.

produce()

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

final

Definition at line 85 of file ZdcHitReconstructor.cc.

                                                                              {
  // get conditions
  edm::ESHandle<HcalDbService> conditions;
  eventSetup.get<HcalDbRecord>().get(conditions);
 
  edm::ESHandle<HcalChannelQuality> p;
  eventSetup.get<HcalChannelQualityRcd>().get("withTopo", p);
  const HcalChannelQuality* myqual = p.product();
 
  edm::ESHandle<HcalSeverityLevelComputer> mycomputer;
  eventSetup.get<HcalSeverityLevelComputerRcd>().get(mycomputer);
  const HcalSeverityLevelComputer* mySeverity = mycomputer.product();
 
  // 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 = myobject->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, det_, HcalSeverityLevelComputer::dropChannel(), dropZSmarkedPassed_, MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), edm::EventSetup::get(), get, HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), mps_fire::i, edm::HandleBase::isValid(), eostools::move(), myobject, HcalLongRecoParam::noiseTS(), AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< T >::product(), DetId::rawId(), reco_, ZdcSimpleRecAlgo::reconstruct(), saturationFlagSetter_, HcalADCSaturationFlag::setSaturationFlag(), setSaturationFlags_, 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 65 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

correctTiming_

bool ZdcHitReconstructor::correctTiming_

private

Definition at line 58 of file ZdcHitReconstructor.h.

det_

DetId::Detector ZdcHitReconstructor::det_

private

Definition at line 52 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

dropZSmarkedPassed_

bool ZdcHitReconstructor::dropZSmarkedPassed_

private

Definition at line 64 of file ZdcHitReconstructor.h.

Referenced by produce().

hbheHSCPFlagSetter_

HBHETimeProfileStatusBitSetter* ZdcHitReconstructor::hbheHSCPFlagSetter_

private

Definition at line 47 of file ZdcHitReconstructor.h.

hbheTimingShapedFlagSetter_

HBHETimingShapedFlagSetter* ZdcHitReconstructor::hbheTimingShapedFlagSetter_

private

Definition at line 48 of file ZdcHitReconstructor.h.

hfdigibit_

HcalHFStatusBitFromDigis* ZdcHitReconstructor::hfdigibit_

private

Definition at line 50 of file ZdcHitReconstructor.h.

hfrechitbit_

HcalHFStatusBitFromRecHits* ZdcHitReconstructor::hfrechitbit_

private

Definition at line 49 of file ZdcHitReconstructor.h.

HFTimingTrustFlagSetter_

HFTimingTrustFlag* ZdcHitReconstructor::HFTimingTrustFlagSetter_

private

Definition at line 46 of file ZdcHitReconstructor.h.

lowGainFrac_

double ZdcHitReconstructor::lowGainFrac_

private

Definition at line 68 of file ZdcHitReconstructor.h.

lowGainOffset_

int ZdcHitReconstructor::lowGainOffset_

private

Definition at line 67 of file ZdcHitReconstructor.h.

myobject

HcalLongRecoParams* ZdcHitReconstructor::myobject

private

Definition at line 70 of file ZdcHitReconstructor.h.

Referenced by beginRun(), endRun(), and produce().

reco_

ZdcSimpleRecAlgo ZdcHitReconstructor::reco_

private

Definition at line 44 of file ZdcHitReconstructor.h.

Referenced by produce().

saturationFlagSetter_

HcalADCSaturationFlag* ZdcHitReconstructor::saturationFlagSetter_

private

Definition at line 45 of file ZdcHitReconstructor.h.

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

setHSCPFlags_

bool ZdcHitReconstructor::setHSCPFlags_

private

Definition at line 60 of file ZdcHitReconstructor.h.

setNoiseFlags_

bool ZdcHitReconstructor::setNoiseFlags_

private

Definition at line 59 of file ZdcHitReconstructor.h.

setSaturationFlags_

bool ZdcHitReconstructor::setSaturationFlags_

private

Definition at line 61 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

setTimingTrustFlags_

bool ZdcHitReconstructor::setTimingTrustFlags_

private

Definition at line 62 of file ZdcHitReconstructor.h.

subdet_

int ZdcHitReconstructor::subdet_

private

Definition at line 53 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

subdetOther_

HcalOtherSubdetector ZdcHitReconstructor::subdetOther_

private

Definition at line 54 of file ZdcHitReconstructor.h.

Referenced by ZdcHitReconstructor().

theTopology

HcalTopology* ZdcHitReconstructor::theTopology

private

Definition at line 71 of file ZdcHitReconstructor.h.

Referenced by beginRun(), and endRun().

tok_input_castor

edm::EDGetTokenT<ZDCDigiCollection> ZdcHitReconstructor::tok_input_castor

private

Definition at line 56 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

tok_input_hcal

edm::EDGetTokenT<ZDCDigiCollection> ZdcHitReconstructor::tok_input_hcal

private

Definition at line 55 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().