HcalTBDigiProducer Class Reference

#include <SimCalorimetry/HcalTestBeam/interface/HcalTBDigiProducer.h>

Inheritance diagram for HcalTBDigiProducer:

List of all members.

Public Member Functions

HcalTBDigiProducer (const edm::ParameterSet &ps)

virtual void produce (edm::Event &e, const edm::EventSetup &c)

Produces the EDM products,.

virtual ~HcalTBDigiProducer ()

Private Types

typedef CaloTDigitizer
< HBHEDigitizerTraits > HBHEDigitizer

Reconstruction algorithm.

typedef CaloTDigitizer
< HODigitizerTraits > HODigitizer

Private Member Functions

void checkGeometry (const edm::EventSetup &eventSetup)

make sure the digitizer has the correct list of all cells that exist in the geometry

void setPhaseShift (const DetId &detId)

void sortHits (const edm::PCaloHitContainer &hits)

fills the vectors for each subdetector

void updateGeometry ()

Private Attributes

bool doPhaseShift

std::string ecalTBInfoLabel

std::vector< DetId > hbheCells

std::vector< DetId > hoCells

HcalAmplifier * theAmplifier

HcalCoderFactory * theCoderFactory

HcalElectronicsSim * theElectronicsSim

const CaloGeometry * theGeometry

HBHEDigitizer * theHBHEDigitizer

HBHEHitFilter theHBHEHitFilter

std::vector< PCaloHit > theHBHEHits

CaloHitResponse * theHBHEResponse

CaloVShape * theHcalIntegratedShape

CaloVShape * theHcalShape

HcalHitCorrection * theHitCorrection

HODigitizer * theHODigitizer

HOHitFilter theHOHitFilter

std::vector< PCaloHit > theHOHits

CaloHitResponse * theHOResponse

HcalTBSimParameterMap * theParameterMap

double thisPhaseShift

double tunePhaseShift

Detailed Description

Definition at line 24 of file HcalTBDigiProducer.h.

Member Typedef Documentation

typedef CaloTDigitizer<HBHEDigitizerTraits> HcalTBDigiProducer::HBHEDigitizer [private]

Reconstruction algorithm.

Definition at line 48 of file HcalTBDigiProducer.h.

typedef CaloTDigitizer<HODigitizerTraits> HcalTBDigiProducer::HODigitizer [private]

Definition at line 49 of file HcalTBDigiProducer.h.

Constructor & Destructor Documentation

HcalTBDigiProducer::HcalTBDigiProducer ( const edm::ParameterSet & ps ) [explicit]

Definition at line 19 of file HcalTBDigiProducer.cc.

References HcalCoderFactory::DB, detId, doPhaseShift, ecalTBInfoLabel, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), DetId::Hcal, CaloHitResponse::setHitCorrection(), CaloHitResponse::setHitFilter(), HcalTBSimParameterMap::simParameters(), theAmplifier, theCoderFactory, theElectronicsSim, theHBHEDigitizer, theHBHEHitFilter, theHBHEResponse, theHitCorrection, theHODigitizer, theHOHitFilter, theHOResponse, theParameterMap, and tunePhaseShift.

00019                                                                 :
00020   theParameterMap(new HcalTBSimParameterMap(ps)), 
00021   theHcalShape(new HcalShape()),
00022   theHcalIntegratedShape(new CaloShapeIntegrator(theHcalShape)),
00023   theHBHEResponse(new CaloHitResponse(theParameterMap, theHcalIntegratedShape)),
00024   theHOResponse(new CaloHitResponse(theParameterMap, theHcalIntegratedShape)),
00025   theAmplifier(0), theCoderFactory(0), theElectronicsSim(0), 
00026   theHitCorrection(0), theHBHEDigitizer(0), theHODigitizer(0), theHBHEHits(),
00027   theHOHits(), thisPhaseShift(0) {
00028   
00029   produces<HBHEDigiCollection>();
00030   produces<HODigiCollection>();
00031 
00032   DetId detId(DetId::Hcal, 1);
00033   bool syncPhase = (theParameterMap->simParameters(detId)).syncPhase();
00034   doPhaseShift   = !syncPhase;
00035 
00036   theHBHEResponse->setHitFilter(&theHBHEHitFilter);
00037   theHOResponse->setHitFilter(&theHOHitFilter);
00038 
00039   bool doTimeSlew = ps.getParameter<bool>("doTimeSlew");
00040   if(doTimeSlew) {
00041     // no time slewing for HF
00042     theHitCorrection = new HcalHitCorrection(theParameterMap);
00043     theHBHEResponse->setHitCorrection(theHitCorrection);
00044     theHOResponse->setHitCorrection(theHitCorrection);
00045   }
00046 
00047   bool doNoise = ps.getParameter<bool>("doNoise");
00048   theAmplifier = new HcalAmplifier(theParameterMap, doNoise);
00049   theCoderFactory = new HcalCoderFactory(HcalCoderFactory::DB);
00050   theElectronicsSim = new HcalElectronicsSim(theAmplifier, theCoderFactory);
00051 
00052   theHBHEDigitizer = new HBHEDigitizer(theHBHEResponse, theElectronicsSim, doNoise);
00053   theHODigitizer = new HODigitizer(theHOResponse, theElectronicsSim, doNoise);
00054 
00055   tunePhaseShift =  ps.getUntrackedParameter<double>("tunePhaseShiftTB", 1.);
00056   ecalTBInfoLabel = ps.getUntrackedParameter<std::string>("EcalTBInfoLabel","SimEcalTBG4Object");
00057   edm::LogInfo("HcalSim") << "HcalTBDigiProducer initialized with doNoise = "
00058                           << doNoise << ", doTimeSlew = " << doTimeSlew
00059                           << " and doPhaseShift = " << doPhaseShift
00060                           << " tunePhasShift = " << tunePhaseShift;
00061 
00062 }

HcalTBDigiProducer::~HcalTBDigiProducer ( ) [virtual]

Definition at line 64 of file HcalTBDigiProducer.cc.

References theAmplifier, theCoderFactory, theElectronicsSim, theHBHEDigitizer, theHBHEResponse, theHcalIntegratedShape, theHcalShape, theHitCorrection, theHODigitizer, theHOResponse, and theParameterMap.

00064                                         {
00065 
00066   if (theHBHEDigitizer)       delete theHBHEDigitizer;
00067   if (theHODigitizer)         delete theHODigitizer;
00068   if (theParameterMap)        delete theParameterMap;
00069   if (theHcalShape)           delete theHcalShape;
00070   if (theHcalIntegratedShape) delete theHcalIntegratedShape;
00071   if (theHBHEResponse)        delete theHBHEResponse;
00072   if (theHOResponse)          delete theHOResponse;
00073   if (theElectronicsSim)      delete theElectronicsSim;
00074   if (theAmplifier)           delete theAmplifier;
00075   if (theCoderFactory)        delete theCoderFactory;
00076   if (theHitCorrection)       delete theHitCorrection;
00077 }

Member Function Documentation

void HcalTBDigiProducer::checkGeometry ( const edm::EventSetup & eventSetup ) [private]

make sure the digitizer has the correct list of all cells that exist in the geometry

Definition at line 154 of file HcalTBDigiProducer.cc.

References edm::EventSetup::get(), theGeometry, and updateGeometry().

Referenced by produce().

00154                                                                        {
00155 
00156   // TODO find a way to avoid doing this every event
00157   edm::ESHandle<CaloGeometry> geometry;
00158   eventSetup.get<CaloGeometryRecord>().get(geometry);
00159  
00160   const CaloGeometry * pGeometry = &*geometry;
00161 
00162   // see if we need to update
00163   if(pGeometry != theGeometry) {
00164     theGeometry = pGeometry;
00165     updateGeometry();
00166   }
00167 }

void HcalTBDigiProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)			`[virtual]`

Produces the EDM products,.

Implements edm::EDProducer.

Definition at line 80 of file HcalTBDigiProducer.cc.

References checkGeometry(), parsecf::pyparsing::col(), doPhaseShift, ecalTBInfoLabel, HcalHitCorrection::fillChargeSums(), edm::EventSetup::get(), edm::Event::getByLabel(), DetId::Hcal, LogDebug, edm::ESHandle< T >::product(), edm::Handle< T >::product(), edm::Event::put(), CaloTDigitizer< Traits >::run(), HcalCoderFactory::setDbService(), HcalAmplifier::setDbService(), setPhaseShift(), theAmplifier, theCoderFactory, theEcalTBInfo, theHBHEDigitizer, theHBHEHits, theHitCorrection, theHODigitizer, theHOHits, and thisPhaseShift.

00080                                                                              {
00081 
00082   // get the appropriate gains, noises, & widths for this event
00083   edm::ESHandle<HcalDbService> conditions;
00084   eventSetup.get<HcalDbRecord>().get(conditions);
00085   theAmplifier->setDbService(conditions.product());
00086   theCoderFactory->setDbService(conditions.product());
00087 
00088   // get the correct geometry
00089   checkGeometry(eventSetup);
00090 
00091   theHBHEHits.clear();
00092   theHOHits.clear();
00093 
00094   // Step A: Get Inputs
00095   edm::Handle<CrossingFrame<PCaloHit> > cf;
00096   // e.getByType(cf);
00097 
00098   LogDebug("HcalSim") << "HcalTBDigiProducer::produce trying to ger SimHit";
00099   // test access to SimHits
00100   const std::string subdet("HcalHits");
00101   // std::auto_ptr<MixCollection<PCaloHit> > col(new MixCollection<PCaloHit>(cf.product(), subdet));
00102   e.getByLabel("mix", subdet, cf);
00103   std::auto_ptr<MixCollection<PCaloHit> > col(new MixCollection<PCaloHit>(cf.product() ));
00104 
00105   LogDebug("HcalSim") << "HcalTBDigiProducer::produce Collection of SimHit found";
00106   if(theHitCorrection != 0) {
00107     theHitCorrection->fillChargeSums(*col);
00108   }
00109   LogDebug("HcalSim") << "HcalTBDigiProducer::produce Hits corrected";
00110 
00111   // Step B: Create empty output
00112   std::auto_ptr<HBHEDigiCollection> hbheResult(new HBHEDigiCollection());
00113   std::auto_ptr<HODigiCollection> hoResult(new HODigiCollection());
00114   LogDebug("HcalSim") << "HcalTBDigiProducer::produce Empty collection created";
00115   if (doPhaseShift) {
00116     
00117     edm::Handle<PEcalTBInfo> theEcalTBInfo;
00118     e.getByLabel(ecalTBInfoLabel,theEcalTBInfo);
00119     thisPhaseShift = theEcalTBInfo->phaseShift();
00120 
00121     DetId detIdHB(DetId::Hcal, 1);
00122     setPhaseShift(detIdHB);
00123     DetId detIdHO(DetId::Hcal, 3);
00124     setPhaseShift(detIdHO);
00125   }
00126 
00127   // Step C: Invoke the algorithm, passing in inputs and getting back outputs.
00128   theHBHEDigitizer->run(*col, *hbheResult);
00129   edm::LogInfo("HcalSim") << "HcalTBDigiProducer: HBHE digis : " << hbheResult->size();
00130   theHODigitizer->run(*col, *hoResult);
00131   edm::LogInfo("HcalSim") << "HcalTBDigiProducer: HO digis   : " << hoResult->size();
00132 
00133   // Step D: Put outputs into event
00134   e.put(hbheResult);
00135   e.put(hoResult);
00136 
00137 }

void HcalTBDigiProducer::setPhaseShift ( const DetId & detId ) [private]

Definition at line 195 of file HcalTBDigiProducer.cc.

References python::trackProbabilityAnalysis_cff::parameters, CaloHitResponse::setPhaseShift(), HcalTBSimParameterMap::simParameters(), DetId::subdetId(), CaloSimParameters::syncPhase(), theHBHEResponse, theHOResponse, theParameterMap, thisPhaseShift, and tunePhaseShift.

Referenced by produce().

00195                                                           {
00196 
00197   const CaloSimParameters & parameters = theParameterMap->simParameters(detId);
00198   if ( !parameters.syncPhase() ) {
00199     int    myDet          = detId.subdetId();
00200     double passPhaseShift = thisPhaseShift + tunePhaseShift;
00201     if (myDet <= 2) {
00202       theHBHEResponse->setPhaseShift(passPhaseShift);
00203     } else {
00204       theHOResponse->setPhaseShift(passPhaseShift);
00205     }
00206   }
00207 }

void HcalTBDigiProducer::sortHits ( const edm::PCaloHitContainer & hits ) [private]

fills the vectors for each subdetector

Definition at line 139 of file HcalTBDigiProducer.cc.

References HcalBarrel, HcalEndcap, HcalOuter, theHBHEHits, and theHOHits.

00139                                                                    {
00140 
00141   for (edm::PCaloHitContainer::const_iterator hitItr = hits.begin();
00142        hitItr != hits.end(); ++hitItr) {
00143     HcalSubdetector subdet = HcalDetId(hitItr->id()).subdet();
00144     if(subdet == HcalBarrel || subdet == HcalEndcap) {
00145       theHBHEHits.push_back(*hitItr);
00146     } else if(subdet == HcalOuter) {
00147       theHOHits.push_back(*hitItr);
00148     } else {
00149       edm::LogError("HcalSim") << "Bad HcalHit subdetector " << subdet;
00150     }
00151   }
00152 }

void HcalTBDigiProducer::updateGeometry ( ) [private]

Definition at line 169 of file HcalTBDigiProducer.cc.

References CaloGeometry::getValidDetIds(), hbheCells, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, hoCells, LogDebug, CaloTDigitizer< Traits >::setDetIds(), CaloHitResponse::setGeometry(), theGeometry, theHBHEDigitizer, theHBHEResponse, theHODigitizer, and theHOResponse.

Referenced by checkGeometry().

00169                                         {
00170 
00171   theHBHEResponse->setGeometry(theGeometry);
00172   theHOResponse->setGeometry(theGeometry);
00173 
00174   // Get cells for HB and HE
00175   hbheCells.clear();
00176   hbheCells = theGeometry->getValidDetIds(DetId::Hcal, HcalBarrel);
00177   std::vector<DetId> heCells = theGeometry->getValidDetIds(DetId::Hcal, HcalEndcap);
00178   // combine HB & HE
00179   hbheCells.insert(hbheCells.end(), heCells.begin(), heCells.end());
00180 
00181   // Get cells for HO
00182   hoCells.clear();
00183   hoCells = theGeometry->getValidDetIds(DetId::Hcal, HcalOuter);
00184 
00185   edm::LogInfo("HcalSim") << "HcalTBDigiProducer update Geometry with "
00186                           << hbheCells.size() << " cells in HB/HE and "
00187                           << hoCells.size() << " cells in HO";
00188 
00189   theHBHEDigitizer->setDetIds(hbheCells);
00190   LogDebug("HcalSim") << "HcalTBDigiProducer: Set DetID's for HB/HE";
00191   theHODigitizer->setDetIds(hoCells);
00192   LogDebug("HcalSim") << "HcalTBDigiProducer: Set DetID's for HO";
00193 }

Member Data Documentation

bool HcalTBDigiProducer::doPhaseShift [private]

Definition at line 79 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), and produce().

std::string HcalTBDigiProducer::ecalTBInfoLabel [private]

Definition at line 76 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), and produce().

std::vector<DetId> HcalTBDigiProducer::hbheCells [private]

Definition at line 71 of file HcalTBDigiProducer.h.

Referenced by updateGeometry().

std::vector<DetId> HcalTBDigiProducer::hoCells [private]

Definition at line 72 of file HcalTBDigiProducer.h.

Referenced by updateGeometry().

HcalAmplifier* HcalTBDigiProducer::theAmplifier [private]

Definition at line 58 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), produce(), and ~HcalTBDigiProducer().

HcalCoderFactory* HcalTBDigiProducer::theCoderFactory [private]

Definition at line 59 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), produce(), and ~HcalTBDigiProducer().

HcalElectronicsSim* HcalTBDigiProducer::theElectronicsSim [private]

Definition at line 60 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), and ~HcalTBDigiProducer().

const CaloGeometry* HcalTBDigiProducer::theGeometry [private]

Definition at line 70 of file HcalTBDigiProducer.h.

Referenced by checkGeometry(), and updateGeometry().

HBHEDigitizer* HcalTBDigiProducer::theHBHEDigitizer [private]

Definition at line 67 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), produce(), updateGeometry(), and ~HcalTBDigiProducer().

HBHEHitFilter HcalTBDigiProducer::theHBHEHitFilter [private]

Definition at line 62 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer().

std::vector<PCaloHit> HcalTBDigiProducer::theHBHEHits [private]

Definition at line 74 of file HcalTBDigiProducer.h.

Referenced by produce(), and sortHits().

CaloHitResponse* HcalTBDigiProducer::theHBHEResponse [private]

Definition at line 55 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), setPhaseShift(), updateGeometry(), and ~HcalTBDigiProducer().

CaloVShape* HcalTBDigiProducer::theHcalIntegratedShape [private]

Definition at line 53 of file HcalTBDigiProducer.h.

Referenced by ~HcalTBDigiProducer().

CaloVShape* HcalTBDigiProducer::theHcalShape [private]

Definition at line 52 of file HcalTBDigiProducer.h.

Referenced by ~HcalTBDigiProducer().

HcalHitCorrection* HcalTBDigiProducer::theHitCorrection [private]

Definition at line 65 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), produce(), and ~HcalTBDigiProducer().

HODigitizer* HcalTBDigiProducer::theHODigitizer [private]

Definition at line 68 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), produce(), updateGeometry(), and ~HcalTBDigiProducer().

HOHitFilter HcalTBDigiProducer::theHOHitFilter [private]

Definition at line 63 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer().

std::vector<PCaloHit> HcalTBDigiProducer::theHOHits [private]

Definition at line 74 of file HcalTBDigiProducer.h.

Referenced by produce(), and sortHits().

CaloHitResponse* HcalTBDigiProducer::theHOResponse [private]

Definition at line 56 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), setPhaseShift(), updateGeometry(), and ~HcalTBDigiProducer().

HcalTBSimParameterMap* HcalTBDigiProducer::theParameterMap [private]

Definition at line 51 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), setPhaseShift(), and ~HcalTBDigiProducer().

double HcalTBDigiProducer::thisPhaseShift [private]

Definition at line 77 of file HcalTBDigiProducer.h.

Referenced by produce(), and setPhaseShift().

double HcalTBDigiProducer::tunePhaseShift [private]

Definition at line 80 of file HcalTBDigiProducer.h.

Referenced by HcalTBDigiProducer(), and setPhaseShift().

The documentation for this class was generated from the following files:

SimCalorimetry/HcalTestBeam/interface/HcalTBDigiProducer.h
SimCalorimetry/HcalTestBeam/src/HcalTBDigiProducer.cc

Generated on Tue Jun 9 18:23:58 2009 for CMSSW by

1.5.4


Public Member Functions
	HcalTBDigiProducer (const edm::ParameterSet &ps)
virtual void	produce (edm::Event &e, const edm::EventSetup &c)
	Produces the EDM products,.
virtual	~HcalTBDigiProducer ()
Private Types
typedef CaloTDigitizer < HBHEDigitizerTraits >	HBHEDigitizer
	Reconstruction algorithm.
typedef CaloTDigitizer < HODigitizerTraits >	HODigitizer
Private Member Functions
void	checkGeometry (const edm::EventSetup &eventSetup)
	make sure the digitizer has the correct list of all cells that exist in the geometry
void	setPhaseShift (const DetId &detId)
void	sortHits (const edm::PCaloHitContainer &hits)
	fills the vectors for each subdetector
void	updateGeometry ()
Private Attributes
bool	doPhaseShift
std::string	ecalTBInfoLabel
std::vector< DetId >	hbheCells
std::vector< DetId >	hoCells
HcalAmplifier *	theAmplifier
HcalCoderFactory *	theCoderFactory
HcalElectronicsSim *	theElectronicsSim
const CaloGeometry *	theGeometry
HBHEDigitizer *	theHBHEDigitizer
HBHEHitFilter	theHBHEHitFilter
std::vector< PCaloHit >	theHBHEHits
CaloHitResponse *	theHBHEResponse
CaloVShape *	theHcalIntegratedShape
CaloVShape *	theHcalShape
HcalHitCorrection *	theHitCorrection
HODigitizer *	theHODigitizer
HOHitFilter	theHOHitFilter
std::vector< PCaloHit >	theHOHits
CaloHitResponse *	theHOResponse
HcalTBSimParameterMap *	theParameterMap
double	thisPhaseShift
double	tunePhaseShift