#include <ZdcSimpleReconstructor.h>

Inheritance diagram for ZdcSimpleReconstructor:

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)

	ZdcSimpleReconstructor (const edm::ParameterSet &ps)

virtual	~ZdcSimpleReconstructor ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
DetId::Detector	det_

bool	dropZSmarkedPassed_

HcalLongRecoParams *	myobject

ZdcSimpleRecAlgo	reco_

int	subdet_

HcalOtherSubdetector	subdetOther_

edm::EDGetTokenT < ZDCDigiCollection >	tok_input_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

Definition at line 22 of file ZdcSimpleReconstructor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file ZdcSimpleReconstructor.cc.

References DetId::Calo, gather_cfg::cout, det_, edm::ParameterSet::getParameter(), HcalCalibration, HcalOther, 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),
   dropZSmarkedPassed_(conf.getParameter<bool>("dropZSmarkedPassed"))
 {
   tok_input_ = consumes<ZDCDigiCollection>(conf.getParameter<edm::InputTag>("digiLabel"));
 
   std::string subd=conf.getParameter<std::string>("Subdetector");
   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 << "ZdcSimpleReconstructor is not associated with a specific subdetector!" << std::endl;
   }       
   
 }

ZdcSimpleReconstructor::~ZdcSimpleReconstructor ( )

virtual

Definition at line 40 of file ZdcSimpleReconstructor.cc.

40 {

41 }

Member Function Documentation

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

finaloverridevirtual

Reimplemented from edm::EDProducer.

Definition at line 42 of file ZdcSimpleReconstructor.cc.

References edm::EventSetup::get(), myobject, AlCaHLTBitMon_ParallelJobs::p, and edm::ESHandle< class >::product().

                                                                              {
 
    edm::ESHandle<HcalLongRecoParams> p;
    es.get<HcalLongRecoParamsRcd>().get(p);
    myobject = new HcalLongRecoParams(*p.product());
 }

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

finaloverridevirtual

Reimplemented from edm::EDProducer.

Definition at line 49 of file ZdcSimpleReconstructor.cc.

References myobject.

                                                                            {
   delete myobject; myobject = 0;
 }

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

virtual

Implements edm::EDProducer.

Definition at line 52 of file ZdcSimpleReconstructor.cc.

References DetId::Calo, det_, dropZSmarkedPassed_, edm::EventSetup::get(), edm::Event::getByToken(), HcalCondObjectContainer< Item >::getValues(), i, ZdcSimpleRecAlgo::initPulseCorr(), myobject, HcalLongRecoParam::noiseTS(), edm::Event::put(), reco_, ZdcSimpleRecAlgo::reconstruct(), HcalLongRecoParam::signalTS(), subdet_, HcalZDCDetId::SubdetectorId, and tok_input_.

 {
   // get conditions
   edm::ESHandle<HcalDbService> conditions;
   eventSetup.get<HcalDbRecord>().get(conditions);
   // 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
     unsigned int toaddMem = 0;
     
     ZDCDigiCollection::const_iterator i;
     for (i=digi->begin(); i!=digi->end(); i++) {
       HcalZDCDetId cell = i->id();    
       DetId detcell=(DetId)cell;
     // rof 27.03.09: drop ZS marked and passed digis:
     if (dropZSmarkedPassed_)
       if (i->zsMarkAndPass()) continue;
 
 // get db values for signalTSs and noiseTSs
       const HcalLongRecoParam* myParams = myobject->getValues(detcell);
       mySignalTS.clear();
       myNoiseTS.clear();
       mySignalTS = myParams->signalTS();
       myNoiseTS = myParams->noiseTS(); 
 // warning: the PulseCorrection is not used by ZDC. If it gets a non-contingious set of 
 // signal TS, it may not work properly. Assume contiguous here....
         unsigned int toadd = mySignalTS.size();    
         if(toaddMem != toadd) {
       reco_.initPulseCorr(toadd);
           toaddMem = toadd;
     }   
       const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
       const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
       const HcalQIEShape* shape = conditions->getHcalShape (channelCoder); 
       HcalCoderDb coder (*channelCoder, *shape);
       rec->push_back(reco_.reconstruct(*i,myNoiseTS,mySignalTS,coder,calibrations));
     }
     // return result
     e.put(rec);     
   }
 }