CMake: removing qpythia from the depedencies
[u/mrichter/AliRoot.git] / EMCAL / AliCaloRawAnalyzerComparison.cxx
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 "AliCaloRawAnalyzerKStandard.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)
40 {
41   // ctor
42   
43   fReferenceAnalyzer                       = new  AliCaloRawAnalyzerKStandard();
44
45   AliCaloRawAnalyzerCrude *crude           = new AliCaloRawAnalyzerCrude();
46   AliCaloRawAnalyzerKStandard *lms         = new AliCaloRawAnalyzerKStandard();
47   AliCaloRawAnalyzerFastFit *fastfit       = new AliCaloRawAnalyzerFastFit();
48   AliCaloRawAnalyzerNN *neuralnet          = new AliCaloRawAnalyzerNN();
49   AliCaloRawAnalyzerPeakFinder *peakfinder = new AliCaloRawAnalyzerPeakFinder();
50
51   fReferenceAnalyzer->SetIsZeroSuppressed();
52   peakfinder        ->SetIsZeroSuppressed();
53   lms               ->SetIsZeroSuppressed();
54   fastfit           ->SetIsZeroSuppressed();
55   neuralnet         ->SetIsZeroSuppressed();
56   crude             ->SetIsZeroSuppressed();
57   
58   fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)crude);
59   fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)lms);
60   fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)fastfit);
61   fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)neuralnet);
62   
63   fRawAnalyzers.push_back( (AliCaloRawAnalyzer*)peakfinder);
64   InitHistograms( fRawAnalyzers, fReferenceAnalyzer );
65
66
67 void 
68 AliCaloRawAnalyzerComparison::Evaluate( const vector<AliCaloBunchInfo> &bunchvector, const UInt_t altrocfg1,  
69                                         const UInt_t altrocfg2, const int event, const int col, const int row )
70 { // evaluate methods
71   //  cout << __FILE__ << __LINE__ << endl;
72   
73   AliCaloFitResults ref = fReferenceAnalyzer->Evaluate( bunchvector, altrocfg1, altrocfg2 );
74   //  AliCaloFitResults ref;
75   
76   for(int i=0; i < fRawAnalyzers.size(); i++ )
77   {
78     AliCaloFitResults an = fRawAnalyzers.at(i)->Evaluate( bunchvector, altrocfg1,  altrocfg2 );
79     
80     //     cout << __FILE__ << __LINE__ << ":" << fRawAnalyzers.at(i)->GetAlgoAbbr() << "  col = " << col << " row " << row << "  amp= " <<  an.GetAmp() << endl ;
81     
82     //   fAmplitudeVsEvent[i]->Fill( event, an.GetAmp() );
83     //   fTofVsEvent[i]->Fill( event, an.GetTof() );
84     
85     fRefAmpVsAnalyzers[i]->Fill(ref.GetAmp(), an.GetAmp() );
86     fRefTofVsAnalyzers[i]->Fill(ref.GetTof(), an.GetTof() );
87     fAmpDiff[i]->Fill(ref.GetAmp() - an.GetAmp() );
88     fTofDiff[i]->Fill(ref.GetTof() - an.GetTof() );
89     
90     if(col  == fMonCol1 && fMonRow1 == row )
91     {
92       fAmplitudeVsEvent[i]->Fill( event, an.GetAmp() );
93       fTofVsEvent[i]->Fill( event, an.GetTof() );
94       fMon1[i] = an;
95     }
96     
97     if(col  == fMonCol2 &&  row == fMonRow2  )
98     {
99       fMon2[i] = an;
100     }
101     
102     if(row >= 0 && col >= 0 )
103     {
104       fAmpHistograms[i][col][row]->Fill( an.GetAmp() );
105     }
106   }
107 }  
108
109 void 
110 AliCaloRawAnalyzerComparison::EventChanged()
111 { // new event
112   for(int i=0; i < fRawAnalyzers.size(); i++ )
113   {
114     /*
115      if( ( fMon1[i].GetAmp() > 50 &&   fMon2[i].GetAmp() > 50 ) &&  (  fMon1[i].GetAmp() < 1023  &&   fMon2[i].GetAmp() < 1023 )  )
116      {
117      cout << __FILE__ << __LINE__ << fRawAnalyzers.at(i)->GetAlgoAbbr() << "\tamp1=" <<  fMon1[i].GetAmp()  << "\tamp2 = " <<  fMon2[i].GetAmp() << endl;
118      }
119      */
120     
121     if( fMon1[i].GetAmp()  > 50  &&  fMon2[i].GetAmp() > 50  )
122     {
123       //  if(  fMon1[i].GetTof() != fMon2[i].GetTof())
124       {
125         //    fTofResDifferential[i]->Fill( ( fMon1[i].GetTof() - fMon2[i].GetTof())/TMath::Sqrt(2) );
126         fTofResDifferential[i]->Fill( ( fMon1[i].GetTof() - fMon2[i].GetTof()) );
127         fTofResAbsolute[i]->Fill(  fMon1[i].GetTof() );
128       }
129     }
130   }
131 }
132
133
134 void 
135 AliCaloRawAnalyzerComparison::WriteHistograms()
136 { // write histograms
137   TFile *f = new TFile("comparison2.root", "recreate");
138   
139   /*
140    for(int col=0; col < NZCOLSSMOD; col ++ )
141    {
142    for(int row=0; row < NXROWSSMOD; row ++ )
143    {
144    fAmpHistograms[col][row]->Write();
145    }
146    }
147    */
148   
149   for(int i=0; i < fRawAnalyzers.size(); i++ )
150   {
151     for(int col=0; col < NZCOLSSMOD; col ++ )
152     {
153       for(int row=0; row < NXROWSSMOD; row ++ )
154             {
155               fAmpHistograms[i][col][row]->Write();
156             }
157     }
158     
159     fAmplitudeVsEvent[i]->Write();
160     fTofVsEvent[i]->Write();
161     fRefAmpVsAnalyzers[i]->Write();
162     fRefTofVsAnalyzers[i]->Write();
163     fAmpDiff[i]->Write();
164     fTofDiff[i]->Write();
165     fTofResDifferential[i]->Write();
166     fTofResAbsolute[i]->Write();
167     
168   }
169   
170   f->Close();
171 }
172
173
174 void
175 AliCaloRawAnalyzerComparison::InitHistograms( vector <AliCaloRawAnalyzer*> analyzers, AliCaloRawAnalyzer* ref )
176 { // init histograms
177   char tmpname[256];
178   
179   /*
180    for(int col=0; col < NZCOLSSMOD; col ++ )
181    {
182    for(int row=0; row < NXROWSSMOD; row ++ )
183    {
184    sprintf(tmpname, "z(col)%d_x(row)%d_amplitude;Counts;Amplitude/ADC counts", col, row);
185    fAmpHistograms[col][row] = new TH1D(tmpname, tmpname, 1024, 0, 1023 );
186    }
187    }
188    */
189   
190   // TH1D *fAmpHistograms[NZCOLSSMOD][NXROWSSMOD];
191   
192   for(int i=0; i < analyzers.size(); i++ )
193   {
194     for(int col=0; col < NZCOLSSMOD; col ++ )
195     {
196       for(int row=0; row < NXROWSSMOD; row ++ )
197             {
198               sprintf(tmpname, "z(col)%d_x(row)%d_amplitude_%s;Amplitude/ADC counts;Counts", col, row,  analyzers.at(i)->GetAlgoAbbr() );
199               fAmpHistograms[i][col][row] = new TH1D(tmpname, tmpname, 1024, 0, 1023 );
200             }
201     }
202     
203     sprintf(tmpname, "%s_amplitude_vs_event_row%d_col%d;Event;Amplitude/ADC counts", analyzers.at(i)->GetAlgoAbbr(), fMonRow1, fMonCol1 );
204     fAmplitudeVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 1024, 0, 1023 );
205     
206     sprintf(tmpname, "%s_tof_vs_event_row%d_col%d;Event;Amplitude/ADC counts", analyzers.at(i)->GetAlgoAbbr(), fMonRow1, fMonCol1 );
207     
208     //     fTofVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 2000, 3000, 5000);
209     fTofVsEvent[i]=new TH2D(tmpname, tmpname, 8000, 0, 7999, 2000, 1000, 4000);
210     
211     
212     sprintf(tmpname, "%s_vs_%s_ampltude; Amplitude_{%s}/ADC counts; Amplitude_{%s}/ADC counts", ref->GetAlgoAbbr(),
213             analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr() );
214     fRefAmpVsAnalyzers[i] = new TH2D(tmpname, tmpname, 1024, 0, 1023, 1024, 0, 1023);
215     sprintf(tmpname, "%s_vs_%s_tof; tof_{%s}/ns; tof_{%s}/ns", ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr() );
216     fRefTofVsAnalyzers[i] = new TH2D(tmpname, tmpname, 500, 2000, 4999, 500, 2000, 4999 );
217     sprintf( tmpname, "%s-%s amplitude;counts;A_{%s} - A_{%s}",   ref->GetAlgoAbbr(),
218             analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr());
219     fAmpDiff[i] = new TH1D(tmpname, tmpname, 100, -10, 10 );
220     sprintf( tmpname, "%s-%s tof;counts;A_{%s} - A_{%s}",   ref->GetAlgoAbbr(),
221             analyzers.at(i)->GetAlgoAbbr(), ref->GetAlgoAbbr(), analyzers.at(i)->GetAlgoAbbr());
222     fTofDiff[i] = new TH1D(tmpname, tmpname, 1000, -5000, 5000 );
223     sprintf( tmpname, "%s Differential tof resolution (%d, %d) vs (%d, %d);#sigma_{tof}^{%s}/ns;Counts",  analyzers.at(i)->GetAlgoAbbr(), fMonCol1, fMonRow1, fMonCol2, fMonRow2,
224             analyzers.at(i)->GetAlgoAbbr() );
225     fTofResDifferential[i] = new TH1D(tmpname, tmpname, 1000, -250, 250 );
226     sprintf( tmpname, "%s Absolute tof distribution (%d, %d);#sigma_{tof}^{%s}/ns;Counts",  analyzers.at(i)->GetAlgoAbbr(), fMonCol1, fMonRow1, analyzers.at(i)->GetAlgoAbbr() );
227     
228     fTofResAbsolute[i] = new TH1D(tmpname, tmpname, 2000 , -1000, 7000 );
229     
230     
231   }
232 }
233