#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	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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 51 of file CastorSimpleReconstructor.cc.

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