#include <PndSdsFE.h>

Inheritance diagram for PndSdsFE:

Public Member Functions
	PndSdsFE ()

	PndSdsFE (PndSdsFE &other)

virtual	~PndSdsFE ()

virtual double	GetTotFromCharge (Double_t charge)

virtual double	GetChargeFromTot (double tot)

virtual double	GetTimeWalkFromCharge (double charge)

virtual double	GetTimeWalkFromTot (double tot)

virtual double	GetTimeBackToBaseline (double charge)

virtual double	GetTimeStamp (double eventtime, double tof, double charge)

virtual double	GetTimeStep ()

virtual void	SetParameter (TString parName, double parValue)

virtual void	SetThreshold (double threshold)

PndSdsFE &	operator= (PndSdsFE &other)

Protected Member Functions
double	DigitizeTime (double time)

void	CalcTimeOffSet ()

void	CreateInterpolatorList ()

void	SaveInterpolatorList (std::vector< double > charge, std::vector< double > tot)

void	LoadInterpolatorList ()

void	GetInterpolatorList ()

Protected Attributes
PndSdsFEAmpModelSimple *	fFrontEndModel

int	fNumberOfSupportPoints

int	fNumberOfMaxElectons

TF1 *	fFunction

double	fMaximumAmplitude

double	fThreshold

double	fBaselineEpsilon

double	fToF

double	fEventTime

double	fTimeOffSet

double	fTimeStep

double	fFunctionRange

TRandom2	fRand

int	i

double	stepsize

std::vector< double >	fCharge_list

std::vector< double >	fTot_list

ROOT::Math::Interpolator *	inter

Private Member Functions
	ClassDef (PndSdsFE, 1)

Detailed Description

Definition at line 24 of file PndSdsFE.h.

Constructor & Destructor Documentation

PndSdsFE::PndSdsFE ( )

inline

Definition at line 26 of file PndSdsFE.h.

References CalcTimeOffSet().

                   : fFrontEndModel(0), fNumberOfSupportPoints(0), fNumberOfMaxElectons(0), fFunction(0), fMaximumAmplitude(0),fThreshold(0), fBaselineEpsilon(1), fToF(0), fEventTime(0),       fTimeOffSet(0), fTimeStep(0),fFunctionRange(0), fRand(0), i(0), stepsize(0), fCharge_list(0), fTot_list(0), inter(0) 
         {
                 CalcTimeOffSet();
         };

PndSdsFE::PndSdsFE ( PndSdsFE & other )

inline

Definition at line 31 of file PndSdsFE.h.

                            :TObject(other),
       fFrontEndModel(other.fFrontEndModel),
       fNumberOfSupportPoints(other.fNumberOfSupportPoints),
       fNumberOfMaxElectons(other.fNumberOfMaxElectons),
       fFunction(other.fFunction),
       fMaximumAmplitude(other.fMaximumAmplitude),
       fThreshold(other.fThreshold),
       fBaselineEpsilon(other.fBaselineEpsilon),
       fToF(other.fToF),
       fEventTime(other.fEventTime),
           fTimeOffSet(other.fTimeOffSet),
       fTimeStep(other.fTimeStep),
       fFunctionRange(other.fFunctionRange),
       fRand(other.fRand),
       i(other.i),
       stepsize(other.stepsize),
       fCharge_list(other.fCharge_list),
       fTot_list(other.fTot_list),
       inter(other.inter)
   {};

PndSdsFE::~PndSdsFE ( )

virtual

Definition at line 16 of file PndSdsFE.cxx.

References fFrontEndModel, and fFunction.

                     {
         delete fFrontEndModel;
         delete fFunction;
 }

Member Function Documentation

void PndSdsFE::CalcTimeOffSet ( )

protected

Definition at line 100 of file PndSdsFE.cxx.

References fRand, fTimeOffSet, and fTimeStep.

Referenced by PndSdsFE().

                              {
         fTimeOffSet = fRand.Uniform(fTimeStep);
 }

PndSdsFE::ClassDef	(	PndSdsFE	,
		1
	)

private

void PndSdsFE::CreateInterpolatorList ( )

protected

Definition at line 104 of file PndSdsFE.cxx.

References fCharge_list, fNumberOfMaxElectons, fNumberOfSupportPoints, fThreshold, fTot_list, GetTotFromCharge(), i, printf(), SaveInterpolatorList(), and stepsize.

Referenced by GetInterpolatorList().

                                      {
 
         fNumberOfSupportPoints = 2000;
         fNumberOfMaxElectons = 630000;
         stepsize = fNumberOfMaxElectons/fNumberOfSupportPoints;
 
         //printf("stepsize %f \n",stepsize);
 
         for( i=0; i < fNumberOfSupportPoints; i++){
                 if(i*stepsize<fThreshold){
                                 continue;
                 }
                 else{
 
                 fTot_list.push_back(i*stepsize);
                 fCharge_list.push_back(GetTotFromCharge(i*stepsize));
 
                 //printf("x-value %f ,y-value %f \n",fCharge_list.back(),fTot_list.back());
                 //printf("x-value %f ,y-value %f \n",fCharge_list[],fTot_list.back());
 
                         if(fTot_list.back()< fTot_list.back()-1) {
                                 printf("++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++%d \n",i );
                                 printf("x-value[i-3] %f ,x-value[i-2] %f \n",fTot_list.back()-3,fTot_list.back()-2);
 
                                 printf("x-value[i] %f ,x-value[i-1] %f \n",fTot_list.back(),fTot_list.back());
                         }
                 }
         }
 
         SaveInterpolatorList(fCharge_list,fTot_list);
 }

double PndSdsFE::DigitizeTime ( double time )

protected

Definition at line 90 of file PndSdsFE.cxx.

References fTimeStep.

Referenced by GetTimeStamp(), and GetTotFromCharge().

                                         {
 
                 return TMath::Floor(( time / fTimeStep )) * fTimeStep + fTimeStep;
 }

double PndSdsFE::GetChargeFromTot ( double tot )

virtual

Definition at line 65 of file PndSdsFE.cxx.

References inter.

Referenced by GetTimeWalkFromTot().

                                            {
         return   inter->Eval(tot);
 }

void PndSdsFE::GetInterpolatorList ( )

protected

Definition at line 175 of file PndSdsFE.cxx.

References CreateInterpolatorList(), fFrontEndModel, fThreshold, PndSdsFEAmpModelSimple::GetModelNumber(), LoadInterpolatorList(), and s.

Referenced by PndSdsFESimple::PndSdsFESimple().

                                   {
         // is list available?
 
         char s[20];
         sprintf(s,"interpolatorlist-modelnumber_%i-thr_%g.root",fFrontEndModel->GetModelNumber(),fThreshold);
 
          TFile *fInFile = new TFile(s,"read");
 
          if(fInFile->IsZombie()) {
                  std::cout << "No interpolator list to load, create one... " << std::endl;
                  CreateInterpolatorList();
                  LoadInterpolatorList();
                  return;
          }
          else
          {
                  std::cout << "Loading Interpolator List " << std::endl;
                  LoadInterpolatorList();
                  return;
          }
 
 }

double PndSdsFE::GetTimeBackToBaseline ( double charge )

virtual

Definition at line 84 of file PndSdsFE.cxx.

References fBaselineEpsilon, fFunction, and fMaximumAmplitude.

                                                    {
         fFunction->SetParameter("charge",Charge);
         fMaximumAmplitude = fFunction->GetMaximumX();
         return fFunction->GetX(fBaselineEpsilon,fMaximumAmplitude,fFunction->GetXmax());
 }

double PndSdsFE::GetTimeStamp	(	double	eventtime,
		double	tof,
		double	charge
	)

virtual

Definition at line 95 of file PndSdsFE.cxx.

References DigitizeTime(), and GetTimeWalkFromCharge().

                                                                         {
         //GetChargeToTimeWalk(charge);
         return DigitizeTime(GetTimeWalkFromCharge(charge)+eventtime+tof);
 
 }

virtual double PndSdsFE::GetTimeStep ( )

inlinevirtual

Definition at line 60 of file PndSdsFE.h.

References fTimeStep.

Referenced by PndSdsTimeWalkCorrTask::Exec(), and PndSdsTimeWalkCorrTask::FinishTask().

60 {return fTimeStep;}

PndSdsFE::fTimeStep

double fTimeStep

Definition: PndSdsFE.h:111

double PndSdsFE::GetTimeWalkFromCharge ( double charge )

virtual

Definition at line 77 of file PndSdsFE.cxx.

References fFunction, fMaximumAmplitude, and fThreshold.

Referenced by GetTimeStamp().

                                                    {
         fFunction->SetParameter("charge",Charge);
         fMaximumAmplitude = fFunction->GetMaximumX();
         double result = fFunction->GetX(fThreshold,0,fMaximumAmplitude); //(xmin<x<xmax)
         return result;
 }

double PndSdsFE::GetTimeWalkFromTot ( double tot )

virtual

Definition at line 69 of file PndSdsFE.cxx.

References fFunction, fMaximumAmplitude, fThreshold, and GetChargeFromTot().

Referenced by PndSdsTimeWalkCorrTask::Exec(), and PndSdsTimeWalkCorrTask::FinishTask().

                                              {
         fFunction->SetParameter("charge",GetChargeFromTot(tot));
         fMaximumAmplitude = fFunction->GetMaximumX();
         return fFunction->GetX(fThreshold,0,fMaximumAmplitude); //(xmin<x<xmax)
 
 }

double PndSdsFE::GetTotFromCharge ( Double_t charge )

virtual

Definition at line 26 of file PndSdsFE.cxx.

References DigitizeTime(), fFunction, fMaximumAmplitude, fThreshold, and fTimeOffSet.

Referenced by CreateInterpolatorList().

                                                 {
 
         fFunction->SetParameter("charge",Charge);
         fMaximumAmplitude = fFunction->GetMaximumX();
 
         if (fFunction->GetParameter("constantcurrent") <= 0){
                 Error("PndSdsFE::GetTotFromCharge","const. current is less or equal zero -> now set to 60 e/ns");
                 fFunction->SetParameter("constantcurrent",60);
         }
 
         if (fFunction->GetParameter("chargetime")<= 0){
                 Error("PndSdsFE::GetTotFromCharge","charge time is less than zero -> now set to 100 ns");
                 fFunction->SetParameter("chargetime",100);
         }
 
         if (fThreshold < 0){
                 Error("PndSdsFE::GetTotFromCharge","Threshold is less than zero -> now set to 0 e");
                 fThreshold = 0;
         }
 
         if (fFunction->GetParameter("charge")<= fThreshold){
                 Warning("PndSdsFE::GetTotFromCharge","charge is equal or less than threshold -> zero TOT");
                 return 0;
         }
 
         double start_exakt=fFunction->GetX(fThreshold,0,fMaximumAmplitude );
         double stop_exakt=fFunction->GetX(fThreshold,fMaximumAmplitude,fFunction->GetXmax());
 
         double start=DigitizeTime(start_exakt+fTimeOffSet);
         double stop=DigitizeTime(stop_exakt+fTimeOffSet);
 //      printf("start exakt:%f   \n",start_exakt);
 //      printf("start      :%f   \n",start);
 //      printf("stop exakt:%f    \n",stop_exakt);
 //      printf("stop       :%f    \n",stop);
 
                 return stop-start;              // with digitization
                 //return stop_exakt - start_exakt;  // without digitization
 }

void PndSdsFE::LoadInterpolatorList ( )

protected

Definition at line 198 of file PndSdsFE.cxx.

References fCharge_list, fFrontEndModel, fThreshold, fTot_list, PndSdsFEAmpModelSimple::GetModelNumber(), and inter.

Referenced by GetInterpolatorList().

                                    {
 
         char dd[20];
         sprintf(dd,"interpolatorlist-modelnumber_%i-thr_%g.root",fFrontEndModel->GetModelNumber(),fThreshold);
 
          TFile *fInFile = new TFile(dd);
 
          if(fInFile->IsZombie()) {
                  std::cout << "PndSdsFE::LoadInterpolatorList(): error loading root file " << std::endl;
 
                  return;
          }
          else
          {
                  std::cout << "PndSdsFE::LoadInterpolatorList(): Loading list done " << std::endl;
 
          }
 
          TVectorD *charge = (TVectorD*)fInFile->Get("charge");
          TVectorD *tot = (TVectorD*)fInFile->Get("tot");
 
          int size_charge = charge->GetNoElements();
          int size_tot = tot->GetNoElements();
 
          fCharge_list.clear();
          fCharge_list.clear();
 
                 for(int bxx=0; bxx < size_charge;bxx++)
                 {
                         //fCharge_list.push_back(charge[bxx]);
                         fCharge_list.push_back((*charge)[bxx]);
 
                 }
 
                 for(int byy=0; byy < size_tot;byy++)
                 {
                 //      fTot_list.push_back(tot);
                                 fTot_list.push_back((*tot)[byy]);
                 }
 
                 inter =  new ROOT::Math::Interpolator(fCharge_list,fTot_list, ROOT::Math::Interpolation::kCSPLINE);
 }

PndSdsFE& PndSdsFE::operator= ( PndSdsFE & other )

inline

Definition at line 64 of file PndSdsFE.h.

References fBaselineEpsilon, fCharge_list, fEventTime, fFrontEndModel, fFunction, fFunctionRange, fMaximumAmplitude, fNumberOfMaxElectons, fNumberOfSupportPoints, fRand, fThreshold, fTimeOffSet, fTimeStep, fToF, fTot_list, i, inter, and stepsize.

   {
     if (this != &other) // protect against invalid self-assignment
     {   
       fFrontEndModel=other.fFrontEndModel;
       fNumberOfSupportPoints=other.fNumberOfSupportPoints;
       fNumberOfMaxElectons=other.fNumberOfMaxElectons;
       fFunction=other.fFunction;
       fMaximumAmplitude=other.fMaximumAmplitude;
       fThreshold=other.fThreshold;
       fBaselineEpsilon=other.fBaselineEpsilon;
       fToF=other.fToF;
       fEventTime=other.fEventTime;
           fTimeOffSet=other.fTimeOffSet;
       fTimeStep=other.fTimeStep;
       fFunctionRange=other.fFunctionRange;
       fRand=other.fRand;
       i=other.i;
       stepsize=other.stepsize;
       fCharge_list=other.fCharge_list;
       fTot_list=other.fTot_list;
       inter=other.inter;
     }
     return *this;
   };

void PndSdsFE::SaveInterpolatorList	(	std::vector< double >	charge,
		std::vector< double >	tot
	)

protected

Definition at line 136 of file PndSdsFE.cxx.

References f1, fFrontEndModel, fThreshold, PndSdsFEAmpModelSimple::GetModelNumber(), and s.

Referenced by CreateInterpolatorList().

                                                                                   {
 
         // To save an object to a root file it needs to inherit from TObject.
         // std::vector does not inherit from TObject, but TVectorD does.
         // so i convert the std::vector to a TVectorD and save them to a root file.
 
         int size_x = charge.size();
         int size_y = tot.size();
 
         TVectorD charge_root(size_x);
         TVectorD tot_root(size_y);
 
         for(int bx=0; bx < size_x;bx++)
         {
                 charge_root[bx] = charge.at(bx);
         }
 
         for(int by=0; by < size_x;by++)
         {
                 tot_root[by] = tot.at(by);
         }
 
         // in den namen der interpolatorliste sollte die threshold und das model hinein. sind diese beiden angaben identisch
         // ist die interpolatorliste eindeutig.
 
         char s[20];
 
         sprintf(s,"interpolatorlist-modelnumber_%i-thr_%g.root",fFrontEndModel->GetModelNumber(),fThreshold);
 
         TFile f1(s,"RECREATE");
 
         TGraph* graph = new TGraph(charge.size(), &charge.at(0),&tot.at(0));
         graph->Draw("AP");
         graph->Write();
 
         charge_root.Write("charge");
         tot_root.Write("tot");
 }

void PndSdsFE::SetParameter	(	TString	parName,
		double	parValue
	)

virtual

Definition at line 22 of file PndSdsFE.cxx.

References fFunction.

                                                            {
         fFunction->SetParameter(parName.Data(), parValue);
 }

virtual void PndSdsFE::SetThreshold ( double threshold )

inlinevirtual

Definition at line 62 of file PndSdsFE.h.

References fThreshold, and threshold.

62 {fThreshold = threshold;}

threshold

double threshold

Definition: outdated/dc1/stt/runreco.C:105

PndSdsFE::fThreshold

double fThreshold

Definition: PndSdsFE.h:106

Member Data Documentation

double PndSdsFE::fBaselineEpsilon

protected

Definition at line 107 of file PndSdsFE.h.

Referenced by GetTimeBackToBaseline(), and operator=().

std::vector<double> PndSdsFE::fCharge_list

protected

Definition at line 117 of file PndSdsFE.h.

Referenced by CreateInterpolatorList(), LoadInterpolatorList(), and operator=().

double PndSdsFE::fEventTime

protected

Definition at line 109 of file PndSdsFE.h.

Referenced by operator=().

PndSdsFEAmpModelSimple* PndSdsFE::fFrontEndModel

protected

Definition at line 99 of file PndSdsFE.h.

Referenced by GetInterpolatorList(), LoadInterpolatorList(), operator=(), PndSdsFESimple::PndSdsFESimple(), SaveInterpolatorList(), and ~PndSdsFE().

TF1* PndSdsFE::fFunction

protected

Definition at line 104 of file PndSdsFE.h.

Referenced by GetTimeBackToBaseline(), GetTimeWalkFromCharge(), GetTimeWalkFromTot(), GetTotFromCharge(), operator=(), PndSdsFESimple::PndSdsFESimple(), PndSdsFESimple::SetParameter(), SetParameter(), and ~PndSdsFE().

double PndSdsFE::fFunctionRange

protected

Definition at line 112 of file PndSdsFE.h.

Referenced by operator=(), and PndSdsFESimple::PndSdsFESimple().

double PndSdsFE::fMaximumAmplitude

protected

Definition at line 105 of file PndSdsFE.h.

Referenced by GetTimeBackToBaseline(), GetTimeWalkFromCharge(), GetTimeWalkFromTot(), GetTotFromCharge(), and operator=().

int PndSdsFE::fNumberOfMaxElectons

protected

Definition at line 102 of file PndSdsFE.h.

Referenced by CreateInterpolatorList(), and operator=().

int PndSdsFE::fNumberOfSupportPoints

protected

Definition at line 101 of file PndSdsFE.h.

Referenced by CreateInterpolatorList(), and operator=().

TRandom2 PndSdsFE::fRand

protected

Definition at line 113 of file PndSdsFE.h.

Referenced by CalcTimeOffSet(), and operator=().

double PndSdsFE::fThreshold

protected

Definition at line 106 of file PndSdsFE.h.

Referenced by CreateInterpolatorList(), GetInterpolatorList(), GetTimeWalkFromCharge(), GetTimeWalkFromTot(), GetTotFromCharge(), LoadInterpolatorList(), operator=(), PndSdsFESimple::PndSdsFESimple(), SaveInterpolatorList(), PndSdsFESimple::SetParameter(), and SetThreshold().

double PndSdsFE::fTimeOffSet

protected

Definition at line 110 of file PndSdsFE.h.

Referenced by CalcTimeOffSet(), GetTotFromCharge(), and operator=().

double PndSdsFE::fTimeStep

protected

Definition at line 111 of file PndSdsFE.h.

Referenced by CalcTimeOffSet(), DigitizeTime(), GetTimeStep(), operator=(), PndSdsFESimple::PndSdsFESimple(), and PndSdsFESimple::SetParameter().

double PndSdsFE::fToF

protected

Definition at line 108 of file PndSdsFE.h.

Referenced by operator=().

std::vector<double> PndSdsFE::fTot_list

protected

Definition at line 118 of file PndSdsFE.h.

Referenced by CreateInterpolatorList(), LoadInterpolatorList(), and operator=().

int PndSdsFE::i

protected

Definition at line 114 of file PndSdsFE.h.

Referenced by CreateInterpolatorList(), and operator=().

ROOT::Math::Interpolator* PndSdsFE::inter

protected

Definition at line 119 of file PndSdsFE.h.

Referenced by GetChargeFromTot(), LoadInterpolatorList(), and operator=().

double PndSdsFE::stepsize

protected

Definition at line 115 of file PndSdsFE.h.

Referenced by CreateInterpolatorList(), and operator=().

The documentation for this class was generated from the following files:

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation