ClassImp(AliPHOSPulseGenerator)
-Double_t AliPHOSPulseGenerator::fgCapa = 1.1; // 1pF
Int_t AliPHOSPulseGenerator::fgOrder = 2 ; // order of the Gamma function
Double_t AliPHOSPulseGenerator::fgTimePeak = 2.1E-6 ; // tau=2.1 micro seconds
Double_t AliPHOSPulseGenerator::fgTimeTrigger = 100E-9 ; // one tick 100 ns
-Double_t AliPHOSPulseGenerator::fgHighCharge = 8.8; // adjusted for a high gain range of 5 GeV (10 bits)
-Double_t AliPHOSPulseGenerator::fgHighGain = 6.85;
-Double_t AliPHOSPulseGenerator::fgHighLowGainFactor = 16.; // adjusted for a low gain range of 80 GeV (10 bits)
//-----------------------------------------------------------------------------
AliPHOSPulseGenerator::AliPHOSPulseGenerator(Double_t a, Double_t t0)
- : TObject(), fAmplitude(a), fTZero(t0), fDataHG(0), fDataLG(0), fDigitize(kTRUE)
+ : TObject(), fAmplitude(a), fTZero(t0), fHG2LGratio(16.), fDataHG(0), fDataLG(0), fDigitize(kTRUE)
{
// Contruct a pulsegenrator object and initializes all necessary parameters
// @param a digit amplitude in GeV
//-----------------------------------------------------------------------------
AliPHOSPulseGenerator::AliPHOSPulseGenerator(const AliPHOSPulseGenerator & pulse)
- : TObject(), fAmplitude(pulse.fAmplitude), fTZero(pulse.fTZero), fDataHG(0), fDataLG(0), fDigitize(kTRUE)
+ : TObject(), fAmplitude(pulse.fAmplitude), fTZero(pulse.fTZero),fHG2LGratio(pulse.fHG2LGratio),
+ fDataHG(0), fDataLG(0), fDigitize(kTRUE)
{
fDataHG = new Double_t[pulse.fkTimeBins];
fDataLG = new Double_t[pulse.fkTimeBins];
//-----------------------------------------------------------------------------
void AliPHOSPulseGenerator::Digitize()
{
- // Emulates ADC: rounds down to nearest integer value all amplitudes
+ // Emulates ADC: rounds up to nearest integer value all amplitudes
for (Int_t i=0; i<fkTimeBins; i++) {
- fDataHG[i] = (Double_t) ((Int_t)(fDataHG[i] + 0.5));
- fDataLG[i] = (Double_t) ((Int_t)(fDataLG[i] + 0.5));
+ fDataHG[i] = (Int_t)(fDataHG[i]);
+ fDataLG[i] = (Int_t)(fDataLG[i]);
}
}
{
// Shape of the electronics raw reponse:
// It is a semi-gaussian, 2nd order Gamma function of the general form
- // v(t) = n**n * Q * A**n / C *(t/tp)**n * exp(-n * t/tp) with
- // tp : peaking time par[0]
+ // v(t) = A *(t/tp)**n * exp(-n * t/tp-n) with
+ // tp : peaking time fgTimePeak
// n : order of the function
- // C : integrating capacitor in the preamplifier
- // A : open loop gain of the preamplifier
- // Q : the total APD charge to be measured Q = C * energy
Double_t signal ;
- Double_t xx = x[0] - ( fgTimeTrigger + par[3] ) ;
+ Double_t xx = x[0] - ( fgTimeTrigger + par[1] ) ;
if (xx < 0 || xx > GetRawFormatTimeMax())
signal = 0. ;
else {
- Double_t fac = par[0] * TMath::Power(fgOrder, fgOrder) * TMath::Power(par[1], fgOrder)/fgCapa ;
- signal = fac * par[2] * TMath::Power(xx/fgTimePeak, fgOrder) * TMath::Exp(-fgOrder*(xx/fgTimePeak)) ;
+ signal = par[0] * TMath::Power(xx/fgTimePeak, fgOrder) * TMath::Exp(-fgOrder*(xx/fgTimePeak-1.)) ; //normalized to par[2] at maximum
}
return signal ;
}
-//__________________________________________________________________
-Double_t AliPHOSPulseGenerator::RawResponseFunctionMax(Double_t charge, Double_t gain)
-{
- // Maximum value of the shaper response function
- return ( charge * TMath::Power(fgOrder, fgOrder) * TMath::Power(gain, fgOrder)
- / ( fgCapa * TMath::Exp(fgOrder) ) );
-}
-
//__________________________________________________________________
Bool_t AliPHOSPulseGenerator::MakeSamples()
{
TF1 signalF("signal", RawResponseFunction, 0, GetRawFormatTimeMax(), 4);
for (Int_t iTime = 0; iTime < GetRawFormatTimeBins(); iTime++) {
- signalF.SetParameter(0, fgHighCharge) ;
- signalF.SetParameter(1, fgHighGain) ;
- signalF.SetParameter(2, fAmplitude) ;
- signalF.SetParameter(3, fTZero) ;
+ signalF.SetParameter(0, fAmplitude) ;
+ signalF.SetParameter(1, fTZero) ;
Double_t time = iTime * GetRawFormatTimeMax() / GetRawFormatTimeBins() ;
Double_t signal = signalF.Eval(time) ;
fDataHG[iTime] += signal;
lowGain = kTRUE ;
}
- signalF.SetParameter(0, GetRawFormatLowCharge() ) ;
- signalF.SetParameter(1, GetRawFormatLowGain() ) ;
+ Double_t aLGamp = fAmplitude/fHG2LGratio ;
+ signalF.SetParameter(0, aLGamp) ;
signal = signalF.Eval(time) ;
fDataLG[iTime] += signal;
if ( static_cast<Int_t>(fDataLG[iTime]+0.5) > kRawSignalOverflow) // larger than 10 bits
fDataLG[iTime] = kRawSignalOverflow ;
-
}
// Digitize floating point amplitudes to integers
if (fDigitize) Digitize();
}
//__________________________________________________________________
-void AliPHOSPulseGenerator::Draw(Option_t*)
+void AliPHOSPulseGenerator::Draw(Option_t* opt)
{
// Draw graphs with high and low gain samples
+ // Option_t* opt="all": draw both HG and LG in one canvas
+ // "HG" : draw HG only
+ // "LG" : draw LG only
Double_t *time = new Double_t[fkTimeBins];
for (Int_t iTime = 0; iTime < GetRawFormatTimeBins(); iTime++) {
TCanvas *c1 = new TCanvas("c1","Raw ALTRO samples",10,10,700,500);
c1->SetFillColor(0);
- c1->Divide(2,1);
- c1->cd(1);
- gPad->SetLeftMargin(0.15);
- graphHG->Draw("AP");
- graphHG->GetHistogram()->SetTitleOffset(1.6,"Y");
- graphHG->GetHistogram()->SetXTitle("time, #musec");
- graphHG->GetHistogram()->SetYTitle("Amplitude, ADC counts");
- c1->cd(2);
- gPad->SetLeftMargin(0.15);
- graphLG->Draw("AP");
- graphLG->GetHistogram()->SetTitleOffset(1.6,"Y");
- graphLG->GetHistogram()->SetXTitle("time, #musec");
- graphLG->GetHistogram()->SetYTitle("Amplitude, ADC counts");
+
+ if (strstr(opt,"all")){
+ c1->Divide(2,1);
+ c1->cd(1);
+ gPad->SetLeftMargin(0.15);
+ }
+ if (strstr(opt,"LG") == 0){
+ graphHG->Draw("AP");
+ graphHG->GetHistogram()->SetTitleOffset(1.0,"Y");
+ graphHG->GetHistogram()->SetXTitle("time, sec");
+ graphHG->GetHistogram()->SetYTitle("Amplitude, ADC counts");
+ }
+ if (strstr(opt,"all")){
+ c1->cd(2);
+ gPad->SetLeftMargin(0.15);
+ }
+ if (strstr(opt,"HG") == 0){
+ graphLG->Draw("AP");
+ graphLG->GetHistogram()->SetTitleOffset(1.0,"Y");
+ graphLG->GetHistogram()->SetXTitle("time, sec");
+ graphLG->GetHistogram()->SetYTitle("Amplitude, ADC counts");
+ }
c1->Update();
}
-//__________________________________________________________________
-Double_t AliPHOSPulseGenerator::GeV2ADC()
-{
- //Return GeV to ADC counts conversion factor.
- //adc_counts = energy[GeV]*AliPHOSPulseGenerator::GeV2ADC().
-
- return RawResponseFunctionMax(fgHighCharge,fgHighGain);
-}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff