HGCDigitizerBase< DFr > Class Template Reference

#include <HGCDigitizerBase.h>

Public Types
typedef edm::SortedCollection < DFr >	DColl

Public Member Functions
	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 DFr &rawDataFrame)
	prepares the output according to the number of time samples to produce More...
	~HGCDigitizerBase ()
	DTOR. More...

Protected Attributes
double	bxTime_
bool	doTimeSamples_
float	keV2fC_
edm::ParameterSet	myCfg_
std::unique_ptr < HGCFEElectronics< DFr > >	myFEelectronics_
float	noise_fC_

Detailed Description

template<class DFr>
class HGCDigitizerBase< DFr >

Definition at line 22 of file HGCDigitizerBase.h.

Member Typedef Documentation

template<class DFr>

typedef edm::SortedCollection<DFr> HGCDigitizerBase< DFr >::DColl

Definition at line 25 of file HGCDigitizerBase.h.

Constructor & Destructor Documentation

template<class DFr>

HGCDigitizerBase< DFr >::HGCDigitizerBase ( const edm::ParameterSet &  ps )

inline

CTOR.

Definition at line 30 of file HGCDigitizerBase.h.

                                               {
    bxTime_        = ps.getParameter<double>("bxTime");
    myCfg_         = ps.getParameter<edm::ParameterSet>("digiCfg"); 
    doTimeSamples_ = myCfg_.getParameter< bool >("doTimeSamples");
    if(myCfg_.exists("keV2fC"))   keV2fC_   = myCfg_.getParameter<double>("keV2fC");    
    else                          keV2fC_   = 1.0;
    if(myCfg_.exists("noise_fC")) noise_fC_ = myCfg_.getParameter<double>("noise_fC");
    else                          noise_fC_ = 1.0;
    edm::ParameterSet feCfg = myCfg_.getParameter<edm::ParameterSet>("feCfg");
    myFEelectronics_        = std::unique_ptr<HGCFEElectronics<DFr> >( new HGCFEElectronics<DFr>(feCfg) );
  }

template<class DFr>

HGCDigitizerBase< DFr >::~HGCDigitizerBase ( )

inline

DTOR.

Definition at line 75 of file HGCDigitizerBase.h.

    { };

Member Function Documentation

template<class DFr>

float HGCDigitizerBase< DFr >::keV2fC ( ) const

inline

getters

Definition at line 50 of file HGCDigitizerBase.h.

{ return keV2fC_; }

template<class DFr>

void HGCDigitizerBase< DFr >::run	(	std::auto_ptr< DColl > &	digiColl,
		hgc::HGCSimHitDataAccumulator &	simData,
		uint32_t	digitizationType,
		CLHEP::HepRandomEngine *	engine
	)

steer digitization mode

Definition at line 6 of file HGCDigitizerBase.cc.

                                                                {
  if(digitizationType==0) runSimple(digiColl,simData,engine);
  else                    runDigitizer(digiColl,simData,digitizationType,engine);
}

template<class DFr>

virtual void HGCDigitizerBase< DFr >::runDigitizer ( std::auto_ptr< DColl > &  coll, hgc::HGCSimHitDataAccumulator &  simData, uint32_t  digitizerType, CLHEP::HepRandomEngine *  engine  )

inlinevirtual

to be specialized by top class

Definition at line 67 of file HGCDigitizerBase.h.

  {
    throw cms::Exception("HGCDigitizerBaseException") << " Failed to find specialization of runDigitizer";
  }

template<class DFr>

void HGCDigitizerBase< DFr >::runSimple	(	std::auto_ptr< DColl > &	coll,
		hgc::HGCSimHitDataAccumulator &	simData,
		CLHEP::HepRandomEngine *	engine
	)

a trivial digitization: sum energies and digitize without noise

Definition at line 15 of file HGCDigitizerBase.cc.

References f, i, and bookConverter::max.

                                                                     {
  HGCSimHitData chargeColl,toa;
  for(HGCSimHitDataAccumulator::iterator it=simData.begin();
      it!=simData.end();
      it++) {
    chargeColl.fill(0.f); 
    toa.fill(0.f);
    for(size_t i=0; i<it->second[0].size(); i++) {
      double rawCharge((it->second)[0][i]);
      
      //time of arrival
      toa[i]=(it->second)[1][i];
      if(myFEelectronics_->toaMode()==HGCFEElectronics<DFr>::WEIGHTEDBYE && rawCharge>0) 
        toa[i]=(it->second)[1][i]/rawCharge;
      
      //convert total energy in GeV to charge (fC)
      //double totalEn=rawEn*1e6*keV2fC_;
      double totalCharge=rawCharge;
      
      //add noise (in fC)
      //we assume it's randomly distributed and won't impact ToA measurement
      totalCharge += std::max( CLHEP::RandGaussQ::shoot(engine,0,noise_fC_) , 0. );
      if(totalCharge<0) totalCharge=0;
      
      chargeColl[i]= totalCharge;
    }
    
    //run the shaper to create a new data frame
    DFr rawDataFrame( it->first );
    myFEelectronics_->runShaper(rawDataFrame, chargeColl, toa, engine);
    
    //update the output according to the final shape
    updateOutput(coll,rawDataFrame);
  }  
}

template<class DFr>

float HGCDigitizerBase< DFr >::tdcOnset ( ) const

inline

Definition at line 52 of file HGCDigitizerBase.h.

{ return myFEelectronics_->getTDCOnset(); }

template<class DFr>

bool HGCDigitizerBase< DFr >::toaModeByEnergy ( ) const

inline

Definition at line 51 of file HGCDigitizerBase.h.

{ return (myFEelectronics_->toaMode()==HGCFEElectronics<DFr>::WEIGHTEDBYE); }

template<class DFr>

void HGCDigitizerBase< DFr >::updateOutput	(	std::auto_ptr< DColl > &	coll,
		const DFr &	rawDataFrame
	)

prepares the output according to the number of time samples to produce

Definition at line 54 of file HGCDigitizerBase.cc.

References data, universalConfigTemplate::mode, HGCSample::set(), and dtDQMClient_cfg::threshold.

                                                                   {
  int itIdx(9);
  if(rawDataFrame.size()<=itIdx+2) return;
  
  DFr dataFrame( rawDataFrame.id() );
  dataFrame.resize(5);
  bool putInEvent(false);
  HGCSample singleSample;
  for(int it=0;it<5; it++) {    
    singleSample.set(rawDataFrame[itIdx-2+it].threshold(),
                     rawDataFrame[itIdx-2+it].mode(),
                     rawDataFrame[itIdx-2+it].toa(),
                     rawDataFrame[itIdx-2+it].data());
    dataFrame.setSample(it, singleSample);
    if(it==2) { putInEvent=rawDataFrame[itIdx-2+it].threshold(); }
  }
  if(putInEvent) coll->push_back(dataFrame);
}

Member Data Documentation

template<class DFr>

double HGCDigitizerBase< DFr >::bxTime_

protected

Definition at line 95 of file HGCDigitizerBase.h.

Referenced by HGCDigitizerBase< HGCEEDataFrame >::HGCDigitizerBase().

template<class DFr>

bool HGCDigitizerBase< DFr >::doTimeSamples_

protected

Definition at line 98 of file HGCDigitizerBase.h.

Referenced by HGCDigitizerBase< HGCEEDataFrame >::HGCDigitizerBase().

template<class DFr>

float HGCDigitizerBase< DFr >::keV2fC_

protected

Definition at line 86 of file HGCDigitizerBase.h.

Referenced by HGCDigitizerBase< HGCEEDataFrame >::HGCDigitizerBase(), and HGCDigitizerBase< HGCEEDataFrame >::keV2fC().

template<class DFr>

edm::ParameterSet HGCDigitizerBase< DFr >::myCfg_

protected

Definition at line 76 of file HGCDigitizerBase.h.

Referenced by HGCDigitizerBase< HGCEEDataFrame >::HGCDigitizerBase().

template<class DFr>

std::unique_ptr<HGCFEElectronics<DFr> > HGCDigitizerBase< DFr >::myFEelectronics_

protected

Definition at line 92 of file HGCDigitizerBase.h.

Referenced by HGCDigitizerBase< HGCEEDataFrame >::HGCDigitizerBase(), HGCDigitizerBase< HGCEEDataFrame >::tdcOnset(), and HGCDigitizerBase< HGCEEDataFrame >::toaModeByEnergy().

template<class DFr>

float HGCDigitizerBase< DFr >::noise_fC_

protected

Definition at line 89 of file HGCDigitizerBase.h.

Referenced by HGCDigitizerBase< HGCEEDataFrame >::HGCDigitizerBase().