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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 ( const edm::ParameterSet &  ps )

explicit

Definition at line 18 of file ZdcHitReconstructor.cc.

References AuxTSvec_, DetId::Calo, gather_cfg::cout, det_, edm::ParameterSet::getParameter(), HcalCalibration, HcalOther, saturationFlagSetter_, setSaturationFlags_, 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),
  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;
  }       
  
}

ZdcHitReconstructor::~ZdcHitReconstructor ( )

override

Definition at line 68 of file ZdcHitReconstructor.cc.

References saturationFlagSetter_.

{
    delete saturationFlagSetter_;
}

Member Function Documentation

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

final

Definition at line 73 of file ZdcHitReconstructor.cc.

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

                                                                          {

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

}

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

final

Definition at line 86 of file ZdcHitReconstructor.cc.

References myobject, and theTopology.

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

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

final

Definition at line 90 of file ZdcHitReconstructor.cc.

References AuxTSvec_, edm::SortedCollection< T, SORT >::begin(), AlignmentProducer_cff::calibrations, DetId::Calo, det_, HcalSeverityLevelComputer::dropChannel(), dropZSmarkedPassed_, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), edm::EventSetup::get(), edm::Event::getByToken(), 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(), edm::Event::put(), 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.

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

Member Data Documentation

std::vector<int> ZdcHitReconstructor::AuxTSvec_

private

Definition at line 64 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

bool ZdcHitReconstructor::correctTiming_

private

Definition at line 57 of file ZdcHitReconstructor.h.

DetId::Detector ZdcHitReconstructor::det_

private

Definition at line 51 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

bool ZdcHitReconstructor::dropZSmarkedPassed_

private

Definition at line 63 of file ZdcHitReconstructor.h.

Referenced by produce().

HBHETimeProfileStatusBitSetter* ZdcHitReconstructor::hbheHSCPFlagSetter_

private

Definition at line 46 of file ZdcHitReconstructor.h.

HBHETimingShapedFlagSetter* ZdcHitReconstructor::hbheTimingShapedFlagSetter_

private

Definition at line 47 of file ZdcHitReconstructor.h.

HcalHFStatusBitFromDigis* ZdcHitReconstructor::hfdigibit_

private

Definition at line 49 of file ZdcHitReconstructor.h.

HcalHFStatusBitFromRecHits* ZdcHitReconstructor::hfrechitbit_

private

Definition at line 48 of file ZdcHitReconstructor.h.

HFTimingTrustFlag* ZdcHitReconstructor::HFTimingTrustFlagSetter_

private

Definition at line 45 of file ZdcHitReconstructor.h.

double ZdcHitReconstructor::lowGainFrac_

private

Definition at line 67 of file ZdcHitReconstructor.h.

int ZdcHitReconstructor::lowGainOffset_

private

Definition at line 66 of file ZdcHitReconstructor.h.

HcalLongRecoParams* ZdcHitReconstructor::myobject

private

Definition at line 69 of file ZdcHitReconstructor.h.

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

ZdcSimpleRecAlgo ZdcHitReconstructor::reco_

private

Definition at line 43 of file ZdcHitReconstructor.h.

Referenced by produce().

HcalADCSaturationFlag* ZdcHitReconstructor::saturationFlagSetter_

private

Definition at line 44 of file ZdcHitReconstructor.h.

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

bool ZdcHitReconstructor::setHSCPFlags_

private

Definition at line 59 of file ZdcHitReconstructor.h.

bool ZdcHitReconstructor::setNoiseFlags_

private

Definition at line 58 of file ZdcHitReconstructor.h.

bool ZdcHitReconstructor::setSaturationFlags_

private

Definition at line 60 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

bool ZdcHitReconstructor::setTimingTrustFlags_

private

Definition at line 61 of file ZdcHitReconstructor.h.

int ZdcHitReconstructor::subdet_

private

Definition at line 52 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

HcalOtherSubdetector ZdcHitReconstructor::subdetOther_

private

Definition at line 53 of file ZdcHitReconstructor.h.

Referenced by ZdcHitReconstructor().

HcalTopology* ZdcHitReconstructor::theTopology

private

Definition at line 70 of file ZdcHitReconstructor.h.

Referenced by beginRun(), and endRun().

edm::EDGetTokenT<ZDCDigiCollection> ZdcHitReconstructor::tok_input_castor

private

Definition at line 55 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

edm::EDGetTokenT<ZDCDigiCollection> ZdcHitReconstructor::tok_input_hcal

private

Definition at line 54 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().