[u/mrichter/AliRoot.git] / EMCAL / AliCaloRawAnalyzerComparison.cxx

/**************************************************************************
 * This file is property of and copyright by the Experimental Nuclear     *
 * Physics Group, Dep. of Physics                                         *
 * University of Oslo, Norway, 2007                                       *
 *                                                                        *
 * Author: Per Thomas Hille <perthi@fys.uio.no> for the ALICE HLT Project.*
 * Contributors are mentioned in the code where appropriate.              *
 * Please report bugs to perthi@fys.uio.no                                *
 *                                                                        *
 * 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.                  *
 **************************************************************************/
#include "AliCaloRawAnalyzerComparison.h"

#include "AliCaloRawAnalyzerCrude.h"
#include "AliCaloRawAnalyzerLMS.h"
#include "AliCaloRawAnalyzerFastFit.h"
#include "AliCaloRawAnalyzerNN.h"
#include "AliCaloRawAnalyzerPeakFinder.h"

#include "AliCaloFitResults.h"
#include "TH2D.h"
#include "TMath.h"
#include "TFile.h"

#include <iostream>

using namespace std;

AliCaloRawAnalyzerComparison::AliCaloRawAnalyzerComparison() : fMod(0),
							       fMonCol1(1),
							       fMonRow1(15),
							       fMonCol2(1),
							       fMonRow2(16)
{ // ctor
  fReferenceAnalyzer = new  AliCaloRawAnalyzerLMS();

  AliCaloRawAnalyzerCrude *crude = new AliCaloRawAnalyzerCrude();
  AliCaloRawAnalyzerLMS *lms = new AliCaloRawAnalyzerLMS();
  AliCaloRawAnalyzerFastFit *fastfit = new AliCaloRawAnalyzerFastFit();
  AliCaloRawAnalyzerNN *neuralnet = new AliCaloRawAnalyzerNN();
  AliCaloRawAnalyzerPeakFinder *peakfinder = new AliCaloRawAnalyzerPeakFinder();

  fReferenceAnalyzer->SetIsZeroSuppressed();
  peakfinder->SetIsZeroSuppressed();
  lms->SetIsZeroSuppressed();
  fastfit->SetIsZeroSuppressed();
  neuralnet->SetIsZeroSuppressed();
  crude->SetIsZeroSuppressed(); 
  
  fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)crude);
  fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)lms);
  fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)fastfit);
  fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)neuralnet);
  
  fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)peakfinder);
  InitHistograms( fRawAnalyzers, fReferenceAnalyzer );
}


AliCaloRawAnalyzerComparison::~AliCaloRawAnalyzerComparison()
{ // dtor

}
 

void 
AliCaloRawAnalyzerComparison::Evaluate( const vector<AliCaloBunchInfo> &bunchvector, const UInt_t altrocfg1,  
					const UInt_t altrocfg2, const int event, const int col, const int row )
{ // evaluate methods
  //  cout << __FILE__ << __LINE__ << endl;

  AliCaloFitResults ref = fReferenceAnalyzer->Evaluate( bunchvector, altrocfg1, altrocfg2 );
  //  AliCaloFitResults ref;

  for(int i=0; i < fRawAnalyzers.size(); i++ )
    {
      AliCaloFitResults an = fRawAnalyzers.at(i)->Evaluate( bunchvector, altrocfg1,  altrocfg2 ); 
      
      //     cout << __FILE__ << __LINE__ << ":" << fRawAnalyzers.at(i)->GetAlgoAbbr() << "  col = " << col << " row " << row << "  amp= " <<  an.GetAmp() << endl ;

      //   fAmplitudeVsEvent[i]->Fill( event, an.GetAmp() );
      //   fTofVsEvent[i]->Fill( event, an.GetTof() );
    
      fRefAmpVsAnalyzers[i]->Fill(ref.GetAmp(), an.GetAmp() );
      fRefTofVsAnalyzers[i]->Fill(ref.GetTof(), an.GetTof() ); 
      fAmpDiff[i]->Fill(ref.GetAmp() - an.GetAmp() ); 
      fTofDiff[i]->Fill(ref.GetTof() - an.GetTof() );  
     
      if(col  == fMonCol1 && fMonRow1 == row )
	{
	  fAmplitudeVsEvent[i]->Fill( event, an.GetAmp() );
	  fTofVsEvent[i]->Fill( event, an.GetTof() );
	  fMon1[i] = an;
	}
      
      if(col  == fMonCol2 &&  row == fMonRow2  )
	{
	  fMon2[i] = an;
	}
 
      if(row >= 0 && col >= 0 )
	{
	  fAmpHistograms[i][col][row]->Fill( an.GetAmp() );
	}
    }
}  


void 
AliCaloRawAnalyzerComparison::EventChanged()
{ // new event
  for(int i=0; i < fRawAnalyzers.size(); i++ )
    {
      /*
      if( ( fMon1[i].GetAmp() > 50 &&   fMon2[i].GetAmp() > 50 ) &&  (  fMon1[i].GetAmp() < 1023  &&   fMon2[i].GetAmp() < 1023 )  )
	{
	  cout << __FILE__ << __LINE__ << fRawAnalyzers.at(i)->GetAlgoAbbr() << "\tamp1=" <<  fMon1[i].GetAmp()  << "\tamp2 = " <<  fMon2[i].GetAmp() << endl;
	}
      */

      if( fMon1[i].GetAmp()  > 50  &&  fMon2[i].GetAmp() > 50  )
	{
	  //  if(  fMon1[i].GetTof() != fMon2[i].GetTof())
	  {
	    //    fTofResDifferential[i]->Fill( ( fMon1[i].GetTof() - fMon2[i].GetTof())/TMath::Sqrt(2) );
	    fTofResDifferential[i]->Fill( ( fMon1[i].GetTof() - fMon2[i].GetTof()) );
	    fTofResAbsolute[i]->Fill(  fMon1[i].GetTof() ); 
	  }
	}
    }
}


void 
AliCaloRawAnalyzerComparison::WriteHistograms()
{ // write histograms 
  TFile *f = new TFile("comparison2.root", "recreate");
  
  /*
  for(int col=0; col < NZCOLSSMOD; col ++ )
    {
      for(int row=0; row < NXROWSSMOD; row ++ )
	{
	  fAmpHistograms[col][row]->Write();
	}
    }
  */  

  for(int i=0; i < fRawAnalyzers.size(); i++ )
    {
      for(int col=0; col < NZCOLSSMOD; col ++ )
	{
	  for(int row=0; row < NXROWSSMOD; row ++ )
	    {
	      fAmpHistograms[i][col][row]->Write();
	    }
	}
      
      fAmplitudeVsEvent[i]->Write();
      fTofVsEvent[i]->Write();
      fRefAmpVsAnalyzers[i]->Write();
      fRefTofVsAnalyzers[i]->Write(); 
      fAmpDiff[i]->Write(); 
      fTofDiff[i]->Write();  
      fTofResDifferential[i]->Write();
      fTofResAbsolute[i]->Write(); 

   }
 
  f->Close();
}


void
AliCaloRawAnalyzerComparison::InitHistograms( vector <AliCaloRawAnalyzer*> analyzers, AliCaloRawAnalyzer* ref )
{ // init histograms 
  char tmpname[256];
 
  /* 
  for(int col=0; col < NZCOLSSMOD; col ++ )
    {
      for(int row=0; row < NXROWSSMOD; row ++ )
	{
	  sprintf(tmpname, "z(col)%d_x(row)%d_amplitude;Counts;Amplitude/ADC counts", col, row);
	  fAmpHistograms[col][row] = new TH1D(tmpname, tmpname, 1024, 0, 1023 );
	}
    }
  */

  // TH1D *fAmpHistograms[NZCOLSSMOD][NXROWSSMOD];

  for(int i=0; i < analyzers.size(); i++ )
    {
      for(int col=0; col < NZCOLSSMOD; col ++ )
	{
	  for(int row=0; row < NXROWSSMOD; row ++ )
	    {
	      sprintf(tmpname, "z(col)%d_x(row)%d_amplitude_%s;Amplitude/ADC counts;Counts", col, row,  analyzers.at(i)->GetAlgoAbbr() );
	      fAmpHistograms[i][col][row] = new TH1D(tmpname, tmpname, 1024, 0, 1023 );
	    }
	}
      
      sprintf(tmpname, "%s_amplitude_vs_event_row%d_col%d;Event;Amplitude/ADC counts", analyzers.at(i)->GetAlgoAbbr(), fMonRow1, fMonCol1 );
      fAmplitudeVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 1024, 0, 1023 ); 
     
      sprintf(tmpname, "%s_tof_vs_event_row%d_col%d;Event;Amplitude/ADC counts", analyzers.at(i)->GetAlgoAbbr(), fMonRow1, fMonCol1 );  
      
      //     fTofVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 2000, 3000, 5000); 
      fTofVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 2000, 1000, 4000);  


      sprintf(tmpname, "%s_vs_%s_ampltude; Amplitude_{%s}/ADC counts; Amplitude_{%s}/ADC counts", ref->GetAlgoAbbr(), 
	      analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr() );
      fRefAmpVsAnalyzers[i] = new TH2D(tmpname, tmpname, 1024, 0, 1023, 1024, 0, 1023); 
      sprintf(tmpname, "%s_vs_%s_tof; tof_{%s}/ns; tof_{%s}/ns", ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr() );
      fRefTofVsAnalyzers[i] = new TH2D(tmpname, tmpname, 500, 2000, 4999, 500, 2000, 4999 ); 
      sprintf( tmpname, "%s-%s amplitude;counts;A_{%s} - A_{%s}",   ref->GetAlgoAbbr(), 
	       analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr());
      fAmpDiff[i] = new TH1D(tmpname, tmpname, 100, -10, 10 );
      sprintf( tmpname, "%s-%s tof;counts;A_{%s} - A_{%s}",   ref->GetAlgoAbbr(), 
	       analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr());
      fTofDiff[i] = new TH1D(tmpname, tmpname, 1000, -5000, 5000 );
      sprintf( tmpname, "%s Differential tof resolution (%d, %d) vs (%d, %d);#sigma_{tof}^{%s}/ns;Counts",  analyzers.at(i)->GetAlgoAbbr(), fMonCol1, fMonRow1, fMonCol2, fMonRow2, 
	       analyzers.at(i)->GetAlgoAbbr() );
      fTofResDifferential[i] = new TH1D(tmpname, tmpname, 1000, -250, 250 );
      sprintf( tmpname, "%s Absolute tof distribution (%d, %d);#sigma_{tof}^{%s}/ns;Counts",  analyzers.at(i)->GetAlgoAbbr(), fMonCol1, fMonRow1, analyzers.at(i)->GetAlgoAbbr() );
      
      fTofResAbsolute[i] = new TH1D(tmpname, tmpname, 2000 , -1000, 7000 );


    }
}
Commit	Line	Data
48a2e3eb	1	/**************************************************************************
	2	* This file is property of and copyright by the Experimental Nuclear *
	3	* Physics Group, Dep. of Physics *
	4	* University of Oslo, Norway, 2007 *
	5	* *
	6	* Author: Per Thomas Hille <perthi@fys.uio.no> for the ALICE HLT Project.*
	7	* Contributors are mentioned in the code where appropriate. *
	8	* Please report bugs to perthi@fys.uio.no *
	9	* *
	10	* Permission to use, copy, modify and distribute this software and its *
	11	* documentation strictly for non-commercial purposes is hereby granted *
	12	* without fee, provided that the above copyright notice appears in all *
	13	* copies and that both the copyright notice and this permission notice *
	14	* appear in the supporting documentation. The authors make no claims *
	15	* about the suitability of this software for any purpose. It is *
	16	* provided "as is" without express or implied warranty. *
	17	**************************************************************************/
	18	#include "AliCaloRawAnalyzerComparison.h"
	19
	20	#include "AliCaloRawAnalyzerCrude.h"
	21	#include "AliCaloRawAnalyzerLMS.h"
	22	#include "AliCaloRawAnalyzerFastFit.h"
	23	#include "AliCaloRawAnalyzerNN.h"
	24	#include "AliCaloRawAnalyzerPeakFinder.h"
	25
	26	#include "AliCaloFitResults.h"
	27	#include "TH2D.h"
	28	#include "TMath.h"
	29	#include "TFile.h"
	30
	31	#include <iostream>
	32
	33	using namespace std;
	34
	35	AliCaloRawAnalyzerComparison::AliCaloRawAnalyzerComparison() : fMod(0),
	36	fMonCol1(1),
	37	fMonRow1(15),
	38	fMonCol2(1),
	39	fMonRow2(16)
85a699d8	40	{ // ctor
48a2e3eb	41	fReferenceAnalyzer = new AliCaloRawAnalyzerLMS();
	42
	43	AliCaloRawAnalyzerCrude *crude = new AliCaloRawAnalyzerCrude();
	44	AliCaloRawAnalyzerLMS *lms = new AliCaloRawAnalyzerLMS();
	45	AliCaloRawAnalyzerFastFit *fastfit = new AliCaloRawAnalyzerFastFit();
	46	AliCaloRawAnalyzerNN *neuralnet = new AliCaloRawAnalyzerNN();
	47	AliCaloRawAnalyzerPeakFinder *peakfinder = new AliCaloRawAnalyzerPeakFinder();
	48
	49	fReferenceAnalyzer->SetIsZeroSuppressed();
	50	peakfinder->SetIsZeroSuppressed();
	51	lms->SetIsZeroSuppressed();
	52	fastfit->SetIsZeroSuppressed();
	53	neuralnet->SetIsZeroSuppressed();
	54	crude->SetIsZeroSuppressed();
	55
	56	fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)crude);
	57	fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)lms);
	58	fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)fastfit);
	59	fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)neuralnet);
	60
	61	fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)peakfinder);
85a699d8	62	InitHistograms( fRawAnalyzers, fReferenceAnalyzer );
48a2e3eb	63	}
	64
	65
	66	AliCaloRawAnalyzerComparison::~AliCaloRawAnalyzerComparison()
85a699d8	67	{ // dtor
48a2e3eb	68
	69	}
	70
	71
	72	void
	73	AliCaloRawAnalyzerComparison::Evaluate( const vector<AliCaloBunchInfo> &bunchvector, const UInt_t altrocfg1,
	74	const UInt_t altrocfg2, const int event, const int col, const int row )
85a699d8	75	{ // evaluate methods
48a2e3eb	76	// cout << __FILE__ << __LINE__ << endl;
	77
	78	AliCaloFitResults ref = fReferenceAnalyzer->Evaluate( bunchvector, altrocfg1, altrocfg2 );
	79	// AliCaloFitResults ref;
	80
	81	for(int i=0; i < fRawAnalyzers.size(); i++ )
	82	{
	83	AliCaloFitResults an = fRawAnalyzers.at(i)->Evaluate( bunchvector, altrocfg1, altrocfg2 );
	84
	85	// cout << __FILE__ << __LINE__ << ":" << fRawAnalyzers.at(i)->GetAlgoAbbr() << " col = " << col << " row " << row << " amp= " << an.GetAmp() << endl ;
	86
	87	// fAmplitudeVsEvent[i]->Fill( event, an.GetAmp() );
	88	// fTofVsEvent[i]->Fill( event, an.GetTof() );
	89
	90	fRefAmpVsAnalyzers[i]->Fill(ref.GetAmp(), an.GetAmp() );
	91	fRefTofVsAnalyzers[i]->Fill(ref.GetTof(), an.GetTof() );
	92	fAmpDiff[i]->Fill(ref.GetAmp() - an.GetAmp() );
	93	fTofDiff[i]->Fill(ref.GetTof() - an.GetTof() );
	94
	95	if(col == fMonCol1 && fMonRow1 == row )
	96	{
	97	fAmplitudeVsEvent[i]->Fill( event, an.GetAmp() );
	98	fTofVsEvent[i]->Fill( event, an.GetTof() );
	99	fMon1[i] = an;
	100	}
	101
	102	if(col == fMonCol2 && row == fMonRow2 )
	103	{
	104	fMon2[i] = an;
	105	}
	106
	107	if(row >= 0 && col >= 0 )
	108	{
	109	fAmpHistograms[i][col][row]->Fill( an.GetAmp() );
	110	}
	111	}
	112	}
	113
	114
	115
	116	void
	117	AliCaloRawAnalyzerComparison::EventChanged()
85a699d8	118	{ // new event
48a2e3eb	119	for(int i=0; i < fRawAnalyzers.size(); i++ )
	120	{
	121	/*
	122	if( ( fMon1[i].GetAmp() > 50 && fMon2[i].GetAmp() > 50 ) && ( fMon1[i].GetAmp() < 1023 && fMon2[i].GetAmp() < 1023 ) )
	123	{
	124	cout << __FILE__ << __LINE__ << fRawAnalyzers.at(i)->GetAlgoAbbr() << "\tamp1=" << fMon1[i].GetAmp() << "\tamp2 = " << fMon2[i].GetAmp() << endl;
	125	}
	126	*/
	127
	128	if( fMon1[i].GetAmp() > 50 && fMon2[i].GetAmp() > 50 )
	129	{
	130	// if( fMon1[i].GetTof() != fMon2[i].GetTof())
	131	{
	132	// fTofResDifferential[i]->Fill( ( fMon1[i].GetTof() - fMon2[i].GetTof())/TMath::Sqrt(2) );
	133	fTofResDifferential[i]->Fill( ( fMon1[i].GetTof() - fMon2[i].GetTof()) );
	134	fTofResAbsolute[i]->Fill( fMon1[i].GetTof() );
	135	}
	136	}
	137	}
	138	}
	139
	140
	141	void
	142	AliCaloRawAnalyzerComparison::WriteHistograms()
85a699d8	143	{ // write histograms
48a2e3eb	144	TFile *f = new TFile("comparison2.root", "recreate");
	145
	146	/*
	147	for(int col=0; col < NZCOLSSMOD; col ++ )
	148	{
	149	for(int row=0; row < NXROWSSMOD; row ++ )
	150	{
	151	fAmpHistograms[col][row]->Write();
	152	}
	153	}
	154	*/
	155
	156	for(int i=0; i < fRawAnalyzers.size(); i++ )
	157	{
	158	for(int col=0; col < NZCOLSSMOD; col ++ )
	159	{
	160	for(int row=0; row < NXROWSSMOD; row ++ )
	161	{
	162	fAmpHistograms[i][col][row]->Write();
	163	}
	164	}
	165
	166	fAmplitudeVsEvent[i]->Write();
	167	fTofVsEvent[i]->Write();
	168	fRefAmpVsAnalyzers[i]->Write();
	169	fRefTofVsAnalyzers[i]->Write();
	170	fAmpDiff[i]->Write();
	171	fTofDiff[i]->Write();
	172	fTofResDifferential[i]->Write();
	173	fTofResAbsolute[i]->Write();
	174
	175	}
	176
	177	f->Close();
	178	}
	179
	180
	181	void
85a699d8	182	AliCaloRawAnalyzerComparison::InitHistograms( vector <AliCaloRawAnalyzer> analyzers, AliCaloRawAnalyzer ref )
85a699d8	183	{ // init histograms
48a2e3eb	184	char tmpname[256];
	185
	186	/*
	187	for(int col=0; col < NZCOLSSMOD; col ++ )
	188	{
	189	for(int row=0; row < NXROWSSMOD; row ++ )
	190	{
	191	sprintf(tmpname, "z(col)%d_x(row)%d_amplitude;Counts;Amplitude/ADC counts", col, row);
	192	fAmpHistograms[col][row] = new TH1D(tmpname, tmpname, 1024, 0, 1023 );
	193	}
	194	}
	195	*/
	196
	197	// TH1D *fAmpHistograms[NZCOLSSMOD][NXROWSSMOD];
	198
	199	for(int i=0; i < analyzers.size(); i++ )
	200	{
	201	for(int col=0; col < NZCOLSSMOD; col ++ )
	202	{
	203	for(int row=0; row < NXROWSSMOD; row ++ )
	204	{
	205	sprintf(tmpname, "z(col)%d_x(row)%d_amplitude_%s;Amplitude/ADC counts;Counts", col, row, analyzers.at(i)->GetAlgoAbbr() );
	206	fAmpHistograms[i][col][row] = new TH1D(tmpname, tmpname, 1024, 0, 1023 );
	207	}
	208	}
	209
	210	sprintf(tmpname, "%s_amplitude_vs_event_row%d_col%d;Event;Amplitude/ADC counts", analyzers.at(i)->GetAlgoAbbr(), fMonRow1, fMonCol1 );
	211	fAmplitudeVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 1024, 0, 1023 );
	212
	213	sprintf(tmpname, "%s_tof_vs_event_row%d_col%d;Event;Amplitude/ADC counts", analyzers.at(i)->GetAlgoAbbr(), fMonRow1, fMonCol1 );
	214
	215	// fTofVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 2000, 3000, 5000);
	216	fTofVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 2000, 1000, 4000);
	217
	218
	219	sprintf(tmpname, "%s_vs_%s_ampltude; Amplitude_{%s}/ADC counts; Amplitude_{%s}/ADC counts", ref->GetAlgoAbbr(),
	220	analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr() );
	221	fRefAmpVsAnalyzers[i] = new TH2D(tmpname, tmpname, 1024, 0, 1023, 1024, 0, 1023);
	222	sprintf(tmpname, "%s_vs_%s_tof; tof_{%s}/ns; tof_{%s}/ns", ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr() );
	223	fRefTofVsAnalyzers[i] = new TH2D(tmpname, tmpname, 500, 2000, 4999, 500, 2000, 4999 );
	224	sprintf( tmpname, "%s-%s amplitude;counts;A_{%s} - A_{%s}", ref->GetAlgoAbbr(),
	225	analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr());
	226	fAmpDiff[i] = new TH1D(tmpname, tmpname, 100, -10, 10 );
	227	sprintf( tmpname, "%s-%s tof;counts;A_{%s} - A_{%s}", ref->GetAlgoAbbr(),
	228	analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr());
	229	fTofDiff[i] = new TH1D(tmpname, tmpname, 1000, -5000, 5000 );
	230	sprintf( tmpname, "%s Differential tof resolution (%d, %d) vs (%d, %d);#sigma_{tof}^{%s}/ns;Counts", analyzers.at(i)->GetAlgoAbbr(), fMonCol1, fMonRow1, fMonCol2, fMonRow2,
	231	analyzers.at(i)->GetAlgoAbbr() );
	232	fTofResDifferential[i] = new TH1D(tmpname, tmpname, 1000, -250, 250 );
	233	sprintf( tmpname, "%s Absolute tof distribution (%d, %d);#sigma_{tof}^{%s}/ns;Counts", analyzers.at(i)->GetAlgoAbbr(), fMonCol1, fMonRow1, analyzers.at(i)->GetAlgoAbbr() );
	234
	235	fTofResAbsolute[i] = new TH1D(tmpname, tmpname, 2000 , -1000, 7000 );
	236
	237
	238	}
	239	}
	240