#include <HcalTDC.h>

Public Member Functions
unsigned int	getThresholdDAC ()

	HcalTDC (unsigned int thresholdDAC=12)

void	setDbService (const HcalDbService *service)
	the Producer will probably update this every event More...

void	setRandomEngine (CLHEP::HepRandomEngine &engine)

void	setThresholdDAC (unsigned int DAC)

void	timing (const CaloSamples &lf, HcalUpgradeDataFrame &digi) const
	adds timing information to the digi More...

	~HcalTDC ()

Private Member Functions
double	getThreshold (const HcalGenericDetId &detId) const

Private Attributes
unsigned int	theDAC

const HcalDbService *	theDbService

CLHEP::RandGaussQ *	theRandGaussQ

HcalTDCParameters	theTDCParameters

Detailed Description

Definition at line 11 of file HcalTDC.h.

Constructor & Destructor Documentation

HcalTDC::HcalTDC ( unsigned int thresholdDAC = 12 )

Definition at line 6 of file HcalTDC.cc.

                                           : theTDCParameters(), 
                           theDbService(0),
                           theDAC(thresholdDAC),
                           theRandGaussQ(0) {}

HcalTDC::~HcalTDC ( )

Definition at line 11 of file HcalTDC.cc.

References theRandGaussQ.

                   {
   if (theRandGaussQ) delete theRandGaussQ;
 }

Member Function Documentation

double HcalTDC::getThreshold ( const HcalGenericDetId & detId ) const

private

Definition at line 88 of file HcalTDC.cc.

References HcalTDCParameters::deltaT(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCalibrationWidths(), HcalTDCParameters::nbins(), HcalCalibrationWidths::pedestal(), HcalCalibrations::pedestal(), theDAC, theDbService, theRandGaussQ, and theTDCParameters.

Referenced by timing().

                                                                  {
   // subtract off pedestal and noise once
   double pedestal = theDbService->getHcalCalibrations(detId).pedestal(0);
   double pedestalWidth = theDbService->getHcalCalibrationWidths(detId).pedestal(0);
   // here the TDCthreshold_ is a multiple of the least significant bit 
   // for the TDC portion of the QIE.  The nominal reference is 40 uA which is
   // divided by an 8 bit DAC.  This means the least significant bit is 0.135 uA
   // in the TDC circuit or 3.74 uA at the input current.
 
   // the pedestal is assumed to be evenly distributed in time with some
   // random variation.  No hysteresis is implemented here.
 
   double const lsb(3.74);
   return lsb*theDAC - theRandGaussQ->shoot(pedestal,  pedestalWidth)/theTDCParameters.deltaT()/theTDCParameters.nbins();
 }

unsigned int HcalTDC::getThresholdDAC ( )

inline

Definition at line 25 of file HcalTDC.h.

References theDAC.

25 { return theDAC; }

HcalTDC::theDAC

unsigned int theDAC

Definition: HcalTDC.h:32

void HcalTDC::setDbService ( const HcalDbService * service )

the Producer will probably update this every event

Definition at line 108 of file HcalTDC.cc.

References theDbService.

Referenced by HcalElectronicsSim::setDbService().

                                                         {
   theDbService = service;
 }

void HcalTDC::setRandomEngine ( CLHEP::HepRandomEngine & engine )

Definition at line 104 of file HcalTDC.cc.

References theRandGaussQ.

Referenced by HcalElectronicsSim::setRandomEngine().

                                                            {
   theRandGaussQ = new CLHEP::RandGaussQ(engine);
 }

void HcalTDC::setThresholdDAC ( unsigned int DAC )

inline

Definition at line 24 of file HcalTDC.h.

References theDAC.

24 { theDAC = DAC; }

HcalTDC::theDAC

unsigned int theDAC

Definition: HcalTDC.h:32

void HcalTDC::timing	(	const CaloSamples &	lf,
		HcalUpgradeDataFrame &	digi
	)		const

adds timing information to the digi

Definition at line 15 of file HcalTDC.cc.

References HcalUpgradeDataFrame::adc(), HcalTDCParameters::alreadyTransitionCode(), getThreshold(), i, HcalUpgradeDataFrame::id(), HcalTDCParameters::nbins(), HcalTDCParameters::noTransitionCode(), CaloSamples::preciseAt(), CaloSamples::preciseSize(), HcalUpgradeDataFrame::setSample(), CaloSamples::size(), and theTDCParameters.

Referenced by HcalElectronicsSim::analogToDigital().

                                                                             {
 
   float TDC_Threshold = getThreshold(digi.id());
   bool alreadyOn = false;
   int tdcBins = theTDCParameters.nbins();
   // start with a loop over 10 samples
   bool hasTDCValues=true;
   if (lf.preciseSize()==0 ) hasTDCValues=false;
 
   // if (hasTDCValues)
   //   std::cout << digi.id() 
   //          << " threshold: " << TDC_Threshold
   //          << " hasTDCValues: " << hasTDCValues 
   //          << '\n';
   for (int ibin = 0; ibin < lf.size(); ++ibin) {
     /*
     If in a given 25ns bunch/time sample, the pulse is above
     TDC_Thresh already, then TDC_RisingEdge=0 if it was low in the
     last precision bin on the previous bunch crossing, otherwise,
     TDC_RisingEdge=63 if the pulse never crosses the threshold
     having started off, then the special code is 62 and then
     one can still have a TDC_FallingEdge that is valid.  If the pulse
     never falls below threshold having started above threshold (or
     goes above threshold in the bunch crossing and doesn't come down),
     then TDC_FallingEdge=.  If the pulse never went above
     threshold, then TDC_RisingEdge=63 and
     TDC_FallingEdge=62.
     */
     // special codes
     int TDC_RisingEdge = (alreadyOn) ? theTDCParameters.alreadyTransitionCode() 
       : theTDCParameters.noTransitionCode();
     int TDC_FallingEdge = (alreadyOn) ? theTDCParameters.noTransitionCode() : 
       theTDCParameters.alreadyTransitionCode();
     int preciseBegin = ibin * tdcBins;
     int preciseEnd = preciseBegin + tdcBins;
     if ( hasTDCValues) {
       // std::cout << " alreadyOn: " << alreadyOn << '\n';
       for(int i = preciseBegin; i < preciseEnd; ++i) { 
       //    std::cout << " preciseBin: " << i
       //                << " preciseAt(i): " << lf.preciseAt(i);
     if (alreadyOn) {
       if( (i%3 == 0) && (lf.preciseAt(i) < TDC_Threshold) ) {
         alreadyOn = false;
         TDC_FallingEdge = i-preciseBegin;
         // std::cout << " falling ";
       }
     } else {
       if (lf.preciseAt(i) > TDC_Threshold) {
         alreadyOn = true;
         TDC_RisingEdge = i-preciseBegin;
         // the flag for hasn't gone low yet
         TDC_FallingEdge = theTDCParameters.noTransitionCode();
         // std::cout << " rising ";
       }
     }
     // std::cout << '\n';
       }
     }
     // change packing to allow for special codes
     int packedTDC = TDC_RisingEdge + (tdcBins*2) * TDC_FallingEdge;
     digi.setSample(ibin, digi.adc(ibin), packedTDC, true);
     // if ( hasTDCValues) {
     //   std::cout << " sample: " << ibin 
     //      << " adc: " << digi.adc(ibin)
     //      << " fC: " << digi[ibin].nominal_fC()
     //      << " risingEdge: " << TDC_RisingEdge
     //      << " fallingEdge: " << TDC_FallingEdge
     //      << " packedTDC: " << packedTDC
     //      << std::endl;
     // }
   } // loop over bunch crossing bins
 }