#include <PndSdsFESimple.h>

Inheritance diagram for PndSdsFESimple:

Public Member Functions
	PndSdsFESimple ()

	PndSdsFESimple (double charsingtime, double constcurrent, double threshold, double frequency)

	~PndSdsFESimple ()

void	SetParameter (double chargingtime, double constcurrent, double threshold, double frequency)

	ClassDef (PndSdsFESimple, 1)

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)

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

Detailed Description

Definition at line 25 of file PndSdsFESimple.h.

Constructor & Destructor Documentation

PndSdsFESimple::PndSdsFESimple ( )

Definition at line 15 of file PndSdsFESimple.cxx.

References PndSdsFEAmpModelSimple::Definition(), PndSdsFE::fFrontEndModel, PndSdsFE::fFunction, PndSdsFE::fFunctionRange, PndSdsFE::fThreshold, PndSdsFE::fTimeStep, and PndSdsFE::GetInterpolatorList().

                                {
         fFrontEndModel = new PndSdsFEAmpModelSimple();
         fFunctionRange = 22000;
         fFunction = new TF1("fFunction",fFrontEndModel,&PndSdsFEAmpModelSimple::Definition,0,fFunctionRange,3);
         fFunction->SetNpx(20000);
 
         fFunction->SetParName(0,"chargetime");
         fFunction->SetParName(1,"constantcurrent");
         fFunction->SetParName(2,"charge");
 
         fFunction->SetParameter("chargetime",100);
         fFunction->SetParameter("constcurrent",60);
 
         fTimeStep = 1. / 150 * 1000.; // Dont forget to change the frequency here if you change it above
         fThreshold = 1000;
 
         GetInterpolatorList();
 
 }

PndSdsFESimple::PndSdsFESimple	(	double	charsingtime,
		double	constcurrent,
		double	threshold,
		double	frequency
	)

Definition at line 35 of file PndSdsFESimple.cxx.

References PndSdsFEAmpModelSimple::Definition(), PndSdsFE::fFrontEndModel, PndSdsFE::fFunction, PndSdsFE::fFunctionRange, PndSdsFE::fThreshold, PndSdsFE::fTimeStep, PndSdsFE::GetInterpolatorList(), and threshold.

                                                                                                          {
         fFrontEndModel = new PndSdsFEAmpModelSimple();
         fFunctionRange = 22000;
         fFunction = new TF1("fFunction",fFrontEndModel,&PndSdsFEAmpModelSimple::Definition,0,fFunctionRange,3);
         fFunction->SetNpx(20000);
 
         fFunction->SetParName(0,"chargetime");
         fFunction->SetParName(1,"constantcurrent");
         fFunction->SetParName(2,"charge");
 
         fFunction->SetParameter(0,chargetime);
         fFunction->SetParameter(1,constcurrent);
 
         fTimeStep = 1. / frequency * 1000.;
         fThreshold = threshold;
 
         GetInterpolatorList();
 
 }

PndSdsFESimple::~PndSdsFESimple ( )

Definition at line 67 of file PndSdsFESimple.cxx.

67 {

68 }

Member Function Documentation

void PndSdsFE::CalcTimeOffSet ( )

protectedinherited

Definition at line 100 of file PndSdsFE.cxx.

References PndSdsFE::fRand, PndSdsFE::fTimeOffSet, and PndSdsFE::fTimeStep.

Referenced by PndSdsFE::PndSdsFE().

                              {
         fTimeOffSet = fRand.Uniform(fTimeStep);
 }

PndSdsFESimple::ClassDef	(	PndSdsFESimple	,
		1
	)

void PndSdsFE::CreateInterpolatorList ( )

protectedinherited

Definition at line 104 of file PndSdsFE.cxx.

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

Referenced by PndSdsFE::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 )

protectedinherited

Definition at line 90 of file PndSdsFE.cxx.

References PndSdsFE::fTimeStep.

Referenced by PndSdsFE::GetTimeStamp(), and PndSdsFE::GetTotFromCharge().

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

double PndSdsFE::GetChargeFromTot ( double tot )

virtualinherited

Definition at line 65 of file PndSdsFE.cxx.

References PndSdsFE::inter.

Referenced by PndSdsFE::GetTimeWalkFromTot().

                                            {
         return   inter->Eval(tot);
 }

void PndSdsFE::GetInterpolatorList ( )

protectedinherited

Definition at line 175 of file PndSdsFE.cxx.

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

Referenced by 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 )

virtualinherited

Definition at line 84 of file PndSdsFE.cxx.

References PndSdsFE::fBaselineEpsilon, PndSdsFE::fFunction, and PndSdsFE::fMaximumAmplitude.

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

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

virtualinherited

Definition at line 95 of file PndSdsFE.cxx.

References PndSdsFE::DigitizeTime(), and PndSdsFE::GetTimeWalkFromCharge().

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

virtual double PndSdsFE::GetTimeStep ( )

inlinevirtualinherited

Definition at line 60 of file PndSdsFE.h.

References PndSdsFE::fTimeStep.

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

60 {return fTimeStep;}

PndSdsFE::fTimeStep

double fTimeStep

Definition: PndSdsFE.h:111

double PndSdsFE::GetTimeWalkFromCharge ( double charge )

virtualinherited

Definition at line 77 of file PndSdsFE.cxx.

References PndSdsFE::fFunction, PndSdsFE::fMaximumAmplitude, and PndSdsFE::fThreshold.

Referenced by PndSdsFE::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 )

virtualinherited

Definition at line 69 of file PndSdsFE.cxx.

References PndSdsFE::fFunction, PndSdsFE::fMaximumAmplitude, PndSdsFE::fThreshold, and PndSdsFE::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 )

virtualinherited

Definition at line 26 of file PndSdsFE.cxx.

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

Referenced by PndSdsFE::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 ( )

protectedinherited

Definition at line 198 of file PndSdsFE.cxx.

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

Referenced by PndSdsFE::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);
 }

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

protectedinherited

Definition at line 136 of file PndSdsFE.cxx.

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

Referenced by PndSdsFE::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 PndSdsFESimple::SetParameter	(	double	chargingtime,
		double	constcurrent,
		double	threshold,
		double	frequency
	)

Definition at line 55 of file PndSdsFESimple.cxx.

References PndSdsFE::fFunction, PndSdsFE::fThreshold, PndSdsFE::fTimeStep, and threshold.

 {
   if(fFunction)
   {
     fFunction->SetParameter(0,chargingtime);
     fFunction->SetParameter(1,constcurrent);
   }
   fThreshold = threshold;
   fTimeStep = 1. / frequency * 1000.;
   return;
 }

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

virtualinherited

Definition at line 22 of file PndSdsFE.cxx.

References PndSdsFE::fFunction.

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

virtual void PndSdsFE::SetThreshold ( double threshold )

inlinevirtualinherited

Definition at line 62 of file PndSdsFE.h.

References PndSdsFE::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

protectedinherited

Definition at line 107 of file PndSdsFE.h.

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

std::vector<double> PndSdsFE::fCharge_list

protectedinherited

Definition at line 117 of file PndSdsFE.h.

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

double PndSdsFE::fEventTime

protectedinherited

Definition at line 109 of file PndSdsFE.h.

Referenced by PndSdsFE::operator=().

PndSdsFEAmpModelSimple* PndSdsFE::fFrontEndModel

protectedinherited

Definition at line 99 of file PndSdsFE.h.

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

TF1* PndSdsFE::fFunction

protectedinherited

Definition at line 104 of file PndSdsFE.h.

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

double PndSdsFE::fFunctionRange

protectedinherited

Definition at line 112 of file PndSdsFE.h.

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

double PndSdsFE::fMaximumAmplitude

protectedinherited

Definition at line 105 of file PndSdsFE.h.

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

int PndSdsFE::fNumberOfMaxElectons

protectedinherited

Definition at line 102 of file PndSdsFE.h.

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

int PndSdsFE::fNumberOfSupportPoints

protectedinherited

Definition at line 101 of file PndSdsFE.h.

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

TRandom2 PndSdsFE::fRand

protectedinherited

Definition at line 113 of file PndSdsFE.h.

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

double PndSdsFE::fThreshold

protectedinherited

Definition at line 106 of file PndSdsFE.h.

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

double PndSdsFE::fTimeOffSet

protectedinherited

Definition at line 110 of file PndSdsFE.h.

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

double PndSdsFE::fTimeStep

protectedinherited

Definition at line 111 of file PndSdsFE.h.

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

double PndSdsFE::fToF

protectedinherited

Definition at line 108 of file PndSdsFE.h.

Referenced by PndSdsFE::operator=().

std::vector<double> PndSdsFE::fTot_list

protectedinherited

Definition at line 118 of file PndSdsFE.h.

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

int PndSdsFE::i

protectedinherited

Definition at line 114 of file PndSdsFE.h.

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

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

protectedinherited

Definition at line 119 of file PndSdsFE.h.

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

double PndSdsFE::stepsize

protectedinherited

Definition at line 115 of file PndSdsFE.h.

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation