#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
DetId::Detector	det_

bool	doSaturationCorr_

int	firstSample_

int	maxADCvalue_

CastorSimpleRecAlgo	reco_

int	samplesToAdd_

bool	setSaturationFlag_

int	subdet_

edm::ESGetToken< CastorDbService, CastorDbRecord >	tok_conditions_

edm::EDGetTokenT< CastorDigiCollection >	tok_input_

edm::ESGetToken< CastorRecoParams, CastorRecoParamsRcd >	tok_recoParams_

edm::ESGetToken< CastorSaturationCorrs, CastorSaturationCorrsRcd >	tok_satCorr_

bool	tsFromDB_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 21 of file CastorSimpleReconstructor.h.

Constructor & Destructor Documentation

◆ CastorSimpleReconstructor()

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

explicit

Definition at line 19 of file CastorSimpleReconstructor.cc.

     : 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"));
   tok_conditions_ = esConsumes<CastorDbService, CastorDbRecord>();
  
   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;
   }
   if (tsFromDB_) {
     tok_recoParams_ = esConsumes<CastorRecoParams, CastorRecoParamsRcd>();
   }
   if (doSaturationCorr_) {
     tok_satCorr_ = esConsumes<CastorSaturationCorrs, CastorSaturationCorrsRcd>();
   }
 }

References DetId::Calo, det_, doSaturationCorr_, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, subdet_, HcalCastorDetId::SubdetectorId, tok_conditions_, tok_input_, tok_recoParams_, tok_satCorr_, and tsFromDB_.

◆ ~CastorSimpleReconstructor()

CastorSimpleReconstructor::~CastorSimpleReconstructor ( )

override

Definition at line 52 of file CastorSimpleReconstructor.cc.

52 {}

Member Function Documentation

◆ produce()

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

override

Definition at line 54 of file CastorSimpleReconstructor.cc.

                                                                                     {
   // get conditions
   edm::ESHandle<CastorDbService> conditions = eventSetup.getHandle(tok_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_) {
     recoparams = eventSetup.getHandle(tok_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_) {
     satcorr = eventSetup.getHandle(tok_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));
   }
 }