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

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () 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::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<>
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 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));
   }
 }