#include <ZdcHitReconstructor.h>

Inheritance diagram for ZdcHitReconstructor:

Public Member Functions
virtual void	beginRun (edm::Run const &r, edm::EventSetup const &es) overridefinal

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

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

	ZdcHitReconstructor (const edm::ParameterSet &ps)

virtual	~ZdcHitReconstructor ()

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
std::vector< int >	AuxTSvec_

bool	correctTiming_

DetId::Detector	det_

bool	dropZSmarkedPassed_

HBHEStatusBitSetter *	hbheFlagSetter_

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_

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

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

Definition at line 36 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_, AlCaHLTBitMon_QueryRunRegistry::string, subdet_, HcalZDCDetId::SubdetectorId, subdetOther_, and tok_input_.

                                                                    :
   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")),
   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(0),
   theTopology(0)
   
 { 
   tok_input_ = consumes<ZDCDigiCollection>(conf.getParameter<edm::InputTag>("digiLabel"));
 
   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 ( )

virtual

Definition at line 61 of file ZdcHitReconstructor.cc.

61 {;

62 }

Member Function Documentation

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

finaloverridevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 63 of file ZdcHitReconstructor.cc.

References edm::EventSetup::get(), myobject, AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< class >::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
	)

finaloverridevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 76 of file ZdcHitReconstructor.cc.

References myobject, and theTopology.

                                                                         {
   delete myobject; myobject=0;
   delete theTopology; theTopology=0;
 }

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

finaloverridevirtual

Implements edm::stream::EDProducerBase.

Definition at line 80 of file ZdcHitReconstructor.cc.

References AuxTSvec_, DetId::Calo, det_, HcalSeverityLevelComputer::dropChannel(), dropZSmarkedPassed_, edm::EventSetup::get(), edm::Event::getByToken(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), i, myobject, HcalLongRecoParam::noiseTS(), AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< class >::product(), edm::Event::put(), DetId::rawId(), reco_, ZdcSimpleRecAlgo::reconstruct(), saturationFlagSetter_, HcalADCSaturationFlag::setSaturationFlag(), setSaturationFlags_, HcalLongRecoParam::signalTS(), subdet_, HcalZDCDetId::SubdetectorId, and tok_input_.

 {
   // 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_,digi);
      
      // create empty output
      std::auto_ptr<ZDCRecHitCollection> rec(new 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_.size()>0)
       auxflag+=((i->sample(AuxTSvec_[0]).capid())<<28);
     (rec->back()).setAux(auxflag);
      }
      // return result
      e.put(rec);     
    } // else if (det_==DetId::Calo...)
 
 } // void HcalHitReconstructor::produce(...)

Member Data Documentation

std::vector<int> ZdcHitReconstructor::AuxTSvec_

private

Definition at line 65 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

bool ZdcHitReconstructor::correctTiming_

private

Definition at line 58 of file ZdcHitReconstructor.h.

DetId::Detector ZdcHitReconstructor::det_

private

Definition at line 53 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

bool ZdcHitReconstructor::dropZSmarkedPassed_

private

Definition at line 64 of file ZdcHitReconstructor.h.

Referenced by produce().

HBHEStatusBitSetter* ZdcHitReconstructor::hbheFlagSetter_

private

Definition at line 47 of file ZdcHitReconstructor.h.

HBHETimeProfileStatusBitSetter* ZdcHitReconstructor::hbheHSCPFlagSetter_

private

Definition at line 48 of file ZdcHitReconstructor.h.

HBHETimingShapedFlagSetter* ZdcHitReconstructor::hbheTimingShapedFlagSetter_

private

Definition at line 49 of file ZdcHitReconstructor.h.

HcalHFStatusBitFromDigis* ZdcHitReconstructor::hfdigibit_

private

Definition at line 51 of file ZdcHitReconstructor.h.

HcalHFStatusBitFromRecHits* ZdcHitReconstructor::hfrechitbit_

private

Definition at line 50 of file ZdcHitReconstructor.h.

HFTimingTrustFlag* ZdcHitReconstructor::HFTimingTrustFlagSetter_

private

Definition at line 46 of file ZdcHitReconstructor.h.

double ZdcHitReconstructor::lowGainFrac_

private

Definition at line 68 of file ZdcHitReconstructor.h.

int ZdcHitReconstructor::lowGainOffset_

private

Definition at line 67 of file ZdcHitReconstructor.h.

HcalLongRecoParams* ZdcHitReconstructor::myobject

private

Definition at line 70 of file ZdcHitReconstructor.h.

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

ZdcSimpleRecAlgo ZdcHitReconstructor::reco_

private

Definition at line 44 of file ZdcHitReconstructor.h.

Referenced by produce().

HcalADCSaturationFlag* ZdcHitReconstructor::saturationFlagSetter_

private

Definition at line 45 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

bool ZdcHitReconstructor::setHSCPFlags_

private

Definition at line 60 of file ZdcHitReconstructor.h.

bool ZdcHitReconstructor::setNoiseFlags_

private

Definition at line 59 of file ZdcHitReconstructor.h.

bool ZdcHitReconstructor::setSaturationFlags_

private

Definition at line 61 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

bool ZdcHitReconstructor::setTimingTrustFlags_

private

Definition at line 62 of file ZdcHitReconstructor.h.

int ZdcHitReconstructor::subdet_

private

Definition at line 54 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

HcalOtherSubdetector ZdcHitReconstructor::subdetOther_

private

Definition at line 55 of file ZdcHitReconstructor.h.

Referenced by ZdcHitReconstructor().

HcalTopology* ZdcHitReconstructor::theTopology

private

Definition at line 71 of file ZdcHitReconstructor.h.

Referenced by beginRun(), and endRun().

edm::EDGetTokenT<ZDCDigiCollection> ZdcHitReconstructor::tok_input_

private

Definition at line 56 of file ZdcHitReconstructor.h.

Referenced by produce(), and ZdcHitReconstructor().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation