#include <CastorSimpleReconstructor.h>

Inheritance diagram for CastorSimpleReconstructor:

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

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

	~CastorSimpleReconstructor () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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::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

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

override

Definition at line 50 of file CastorSimpleReconstructor.cc.

50 {

51 }

Member Function Documentation

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

override

Definition at line 53 of file CastorSimpleReconstructor.cc.

References HcalCaloFlagLabels::ADCSaturationBit, edm::SortedCollection< T, SORT >::begin(), AlignmentProducer_cff::calibrations, DetId::Calo, CastorSimpleRecAlgo::checkADCSaturation(), det_, doSaturationCorr_, edm::SortedCollection< T, SORT >::end(), CastorRecoParam::firstSample(), edm::EventSetup::get(), edm::Event::getByToken(), CastorDbService::getCastorCalibrations(), CastorDbService::getCastorCoder(), CastorDbService::getCastorShape(), CastorSaturationCorr::getValue(), CastorCondObjectContainer< Item >::getValues(), mps_fire::i, edm::ESHandleBase::isValid(), maxADCvalue_, eostools::move(), 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
         auto rec = std::make_unique<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(std::move(rec));
     }
 }