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::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override
void	endJob () override
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::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 75 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 86 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::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 175 of file RecHitCorrector.cc.

{
}

void RecHitCorrector::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 181 of file RecHitCorrector.cc.

References DEFINE_FWK_MODULE.

                        {
}

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

overrideprivate

Definition at line 101 of file RecHitCorrector.cc.

References AlignmentProducer_cff::calibrations, gather_cfg::cout, 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(), HcalCastorDetId::module(), 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()) std::cout << "No valid CastorRecHitCollection found, please check the InputLabel..." << std::endl;
   
   CastorCalibrations calibrations;
   
   auto rec = std::make_unique<CastorRecHitCollection>();
   
   for (unsigned int i=0;i<rechits->size();i++) {
    CastorRecHit rechit = (*rechits)[i];
    //std::cout << "rechit energy = " << rechit.energy() << std::endl;
    double fC = factor_*rechit.energy();
    double time = rechit.time();
    //std::cout << "rechit energy(fC) = " << fC << " time = " << time << std::endl;
    
    // 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
    
    double correctedenergy = 0;
    if (doInterCalib_) {
        if (rechit.id().module() <= 2) {
            correctedenergy = 0.5*fC*calibrations.gain(capid);
            //std::cout << " correctedenergy = " << correctedenergy << " gain = " << calibrations.gain(capid) << std::endl;
        } else {
            correctedenergy = fC*calibrations.gain(capid);
        }
    } else {
        if (rechit.id().module() <= 2) {
            correctedenergy = 0.5*fC;
        } else {
            correctedenergy = fC;
        }
    }
    
    // now check the channelquality of this rechit
    bool ok = true;
    DetId detcell=(DetId)rechit.id();
    std::vector<DetId> channels = myqual->getAllChannels();
    //std::cout << "number of specified quality flags = " << channels.size() << std::endl;
    for (std::vector<DetId>::iterator channel = channels.begin();channel !=  channels.end();channel++) {    
        if (channel->rawId() == detcell.rawId()) {
            const CastorChannelStatus* mydigistatus=myqual->getValues(*channel);
            //std::cout << "CastorChannelStatus = " << mydigistatus->getValue() << std::endl;
            if (mydigistatus->getValue() == 2989) ok = false; // 2989 = BAD
        }
    }
    
    if (ok) {
        CastorRecHit *correctedhit = new CastorRecHit(rechit.id(),correctedenergy,time);
        rec->push_back(*correctedhit);
    }
   }
   
   iEvent.put(std::move(rec));
 
}

Member Data Documentation

bool RecHitCorrector::doInterCalib_

private

Definition at line 60 of file RecHitCorrector.cc.

Referenced by produce().

double RecHitCorrector::factor_

private

Definition at line 59 of file RecHitCorrector.cc.

Referenced by produce().

edm::EDGetTokenT<CastorRecHitCollection> RecHitCorrector::tok_input_

private

Definition at line 58 of file RecHitCorrector.cc.

Referenced by produce(), and RecHitCorrector().