RecHitCorrector Class Reference

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

Inheritance diagram for RecHitCorrector:

Public Member Functions
	RecHitCorrector (const edm::ParameterSet &)
	~RecHitCorrector ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob () override
virtual void	endJob () override
virtual 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
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

(

)

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.

                        {
}

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 101 of file RecHitCorrector.cc.

References gather_cfg::cout, doInterCalib_, CaloRecHit::energy(), factor_, CastorCalibrations::gain(), edm::EventSetup::get(), edm::Event::getByToken(), CastorChannelStatus::getValue(), i, CastorRecHit::id(), HcalCastorDetId::module(), convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, edm::Event::put(), DetId::rawId(), HI_PhotonSkim_cff::rechits, CaloRecHit::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;
   
   std::auto_ptr<CastorRecHitCollection> rec(new 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(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().