#include <PndSdsTotChargeConversion.h>

Inheritance diagram for PndSdsTotChargeConversion:

Public Member Functions
	PndSdsTotChargeConversion (Int_t VerboseLevel=0)
	Default constructor. More...

	PndSdsTotChargeConversion (Double_t tr, Double_t a, Double_t threshold, Double_t clockfrequency=0, Int_t VerboseLevel=0)
	Main constructor. More...

	~PndSdsTotChargeConversion ()
	Destructor. More...

void	StartExecute ()

virtual Double_t	ChargeToDigiValue (Double_t charge)
	Converts a given charge in electrons into the electronics answer e.g. ToT [ns]. More...

virtual Double_t	DigiValueToCharge (Double_t digivalue)
	Converts a given digitized charge into charge in electrons. More...

virtual Double_t	GetPileUpTime (Double_t charge)
	returns the time the capacitor is loaded and therefore the time this MVD pixel/strip is blind for other events More...

virtual Double_t	GetRelativeError (Double_t Charge)

virtual Double_t	GetTimeStamp (Double_t tof, Double_t charge, Double_t MCEventTime)
	absolute time stamp of a hit in ns (clock is taken into account) More...

virtual Double_t	GetTimeWalk (Double_t Charge)

virtual Double_t	GetTimeStampErrorAfterCorrection ()

virtual Double_t	GetTimeStep ()

virtual void	EndExecute ()

virtual Double_t	DigiValueToCharge (PndSdsDigi &digi)

Double_t	GetParameter (TString param)

void	SetParameter (TString param, Double_t value)

Private Member Functions
Double_t	GetTotWC ()

Double_t	DigitizeTime (Double_t time)
	returns the time in [ns] but binned to clock units More...

	ClassDef (PndSdsTotChargeConversion, 1)

Private Attributes
Double_t	fthreshold

Double_t	Qt

Double_t	t1e

Double_t	t2e

Double_t	ftimestep

Double_t	fstarttime

Double_t	fstoptime

Double_t	ftimeoffset

Double_t	ftimewalk

Int_t	fVerboseLevel

TRandom2	fRand

Detailed Description

Class to calculate the deposited charge in one digi <-> time over threshold and vice versa

Author: D.-L.Pohl d.poh.nosp@m.l@fz.nosp@m.-juel.nosp@m.ich..nosp@m.de

This class calculates the time over threshold by using the charge of a pixel/strip. To calculate the TOT a simple capacitor model (integrator) is used. The charge time tr [ns] is assumed to be constant. This leads to a rising straight line. A constant current unloads the capacitor with a ratio of: a [e/ns]. The distance between the points of intersection between the resulting triangle and the threshold line is the TOT.

q | | /\ a | / \ |—/-------—------------— threshold | / \ ----------------------------— q = 0 (charge) | tr| t[ns] <- TOT ->

Definition at line 32 of file PndSdsTotChargeConversion.h.

Constructor & Destructor Documentation

PndSdsTotChargeConversion::PndSdsTotChargeConversion ( Int_t VerboseLevel = 0 )

Default constructor.

Definition at line 5 of file PndSdsTotChargeConversion.cxx.

References ftimestep, fVerboseLevel, PndSdsChargeConversion::GetParameter(), PndSdsChargeConversion::SetParameter(), and StartExecute().

                                                                        : PndSdsChargeConversion(kToT),
 fthreshold(0.),
 Qt(0.),
 t1e(0.),
 t2e(0.),
 ftimestep(0.),
 fstarttime(0.),
 fstoptime(0.),
 ftimeoffset(0.),
 ftimewalk(0.),
 fVerboseLevel(0),
 fRand()
 {
   SetParameter("ftr", 100.); //todo take data from Database
   SetParameter("fa", 60.);
   SetParameter("fth", 3000.);
   SetParameter("fclk", 156.);
   ftimestep = 1. / GetParameter("fclk") * 1000.;
   fVerboseLevel = VerboseLevel;
   StartExecute();
   Warning("constructor","default constructor used");
   if (fVerboseLevel>0){
     std::cout<<"Tot parameter "<<std::endl;
     std::cout<<"  charge time: "<<GetParameter("ftr")<<" ns"<<std::endl;
     std::cout<<"  const. current: "<<GetParameter("fa")<<" e/ns"<<std::endl;
     std::cout<<"  threshold: "<<GetParameter("fth")<<" e"<<std::endl;
     std::cout<<"  clock frequency: "<<GetParameter("fclk")<<" MHz"<<std::endl;
   }
 };

PndSdsTotChargeConversion::PndSdsTotChargeConversion	(	Double_t	tr,
		Double_t	a,
		Double_t	threshold,
		Double_t	clockfrequency = `0`,
		Int_t	VerboseLevel = `0`
	)

Main constructor.

Definition at line 36 of file PndSdsTotChargeConversion.cxx.

References ftimestep, fVerboseLevel, PndSdsChargeConversion::GetParameter(), PndSdsChargeConversion::SetParameter(), and StartExecute().

                                                                                                                                              : PndSdsChargeConversion(kToT),
 fthreshold(0.),
 Qt(0.),
 t1e(0.),
 t2e(0.),
 ftimestep(0.),
 fstarttime(0.),
 fstoptime(0.),
 ftimeoffset(0.),
 ftimewalk(0.),
 fVerboseLevel(0),
 fRand()
 {
   SetParameter("ftr", tr);                              //charge time [ns]
   SetParameter("fa", a);                                        //const. curren [e/ns]
   SetParameter("fth", threshold);
   SetParameter("fclk", clockfrequency);         //the frequency the readout chip works with
   ftimestep = 1. / clockfrequency * 1000.;
   fVerboseLevel = VerboseLevel;
   if (fVerboseLevel>0){
     std::cout<<"Tot parameter "<<std::endl;
     std::cout<<"  charge time: "<<GetParameter("ftr")<<" ns"<<std::endl;
     std::cout<<"  const. current: "<<GetParameter("fa")<<" e/ns"<<std::endl;
     std::cout<<"  threshold: "<<GetParameter("fth")<<" e"<<std::endl;
     std::cout<<"  clock frequency: "<<GetParameter("fclk")<<" MHz"<<std::endl;
   }
   StartExecute();
 };

PndSdsTotChargeConversion::~PndSdsTotChargeConversion ( )

inline

Destructor.

Definition at line 42 of file PndSdsTotChargeConversion.h.

42 {};

Member Function Documentation

Double_t PndSdsTotChargeConversion::ChargeToDigiValue ( Double_t Charge )

virtual

Converts a given charge in electrons into the electronics answer e.g. ToT [ns].

Implements PndSdsChargeConversion.

Definition at line 65 of file PndSdsTotChargeConversion.cxx.

References Double_t, ftimeoffset, fVerboseLevel, PndSdsChargeConversion::GetParameter(), GetTotWC(), Qt, PndSdsChargeConversion::SetParameter(), t1e, and t2e.

                                                                     { //returns the TOT in ns
         if (GetParameter("fa") <= 0){
                 if (fVerboseLevel>0) 
       Error("ConvertChargeToDigiValue(Double_t charge)","const. current is less or equal zero -> now set to 60 e/ns");
                 SetParameter("fa",60.);
         }
   
         Double_t Q = charge;
         Qt = GetParameter("fth");
   //Error handling: if there is a parameter leading to a division by zero Q=Qt=1 is set to prevent this
   
         if (Qt < 0){
                 if (fVerboseLevel>0)
       Error("ConvertChargeToDigiValue(Double_t charge)","threshold is less than zero -> now set to 0 eV");
                 Qt = 0.;
         }
   
         if ((Q <= Qt) or (Q <= 0)){
                 if (fVerboseLevel>0)
       Warning("ConvertChargeToDigiValue(Double_t charge)","charge (%g) is equal or less than threshold (%g) -> zero TOT",Q,Qt);
                 Q = 1.;
                 Qt = 1.;
         }
   
         t1e = (GetParameter("ftr")*Qt/Q+ftimeoffset);                                                   //exact time when signal is over threshold
         t2e = (Q-Qt)/GetParameter("fa")+ftimeoffset+GetParameter("ftr");                //exact time when signal is again below threshold
   
         //return (t2e-t1e);                                             //turn off clock
         return GetTotWC();
 }

PndSdsTotChargeConversion::ClassDef	(	PndSdsTotChargeConversion	,
		1
	)

private

Double_t PndSdsTotChargeConversion::DigitizeTime ( Double_t time )

private

returns the time in [ns] but binned to clock units

Definition at line 154 of file PndSdsTotChargeConversion.cxx.

References ftimestep, and fVerboseLevel.

Referenced by GetTimeStamp().

 {
         Int_t temp = (Int_t)((time) / ftimestep);
         time =  temp *ftimestep + ftimestep;
   
   if (fVerboseLevel>2) std::cout << "temp " << temp << " time " << time << std::endl;
   
         return time;
 }

Double_t PndSdsChargeConversion::DigiValueToCharge ( PndSdsDigi & digi )

virtualinherited

Definition at line 3 of file PndSdsChargeConversion.cxx.

References PndSdsChargeConversion::DigiValueToCharge(), and PndSdsDigi::GetCharge().

                                                                           {
         return DigiValueToCharge(digi.GetCharge());
 }

Double_t PndSdsTotChargeConversion::DigiValueToCharge ( Double_t digi )

virtual

Converts a given digitized charge into charge in electrons.

Implements PndSdsChargeConversion.

Definition at line 218 of file PndSdsTotChargeConversion.cxx.

References fVerboseLevel, PndSdsChargeConversion::GetParameter(), and sqrt().

                                                                        { //returns the charge for the given tot value
         if (digivalue<0){
                 if (fVerboseLevel>0)
       Error("ConvertDigiValueToELoss(Double_t digi)","charge digitization value not calculated properly");
                 return -1;
         }
         return (-GetParameter("fa")*GetParameter("ftr")+GetParameter("fth")+digivalue*GetParameter("fa"))/2.+sqrt( pow( (GetParameter("fa")*GetParameter("ftr")-GetParameter("fth")-digivalue*GetParameter("fa")),2) / 4. + GetParameter("fa")* GetParameter("fth") * GetParameter("ftr"));
 }

virtual void PndSdsChargeConversion::EndExecute ( )

inlinevirtualinherited

Definition at line 21 of file PndSdsChargeConversion.h.

Referenced by PndSdsHybridHitProducer::Exec().

21 {};

Double_t PndSdsChargeConversion::GetParameter ( TString param )

inlineinherited

Definition at line 37 of file PndSdsChargeConversion.h.

References PndSdsChargeConversion::fParams, and PndSdsChargeConversion::it.

Referenced by ChargeToDigiValue(), DigiValueToCharge(), GetPileUpTime(), GetRelativeError(), GetTimeWalk(), and PndSdsTotChargeConversion().

                                         {
       it=fParams.find(param);
       if (it == fParams.end()){
         Error("GetParameter(TString param)","No parameter named: %s",param.Data());
         return -1;
       }
       return it->second;
     };

Double_t PndSdsTotChargeConversion::GetPileUpTime ( Double_t charge )

virtual

returns the time the capacitor is loaded and therefore the time this MVD pixel/strip is blind for other events

Definition at line 214 of file PndSdsTotChargeConversion.cxx.

References PndSdsChargeConversion::GetParameter().

                                                                  {      //returns the time the detector is not sensitive for other events
         return (GetParameter("ftr") + fcharge / GetParameter("fa"));
 }

Double_t PndSdsTotChargeConversion::GetRelativeError ( Double_t Charge )

virtual

Implements PndSdsChargeConversion.

Definition at line 96 of file PndSdsTotChargeConversion.cxx.

References a, Double_t, PndSdsChargeConversion::GetParameter(), Qt, and sqrt().

 {
   // formula from D.Pohl/FZJ his sources?
   // dQ_rec     a             /               (Q_t - a*t_c)/2  +  (Q - Q_t)(t_c/Q - 1/a)*a/2              \.
   // ------  =  - * dt_max * | 1 + ----------------------------------------------------------------------  |
   //  Q_rec     2             \    sqrt{ [ (Q_t - a*t_c)/2 + (Q - Q_t)(t_c/Q - 1/a)*a/2 ]^2 + a*Q_t*t_c } /
   //
   //            a     1      /           2Qt - Qt*Qc/Q - Q                \.
   //         =  - * ----- * | 1 + ---------------------------------------  |
   //            2   f_clk    \    sqrt{ [2Qt - Qt*Qc/Q - Q]^2 + 4*Qc*Qt } /
   //
   
         Double_t Q = Charge;
         Qt = GetParameter("fth"); // threshold
   Double_t a = GetParameter("fa"); // const current
   Double_t tc = GetParameter("ftr");// time to load capacitor fully
   Double_t Qc = a*tc; // storable capacitor charge
   Double_t temp = 2.*Qt - Qt*Qc/Q - Q;
   temp = temp/(sqrt(temp*temp + 4*Qc*Qt));
   temp = 0.5*a*(1+temp)/GetParameter("fclk");
   return temp;
 }

Double_t PndSdsTotChargeConversion::GetTimeStamp	(	Double_t	tof,
		Double_t	charge,
		Double_t	MCEventTime
	)

virtual

absolute time stamp of a hit in ns (clock is taken into account)

Implements PndSdsChargeConversion.

Definition at line 126 of file PndSdsTotChargeConversion.cxx.

References DigitizeTime(), Double_t, ftimestep, ftimewalk, fVerboseLevel, GetTimeWalk(), and Qt.

 {
         ftimewalk=GetTimeWalk(Charge);
   
         Double_t totaltime = ftimewalk;
         Double_t eventtime =  MCEventTime;
         Double_t flighttime = time;
   
         totaltime += flighttime; // [ns]
         totaltime += eventtime;  // [ns]
   
         Double_t digitizedtime = DigitizeTime(totaltime);
   
   if (fVerboseLevel>2){
     std::cout<<"  +++threshold: "<< Qt <<" "<<std::endl;
     std::cout<<"  charge: "<< Charge <<" "<<std::endl;
     std::cout<<"  time since event: "<< time <<" "<<std::endl;
     std::cout<<"  event time: "<< MCEventTime << std::endl;
     std::cout<<"  timewalk: "<< ftimewalk <<std::endl;
     std::cout<<"  total time: "<< totaltime <<std::endl;
     std::cout<<"  time step: "<< ftimestep<<std::endl;
     std::cout<<"  digitized total time:"<< digitizedtime << "+++"<<std::endl;
   }
   return digitizedtime ; //digitalisiert
   //     return totaltime; // nicht digitalisiert
 }

virtual Double_t PndSdsTotChargeConversion::GetTimeStampErrorAfterCorrection ( )

inlinevirtual

Reimplemented from PndSdsChargeConversion.

Definition at line 53 of file PndSdsTotChargeConversion.h.

References ftimestep, and sqrt().

                                                      {
         return ftimestep/sqrt(12);
   }

virtual Double_t PndSdsTotChargeConversion::GetTimeStep ( )

inlinevirtual

Reimplemented from PndSdsChargeConversion.

Definition at line 57 of file PndSdsTotChargeConversion.h.

References ftimestep.

57 {return ftimestep;}

PndSdsTotChargeConversion::ftimestep

Double_t ftimestep

Definition: PndSdsTotChargeConversion.h:66

Double_t PndSdsTotChargeConversion::GetTimeWalk ( Double_t Charge )

virtual

Reimplemented from PndSdsChargeConversion.

Definition at line 164 of file PndSdsTotChargeConversion.cxx.

References Double_t, ftimewalk, fVerboseLevel, PndSdsChargeConversion::GetParameter(), Qt, and PndSdsChargeConversion::SetParameter().

Referenced by GetTimeStamp().

                                                                { // [ns]
         if (GetParameter("fa") <= 0){
     if (fVerboseLevel>0)
       Error("GetTimeWalk(Double_t charge)","const. current is less or equal zero -> now set to 60 e/ns");
     SetParameter("fa",60.);
   }
 
         Double_t Q = Charge;
 
         if (Qt < 0){
                 if (fVerboseLevel>0)
       Error("GetTimeWalk(Double_t charge)","threshold is less than zero -> now set to 0 eV");
                 Qt = 0.;
         }
   
         if ((Q <= Qt) or (Q <= 0)){
                 if (fVerboseLevel>0)
       Warning("GetTimeWalk(Double_t charge)","charge is equal or less than threshold -> zero TOT -> infinity TimeWalk");
                 Q = 1.;
                 Qt = 100000.;
                 ftimewalk = GetParameter("ftr");
                 return ftimewalk;
         }
   
 
         Qt = GetParameter("fth");
         ftimewalk = (GetParameter("ftr")*Qt/Q);
 //      ftimewalk += ftimestep;
         return ftimewalk;
 }

Double_t PndSdsTotChargeConversion::GetTotWC ( )

private

Definition at line 195 of file PndSdsTotChargeConversion.cxx.

References fstarttime, fstoptime, ftimestep, t1e, and t2e.

Referenced by ChargeToDigiValue().

                                             { //calculates start time, stop time with a clock
         fstarttime = (Int_t) ( t1e / ftimestep ) * ftimestep + ftimestep; //quantization of the start signal
         fstoptime = (Int_t) ( t2e / ftimestep ) * ftimestep + ftimestep; //quantization of the stop signal
   /*
    std::cout<<"  start point exact: "<<t1e<<std::endl;
    std::cout<<"  stop point exact: "<<t2e<<std::endl;
    std::cout<<"  time step is: "<<ftimestep<<" ns"<<std::endl;
    std::cout<<"  fclockfrequency is: "<<GetParameter("fclk")<<" MHz"<<std::endl;
    std::cout<<"  time offset is: "<<ftimeoffset<<" ns"<<std::endl;
    std::cout<<"  start point is: "<<fstarttime<<" ns"<<std::endl;
    std::cout<<"  stop point is: "<<fstoptime<<" ns"<<std::endl;
    */
         return (fstoptime-fstarttime);
 }

void PndSdsChargeConversion::SetParameter	(	TString	param,
		Double_t	value
	)

inlineinherited

Definition at line 45 of file PndSdsChargeConversion.h.

References PndSdsChargeConversion::fParams.

Referenced by ChargeToDigiValue(), GetTimeWalk(), and PndSdsTotChargeConversion().

                                                     {
       if (value < 0 ) Error("SetParameter(TString param, Double_t value)","invalid value for param: %s", param.Data());
       fParams.insert(std::pair<TString, Double_t>(param, value));
     };

void PndSdsTotChargeConversion::StartExecute ( )

virtual

Reimplemented from PndSdsChargeConversion.

Definition at line 210 of file PndSdsTotChargeConversion.cxx.

References fRand, ftimeoffset, and ftimestep.

Referenced by PndSdsTotChargeConversion().

                                             {   //function called for different events to set different time offsets
         ftimeoffset = fRand.Uniform(ftimestep);
 }