RecHitCorrector Class Reference

#include <RecoLocalCalo/Castor/src/RecHitCorrector.cc>

Inheritance diagram for RecHitCorrector:

Public Member Functions
	RecHitCorrector (const edm::ParameterSet &)
	~RecHitCorrector () 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 Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
bool	doInterCalib_
double	factor_
edm::EDGetTokenT< CastorRecHitCollection >	tok_input_

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 46 of file RecHitCorrector.cc.

Constructor & Destructor Documentation

RecHitCorrector::RecHitCorrector ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 71 of file RecHitCorrector.cc.

References edm::ParameterSet::getParameter(), and tok_input_.

    : factor_(iConfig.getParameter<double>("revertFactor")), doInterCalib_(iConfig.getParameter<bool>("doInterCalib")) {
  tok_input_ = consumes<CastorRecHitCollection>(iConfig.getParameter<edm::InputTag>("rechitLabel"));
  //register your products
  produces<CastorRecHitCollection>();
  //now do what ever other initialization is needed
}

RecHitCorrector::~RecHitCorrector ( )

override

Definition at line 79 of file RecHitCorrector.cc.

                                  {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

void RecHitCorrector::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 89 of file RecHitCorrector.cc.

References AlignmentProducer_cff::calibrations, ewkTauDQM_cfi::channels, DEFINE_FWK_MODULE, doInterCalib_, CaloRecHit::energy(), factor_, CastorCalibrations::gain(), edm::EventSetup::get(), CastorCondObjectContainer< Item >::getAllChannels(), edm::Event::getByToken(), CastorDbService::getCastorCalibrations(), CastorChannelStatus::getValue(), CastorCondObjectContainer< Item >::getValues(), mps_fire::i, CastorRecHit::id(), edm::HandleBase::isValid(), eostools::move(), convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< T >::product(), edm::Event::put(), DetId::rawId(), TrackInfoProducer_cfi::rechits, edm::SortedCollection< T, SORT >::size(), CaloRecHit::time(), ntuplemaker::time, and tok_input_.

                                                                           {
  using namespace edm;

  // get original rechits
  edm::Handle<CastorRecHitCollection> rechits;
  iEvent.getByToken(tok_input_, rechits);

  // get conditions
  edm::ESHandle<CastorDbService> conditions;
  iSetup.get<CastorDbRecord>().get(conditions);

  edm::ESHandle<CastorChannelQuality> p;
  iSetup.get<CastorChannelQualityRcd>().get(p);
  CastorChannelQuality* myqual = new CastorChannelQuality(*p.product());

  if (!rechits.isValid())
    edm::LogWarning("CastorRecHitCorrector") << "No valid CastorRecHitCollection found, please check the InputLabel...";

  CastorCalibrations calibrations;

  auto rec = std::make_unique<CastorRecHitCollection>();

  for (unsigned int i = 0; i < rechits->size(); i++) {
    CastorRecHit rechit = (*rechits)[i];
    double time = rechit.time();
    double correctedenergy = factor_ * rechit.energy();

    if (doInterCalib_) {
      // do proper gain calibration reading the latest entries in the condDB
      const CastorCalibrations& calibrations = conditions->getCastorCalibrations(rechit.id());
      int capid = 0;  // take some capid, gains are the same for all capid's
      correctedenergy *= calibrations.gain(capid);
    }

    // now check the channelquality of this rechit
    bool ok = true;
    DetId detcell = (DetId)rechit.id();
    std::vector<DetId> channels = myqual->getAllChannels();
    for (auto channel : channels) {
      if (channel.rawId() == detcell.rawId()) {
        const CastorChannelStatus* mydigistatus = myqual->getValues(channel);
        if (mydigistatus->getValue() == 2989) {
          ok = false;  // 2989 = BAD
          break;
        }
      }
    }

    if (ok) {
      rec->emplace_back(rechit.id(), correctedenergy, time);
    }
  }

  iEvent.put(std::move(rec));

  delete myqual;
}

Member Data Documentation

bool RecHitCorrector::doInterCalib_

private

Definition at line 57 of file RecHitCorrector.cc.

Referenced by produce().

double RecHitCorrector::factor_

private

Definition at line 56 of file RecHitCorrector.cc.

Referenced by produce().

edm::EDGetTokenT<CastorRecHitCollection> RecHitCorrector::tok_input_

private

Definition at line 55 of file RecHitCorrector.cc.

Referenced by produce(), and RecHitCorrector().