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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 47 of file RecHitCorrector.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 73 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 84 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 99 of file RecHitCorrector.cc.

References AlignmentProducer_cff::calibrations, 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(), protons_cff::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 58 of file RecHitCorrector.cc.

Referenced by produce().

double RecHitCorrector::factor_

private

Definition at line 57 of file RecHitCorrector.cc.

Referenced by produce().

edm::EDGetTokenT<CastorRecHitCollection> RecHitCorrector::tok_input_

private

Definition at line 56 of file RecHitCorrector.cc.

Referenced by produce(), and RecHitCorrector().