[u/mrichter/AliRoot.git] / PHOS / AliPHOSPulseGenerator.cxx

/**************************************************************************
 * Copyright(c) 2007, ALICE Experiment at CERN, All rights reserved.      *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

// The class which simulates the pulse shape from the PHOS FEE shaper,
// make sampled amplitudes, digitize them.
// Use case:
//   AliPHOSPulseGenerator *pulse = new AliPHOSPulseGenerator(energy,time);
//   Int_t *adcHG = new Int_t[pulse->GetRawFormatTimeBins()];
//   Int_t *adcLG= new Int_t[pulse->GetRawFormatTimeBins()];
//   pulse->AddNoise(1.);
//   pulse->MakeSamples();
//   pulse->GetSamples(adcHG, adcHG) ; 
//   pulse->Print();
//   pulse->Draw();
//
// Author: Yuri Kharlov, after Yves Schutz and Per Thomas Hille

// --- ROOT system ---

#include <TCanvas.h> 
#include <TF1.h> 
#include <TGraph.h> 
#include <TH1F.h> 
#include <TMath.h> 
#include <TRandom.h>

// --- AliRoot header files ---
#include "AliLog.h"
#include "AliPHOSPulseGenerator.h"

// --- Standard library ---
#include <cmath>
#include <iostream>

using std::cout;
using std::endl; 

ClassImp(AliPHOSPulseGenerator) 

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

//-----------------------------------------------------------------------------
AliPHOSPulseGenerator::AliPHOSPulseGenerator(Double_t a, Double_t t0)
  : 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
  // @param t0 time delay in nanoseconds of signal relative the first sample. 
  // This value should be between 0 and Ts, where Ts is the sample interval

  fDataHG = new Double_t[fkTimeBins];
  fDataLG = new Double_t[fkTimeBins];
  Reset();
}

//-----------------------------------------------------------------------------
AliPHOSPulseGenerator::~AliPHOSPulseGenerator()
{
  // Destructor: delete arrays of samples

  delete [] fDataHG;
  fDataHG=0;
  delete [] fDataLG;
  fDataLG=0;
}

//-----------------------------------------------------------------------------
void AliPHOSPulseGenerator::Reset()
{
  // Reset all sample amplitudes to 0

  for (Int_t i=0; i<fkTimeBins; i++) {
    fDataHG[i] = 0.;
    fDataLG[i] = 0.;
  }
}

//-----------------------------------------------------------------------------
void AliPHOSPulseGenerator::AddBaseline(Double_t baselineLevel)
{
  // Adds a baseline offset to the signal
  // @param baselineLevel The basline level to add
  for (Int_t i=0; i<fkTimeBins; i++) {
    fDataHG[i] += baselineLevel;
    fDataLG[i] += baselineLevel;
  }
  // Digitize floating point amplitudes to integers
  if (fDigitize) Digitize();
}

//-----------------------------------------------------------------------------
void AliPHOSPulseGenerator::AddNoise(Double_t sigma)
{
  // Adds Gaussian uncorrelated to the sample array
  // @param sigma the noise amplitude in entities of ADC levels  
  
  for (Int_t i=0; i<fkTimeBins; i++) {
    fDataHG[i] = gRandom->Gaus(0., sigma) ; 
    fDataLG[i] = gRandom->Gaus(0., sigma) ; 
  }
}

//-----------------------------------------------------------------------------
void AliPHOSPulseGenerator::AddNoise(Double_t * /* sigma */, Double_t /* cutoff */)
{
  //Adds correlated Gaussian noise with cutof frequency "cutoff"
  // @param sigma noise amplitude in entities of ADC levels
  // @param -30DB cutoff frequency of the noise in entities of sampling frequency

  AliError("not implemented yet");
}

//-----------------------------------------------------------------------------
void AliPHOSPulseGenerator::AddPretriggerSamples(Int_t nPresamples)
{
  // Adds pretrigger samples to the sample arrays and replace them
  // with concatinated and truncated arrays

  Double_t *tmpDataHG = new Double_t[fkTimeBins];
  Double_t *tmpDataLG = new Double_t[fkTimeBins];
  Int_t i;
  for (i=0; i<fkTimeBins; i++) {
    tmpDataHG[i] = fDataHG[i];
    tmpDataLG[i] = fDataLG[i];
  }
  for (i=0; i<fkTimeBins; i++) {
    if (i<nPresamples) {
      fDataHG[i] = 0.;
      fDataLG[i] = 0.;
    }
    else {
      fDataHG[i] = tmpDataHG[i-nPresamples];
      fDataLG[i] = tmpDataLG[i-nPresamples];
    }
  }
  delete [] tmpDataHG;
  delete [] tmpDataLG;
}

//-----------------------------------------------------------------------------
void AliPHOSPulseGenerator::Digitize()
{
  // Emulates ADC: rounds up to nearest integer value all amplitudes
  for (Int_t i=0; i<fkTimeBins; i++) {
    fDataHG[i] = (Int_t)(fDataHG[i]);
    fDataLG[i] = (Int_t)(fDataLG[i]);
  }
}

//-----------------------------------------------------------------------------
Double_t AliPHOSPulseGenerator::RawResponseFunction(Double_t *x, Double_t *par) 
{
  // Shape of the electronics raw reponse:
  // It is a semi-gaussian, 2nd order Gamma function of the general form
  // v(t) = A *(t/tp)**n * exp(-n * t/tp-n) with 
  // tp : peaking time  fgTimePeak
  // n  : order of the function
  
  Double_t signal ;
  Double_t xx = x[0] - ( fgTimeTrigger + par[1] ) ; 

  if (xx < 0 || xx > GetRawFormatTimeMax()) 
    signal = 0. ;  
  else {
    signal =  par[0] * TMath::Power(xx/fgTimePeak, fgOrder) * TMath::Exp(-fgOrder*(xx/fgTimePeak-1.)) ; //normalized to par[2] at maximum
  }
  return signal ;  
}

//__________________________________________________________________
Bool_t AliPHOSPulseGenerator::MakeSamples()
{
  // for a start time fTZero and an amplitude fAmplitude given by digit, 
  // calculates the raw sampled response AliPHOSPulseGenerator::RawResponseFunction

  const Int_t kRawSignalOverflow = 0x3FF ; // decimal 1023
  Bool_t lowGain = kFALSE ; 

  TF1 signalF("signal", RawResponseFunction, 0, GetRawFormatTimeMax(), 4);

  for (Int_t iTime = 0; iTime < GetRawFormatTimeBins(); iTime++) {
    signalF.SetParameter(0, fAmplitude) ; 
    signalF.SetParameter(1, fTZero) ; 
    Double_t time = iTime * GetRawFormatTimeMax() / GetRawFormatTimeBins() ;
    Double_t signal = signalF.Eval(time) ;     
    fDataHG[iTime] += signal;
    if ( static_cast<Int_t>(fDataHG[iTime]+0.5) > kRawSignalOverflow ){  // larger than 10 bits 
      fDataHG[iTime] = kRawSignalOverflow ;
      lowGain = kTRUE ; 
    }

    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();
  return lowGain ; 
}

//__________________________________________________________________
void AliPHOSPulseGenerator::GetSamples(Int_t *adcHG, Int_t *adcLG) const
{
  // Return integer sample arrays adcHG and adcLG
  for (Int_t iTime = 0; iTime < GetRawFormatTimeBins(); iTime++) {
    adcHG[iTime] = static_cast<Int_t>(fDataHG[iTime]) ;
    adcLG[iTime] = static_cast<Int_t>(fDataLG[iTime]) ;
  }
}

//__________________________________________________________________
void AliPHOSPulseGenerator::Print(Option_t*) const
{
  // Prints sampled amplitudes to stdout
  Int_t i;
  cout << "High gain: ";
  for (i=0; i<fkTimeBins; i++)
    cout << (Int_t)fDataHG[i] << " ";
  cout << endl;

  cout << "Low  gain: ";
  for (i=0; i<fkTimeBins; i++)
    cout << (Int_t)fDataLG[i] << " ";
  cout << endl;
}

//__________________________________________________________________
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++) {
    time[iTime] = iTime * GetRawFormatTimeMax() / GetRawFormatTimeBins() ;
  }
  Int_t nPoints = fkTimeBins;
  TGraph *graphHG = new TGraph(nPoints,time,fDataHG);
  TGraph *graphLG = new TGraph(nPoints,time,fDataLG);
  graphHG->SetMarkerStyle(20);
  graphLG->SetMarkerStyle(20);
  graphHG->SetMarkerSize(0.4);
  graphLG->SetMarkerSize(0.4);
  graphHG->SetTitle("High gain samples");
  graphLG->SetTitle("Low gain samples");

  TCanvas *c1 = new TCanvas("c1","Raw ALTRO samples",10,10,700,500);
  c1->SetFillColor(0);

  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();
}
Commit	Line	Data
431a9211	1	/**************************************************************************
	2	* Copyright(c) 2007, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/* $Id$ */
	17
	18	// The class which simulates the pulse shape from the PHOS FEE shaper,
	19	// make sampled amplitudes, digitize them.
	20	// Use case:
	21	// AliPHOSPulseGenerator *pulse = new AliPHOSPulseGenerator(energy,time);
	22	// Int_t *adcHG = new Int_t[pulse->GetRawFormatTimeBins()];
	23	// Int_t *adcLG= new Int_t[pulse->GetRawFormatTimeBins()];
2111ab30	24	// pulse->AddNoise(1.);
431a9211	25	// pulse->MakeSamples();
	26	// pulse->GetSamples(adcHG, adcHG) ;
	27	// pulse->Print();
	28	// pulse->Draw();
	29	//
	30	// Author: Yuri Kharlov, after Yves Schutz and Per Thomas Hille
	31
	32	// --- ROOT system ---
a1e17193	33
	34	#include <TCanvas.h>
	35	#include <TF1.h>
	36	#include <TGraph.h>
	37	#include <TH1F.h>
	38	#include <TMath.h>
	39	#include <TRandom.h>
431a9211	40
	41	// --- AliRoot header files ---
	42	#include "AliLog.h"
	43	#include "AliPHOSPulseGenerator.h"
	44
	45	// --- Standard library ---
	46	#include <cmath>
	47	#include <iostream>
	48
	49	using std::cout;
	50	using std::endl;
	51
	52	ClassImp(AliPHOSPulseGenerator)
	53
431a9211	54	Int_t AliPHOSPulseGenerator::fgOrder = 2 ; // order of the Gamma function
ff857cdc	55	Double_t AliPHOSPulseGenerator::fgTimePeak = 2.1E-6 ; // tau=2.1 micro seconds
ff857cdc	56	Double_t AliPHOSPulseGenerator::fgTimeTrigger = 100E-9 ; // one tick 100 ns
431a9211	57
	58	//-----------------------------------------------------------------------------
	59	AliPHOSPulseGenerator::AliPHOSPulseGenerator(Double_t a, Double_t t0)
b28a9811	60	: TObject(), fAmplitude(a), fTZero(t0), fHG2LGratio(16.), fDataHG(0), fDataLG(0), fDigitize(kTRUE)
431a9211	61	{
	62	// Contruct a pulsegenrator object and initializes all necessary parameters
	63	// @param a digit amplitude in GeV
	64	// @param t0 time delay in nanoseconds of signal relative the first sample.
	65	// This value should be between 0 and Ts, where Ts is the sample interval
	66
	67	fDataHG = new Double_t[fkTimeBins];
	68	fDataLG = new Double_t[fkTimeBins];
2111ab30	69	Reset();
431a9211	70	}
	71
	72	//-----------------------------------------------------------------------------
431a9211	73	AliPHOSPulseGenerator::~AliPHOSPulseGenerator()
	74	{
	75	// Destructor: delete arrays of samples
	76
	77	delete [] fDataHG;
	78	fDataHG=0;
	79	delete [] fDataLG;
	80	fDataLG=0;
	81	}
	82
2111ab30	83	//-----------------------------------------------------------------------------
	84	void AliPHOSPulseGenerator::Reset()
	85	{
	86	// Reset all sample amplitudes to 0
	87
	88	for (Int_t i=0; i<fkTimeBins; i++) {
	89	fDataHG[i] = 0.;
	90	fDataLG[i] = 0.;
	91	}
	92	}
	93
431a9211	94	//-----------------------------------------------------------------------------
	95	void AliPHOSPulseGenerator::AddBaseline(Double_t baselineLevel)
	96	{
	97	// Adds a baseline offset to the signal
	98	// @param baselineLevel The basline level to add
	99	for (Int_t i=0; i<fkTimeBins; i++) {
	100	fDataHG[i] += baselineLevel;
	101	fDataLG[i] += baselineLevel;
	102	}
2111ab30	103	// Digitize floating point amplitudes to integers
2111ab30	104	if (fDigitize) Digitize();
431a9211	105	}
	106
	107	//-----------------------------------------------------------------------------
2111ab30	108	void AliPHOSPulseGenerator::AddNoise(Double_t sigma)
431a9211	109	{
2111ab30	110	// Adds Gaussian uncorrelated to the sample array
431a9211	111	// @param sigma the noise amplitude in entities of ADC levels
431a9211	112
2111ab30	113	for (Int_t i=0; i<fkTimeBins; i++) {
	114	fDataHG[i] = gRandom->Gaus(0., sigma) ;
	115	fDataLG[i] = gRandom->Gaus(0., sigma) ;
	116	}
431a9211	117	}
431a9211	118
431a9211	119	//-----------------------------------------------------------------------------
	120	void AliPHOSPulseGenerator::AddNoise(Double_t * /* sigma /, Double_t / cutoff */)
	121	{
	122	//Adds correlated Gaussian noise with cutof frequency "cutoff"
	123	// @param sigma noise amplitude in entities of ADC levels
	124	// @param -30DB cutoff frequency of the noise in entities of sampling frequency
	125
	126	AliError("not implemented yet");
	127	}
	128
	129	//-----------------------------------------------------------------------------
	130	void AliPHOSPulseGenerator::AddPretriggerSamples(Int_t nPresamples)
	131	{
	132	// Adds pretrigger samples to the sample arrays and replace them
	133	// with concatinated and truncated arrays
	134
	135	Double_t *tmpDataHG = new Double_t[fkTimeBins];
	136	Double_t *tmpDataLG = new Double_t[fkTimeBins];
	137	Int_t i;
	138	for (i=0; i<fkTimeBins; i++) {
	139	tmpDataHG[i] = fDataHG[i];
	140	tmpDataLG[i] = fDataLG[i];
	141	}
	142	for (i=0; i<fkTimeBins; i++) {
	143	if (i<nPresamples) {
2111ab30	144	fDataHG[i] = 0.;
2111ab30	145	fDataLG[i] = 0.;
431a9211	146	}
431a9211	147	else {
2111ab30	148	fDataHG[i] = tmpDataHG[i-nPresamples];
2111ab30	149	fDataLG[i] = tmpDataLG[i-nPresamples];
431a9211	150	}
	151	}
	152	delete [] tmpDataHG;
	153	delete [] tmpDataLG;
	154	}
	155
	156	//-----------------------------------------------------------------------------
	157	void AliPHOSPulseGenerator::Digitize()
	158	{
6f47f50d	159	// Emulates ADC: rounds up to nearest integer value all amplitudes
431a9211	160	for (Int_t i=0; i<fkTimeBins; i++) {
54ac223e	161	fDataHG[i] = (Int_t)(fDataHG[i]);
54ac223e	162	fDataLG[i] = (Int_t)(fDataLG[i]);
431a9211	163	}
	164	}
	165
	166	//-----------------------------------------------------------------------------
	167	Double_t AliPHOSPulseGenerator::RawResponseFunction(Double_t x, Double_t par)
	168	{
	169	// Shape of the electronics raw reponse:
	170	// It is a semi-gaussian, 2nd order Gamma function of the general form
b28a9811	171	// v(t) = A (t/tp)n exp(-n * t/tp-n) with
b28a9811	172	// tp : peaking time fgTimePeak
431a9211	173	// n : order of the function
431a9211	174
431a9211	175	Double_t signal ;
b28a9811	176	Double_t xx = x[0] - ( fgTimeTrigger + par[1] ) ;
431a9211	177
ff857cdc	178	if (xx < 0 \|\| xx > GetRawFormatTimeMax())
431a9211	179	signal = 0. ;
431a9211	180	else {
b28a9811	181	signal = par[0] * TMath::Power(xx/fgTimePeak, fgOrder) * TMath::Exp(-fgOrder*(xx/fgTimePeak-1.)) ; //normalized to par[2] at maximum
431a9211	182	}
	183	return signal ;
	184	}
	185
431a9211	186	//__________________________________________________________________
	187	Bool_t AliPHOSPulseGenerator::MakeSamples()
	188	{
	189	// for a start time fTZero and an amplitude fAmplitude given by digit,
	190	// calculates the raw sampled response AliPHOSPulseGenerator::RawResponseFunction
	191
2111ab30	192	const Int_t kRawSignalOverflow = 0x3FF ; // decimal 1023
431a9211	193	Bool_t lowGain = kFALSE ;
	194
	195	TF1 signalF("signal", RawResponseFunction, 0, GetRawFormatTimeMax(), 4);
	196
	197	for (Int_t iTime = 0; iTime < GetRawFormatTimeBins(); iTime++) {
b28a9811	198	signalF.SetParameter(0, fAmplitude) ;
b28a9811	199	signalF.SetParameter(1, fTZero) ;
431a9211	200	Double_t time = iTime * GetRawFormatTimeMax() / GetRawFormatTimeBins() ;
431a9211	201	Double_t signal = signalF.Eval(time) ;
2111ab30	202	fDataHG[iTime] += signal;
	203	if ( static_cast<Int_t>(fDataHG[iTime]+0.5) > kRawSignalOverflow ){ // larger than 10 bits
	204	fDataHG[iTime] = kRawSignalOverflow ;
431a9211	205	lowGain = kTRUE ;
431a9211	206	}
431a9211	207
b28a9811	208	Double_t aLGamp = fAmplitude/fHG2LGratio ;
b28a9811	209	signalF.SetParameter(0, aLGamp) ;
431a9211	210	signal = signalF.Eval(time) ;
2111ab30	211	fDataLG[iTime] += signal;
	212	if ( static_cast<Int_t>(fDataLG[iTime]+0.5) > kRawSignalOverflow) // larger than 10 bits
	213	fDataLG[iTime] = kRawSignalOverflow ;
431a9211	214	}
2111ab30	215	// Digitize floating point amplitudes to integers
2111ab30	216	if (fDigitize) Digitize();
431a9211	217	return lowGain ;
	218	}
	219
	220	//__________________________________________________________________
	221	void AliPHOSPulseGenerator::GetSamples(Int_t adcHG, Int_t adcLG) const
	222	{
	223	// Return integer sample arrays adcHG and adcLG
	224	for (Int_t iTime = 0; iTime < GetRawFormatTimeBins(); iTime++) {
	225	adcHG[iTime] = static_cast<Int_t>(fDataHG[iTime]) ;
	226	adcLG[iTime] = static_cast<Int_t>(fDataLG[iTime]) ;
	227	}
	228	}
	229
	230	//__________________________________________________________________
	231	void AliPHOSPulseGenerator::Print(Option_t*) const
	232	{
	233	// Prints sampled amplitudes to stdout
	234	Int_t i;
	235	cout << "High gain: ";
	236	for (i=0; i<fkTimeBins; i++)
	237	cout << (Int_t)fDataHG[i] << " ";
	238	cout << endl;
	239
	240	cout << "Low gain: ";
	241	for (i=0; i<fkTimeBins; i++)
	242	cout << (Int_t)fDataLG[i] << " ";
	243	cout << endl;
	244	}
	245
	246	//__________________________________________________________________
565550c6	247	void AliPHOSPulseGenerator::Draw(Option_t* opt)
431a9211	248	{
431a9211	249	// Draw graphs with high and low gain samples
565550c6	250	// Option_t* opt="all": draw both HG and LG in one canvas
	251	// "HG" : draw HG only
	252	// "LG" : draw LG only
431a9211	253
	254	Double_t *time = new Double_t[fkTimeBins];
	255	for (Int_t iTime = 0; iTime < GetRawFormatTimeBins(); iTime++) {
	256	time[iTime] = iTime * GetRawFormatTimeMax() / GetRawFormatTimeBins() ;
	257	}
	258	Int_t nPoints = fkTimeBins;
	259	TGraph *graphHG = new TGraph(nPoints,time,fDataHG);
	260	TGraph *graphLG = new TGraph(nPoints,time,fDataLG);
	261	graphHG->SetMarkerStyle(20);
2111ab30	262	graphLG->SetMarkerStyle(20);
	263	graphHG->SetMarkerSize(0.4);
	264	graphLG->SetMarkerSize(0.4);
431a9211	265	graphHG->SetTitle("High gain samples");
	266	graphLG->SetTitle("Low gain samples");
	267
	268	TCanvas *c1 = new TCanvas("c1","Raw ALTRO samples",10,10,700,500);
2111ab30	269	c1->SetFillColor(0);
565550c6	270
	271	if (strstr(opt,"all")){
	272	c1->Divide(2,1);
	273	c1->cd(1);
	274	gPad->SetLeftMargin(0.15);
	275	}
	276	if (strstr(opt,"LG") == 0){
	277	graphHG->Draw("AP");
37f77f86	278	graphHG->GetHistogram()->SetTitleOffset(1.0,"Y");
37f77f86	279	graphHG->GetHistogram()->SetXTitle("time, sec");
565550c6	280	graphHG->GetHistogram()->SetYTitle("Amplitude, ADC counts");
	281	}
	282	if (strstr(opt,"all")){
	283	c1->cd(2);
	284	gPad->SetLeftMargin(0.15);
	285	}
	286	if (strstr(opt,"HG") == 0){
	287	graphLG->Draw("AP");
37f77f86	288	graphLG->GetHistogram()->SetTitleOffset(1.0,"Y");
37f77f86	289	graphLG->GetHistogram()->SetXTitle("time, sec");
565550c6	290	graphLG->GetHistogram()->SetYTitle("Amplitude, ADC counts");
565550c6	291	}
431a9211	292	c1->Update();
431a9211	293	}
7ceef380	294