For Pythia with tune don't switch off MI in ConfigHeavyFlavor
[u/mrichter/AliRoot.git] / ITS / AliITSOnlineSPDfo.cxx
1 /**************************************************************************
2  * Copyright(c) 2008-2010, 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 // Author: A. Mastroserio                                     // 
20 // This class is used within the detector algorithm framework //
21 // to write and read FO scan data.                            //
22 ////////////////////////////////////////////////////////////////
23
24 #include <TFile.h>
25 #include <Riostream.h>
26 #include <TSystem.h>
27 #include <TObjArray.h>
28 #include <TObjString.h>
29 #include <TArrayI.h>
30 #include <TIterator.h>
31 #include <TKey.h>
32 #include <TMath.h>
33 #include <TList.h>
34 #include "AliLog.h"
35 #include "AliITSOnlineSPDfoChipConfig.h"
36 #include "AliITSOnlineSPDfoChip.h"
37 #include "AliITSOnlineSPDfoInfo.h"
38 #include "AliITSOnlineSPDfo.h"
39
40
41 ClassImp(AliITSOnlineSPDfo)
42 //-----------------------------------------------
43 AliITSOnlineSPDfo::AliITSOnlineSPDfo():
44 fRunNr(0),
45 fNdacs(0),
46 fFileName(""),
47 fFile(0x0),
48 fInfo(0x0),
49 fDACnames(0x0),
50 fArray(0x0),
51 fCheckIndex(-1),
52 fIndex(-1),
53 fInitialConfiguration("")
54 {
55 // default constructor
56
57
58 //-------------------------------------------------
59 AliITSOnlineSPDfo::AliITSOnlineSPDfo(TString inputfile, Int_t runNr, Int_t eqId):
60 fRunNr(runNr),
61 fNdacs(0),
62 fFileName(""),
63 fFile(0x0),
64 fInfo(0x0),
65 fDACnames(0x0),
66 fArray(0x0),
67 fCheckIndex(-1),
68 fIndex(-1),
69 fInitialConfiguration("")
70 {
71  //
72  // constructor 
73  //   
74     fFileName=Form("%i_%s%02i.root",runNr,inputfile.Data(),eqId);
75     fFile = TFile::Open(fFileName.Data());
76     fArray = new TObjArray();
77     fDACnames = new THashList();
78 }
79 //--------------------------------------------------
80 AliITSOnlineSPDfo::AliITSOnlineSPDfo(const AliITSOnlineSPDfo &c):
81 fRunNr(c.fRunNr),
82 fNdacs(c.fNdacs),
83 fFileName(c.fFileName),
84 fFile(c.fFile),
85 fInfo(c.fInfo),
86 fDACnames(c.fDACnames),
87 fArray(c.fArray),
88 fCheckIndex(c.fCheckIndex),
89 fIndex(c.fIndex),
90 fInitialConfiguration(c.fInitialConfiguration)
91 {
92   //
93   //copy constructor
94   //
95 }    
96 //--------------------------------------------------
97 void AliITSOnlineSPDfo::SetFile(TString inputfile)
98 {
99   //
100   // open the file where the data are
101   //
102   
103  if(fFile) {
104    fInfo = (AliITSOnlineSPDfoInfo *)fFile->Get("generalinfo");
105    return; 
106  } else {
107    fFile = TFile::Open(inputfile.Data());
108    if(!fFile) {
109      Info("AliITSOnlineSPDfo::SetFile",Form(" %s  not existing.... The scan info are not available....crash is expected \n",inputfile.Data()));
110      return;
111     } else {
112      fInfo = (AliITSOnlineSPDfoInfo *)fFile->Get("generalinfo");
113      fNdacs = fInfo->GetNumDACindex();
114    }
115 }
116
117 }
118
119 //------------------------------------------------------
120 void AliITSOnlineSPDfo::CreateOutputFile()
121 {
122   //
123   // Create the file (needed only in the DA ), will delete the previous one!
124   //  
125   if(fFile) Info("AliITSOnlineSPDfo::CreateOutputFile","removing previous file....\n");
126    fFile = TFile::Open(fFileName.Data(),"RECREATE"); 
127   
128 }
129 //------------------------------------------------------
130 void AliITSOnlineSPDfo::AddMeasurement(const TArrayS dac, Short_t measure[4], Int_t hs, Int_t chipId)
131 {
132   /*
133   // Here a single measurement is added to the chip container. 
134   // A single measurement corresponds to a specific pixel-configuration output
135   // in the Fast-OR chip. If N configurations are considered, then
136   // the data structure is the following:
137   //
138   //                                          -> HS0_CHIP0 -> measure0[4]
139   //                                         /                measure1[4]
140   //                                        /                   ...
141   //  fArray->At(i) =  DAC1-DAC2-DAC3-DAC4                    measureN[4]
142   //                                        \
143   //                                         \ 
144   //                                          -> HS0_CHIP1 -> measure0[4]
145   //                                                          measure1[4]
146   //                                                             ....
147   //                                                          measureN[4]
148   //
149   */
150
151   AliITSOnlineSPDfoChipConfig *counts = new AliITSOnlineSPDfoChipConfig(measure);
152   Int_t arrayelement = CheckDACEntry(dac);
153   
154   if(arrayelement< 0){
155     TString dacname = CreateDACEntry(dac);
156     TObjString *string = new TObjString(dacname.Data());
157     fDACnames->Add(string);
158      
159     TObjArray *array = new TObjArray(60);
160     AliITSOnlineSPDfoChip * chip = new AliITSOnlineSPDfoChip(dac.GetSize());
161     chip->SetActiveHS(hs);
162     chip->SetChipId(chipId); 
163     for(Int_t i=0; i< dac.GetSize() ; i++) chip->SetDACParameter(i,dac.At(i));               
164     chip->AddMeasurement(counts);
165     array->AddAt(chip,hs*10+chipId);  
166     fArray->AddLast(array); 
167      
168   } else {
169     
170     TObjArray *arr = (TObjArray*)fArray->At(arrayelement);  
171     if(!arr->At(hs*10+chipId)){
172       AliITSOnlineSPDfoChip * chip = new AliITSOnlineSPDfoChip(dac.GetSize());
173       chip->SetActiveHS(hs);
174       chip->SetChipId(chipId); 
175       for(Int_t i=0; i< dac.GetSize() ; i++) chip->SetDACParameter(i,dac.At(i));               
176       chip->AddMeasurement(counts);
177       arr->AddAt(chip,hs*10+chipId);     
178       
179       } else {
180       
181       AliITSOnlineSPDfoChip *c = (AliITSOnlineSPDfoChip *)arr->At(hs*10+chipId);
182       if(c)c->AddMeasurement(counts);    
183     }   
184   }
185 }
186 //---------------------------------------
187 Int_t AliITSOnlineSPDfo::CheckDACEntry(const TArrayS dac)
188 {
189   //
190   // Check if the set of dacs has been already added to the array
191   // 
192   
193    
194   TString name = CreateDACEntry(dac);
195   
196   if(!fDACnames) {
197     Info("AliITSOnlineSPDfo::CheckDACEntry"," NO DAC name array is present, exiting.... \n");
198     return -1;
199   }
200   
201    Double_t c = 0;
202   for(Int_t i=0; i< fNdacs; i++)  c+=dac.At(i)*TMath::Power(10,3*i);
203   if(c==fCheckIndex) return fIndex;
204    
205   
206  
207    Int_t idx = -1;
208    
209   if(fDACnames->FindObject( name.Data() )) {
210     idx = fDACnames->IndexOf( fDACnames->FindObject( name.Data() )  );
211     fCheckIndex=0;
212     for(Int_t i=0; i< fNdacs; i++) fCheckIndex+=dac.At(i)*TMath::Power(10,3*i);
213     fIndex=idx;
214
215   }
216   return idx;
217
218 }
219 //_________________________________________________________
220 void AliITSOnlineSPDfo::WriteToFile()
221 {
222   //
223   //The array of DACS and all its content is written to file.
224   // Here the general info on the FO calibration scan are 
225   // written in the same file
226   
227   if(fDACnames->GetEntries() != fArray->GetEntries()) {
228    printf("mismatch names-array. Exiting....");   
229     return;
230   }
231   
232   for(Int_t i=0; i< fDACnames->GetEntries(); i++){
233   fFile->WriteObject(fArray->At(i),fDACnames->At(i)->GetName());
234   }
235   fFile->WriteTObject(fInfo,"generalinfo"); 
236   fFile->Close();
237   
238 }
239 //______________________________________________________
240
241 TString AliITSOnlineSPDfo::CreateDACEntry(const TArrayS dacs) const
242 {
243   //
244   // The string of DACs is build
245   // 
246   
247   TString dacvalues;
248   for(Int_t i=0; i<dacs.GetSize(); i++) dacvalues+=Form("-%i",dacs.At(i));
249   dacvalues.Remove(0,1);
250   return dacvalues;
251  
252 }
253 //_____________________________________________________
254 TArrayI AliITSOnlineSPDfo::GetDACscanParams() const
255
256   //
257   // this method retrieves the DAC value range and its steps
258   //
259   TFile *f =0x0;
260   if(fFile->IsOpen()) f = fFile;
261   else f = TFile::Open(fFileName.Data());
262   
263   TArrayI dacs;
264   if(f->GetNkeys() < 2) return dacs;
265   dacs.Set(fNdacs*3);
266   
267   TArrayI min(fNdacs), max(fNdacs),step(fNdacs), check(fNdacs), refvalues(fNdacs);
268   
269   for(Int_t i=0; i< fNdacs ; i++) {
270     min.AddAt(9999,i);
271     max.AddAt(0,i);
272     step.AddAt(0,i);
273     check.AddAt(0,i);
274   }
275   
276   
277   TKey *key;
278   TIter iter(f->GetListOfKeys());  
279   while ((key = (TKey*)(iter.Next()))) {
280     TString classname = key->GetClassName();
281     if(classname.Contains("OnlineSPD")) break;
282     
283     TString values = key->GetName();
284     Int_t *val = GetDACvalues(values,fNdacs); // the user has to delete it!
285       
286     for(Int_t i=0; i< fNdacs; i++) {
287       if(val[i]<=min.At(i)) min.AddAt(val[i],i);
288       if(val[i] >= max.At(i)) max.AddAt(val[i],i);
289       // procedure to get the step size;
290       if(!check.At(i)) { 
291         refvalues.AddAt(val[i],i);
292         check.AddAt(1,i);
293       }
294       if(step.At(i) ==0 && check.At(i)){
295         if(val[i]!=refvalues.At(i)) step.AddAt(TMath::Abs(refvalues.At(i) - val[i]),i);
296       }  
297       // end procedure to get the step size 
298     }   
299     delete [] val;
300   }
301
302  
303  
304   for(Int_t i=0; i<fNdacs; i++) {
305     dacs.AddAt(min.At(i),3*i);    
306     dacs.AddAt(max.At(i),3*i+1); 
307     dacs.AddAt(step.At(i),3*i+2); 
308   }
309  
310   
311   return dacs;     
312 }
313 //___________________________________________________________________
314 Int_t* AliITSOnlineSPDfo::GetDACvalues(TString s, const Int_t ndacs) const
315 {
316   //
317   // Translates the string of DACS values into an array of integers
318   //
319   
320    Int_t *val = new Int_t[ndacs]; 
321    
322    s.ReplaceAll("-"," ");
323    char *pEnd = Form("%s",s.Data());
324    for(Int_t i=0; i< ndacs; i++) {
325      val[i] = strtol(pEnd,&pEnd,10); // conversion from string to long
326    }  
327    return val;   
328 }
329
330 //___________________________________________________________________
331 Double_t* AliITSOnlineSPDfo::GetDACvaluesD(TString s, const Int_t ndacs) const 
332 {  
333   //
334   // Translates the string of DACS values into an array of doubles
335   // (needed to fill the thnsparse)
336   //
337    Double_t *val = new Double_t[ndacs];
338    Int_t *values =  GetDACvalues(s,ndacs);
339    for(Int_t i=0; i< ndacs; i++) val[i] = (Double_t)values[i];
340    delete [] values;
341    
342    return val;   
343 }
344 //-------------------------------------------------------------------
345 TArrayS AliITSOnlineSPDfo::CreateDACArray(const TArrayS dacs, const TArrayS dacId) const
346 {
347   //
348   // method to order the data according to the DAC index
349   //
350   
351   TArrayS dacarray(dacs.GetSize());
352   for(Int_t i=0; i<dacs.GetSize(); i++) {
353     if(dacId.At(i)==kIdFOPOL) dacarray.AddAt(dacs.At(i),kFOPOL);
354     else if(dacId.At(i)==kIdCONVPOL) dacarray.AddAt(dacs.At(i),kCONVPOL);
355     else if(dacId.At(i)==kIdCOMPREF) dacarray.AddAt(dacs.At(i),kCOMPREF);
356     else if(dacId.At(i)==kIdPreVTH) dacarray.AddAt(dacs.At(i),kPreVTH);
357     else if(dacId.At(i)==kIdCGPOL) dacarray.AddAt(dacs.At(i),kCGPOL);
358     else printf("new DAC included in the scan??\n");
359     
360   }
361  
362   return dacarray;
363   
364 }
365   
366