Updated setting of the storages
[u/mrichter/AliRoot.git] / TOF / TOFPreprocessor.C
1 /* 
2 $Id$ 
3 */
4
5 // This class runs the test TOF preprocessor
6 // It uses AliTestShuttle to simulate a full Shuttle process
7
8 // The input data is created in the functions
9 //   CreateDCSAliasMap() creates input that would in the same way come from DCS
10 //   ReadDCSAliasMap() reads from a file
11 //   CreateInputFilesMap() creates a list of local files, that can be accessed by the shuttle
12
13 extern TBenchmark *gBenchmark;
14 void TOFPreprocessor()
15 {
16   gSystem->Load("$ALICE/SHUTTLE/TestShuttle/libTestShuttle.so");
17
18   // initialize location of CDB
19   AliTestShuttle::SetMainCDB("local://$ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB");
20   AliTestShuttle::SetMainRefStorage("local://$ALICE_ROOT/SHUTTLE/TestShuttle/TestReference");
21
22   // create AliTestShuttle instance
23   AliTestShuttle* shuttle = new AliTestShuttle(0, 0, 1000);
24
25   // Generation of "fake" input DCS data
26   TMap* dcsAliasMap = CreateDCSAliasMap();  
27
28   // now give the alias map to the shuttle
29   shuttle->SetDCSInput(dcsAliasMap);   
30
31   // processing files. for the time being, the files are local.
32   shuttle->AddInputFile(AliTestShuttle::kDAQ, "TOF", "DELAYS", "MON", "$ALICE_ROOT/TOF/ShuttleInput/Total.root");
33   shuttle->AddInputFile(AliTestShuttle::kDAQ, "TOF", "RUNLevel", "MON", "$ALICE_ROOT/TOF/ShuttleInput/Partial.root");
34
35   // instantiation of the preprocessor
36   AliPreprocessor* pp = new AliTOFPreprocessor(shuttle);
37
38   // preprocessing
39   gBenchmark->Start("process");
40   shuttle->Process();
41   gBenchmark->Stop("process");
42   gBenchmark->Print("process");
43
44   // checking the file which should have been created  
45   AliCDBEntry* chkEntry = AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
46                         ->Get("TOF/Calib/OnlineDelay", 0);
47   if (!chkEntry)
48   {
49     printf("The file is not there. Something went wrong.\n");
50     return;
51   }
52
53   AliTOFDataDCS* output = dynamic_cast<AliTOFDataDCS*> (chkEntry->GetObject());
54   // If everything went fine, draw the result
55   if (output)
56     output->Draw();
57   
58 }
59
60 TMap* CreateDCSAliasMap()
61 {
62   // Creates a DCS structure
63   // The structure is the following:
64   //   TMap (key --> value)
65   //     <DCSAlias> --> <valueList>
66   //     <DCSAlias> is a string
67   //     <valueList> is a TObjArray of AliDCSValue
68   //     An AliDCSValue consists of timestamp and a value in form of a AliSimpleValue
69
70   // In this example 6 aliases exists: DCSAlias1 ... DCSAlias6
71   // Each contains 1000 values randomly generated by TRandom::Gaus + 5*nAlias
72
73   TMap* aliasMap = new TMap;
74   aliasMap->SetOwner(1);
75
76   TRandom random;
77   TDatime *datime = new TDatime();
78   Int_t time = datime->GetTime();
79   Int_t date = datime->GetDate();
80   Int_t pid  = gSystem->GetPid();
81   delete datime;
82   Int_t iseed = TMath::Abs(10000 * pid + time - date); 
83
84   //Float_t thrHVv=7.75, thrHVc=3, thrLVv=2.75, thrLVc=2.5,
85   //thrFEEthr=1.5, thrFEEt=10, thrTemp=35, thrPress=1000;
86   //Float_t tentHVv=6.75, tentHVc=2, tentLVv=1.75, tentLVc=1.5,
87   //  tentFEEthr=0.5, te    result=0;
88   //ntFEEt=20, tentTemp=25, tentPress=900;
89   //Float_t sigmaHVv=1, sigmaHVc=0.25, sigmaLVv=0.25, sigmaLVc=0.25,
90   //  sigmaFEEthr=0.05, sigmaFEEt=5, sigmaTemp=1, sigmaPress=10;
91
92   Float_t tentHVv=6500, tentHVi=80, tentLVv=2.7, tentLVi=4.5,
93     tentLVv33=3.3, tentLVv50=5.0, tentLVv48=48,
94     tentLVi33=100, tentLVi50=3.0, tentLVi48=10,
95     tentFEEthr=1.0, tentFEEtfeac=25, tentFEEttrm=45, 
96     tentTemp=25, tentPress=900;
97   Float_t sigmaHVv=10, sigmaHVi=10, sigmaLVv=0.2, sigmaLVi=1.0,
98     sigmaLVv33=0.1, sigmaLVv50=0.1, sigmaLVv48=1,
99     sigmaLVi33=10, sigmaLVi50=0.5, sigmaLVi48=2,
100     sigmaFEEthr=0.1, sigmaFEEtfeac=10, sigmaFEEttrm=4, 
101     sigmaTemp=1, sigmaPress=10;
102
103   Float_t tent=0, sigma=0, thr=0;
104   Int_t NAliases=10514, NHV=90, NLV=576, NLV33=72, NLV50=72, NLV48=72, NFEEthr=1152, NFEEtfeac=576, NFEEttrm=6840, NT=1, NP=1;
105
106   for(int nAlias=0;nAlias<NAliases;nAlias++) {
107
108     TObjArray* valueSet = new TObjArray;
109     valueSet->SetOwner(1);
110
111     TString sindex;
112     TString aliasName;
113     if (nAlias<NHV){
114       aliasName = "tof_hv_vp_";
115       sindex.Form("%02i",nAlias);
116       aliasName += sindex;
117       //aliasName += nAlias;
118       tent=tentHVv;
119       sigma=sigmaHVv;
120       //      thr=thrHVv;
121     }
122     else if (nAlias<NHV*2){
123       //      aliasName = "HVvneg";
124       //aliasName += nAlias-NHV;
125       aliasName = "tof_hv_vn_";
126       sindex.Form("%02i",nAlias-NHV);
127       aliasName += sindex;
128       tent=-tentHVv;
129       sigma=-sigmaHVv;
130       //thr=-thrHVv;
131     }
132     else if (nAlias<NHV*3){
133       //      aliasName = "HVcpos";
134       //aliasName += nAlias-2*NHV;
135       aliasName = "tof_hv_ip_";
136       sindex.Form("%02i",nAlias-2*NHV);
137       aliasName += sindex;
138       tent=tentHVi;
139       sigma=sigmaHVi;
140       //thr=thrHVc;
141     }
142     else if (nAlias<NHV*4){
143       //      aliasName = "HVcneg";
144       //aliasName += nAlias-3*NHV;
145       aliasName = "tof_hv_in_";
146       sindex.Form("%02i",nAlias-3*NHV);
147       aliasName += sindex;
148       tent=-tentHVi;
149       sigma=-sigmaHVi;
150       //thr=-thrHVc;
151     }
152     else if (nAlias<NHV*4+NLV){
153       //      aliasName = "LVv";
154       //aliasName += nAlias-4*NHV;
155       aliasName = "tof_lv_vfea_";
156       sindex.Form("%03i",nAlias-4*NHV);
157       aliasName += sindex;
158       tent=tentLVv;
159       sigma=sigmaLVv;
160       //thr=thrLVv;
161     }
162     else if (nAlias<NHV*4+2*NLV){
163       //      aliasName = "LVc";
164       //aliasName += nAlias-(4*NHV+NLV);
165       aliasName = "tof_lv_ifea_";
166       sindex.Form("%03i",nAlias-(4*NHV+NLV));
167       aliasName += sindex;
168       tent=tentLVi;
169       sigma=sigmaLVi;
170       //thr=thrLVc;
171     }
172     else if (nAlias<NHV*4+2*NLV+NLV33){
173       //      aliasName = "LVc";
174       //aliasName += nAlias-(4*NHV+NLV);
175       aliasName = "tof_lv_v33_";
176       sindex.Form("%02i",nAlias-(4*NHV+2*NLV));
177       aliasName += sindex;
178       tent=tentLVv33;
179       sigma=sigmaLVv33;
180       //thr=thrLVc;
181     }
182     else if (nAlias<NHV*4+2*NLV+2*NLV33){
183       //      aliasName = "LVc";
184       //aliasName += nAlias-(4*NHV+NLV);
185       aliasName = "tof_lv_i33_";
186       sindex.Form("%02i",nAlias-(4*NHV+2*NLV+NLV33));
187       aliasName += sindex;
188       tent=tentLVi33;
189       sigma=sigmaLVi33;
190       //thr=thrLVc;
191     }
192     else if (nAlias<NHV*4+2*NLV+2*NLV33+NLV50){
193       //      aliasName = "LVc";
194       //aliasName += nAlias-(4*NHV+NLV);
195       aliasName = "tof_lv_v50_";
196       sindex.Form("%02i",nAlias-(4*NHV+2*NLV+2*NLV33));
197       aliasName += sindex;
198       tent=tentLVv50;
199       sigma=sigmaLVv50;
200       //thr=thrLVc;
201     }
202     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50){
203       //      aliasName = "LVc";
204       //aliasName += nAlias-(4*NHV+NLV);
205       aliasName = "tof_lv_i50_";
206       sindex.Form("%02i",nAlias-(4*NHV+2*NLV+2*NLV33+NLV50));
207       aliasName += sindex;
208       tent=tentLVi50;
209       sigma=sigmaLVi50;
210       //thr=thrLVc;
211     }
212     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50+NLV48){
213       //      aliasName = "LVc";
214       //aliasName += nAlias-(4*NHV+NLV);
215       aliasName = "tof_lv_v48_";
216       sindex.Form("%02i",nAlias-(4*NHV+2*NLV+2*NLV33+2*NLV50));
217       aliasName += sindex;
218       tent=tentLVv48;
219       sigma=sigmaLVv48;
220       //thr=thrLVc;
221     }
222     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50+2*NLV48){
223       //      aliasName = "LVc";
224       //aliasName += nAlias-(4*NHV+NLV);
225       aliasName = "tof_lv_i48_";
226       sindex.Form("%02i",nAlias-(4*NHV+2*NLV+2*NLV33+2*NLV50+NLV48));
227       aliasName += sindex;
228       tent=tentLVi48;
229       sigma=sigmaLVi48;
230       //thr=thrLVc;
231     }
232     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50+2*NLV48+NFEEthr){
233       //      aliasName = "FEEthr";
234       //aliasName += nAlias-(4*NHV+2*NLV-(4*NHV+2*NLV+2*NLV33+2*NLV50+2*NLV48));
235       aliasName = "tof_fee_th_";
236       sindex.Form("%04i",nAlias-(4*NHV+2*NLV+2*NLV33+2*NLV50+2*NLV48));
237       aliasName += sindex;
238       tent=tentFEEthr;
239       sigma=sigmaFEEthr;
240       //thr=thrFEEthr;
241     }
242     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50+2*NLV48+NFEEthr+NFEEtfeac){
243       //cout << " nalias fee temp = " << nAlias << endl;
244       //      aliasName = "FEEt";
245       //aliasName += nAlias-(4*NHV+2*NLV+NFEEthr);
246       aliasName = "tof_fee_tfeac_";
247       sindex.Form("%03i",nAlias-(4*NHV+2*NLV+2*NLV33+2*NLV50+2*NLV48+NFEEthr));
248       aliasName += sindex;
249       //cout << " nalias fee temp name = " << aliasName << endl;
250       tent=tentFEEtfeac;
251       sigma=sigmaFEEtfeac;
252       //thr=thrFEEthr;
253     }
254     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50+2*NLV48+NFEEthr+NFEEtfeac+NFEEttrm){
255       //cout << " nalias fee temp = " << nAlias << endl;
256       //      aliasName = "FEEt";
257       //aliasName += nAlias-(4*NHV+2*NLV+NFEEthr);
258       aliasName = "tof_fee_ttrm_";
259       sindex.Form("%04i",nAlias-(4*NHV+2*NLV+2*NLV33+2*NLV50+2*NLV48+NFEEthr+NFEEtfeac));
260       aliasName += sindex;
261       //cout << " nalias fee temp name = " << aliasName << endl;
262       tent=tentFEEttrm;
263       sigma=sigmaFEEttrm;
264       //thr=thrFEEthr;
265     }
266     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50+2*NLV48+NFEEthr+NFEEtfeac+NFEEttrm+1){
267       cout << " putting temperature" << endl;
268       aliasName = "temperature";
269       tent=tentTemp;
270       sigma=sigmaTemp;
271       //thr=thrTemp;
272     }
273     else if (nAlias<NHV*4+2*NLV+2*NLV33+2*NLV50+2*NLV48+NFEEthr+NFEEtfeac+NFEEttrm+2){
274       cout << " putting pressure" << endl;
275       aliasName = "pressure";
276       tent=tentPress;
277       sigma=sigmaPress;
278       //thr=thrPress;
279     }
280
281
282     // gauss generation of values 
283     for (int timeStamp=0;timeStamp<1000;timeStamp+=10){
284       Float_t gaussvalue = (Float_t) (random.Gaus(tent,sigma));
285       if (TMath::Abs(gaussvalue-tent)>sigma){
286         AliDCSValue* dcsVal = new AliDCSValue(gaussvalue, timeStamp);
287         valueSet->Add(dcsVal);
288       }
289     }
290
291     aliasMap->Add(new TObjString(aliasName), valueSet);
292   }
293
294   return aliasMap;
295 }
296
297 TMap* ReadDCSAliasMap()
298 {
299   // Open a file that contains DCS input data
300   // The CDB framework is used to open the file, this means the file is located
301   // in $ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB/<detector>/DCS/Data
302   // The file contains an AliCDBEntry that contains a TMap with the DCS structure.
303   // An explanation of the structure can be found in CreateDCSAliasMap()
304
305   AliCDBEntry *entry = AliCDBManager::Instance()->Get("TOF/DCS/Data", 0);
306   return dynamic_cast<TMap*> (entry->GetObject());
307 }
308
309 void WriteDCSAliasMap()
310 {
311   // This writes the output from CreateDCSAliasMap to a CDB file
312
313   TMap* dcsAliasMap = CreateDCSAliasMap();
314
315   AliCDBMetaData metaData;
316         metaData.SetBeamPeriod(0);
317         metaData.SetResponsible("Chiara");
318         metaData.SetComment("Test object for TOFPreprocessor.C");
319
320   AliCDBId id("TOF/DCS/Data", 0, 0);
321
322   // initialize location of CDB
323   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB");
324
325   AliCDBManager::Instance()->Put(dcsAliasMap, id, &metaData);
326 }