HGCHEbackDigitizer Class Reference

#include <HGCHEbackDigitizer.h>

Inheritance diagram for HGCHEbackDigitizer:

Public Member Functions
	HGCHEbackDigitizer (const edm::ParameterSet &ps)
void	runDigitizer (std::auto_ptr< HGCHEDigiCollection > &digiColl, hgc::HGCSimHitDataAccumulator &simData, uint32_t digitizationType, CLHEP::HepRandomEngine *engine)
	~HGCHEbackDigitizer ()
Public Member Functions inherited from HGCDigitizerBase< HGCHEDataFrame >
	HGCDigitizerBase (const edm::ParameterSet &ps)
	CTOR. More...
float	keV2fC () const
	getters More...
void	run (std::auto_ptr< DColl > &digiColl, hgc::HGCSimHitDataAccumulator &simData, uint32_t digitizationType, CLHEP::HepRandomEngine *engine)
	steer digitization mode More...
virtual void	runDigitizer (std::auto_ptr< DColl > &coll, hgc::HGCSimHitDataAccumulator &simData, uint32_t digitizerType, CLHEP::HepRandomEngine *engine)
	to be specialized by top class More...
void	runSimple (std::auto_ptr< DColl > &coll, hgc::HGCSimHitDataAccumulator &simData, CLHEP::HepRandomEngine *engine)
	a trivial digitization: sum energies and digitize without noise More...
float	tdcOnset () const
bool	toaModeByEnergy () const
void	updateOutput (std::auto_ptr< DColl > &coll, const HGCHEDataFrame &rawDataFrame)
	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, hgc::HGCSimHitDataAccumulator &simData, 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< HGCHEDataFrame >
typedef edm::SortedCollection < HGCHEDataFrame >	DColl
Protected Attributes inherited from HGCDigitizerBase< HGCHEDataFrame >
double	bxTime_
bool	doTimeSamples_
float	keV2fC_
edm::ParameterSet	myCfg_
std::unique_ptr < HGCFEElectronics < HGCHEDataFrame > >	myFEelectronics_
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 11 of file HGCHEbackDigitizer.cc.

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

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

HGCHEbackDigitizer::~HGCHEbackDigitizer

(

)

Definition at line 83 of file HGCHEbackDigitizer.cc.

{
}

Member Function Documentation

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

private

Definition at line 30 of file HGCHEbackDigitizer.cc.

References constexpr, gather_cfg::cout, debug, f, myMath::fast_expf(), myMath::fast_logf(), i, keV2MIP_, bookConverter::max, HGCDigitizerBase< HGCHEDataFrame >::myFEelectronics_, noise_MIP_, nPEperMIP_, nTotalPE_, sdPixels_, HGCDigitizerBase< HGCHEDataFrame >::updateOutput(), x, and xTalk_.

Referenced by runDigitizer().

{

  //switch to true if you want to print some details
  constexpr bool debug(false);
  
  HGCSimHitData chargeColl;
  for(HGCSimHitDataAccumulator::iterator it=simData.begin();
      it!=simData.end();
      it++)
    {
      chargeColl.fill(0.f);      
      for(size_t i=0; i<it->second[0].size(); ++i)
    {          
      //convert total energy keV->MIP, since converted to keV in accumulator
      float totalIniMIPs( (it->second)[0][i]*keV2MIP_ );
          //std::cout << "energy in MIP: " << std::scientific << totalIniMIPs << std::endl;

      //generate random number of photon electrons
      uint32_t npe = std::floor(CLHEP::RandPoissonQ::shoot(engine,totalIniMIPs*nPEperMIP_));
          
      //number of pixels    
      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
      nPixel = (uint32_t)std::max( CLHEP::RandGaussQ::shoot(engine,(double)nPixel,sdPixels_), 0. );
                
      //convert to MIP again and saturate
      float totalMIPs(totalIniMIPs);
          const float xtalk = (nTotalPE_-xTalk_*((float)nPixel))/(nTotalPE_-((float)nPixel));
      if( nTotalPE_ != nPixel && xtalk > 0. )
        totalMIPs = (nTotalPE_/nPEperMIP_)*vdt::fast_logf(xtalk);
      else
        totalMIPs = 0.f;
      
      //add noise (in MIPs)
      chargeColl[i] = totalMIPs+std::max( CLHEP::RandGaussQ::shoot(engine,0.,noise_MIP_), 0. );
      if(debug && (it->second)[0][i]>0) 
        std::cout << "[runCaliceLikeDigitizer] xtalk=" << xtalk << " En=" << (it->second)[0][i] << " keV -> " << totalIniMIPs << " raw-MIPs -> " << chargeColl[i] << " digi-MIPs" << std::endl;          
    }   
      
      //init a new data frame and run shaper
      HGCHEDataFrame newDataFrame( it->first );
      myFEelectronics_->runTrivialShaper( newDataFrame, chargeColl );

      //prepare the output
      updateOutput(digiColl,newDataFrame);
    } 
}

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

Definition at line 24 of file HGCHEbackDigitizer.cc.

References runCaliceLikeDigitizer().

{
  runCaliceLikeDigitizer(digiColl,simData,engine);
}

Member Data Documentation

float HGCHEbackDigitizer::keV2MIP_

private

Definition at line 18 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::noise_MIP_

private

Definition at line 18 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::nPEperMIP_

private

Definition at line 19 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::nTotalPE_

private

Definition at line 19 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::sdPixels_

private

Definition at line 19 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().

float HGCHEbackDigitizer::xTalk_

private

Definition at line 19 of file HGCHEbackDigitizer.h.

Referenced by HGCHEbackDigitizer(), and runCaliceLikeDigitizer().