HGCHEbackDigitizer Class Reference

#include <HGCHEbackDigitizer.h>

Inheritance diagram for HGCHEbackDigitizer:

Public Member Functions
	HGCHEbackDigitizer (const edm::ParameterSet &ps)
void	runDigitizer (std::unique_ptr< HGCBHDigiCollection > &digiColl, hgc::HGCSimHitDataAccumulator &simData, const CaloSubdetectorGeometry *theGeom, const std::unordered_set< DetId > &validIds, uint32_t digitizationType, CLHEP::HepRandomEngine *engine) override
	~HGCHEbackDigitizer ()
Public Member Functions inherited from HGCDigitizerBase< HGCBHDataFrame >
	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...
virtual void	runDigitizer (std::unique_ptr< DColl > &coll, hgc::HGCSimHitDataAccumulator &simData, const CaloSubdetectorGeometry *theGeom, const std::unordered_set< DetId > &validIds, uint32_t digitizerType, CLHEP::HepRandomEngine *engine)
	to be specialized by top class 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...
float	tdcOnset () const
bool	toaModeByEnergy () const
void	updateOutput (std::unique_ptr< DColl > &coll, const HGCBHDataFrame &rawDataFrame)
	prepares the output according to the number of time samples to produce More...
	~HGCDigitizerBase ()
	DTOR. More...

Private Member Functions
void	runCaliceLikeDigitizer (std::unique_ptr< HGCBHDigiCollection > &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< HGCBHDataFrame >
typedef edm::SortedCollection< HGCBHDataFrame >	DColl
typedef HGCBHDataFrame	DigiType
Protected Attributes inherited from HGCDigitizerBase< HGCBHDataFrame >
double	bxTime_
bool	doTimeSamples_
float	keV2fC_
edm::ParameterSet	myCfg_
std::unique_ptr< HGCFEElectronics< HGCBHDataFrame > >	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 54 of file HGCHEbackDigitizer.cc.

References looper::cfg, edm::ParameterSet::getParameter(), HGCDigitizerBase< HGCBHDataFrame >::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 137 of file HGCHEbackDigitizer.cc.

{
}

Member Function Documentation

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

private

Definition at line 75 of file HGCHEbackDigitizer.cc.

References constexpr, gather_cfg::cout, debug, f, myMath::fast_expf(), myMath::fast_logf(), objects.autophobj::float, hgc_digi::HGCCellInfo::hit_info, i, keV2MIP_, hpstanc_transforms::max, HGCDigitizerBase< HGCBHDataFrame >::myFEelectronics_, noise_MIP_, nPEperMIP_, nTotalPE_, sdPixels_, HGCDigitizerBase< HGCBHDataFrame >::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
      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+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
      HGCBHDataFrame newDataFrame( id );
      myFEelectronics_->runTrivialShaper( newDataFrame, chargeColl, 1 );

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

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

override

Definition at line 67 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().