#include <CastorSimpleReconstructor.h>

Inheritance diagram for CastorSimpleReconstructor:

Public Member Functions
	CastorSimpleReconstructor (const edm::ParameterSet &ps)

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

virtual	~CastorSimpleReconstructor ()

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	doSaturationCorr_

int	firstSample_

int	maxADCvalue_

CastorSimpleRecAlgo	reco_

int	samplesToAdd_

bool	setSaturationFlag_

int	subdet_

edm::EDGetTokenT < CastorDigiCollection >	tok_input_

bool	tsFromDB_

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 14 of file CastorSimpleReconstructor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 25 of file CastorSimpleReconstructor.cc.

References DetId::Calo, det_, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, subdet_, HcalCastorDetId::SubdetectorId, and tok_input_.

                                                                                :
 reco_(conf.getParameter<int>("firstSample"),conf.getParameter<int>("samplesToAdd"),conf.getParameter<bool>("correctForTimeslew"),
       conf.getParameter<bool>("correctForPhaseContainment"),conf.getParameter<double>("correctionPhaseNS")),
 det_(DetId::Hcal),
 firstSample_(conf.getParameter<int>("firstSample")),
 samplesToAdd_(conf.getParameter<int>("samplesToAdd")),
 maxADCvalue_(conf.getParameter<int>("maxADCvalue")),
 tsFromDB_(conf.getParameter<bool>("tsFromDB")),
 setSaturationFlag_(conf.getParameter<bool>("setSaturationFlag")),
 doSaturationCorr_(conf.getParameter<bool>("doSaturationCorr"))
 {
 
   tok_input_ = consumes<CastorDigiCollection>(conf.getParameter<edm::InputTag>("digiLabel"));
 
     std::string subd=conf.getParameter<std::string>("Subdetector");
     if (!strcasecmp(subd.c_str(),"CASTOR")) {
         det_=DetId::Calo;
         subdet_=HcalCastorDetId::SubdetectorId;
         produces<CastorRecHitCollection>();
     } else {
         edm::LogWarning("CastorSimpleReconstructor") << "CastorSimpleReconstructor is not associated with CASTOR subdetector!" << std::endl;
     }       
     
 }

CastorSimpleReconstructor::~CastorSimpleReconstructor ( )

virtual

Definition at line 50 of file CastorSimpleReconstructor.cc.

50 {

51 }

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 53 of file CastorSimpleReconstructor.cc.

References HcalCaloFlagLabels::ADCSaturationBit, DetId::Calo, CastorSimpleRecAlgo::checkADCSaturation(), det_, doSaturationCorr_, CastorRecoParam::firstSample(), edm::EventSetup::get(), edm::Event::getByToken(), CastorSaturationCorr::getValue(), i, edm::ESHandleBase::isValid(), maxADCvalue_, edm::Event::put(), DetId::rawId(), reco_, CastorSimpleRecAlgo::reconstruct(), CastorSimpleRecAlgo::recoverADCSaturation(), CastorSimpleRecAlgo::resetTimeSamples(), CastorRecoParam::samplesToAdd(), setSaturationFlag_, subdet_, HcalCastorDetId::SubdetectorId, tok_input_, and tsFromDB_.

 {
     // get conditions
     edm::ESHandle<CastorDbService> conditions;
     eventSetup.get<CastorDbRecord>().get(conditions);
     const CastorQIEShape* shape = conditions->getCastorShape (); // this one is generic
     
     CastorCalibrations calibrations;
     
     // try to get the TS windows from the db
     edm::ESHandle<CastorRecoParams> recoparams;
     if (tsFromDB_) {
         eventSetup.get<CastorRecoParamsRcd>().get(recoparams);
         if (!recoparams.isValid()) { 
             tsFromDB_ = false;
             edm::LogWarning("CastorSimpleReconstructor") << "Could not handle the CastorRecoParamsRcd correctly, using parameters from cfg file from this event onwards... These parameters could be wrong for this run... please check" << std::endl;
         }
     }
     
     // try to get the saturation correction constants from the db
     edm::ESHandle<CastorSaturationCorrs> satcorr;
     if (doSaturationCorr_) {
         eventSetup.get<CastorSaturationCorrsRcd>().get(satcorr);
         if (!satcorr.isValid()) {
             doSaturationCorr_ = false;
             edm::LogWarning("CastorSimpleReconstructor") << "Could not handle the CastorSaturationCorrsRcd correctly. We'll not try the saturation correction from this event onwards..." << std::endl;
         }
     }
     
     
     if (det_==DetId::Calo && subdet_==HcalCastorDetId::SubdetectorId) {
         edm::Handle<CastorDigiCollection> digi;
         e.getByToken(tok_input_,digi);
         
         // create empty output
         std::auto_ptr<CastorRecHitCollection> rec(new CastorRecHitCollection);
         // run the algorithm
         CastorDigiCollection::const_iterator i;
         for (i=digi->begin(); i!=digi->end(); i++) {
             HcalCastorDetId cell = i->id();
             DetId detcell=(DetId)cell;    
             const CastorCalibrations& calibrations=conditions->getCastorCalibrations(cell);
             
             if (tsFromDB_) {
                 const CastorRecoParam* param_ts = recoparams->getValues(detcell.rawId());
                 reco_.resetTimeSamples(param_ts->firstSample(),param_ts->samplesToAdd());
             }          
             const CastorQIECoder* channelCoder = conditions->getCastorCoder (cell);
             CastorCoderDb coder (*channelCoder, *shape);
             
             // reconstruct the rechit
             rec->push_back(reco_.reconstruct(*i,coder,calibrations));
             
             // set the saturation flag if needed
             if (setSaturationFlag_) {
                 reco_.checkADCSaturation(rec->back(),*i,maxADCvalue_);
                 
                 //++++ Saturation Correction +++++
                 if (doSaturationCorr_ && rec->back().flagField(HcalCaloFlagLabels::ADCSaturationBit)) {
                     // get saturation correction value
                     const CastorSaturationCorr* saturationCorr = satcorr->getValues(detcell.rawId());
                     double satCorrConst = 1.;
                     satCorrConst = saturationCorr->getValue();
                     reco_.recoverADCSaturation(rec->back(),coder,calibrations,*i,maxADCvalue_,satCorrConst);
                 }
             }
         }
         // return result
         e.put(rec);     
     }
 }