#include <HGCHEbackDigitizer.h>

Inheritance diagram for HGCHEbackDigitizer:

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

void	runDigitizer (std::unique_ptr< HGCalDigiCollection > &digiColl, hgc::HGCSimHitDataAccumulator &simData, const CaloSubdetectorGeometry theGeom, const std::unordered_set< DetId > &validIds, uint32_t digitizationType, CLHEP::HepRandomEngine engine) override
	to be specialized by top class More...

	~HGCHEbackDigitizer () override

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

float	keV2fC () const
	getters More...

void	run (std::unique_ptr< DColl > &digiColl, hgc::HGCSimHitDataAccumulator &simData, const CaloSubdetectorGeometry theGeom, const std::unordered_set< DetId > &validIds, uint32_t digitizationType, CLHEP::HepRandomEngine engine)
	steer digitization mode More...

void	runSimple (std::unique_ptr< DColl > &coll, hgc::HGCSimHitDataAccumulator &simData, const CaloSubdetectorGeometry theGeom, const std::unordered_set< DetId > &validIds, CLHEP::HepRandomEngine engine)
	a trivial digitization: sum energies and digitize without noise More...

std::array< float, 3 >	tdcForToAOnset () const

float	tdcOnset () const

bool	toaModeByEnergy () const

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

virtual	~HGCDigitizerBase ()
	DTOR. More...

Private Member Functions
void	runCaliceLikeDigitizer (std::unique_ptr< HGCalDigiCollection > &digiColl, hgc::HGCSimHitDataAccumulator &simData, const CaloSubdetectorGeometry theGeom, const std::unordered_set< DetId > &validIds, CLHEP::HepRandomEngine engine)

Private Attributes
float	keV2MIP_

float	noise_MIP_

float	nPEperMIP_

float	nTotalPE_

float	sdPixels_

float	xTalk_

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

typedef HGCalDataFrame	DigiType

Protected Attributes inherited from HGCDigitizerBase< HGCalDataFrame >
double	bxTime_

std::vector< double >	cce_

bool	doTimeSamples_

float	keV2fC_

edm::ParameterSet	myCfg_

std::unique_ptr< HGCFEElectronics< HGCalDataFrame > >	myFEelectronics_

std::vector< float >	noise_fC_

Detailed Description

Definition at line 7 of file HGCHEbackDigitizer.h.

Constructor & Destructor Documentation

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

Definition at line 15 of file HGCHEbackDigitizer.cc.

References looper::cfg, edm::ParameterSet::getParameter(), HGCDigitizerBase< HGCalDataFrame >::keV2fC_, keV2MIP_, noise_MIP_, nPEperMIP_, nTotalPE_, sdPixels_, and xTalk_.

                                                                 : HGCDigitizerBase(ps)
 {
   edm::ParameterSet cfg = ps.getParameter<edm::ParameterSet>("digiCfg");
   keV2MIP_   = cfg.getParameter<double>("keV2MIP");
   this->keV2fC_    = 1.0; //keV2MIP_; // hack for HEB
   noise_MIP_ = cfg.getParameter<edm::ParameterSet>("noise_MIP").getParameter<double>("value");
   nPEperMIP_ = cfg.getParameter<double>("nPEperMIP");
   nTotalPE_  = cfg.getParameter<double>("nTotalPE");
   xTalk_     = cfg.getParameter<double>("xTalk");
   sdPixels_  = cfg.getParameter<double>("sdPixels");
 }

HGCHEbackDigitizer::~HGCHEbackDigitizer ( )

override

Definition at line 99 of file HGCHEbackDigitizer.cc.

100 {

101 }

Member Function Documentation

void HGCHEbackDigitizer::runCaliceLikeDigitizer	(	std::unique_ptr< HGCalDigiCollection > &	digiColl,
		hgc::HGCSimHitDataAccumulator &	simData,
		const CaloSubdetectorGeometry *	theGeom,
		const std::unordered_set< DetId > &	validIds,
		CLHEP::HepRandomEngine *	engine
	)

private

Definition at line 36 of file HGCHEbackDigitizer.cc.

References hgc_digi_utils::addCellMetadata(), constexpr, gather_cfg::cout, debug, f, myMath::fast_expf(), myMath::fast_logf(), objects.autophobj::float, hgc_digi::HGCCellInfo::hit_info, mps_fire::i, keV2MIP_, SiStripPI::max, HGCDigitizerBase< HGCalDataFrame >::myFEelectronics_, noise_MIP_, nPEperMIP_, nTotalPE_, sdPixels_, HGCDigitizerBase< HGCalDataFrame >::updateOutput(), x, and xTalk_.

Referenced by runDigitizer().

 {
   //switch to true if you want to print some details
   constexpr bool debug(false);
 
   HGCSimHitData chargeColl;
 
   // this represents a cell with no signal charge
   HGCCellInfo zeroData;
   zeroData.hit_info[0].fill(0.f); //accumulated energy
   zeroData.hit_info[1].fill(0.f); //time-of-flight
 
   for( const auto& id : validIds ) {
     chargeColl.fill(0.f);
     HGCSimHitDataAccumulator::iterator it = simData.find(id);
     HGCCellInfo& cell = ( simData.end() == it ? zeroData : it->second );
     addCellMetadata(cell,theGeom,id);
 
     for(size_t i=0; i<cell.hit_info[0].size(); ++i)
       {
     //convert total energy keV->MIP, since converted to keV in accumulator
     const float totalIniMIPs( cell.hit_info[0][i]*keV2MIP_ );
     //std::cout << "energy in MIP: " << std::scientific << totalIniMIPs << std::endl;
 
       //generate random number of photon electrons
       const uint32_t npe = std::floor(CLHEP::RandPoissonQ::shoot(engine,totalIniMIPs*nPEperMIP_));
 
       //number of pixels
       const float x = vdt::fast_expf( -((float)npe)/nTotalPE_ );
       uint32_t nPixel(0);
       if(xTalk_*x!=1) nPixel=(uint32_t) std::max( nTotalPE_*(1.f-x)/(1.f-xTalk_*x), 0.f );
 
       //update signal
       if(sdPixels_ != 0) nPixel = (uint32_t)std::max( CLHEP::RandGaussQ::shoot(engine,(double)nPixel,sdPixels_), 0. );
 
       //convert to MIP again and saturate
           float totalMIPs(0.f), xtalk = 0.f;
           const float peDiff = nTotalPE_ - (float) nPixel;
           if (peDiff != 0.f) {
             xtalk = (nTotalPE_-xTalk_*((float)nPixel)) / peDiff;
             if( xtalk > 0.f && nPEperMIP_ != 0.f)
               totalMIPs = (nTotalPE_/nPEperMIP_)*vdt::fast_logf(xtalk);
           }
 
       //add noise (in MIPs)
       chargeColl[i] = totalMIPs;
       if(noise_MIP_ != 0) chargeColl[i] += std::max( CLHEP::RandGaussQ::shoot(engine,0.,noise_MIP_), 0. );
       if(debug && cell.hit_info[0][i]>0)
         std::cout << "[runCaliceLikeDigitizer] xtalk=" << xtalk << " En=" << cell.hit_info[0][i] << " keV -> " << totalIniMIPs << " raw-MIPs -> " << chargeColl[i] << " digi-MIPs" << std::endl;
     }
 
       //init a new data frame and run shaper
       HGCalDataFrame newDataFrame( id );
       this->myFEelectronics_->runTrivialShaper( newDataFrame, chargeColl, 1 );
 
       //prepare the output
       this->updateOutput(digiColl,newDataFrame);
     }
 }

void HGCHEbackDigitizer::runDigitizer	(	std::unique_ptr< HGCalDigiCollection > &	coll,
		hgc::HGCSimHitDataAccumulator &	simData,
		const CaloSubdetectorGeometry *	theGeom,
		const std::unordered_set< DetId > &	validIds,
		uint32_t	digitizerType,
		CLHEP::HepRandomEngine *	engine
	)

overridevirtual

to be specialized by top class

Reimplemented from HGCDigitizerBase< HGCalDataFrame >.

Definition at line 28 of file HGCHEbackDigitizer.cc.

References runCaliceLikeDigitizer().

 {
   runCaliceLikeDigitizer(digiColl,simData,theGeom,validIds,engine);
 }

Member Data Documentation

float HGCHEbackDigitizer::keV2MIP_

private

Definition at line 20 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::noise_MIP_

private

Definition at line 20 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::nPEperMIP_

private

Definition at line 21 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::nTotalPE_

private

Definition at line 21 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::sdPixels_

private

Definition at line 21 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::xTalk_

private

Definition at line 21 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation