#include <EcalUncalibRecHitWorkerAnalFit.h>

Inheritance diagram for EcalUncalibRecHitWorkerAnalFit:

Public Member Functions
	EcalUncalibRecHitWorkerAnalFit (const edm::ParameterSet &ps, edm::ConsumesCollector &c)

	EcalUncalibRecHitWorkerAnalFit ()

edm::ParameterSetDescription	getAlgoDescription ()

bool	run (const edm::Event &evt, const EcalDigiCollection::const_iterator &digi, EcalUncalibratedRecHitCollection &result)

void	set (const edm::EventSetup &es)

	~EcalUncalibRecHitWorkerAnalFit ()

Public Member Functions inherited from EcalUncalibRecHitWorkerBaseClass
	EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &, edm::ConsumesCollector &c)

	EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &)

	EcalUncalibRecHitWorkerBaseClass ()

virtual void	set (const edm::Event &evt)

virtual	~EcalUncalibRecHitWorkerBaseClass ()

Private Attributes
EcalUncalibRecHitRecAnalFitAlgo < EBDataFrame >	algoEB_

EcalUncalibRecHitRecAnalFitAlgo < EEDataFrame >	algoEE_

edm::ESHandle< EcalPedestals >	pedHandle

edm::ESHandle< EcalGainRatios >	pRatio

Detailed Description

produce ECAL uncalibrated rechits from dataframes with the analytical fit method

Author: Shahram Rahatlou, University of Rome & INFN, Sept 2005

Definition at line 27 of file EcalUncalibRecHitWorkerAnalFit.h.

Constructor & Destructor Documentation

EcalUncalibRecHitWorkerAnalFit::EcalUncalibRecHitWorkerAnalFit	(	const edm::ParameterSet &	ps,
		edm::ConsumesCollector &	c
	)

Definition at line 41 of file EcalUncalibRecHitWorkerAnalFit.cc.

                                                                                                                  :
   EcalUncalibRecHitWorkerBaseClass( ps ,c)
 {
 }

EcalUncalibRecHitWorkerAnalFit::EcalUncalibRecHitWorkerAnalFit ( )

inline

Definition at line 31 of file EcalUncalibRecHitWorkerAnalFit.h.

31 {};

EcalUncalibRecHitWorkerAnalFit::~EcalUncalibRecHitWorkerAnalFit ( )

inline

Definition at line 32 of file EcalUncalibRecHitWorkerAnalFit.h.

32 {};

Member Function Documentation

edm::ParameterSetDescription EcalUncalibRecHitWorkerAnalFit::getAlgoDescription ( )

virtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 141 of file EcalUncalibRecHitWorkerAnalFit.cc.

                                                    {
 
   edm::ParameterSetDescription psd;
 //psd.addNode(edm::ParameterSet<edm::InputTag>("EBdigiCollection", edm::InputTag("ecalDigis","ebDigis")) and
 //         edm::ParameterSet<std::string>("EEhitCollection", "EcalUncalibRecHitsEE") and
 //         edm::ParameterSet<edm::InputTag>("EEdigiCollection", edm::InputTag("ecalDigis","eeDigis")) and
 //         edm::ParameterSet<std::string>("algo", "EcalUncalibRecHitWorkerAnalFit") and
 //         edm::ParameterSet<std::string>("EBhitCollection", "EcalUncalibRecHitsEB"));
 //
   return psd;
 }

bool EcalUncalibRecHitWorkerAnalFit::run	(	const edm::Event &	evt,
		const EcalDigiCollection::const_iterator &	digi,
		EcalUncalibratedRecHitCollection &	result
	)

virtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 62 of file EcalUncalibRecHitWorkerAnalFit.cc.

References algoEB_, algoEE_, EcalUncalibratedRecHit::amplitude(), EcalBarrel, LogDebug, EcalUncalibRecHitRecAnalFitAlgo< C >::makeRecHit(), pedHandle, pRatio, and edm::SortedCollection< T, SORT >::push_back().

 {
         using namespace edm;
 
         const EcalGainRatioMap & gainMap = pRatio.product()->getMap(); // map of gain ratios
         const EcalPedestalsMap & pedMap = pedHandle.product()->getMap(); // map of pedestals
 
         EcalPedestalsMapIterator pedIter; // pedestal iterator
         EcalPedestals::Item aped; // pedestal object for a single xtal
 
         EcalGainRatioMap::const_iterator gainIter; // gain iterator
         EcalMGPAGainRatio aGain; // gain object for a single xtal
 
         DetId detid( itdg->id() );
 
         // find pedestals for this channel
         //LogDebug("EcalUncalibRecHitDebug") << "looking up pedestal for crystal: " << itdg->id(); // FIXME
         pedIter = pedMap.find( detid );
         if( pedIter != pedMap.end() ) {
                 aped = (*pedIter);
         } else {
                 edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "error!! could not find pedestals for channel: ";
                 if ( detid.subdetId() == EcalBarrel ) {
                         edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EBDetId( detid );
                 } else {
                         edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EEDetId( detid );
                 }
                 edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "\n  no uncalib rechit will be made for this digi!";
                 return false;
         }
         double pedVec[3];
         pedVec[0] = aped.mean_x12;
         pedVec[1] = aped.mean_x6;
         pedVec[2] = aped.mean_x1;
 
 
         // find gain ratios
         //LogDebug("EcalUncalibRecHitDebug") << "looking up gainRatios for crystal: " << itdg->id(); // FIXME
         gainIter = gainMap.find( detid );
         if( gainIter != gainMap.end() ) {
                 aGain = (*gainIter);
         } else {
                 edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "error!! could not find gain ratios for channel: ";
                 if ( detid.subdetId() == EcalBarrel ) {
                         edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EBDetId( detid );
                 } else {
                         edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EEDetId( detid );
                 }
                 edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "\n  no uncalib rechit will be made for this digi!";
                 return false;
         }
         double gainRatios[3];
         gainRatios[0] = 1.;
         gainRatios[1] = aGain.gain12Over6();
         gainRatios[2] = aGain.gain6Over1()*aGain.gain12Over6();
 
         if ( detid.subdetId() == EcalBarrel ) {
                 EcalUncalibratedRecHit aHit = algoEB_.makeRecHit(*itdg, pedVec, gainRatios, 0 ,0);
                 result.push_back( aHit );
                 if(aHit.amplitude()>0.) {
                         LogDebug("EcalUncalibRecHitInfo") << "EcalUncalibRecHitWorkerAnalFit: processed EBDataFrame with id: "
                                 << EBDetId( detid )
                                 << "\n" << "uncalib rechit amplitude: " << aHit.amplitude();
                 }
         } else {
                 EcalUncalibratedRecHit aHit = algoEE_.makeRecHit(*itdg, pedVec, gainRatios, 0, 0);
                 result.push_back( aHit );
                 if(aHit.amplitude()>0.) {
                         LogDebug("EcalUncalibRecHitInfo") << "EcalUncalibRecHitWorkerAnalFit: processed EEDataFrame with id: "
                                 << EEDetId( detid ) << "\n"
                                 << "uncalib rechit amplitude: " << aHit.amplitude();
                 }
         }
         return true;
 }

void EcalUncalibRecHitWorkerAnalFit::set ( const edm::EventSetup & es )

virtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 48 of file EcalUncalibRecHitWorkerAnalFit.cc.

References edm::EventSetup::get(), LogDebug, pedHandle, and pRatio.

 {
         // Gain Ratios
         LogDebug("EcalUncalibRecHitDebug") << "fetching gainRatios....";
         es.get<EcalGainRatiosRcd>().get(pRatio);
         LogDebug("EcalUncalibRecHitDebug") << "done." ;
 
         // fetch the pedestals from the cond DB via EventSetup
         LogDebug("EcalUncalibRecHitDebug") << "fetching pedestals....";
         es.get<EcalPedestalsRcd>().get( pedHandle );
         LogDebug("EcalUncalibRecHitDebug") << "done." ;
 }