]> git.uio.no Git - u/mrichter/AliRoot.git/blob - ITS/AliITSQASSDRefData.cxx
A skeleton of a comparison macro for pileup studies
[u/mrichter/AliRoot.git] / ITS / AliITSQASSDRefData.cxx
1 /**************************************************************************
2  * Copyright(c) 2009-2011, 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
19 //-------------------------------------------------------------------------
20 //                          Class AliITSQASSDRefData
21 //                     ITS SSD reference values for the QA
22 //
23 //         Origin: Panos.Christakoglou@cern.ch, NIKHEF-Utrecht University
24 //-------------------------------------------------------------------------
25
26 #include <Riostream.h>
27 #include <fstream>
28 #include <TArray.h>
29 #include <TString.h>
30 #include <TObjString.h>
31 #include <TObjArray.h>
32
33 #include "AliLog.h"
34 #include "AliITSQASSDRefData.h"
35
36 using std::ifstream;
37 ClassImp(AliITSQASSDRefData)
38
39 //___________________________________________________________________________
40 AliITSQASSDRefData::AliITSQASSDRefData() :
41   TObject(),
42   fRefList(0),
43   fNameList(0) { 
44   //Default constructor
45 }
46
47 //___________________________________________________________________________
48 AliITSQASSDRefData::AliITSQASSDRefData(Int_t specie) :
49   TObject(),
50   fRefList(0),
51   fNameList(0) { 
52   //Default constructor
53   SetDefault(specie);
54 }
55
56 //___________________________________________________________________________
57 AliITSQASSDRefData::AliITSQASSDRefData(const char* path) :
58   TObject(),
59   fRefList(0),
60   fNameList(0) {
61   //Constructor with the path of the ascii file as an argument
62   SetReferenceData(path);
63 }
64
65 //___________________________________________________________________________
66 AliITSQASSDRefData::AliITSQASSDRefData(const AliITSQASSDRefData& refData):
67 TObject(),
68 fRefList(refData.fRefList),
69 fNameList(refData.fNameList) {
70   //Copy constructor
71 }
72
73 //___________________________________________________________________________
74 AliITSQASSDRefData& AliITSQASSDRefData::operator = (const AliITSQASSDRefData& refData) {
75   //assignment operator
76   if(&refData != this) {
77     fRefList = refData.fRefList;
78     fNameList = refData.fNameList;
79   }
80   return *this ;
81 }
82
83 //___________________________________________________________________________
84 AliITSQASSDRefData::~AliITSQASSDRefData() { 
85   //Destructor
86   if(fRefList) delete fRefList;
87   if(fNameList) delete fNameList;
88 }
89
90 //___________________________________________________________________________
91 void AliITSQASSDRefData::AddReference(const char* name="", 
92                                       Int_t id=-1, 
93                                       Double_t value=0) {
94   //Adding a ref. value to the list
95   //Printf("(AliITSQASSDRefData::AddReference) Name: %s - Id: %d - Value: %lf",name,id,value);
96   if(id>-1&&id<fRefList->GetSize()) {
97     AliError(Form("Reference with id %i already exists. Choose other id or use SetReferenceValue(Int_t, Double_t) to overwrite",id));
98     return;
99   }
100   
101   if( (strcmp(name,"")!=0) && GetID(name)!=-1) {
102     AliError(Form("Reference with name %s already exists. Choose other name or use SetReferenceValue(const char*, Double_t) to overwrite",name));
103     return;
104   }
105   
106   if(id==-1) id=fRefList->GetSize();
107   fRefList->Set(id+1);
108   fRefList->AddAt(value,id);
109   fNameList->AddAt(new TObjString(name),id);
110 }
111
112 //___________________________________________________________________________
113 Int_t AliITSQASSDRefData::GetID(const char* name) {
114   //Get the id of the reference value
115   Int_t status = -1;
116   TString refName = "";
117   TString stringName = name;
118   TObjString *dummyString = 0;
119   for (Int_t id=0; id<fNameList->GetEntriesFast(); id++){
120     dummyString = static_cast <TObjString *>(fNameList->At(id));
121     refName = dummyString->GetString();
122     if(refName == stringName) {
123       status = id;
124     }
125   }
126
127   return status;
128 }
129
130 //___________________________________________________________________________
131 Double_t AliITSQASSDRefData::GetReferenceValue(const char* name) {
132   //Returns the ref. value based on the given name
133   TString refName = "";
134   TObjString *dummyString = 0;
135   for (Int_t id=0; id<fNameList->GetEntriesFast(); id++){
136     dummyString = static_cast <TObjString *>(fNameList->At(id));
137     refName = dummyString->GetString();
138
139     if(refName.Data()==name) return fRefList->At(id);
140   }
141   AliError(Form("Reference name %s unknown",name));
142   return -1;
143 }
144
145 //___________________________________________________________________________
146 Double_t AliITSQASSDRefData::GetReferenceValue(Int_t id) {
147   //Returns the ref. value based on the given id
148   if (id<0||id>fRefList->GetSize()-1){
149     AliError("Reference ID out of range");
150     return 0;
151   }
152   return fRefList->At(id);
153
154 }
155
156 //___________________________________________________________________________
157 void AliITSQASSDRefData::PrintTable() {
158   //Prints the list of reference values
159   Printf("___ SSD REFERENCE DATA ___ ");
160   Printf("ID ----- Value ------- Name");
161   Int_t id=0;
162   TString refName = "";
163   TObjString *dummyString = 0;
164   for(id=0; id<fRefList->GetSize()-1; id++) {
165     dummyString = static_cast <TObjString *>(fNameList->At(id));
166     refName = dummyString->GetString();
167     Printf("%i ------ %4.3g -------- %s",id,fRefList->At(id),refName.Data());
168            
169   }
170
171 }
172
173 //___________________________________________________________________________
174 void AliITSQASSDRefData::SetDefault(Int_t specie) {
175   //Sets the default values
176   if(!fNameList) fNameList = new TObjArray();
177   fNameList->Add(new TObjString("minSSDDataSize"));
178   fNameList->Add(new TObjString("maxSSDDataSize"));
179   fNameList->Add(new TObjString("minDDLDataSize"));
180   fNameList->Add(new TObjString("maxDDLDataSize"));
181   fNameList->Add(new TObjString("minLDCDataSize"));
182   fNameList->Add(new TObjString("maxLDCDataSize"));
183   fNameList->Add(new TObjString("minMeanDDLDataSize"));
184   fNameList->Add(new TObjString("maxMeanDDLDataSize"));
185   fNameList->Add(new TObjString("minMeanLDCDataSize"));
186   fNameList->Add(new TObjString("maxMeanLDCDataSize"));
187   fNameList->Add(new TObjString("maxOccupancy"));
188   fNameList->SetOwner(kTRUE);
189   
190   //specie == 0 ==> Default
191   Double_t refValues[11] = {0,0.0,0,0,0,0,0,0,0,0};
192   //specie == 1 ==> Low multiplicity
193   if(specie == 1) {
194     refValues[0] = 0; refValues[1] = 500; refValues[2] = 0; refValues[3] = 50;
195     refValues[4] = 0; refValues[5] = 100; refValues[6] = 0; refValues[7] = 50;
196     refValues[8] = 0; refValues[9] = 100; refValues[10] = 5;
197   }
198   //specie == 2 ==> High multiplicity
199   if(specie == 2) {
200     refValues[0] = 0; refValues[1] = 500; refValues[2] = 0; refValues[3] = 50;
201     refValues[4] = 0; refValues[5] = 100; refValues[6] = 0; refValues[7] = 50;
202     refValues[8] = 0; refValues[9] = 100; refValues[10] = 5;
203   }
204   //specie == 3 ==> Cosmics
205   if(specie == 3) {
206     refValues[0] = 0; refValues[1] = 500; refValues[2] = 0; refValues[3] = 50;
207     refValues[4] = 0; refValues[5] = 100; refValues[6] = 0; refValues[7] = 50;
208     refValues[8] = 0; refValues[9] = 100; refValues[10] = 5;
209   }
210   //specie == 4 ==> Calibration
211   if(specie == 4) {
212     refValues[0] = 0; refValues[1] = 500; refValues[2] = 0; refValues[3] = 50;
213     refValues[4] = 0; refValues[5] = 100; refValues[6] = 0; refValues[7] = 50;
214     refValues[8] = 0; refValues[9] = 100; refValues[10] = 5;
215   }
216
217   if(!fRefList) fRefList = new TArrayD();
218   fRefList->Set(11,refValues);
219 }
220
221 //___________________________________________________________________________
222 void AliITSQASSDRefData::SetReferenceData(const char* path) {
223   //Parses an ascii file with the reference values
224   if(!fRefList) fRefList = new TArrayD();
225   if(!fNameList) fNameList = new TObjArray();
226
227   ifstream file;
228   file.open(path);
229   
230   if (!file) {
231     AliWarning(Form("No file found at %s",path));
232     SetDefault(0);
233     return;
234   }
235   if(file.bad()){
236     AliWarning("Reference data could not be read: Default values are used.");
237     SetDefault(0);
238     return;
239   }
240   
241   fRefList->Set(0);
242   Int_t id = 0, newid = -1;
243   Double_t value = 0;
244   TString name = "";
245   
246   while (!file.eof()) {
247     //file >> newid;
248     file >> name >> id >> value;
249     //Printf("Name: %s - Id: %d - Value: %lf",name.Data(),id,value);
250     
251     if (newid==id) continue; //skip line if id is the same as previous
252     AddReference(name.Data(), id, value);
253     newid = id;
254   }
255
256   file.close();
257 }
258
259 //___________________________________________________________________________
260 void AliITSQASSDRefData::SetReferenceValue(Int_t id, Double_t value) {
261   //Adding a single reference value by id
262   if(id<0||id>fRefList->GetSize()-1) {
263     AliWarning(Form("Reference ID %i out of range: value not set",id));
264   }
265   else fRefList->SetAt(value,id);
266
267 }
268
269 //___________________________________________________________________________
270 void AliITSQASSDRefData::SetReferenceValue(const char* name, Double_t value) {
271   //Adding a single reference value by name
272   Int_t id = GetID(name);
273   //Printf("Name: %s - Id: %d",name,id);
274   if(id == -1) {
275     AliError(Form("Reference name %s unknown: value not set",name));
276     return;
277   }
278    
279   fRefList->SetAt(value,id);
280 }