Code cleanup (F. Prino)
[u/mrichter/AliRoot.git] / ITS / AliITSPreprocessorSDD.cxx
1 /**************************************************************************
2  * Copyright(c) 2007-2009, 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 // Implementation of the class for SDD preprocessing             //
21 // Origin: E.Crescio, Torino, crescio@to.infn.it                 //
22 //         F.Prino, Torino, prino@to.infn.it                     //
23 //                                                               //
24 ///////////////////////////////////////////////////////////////////
25
26 #include "AliITSPreprocessorSDD.h"
27 #include "AliITSCalibrationSDD.h"
28 #include "AliITSDriftSpeedSDD.h"
29 #include "AliITSDriftSpeedArraySDD.h"
30 #include "AliITSDCSAnalyzerSDD.h"
31 #include "AliShuttleInterface.h"
32 #include "AliCDBEntry.h"
33 #include "AliCDBMetaData.h"
34 #include "TObjArray.h"
35 #include "AliLog.h"
36 #include <TObjString.h>
37 #include <TSystem.h>
38 #include <TList.h>
39
40 const TString AliITSPreprocessorSDD::fgkNameHistoPedestals = "hpedestal";
41 const TString AliITSPreprocessorSDD::fgkNameHistoNoise = "hnoise";
42 ClassImp(AliITSPreprocessorSDD)
43
44
45 //______________________________________________________________________
46 UInt_t AliITSPreprocessorSDD::Process(TMap* dcsAliasMap){
47
48   // Get DDL map from OCDB
49   AliCDBEntry* entry = GetFromOCDB("Calib", "DDLMapSDD");
50   if(!entry){
51     Log("DDL map file not found in OCDB.");  
52     return 2;
53   }
54   AliITSDDLModuleMapSDD* ddlmap = (AliITSDDLModuleMapSDD*)entry->GetObject();
55   if(!ddlmap){ 
56     Log("AliITSDDLModuleMapSDD object not in file.");
57     return 2;
58   }
59
60   //preprocessing
61   TString runType = GetRunType();
62   Bool_t retcode;
63
64   if (runType == "PULSER"){
65     retcode=ProcessPulser(ddlmap);
66   }else if(runType== "INJECTOR"){
67     retcode=ProcessInjector(ddlmap);
68   }else{
69     // do nothing for other run types
70     retcode=1;
71   }
72   if(retcode){
73     Bool_t retcodedcs =ProcessDCSDataPoints(dcsAliasMap);
74     if(retcodedcs) return 0;
75   }
76   return 1;
77 }
78 //______________________________________________________________________
79 Bool_t AliITSPreprocessorSDD::ProcessPulser(AliITSDDLModuleMapSDD* ddlmap){
80   // Process FXS files from PULSER run (baseline, noise, gain)
81   TObjArray calSDD(kNumberOfSDD);
82   calSDD.SetOwner(kFALSE);
83   Char_t command[100];
84   Char_t inpFileName[100];
85   Float_t baseline,rawnoise,cmn,corn,gain;
86   Int_t isgoodan,i,im,is,isgoodmod,basmin,basoff;
87   Int_t numOfBadChannels[kNumberOfSDD];
88   
89   TList* sourceList = GetFileSources(kDAQ, "SDD_Calib");
90   if (!sourceList){ 
91     Log("Error: no sources found for SDD_Calib");
92     return 2;
93   }
94
95   Int_t ind = 0;
96   while (sourceList->At(ind)!=NULL) {
97     TObjString* tarId = (TObjString*) sourceList->At(ind);
98     TString tarName = GetFile(kDAQ, "SDD_Calib", tarId->GetString().Data());
99     if(tarName.Length()==0){
100       Log(Form("Baseline tar file from source %d not found.",ind));
101       return 2;
102     }
103     sprintf(command,"tar -xf %s",tarName.Data());
104     gSystem->Exec(command);
105     ind++;
106   }
107   delete sourceList;
108   
109   for(Int_t iddl=0;iddl<kNumberOfDDL;iddl++){
110     for(Int_t imod=0;imod<kModulesPerDDL;imod++){
111       Int_t modID=ddlmap->GetModuleNumber(iddl,imod);
112       if(modID==-1) continue;
113       modID-=240; // to have SDD modules numbering from 0 to 260
114       AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");
115       numOfBadChannels[modID]=0;
116       Int_t badch[kNumberOfChannels];
117       for(Int_t isid=0;isid<=1;isid++){
118         sprintf(inpFileName,"./SDDbase_ddl%02dc%02d_sid%d.data",iddl,imod,isid);
119         FILE* basFil = fopen(inpFileName,"read");
120         if (basFil == 0) {
121           Log(Form("File %s not found.",inpFileName));
122           cal->SetBad();
123           continue;
124         }
125         fscanf(basFil,"%d %d %d\n",&im,&is,&isgoodmod);
126         if(!isgoodmod) cal->SetBad();
127         for(Int_t ian=0;ian<(kNumberOfChannels/2);ian++){
128           fscanf(basFil,"%d %d %f %d %d %f %f %f %f\n",&i,&isgoodan,&baseline,&basmin,&basoff,&rawnoise,&cmn,&corn,&gain);
129           Int_t ich=ian;
130           if(isid==1) ich+=256;
131           if(!isgoodan){ 
132             Int_t ibad=numOfBadChannels[modID];
133             badch[ibad]=ich;
134             numOfBadChannels[modID]++;
135           }
136           cal->SetBaseline(ich,baseline);
137           cal->SetNoiseAfterElectronics(ich,rawnoise);
138           cal->SetGain(ich,gain);
139         }
140         cal->SetDeadChannels(numOfBadChannels[modID]);
141         for(Int_t ibad=0;ibad<numOfBadChannels[modID];ibad++){
142           cal->SetBadChannel(ibad,badch[ibad]);
143         }
144         fclose(basFil);
145       }
146       calSDD.AddAt(cal,modID);
147     }
148   }
149   AliCDBMetaData *md= new AliCDBMetaData();
150   md->SetResponsible("Francesco Prino");
151   md->SetBeamPeriod(0);
152   md->SetComment("AliITSCalibrationSDD from PEDESTAL+PULSER runs");
153   Bool_t retCode = Store("Calib","CalibSDD",&calSDD,md, 0, kTRUE);
154   return retCode;
155 }
156 //______________________________________________________________________
157 Bool_t AliITSPreprocessorSDD::ProcessInjector(AliITSDDLModuleMapSDD* ddlmap){
158   // Process FXS files from injector events (INJECTOR or PHYSICS runs)
159   TObjArray vdrift(2*kNumberOfSDD);
160   vdrift.SetOwner(kFALSE);
161   Char_t command[100];
162   Char_t inpFileName[100];
163   Int_t evNumb,polDeg; 
164   UInt_t timeStamp;
165   Float_t param[4];
166
167   TList* sourceList = GetFileSources(kDAQ, "SDD_Injec");
168   if (!sourceList){ 
169     Log("Error: no sources found for SDD_Injec");
170     return 2;
171   }
172   Int_t ind = 0;
173   while (sourceList->At(ind)!=NULL) {
174     TObjString* tarId = (TObjString*) sourceList->At(ind);
175     TString tarName = GetFile(kDAQ, "SDD_Injec", tarId->GetString().Data());
176     if(tarName.Length()==0){
177       Log(Form("Injector tar file from source %d not found.",ind));
178       return 2;
179     }
180     sprintf(command,"tar -xf %s",tarName.Data());
181     gSystem->Exec(command);
182     ind++;
183   }
184   delete sourceList;
185   
186   for(Int_t iddl=0;iddl<kNumberOfDDL;iddl++){
187     for(Int_t imod=0;imod<kModulesPerDDL;imod++){
188       Int_t modID=ddlmap->GetModuleNumber(iddl,imod);
189       if(modID==-1) continue;
190       modID-=240; // to have SDD modules numbering from 0 to 260
191       for(Int_t isid=0;isid<=1;isid++){
192         AliITSDriftSpeedArraySDD *arr=new AliITSDriftSpeedArraySDD();
193         sprintf(inpFileName,"./SDDinj_ddl%02dc%02d_sid%d.data",iddl,imod,isid);
194         FILE* injFil = fopen(inpFileName,"read");
195         if (injFil == 0) {
196           Log(Form("File %s not found.",inpFileName));
197           AliITSDriftSpeedSDD *dsp=new AliITSDriftSpeedSDD();
198           arr->AddDriftSpeed(dsp);
199           vdrift.AddAt(arr,2*modID+isid);
200           continue;
201         }
202         fscanf(injFil,"%d",&polDeg);
203         while (!feof(injFil)){
204           fscanf(injFil,"%d %d",&evNumb,&timeStamp);
205           if(feof(injFil)) break;
206           for(Int_t ic=0;ic<4;ic++) fscanf(injFil,"%f",&param[ic]);
207           AliITSDriftSpeedSDD *dsp=new AliITSDriftSpeedSDD(evNumb,timeStamp,polDeg,param);
208           arr->AddDriftSpeed(dsp);
209         }
210         vdrift.AddAt(arr,2*modID+isid);
211       }
212     }
213   }
214   AliCDBMetaData *md= new AliCDBMetaData();
215   md->SetResponsible("Francesco Prino");
216   md->SetBeamPeriod(0);
217   md->SetComment("AliITSDriftSpeedSDD from injector events");
218   Bool_t retCode = Store("Calib","DriftSpeedSDD",&vdrift,md,0, kTRUE);
219   return retCode;
220 }
221 //______________________________________________________________________
222 Bool_t AliITSPreprocessorSDD::ProcessDCSDataPoints(TMap* dcsAliasMap){
223   // Process DCS data
224   AliITSDCSAnalyzerSDD *dcs=new AliITSDCSAnalyzerSDD();
225   dcs->AnalyzeData(dcsAliasMap);
226   TObjArray refDCS(kNumberOfSDD);
227   refDCS.SetOwner(kFALSE);
228   for(Int_t imod=0;imod<kNumberOfSDD;imod++){
229     AliITSDCSDataSDD *dcsdata=dcs->GetDCSData(imod);
230     refDCS.Add(dcsdata);
231   }    
232   AliCDBMetaData *md= new AliCDBMetaData();
233   md->SetResponsible("Francesco Prino");
234   md->SetBeamPeriod(0);
235   md->SetComment("AliITSDCSDataSDD objects from DCS DB");
236   Bool_t retCode = StoreReferenceData("DCS","DataSDD",&refDCS,md);
237   return retCode;
238 }