#include <HGCHEbackDigitizer.h>

Inheritance diagram for HGCHEbackDigitizer:

Public Member Functions
	HGCHEbackDigitizer (const edm::ParameterSet &ps)

void	runDigitizer (std::auto_ptr< HGCHEDigiCollection > &digiColl, HGCSimHitDataAccumulator &simData, uint32_t digitizationType, CLHEP::HepRandomEngine *engine)

	~HGCHEbackDigitizer ()

Public Member Functions inherited from HGCDigitizerBase< HGCHEDataFrame >
	HGCDigitizerBase (const edm::ParameterSet &ps)
	CTOR. More...

void	run (std::auto_ptr< DColl > &digiColl, HGCSimHitDataAccumulator &simData, uint32_t digitizationType, CLHEP::HepRandomEngine *engine)
	steer digitization mode More...

virtual void	runDigitizer (std::auto_ptr< DColl > &coll, HGCSimHitDataAccumulator &simData, uint32_t digitizerType, CLHEP::HepRandomEngine *engine)
	to be specialized by top class More...

void	runShaper (HGCHEDataFrame &dataFrame)
	applies a shape to each time sample and propagates the tails to the subsequent time samples More...

void	runTrivial (std::auto_ptr< DColl > &coll, HGCSimHitDataAccumulator &simData, CLHEP::HepRandomEngine *engine)
	a trivial digitization: sum energies and digitize without noise More...

void	updateOutput (std::auto_ptr< DColl > &coll, HGCHEDataFramerawDataFrame)
	prepares the output according to the number of time samples to produce More...

	~HGCDigitizerBase ()
	DTOR. More...

Private Member Functions
void	runCaliceLikeDigitizer (std::auto_ptr< HGCHEDigiCollection > &digiColl, HGCSimHitDataAccumulator &simData, CLHEP::HepRandomEngine *engine)

Private Attributes
float	nPEperMIP_

float	nTotalPE_

float	sdPixels_

float	xTalk_

Additional Inherited Members
Public Types inherited from HGCDigitizerBase< HGCHEDataFrame >
typedef edm::SortedCollection < HGCHEDataFrame >	DColl

Public Attributes inherited from HGCDigitizerBase< HGCHEDataFrame >
uint32_t	adcThreshold_

int	bxTime_

bool	doTimeSamples_

double	lsbInMIP_

double	mip2noise_

double	mipInKeV_

edm::ParameterSet	myCfg_

double	shaperN_

double	shaperTau_

Detailed Description

Definition at line 8 of file HGCHEbackDigitizer.h.

Constructor & Destructor Documentation

HGCHEbackDigitizer::HGCHEbackDigitizer ( const edm::ParameterSet & ps )

Definition at line 6 of file HGCHEbackDigitizer.cc.

References alignCSCRings::e, cppFunctionSkipper::exception, edm::ParameterSet::getParameter(), nPEperMIP_, nTotalPE_, sdPixels_, and xTalk_.

                                                                 : HGCDigitizerBase(ps) {
   try{
     edm::ParameterSet caliceSpec =  ps.getParameter<edm::ParameterSet>("digiCfg").getParameter<edm::ParameterSet>("caliceSpecific");
     nPEperMIP_ = caliceSpec.getParameter<double>("nPEperMIP");
     nTotalPE_  = caliceSpec.getParameter<double>("nTotalPE");
     xTalk_     = caliceSpec.getParameter<double>("xTalk");
     sdPixels_  = caliceSpec.getParameter<double>("sdPixels");
   }catch(std::exception &e){
     //no need to propagate
   }
 }

HGCHEbackDigitizer::~HGCHEbackDigitizer ( )

Definition at line 86 of file HGCHEbackDigitizer.cc.

86 { }

Member Function Documentation

void HGCHEbackDigitizer::runCaliceLikeDigitizer	(	std::auto_ptr< HGCHEDigiCollection > &	digiColl,
		HGCSimHitDataAccumulator &	simData,
		CLHEP::HepRandomEngine *	engine
	)

private

Definition at line 30 of file HGCHEbackDigitizer.cc.

References create_public_lumi_plots::exp, i, cmsBatch::log, HGCDigitizerBase< HGCHEDataFrame >::lsbInMIP_, bookConverter::max, HGCDigitizerBase< HGCHEDataFrame >::mip2noise_, HGCDigitizerBase< HGCHEDataFrame >::mipInKeV_, nPEperMIP_, nTotalPE_, HGCDigitizerBase< HGCHEDataFrame >::runShaper(), sdPixels_, HGCSample::set(), HGCDataFrame< D, S >::setSample(), HGCDigitizerBase< HGCHEDataFrame >::updateOutput(), x, and xTalk_.

Referenced by runDigitizer().

                                                                                                                                                           {
  
   for(HGCSimHitDataAccumulator::iterator it=simData.begin();
       it!=simData.end(); it++) {
     //init a new data frame
     HGCHEDataFrame newDataFrame( it->first );
 
     for(size_t i=0; i<it->second.size(); i++) {
       //convert total energy GeV->keV->ADC counts
       float totalEn( (it->second)[i]*1e6 );
       
       //convert energy to MIP
       float totalIniMIPs = totalEn/mipInKeV_;
       
       //generate random number of photon electrons
       CLHEP::RandPoissonQ randPoissonQ(*engine, totalIniMIPs*nPEperMIP_);
       uint32_t npe = (uint32_t)randPoissonQ.fire();
       
       //number of pixels    
       float x=exp(-(float)(npe)/(float)(nTotalPE_));
       uint32_t nPixel(0);
       if(xTalk_*x!=1) nPixel=(uint32_t) std::max( float(nTotalPE_*(1-x)/(1-xTalk_*x)), float(0.) );
       
       //update signal
       nPixel=(uint32_t)std::max( float(CLHEP::RandGaussQ::shoot(engine,(float)nPixel,(float)sdPixels_)),float(0.) );
       
       //convert to MIP again and saturate
       float totalMIPs(totalIniMIPs);
       if(nTotalPE_!=nPixel && (nTotalPE_-xTalk_*nPixel)/(nTotalPE_-nPixel)>0 )
     totalMIPs = (nTotalPE_/nPEperMIP_)*log((nTotalPE_-xTalk_*nPixel)/(nTotalPE_-nPixel));
       else
     totalMIPs = 0;
       
       //add noise (in MIPs)
       double noiseMIPs=CLHEP::RandGaussQ::shoot(engine,0.,1./mip2noise_);
       totalMIPs=std::max(float(totalMIPs+noiseMIPs),float(0.));
       
       //round to integer (sample will saturate the value according to available bits)
       uint16_t totalEnInt = floor( totalMIPs / lsbInMIP_ );
           
       //0 gain for the moment
       HGCSample singleSample;
       singleSample.set(0, totalEnInt );
       newDataFrame.setSample(i, singleSample);
 
     }   
       
     //run shaper
     runShaper(newDataFrame);
 
     //prepare the output
     updateOutput(digiColl,newDataFrame);
   } 
 }

void HGCHEbackDigitizer::runDigitizer	(	std::auto_ptr< HGCHEDigiCollection > &	digiColl,
		HGCSimHitDataAccumulator &	simData,
		uint32_t	digitizationType,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 19 of file HGCHEbackDigitizer.cc.

References runCaliceLikeDigitizer().

                                                                                                                                                                           {
   switch(digitizationType) {
   case 1: 
     {
       runCaliceLikeDigitizer(digiColl,simData,engine);
       break;
     }
   }
 }