Adding a reminder for coders
[u/mrichter/AliRoot.git] / ZDC / TestZDCPreprocessor.C
CommitLineData
438a188c 1/* $Id$ */
2
3// This class runs the test preprocessor
4// It uses AliTestShuttle to simulate a full Shuttle process
5
6// The input data is created in the functions
7// CreateDCSAliasMap() creates input that would in the same way come from DCS
8// ReadDCSAliasMap() reads from a file
9// CreateInputFilesMap() creates a list of local files, that can be accessed by the shuttle
10
bf867209 11void TestZDCPreprocessor(const char* runType="PHYSICS")
438a188c 12{
13 // load library
14 gSystem->Load("libTestShuttle.so");
15
16 // create AliTestShuttle instance
17 // The parameters are run, startTime, endTime
86dc0caa 18 AliTestShuttle* shuttle = new AliTestShuttle(0, 0, 1);
438a188c 19
438a188c 20 // TODO if needed, change location of OCDB and Reference test folders
21 // by default they are set to $ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB and TestReference
1a1c78d4 22 //AliTestShuttle::SetMainCDB("local://$ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB");
23 AliTestShuttle::SetMainCDB("local://$ALICE_ROOT/OCDB");
da2b6160 24 //AliTestShuttle::SetMainCDB("alien://folder=/alice/data/2009/OCDB/");
e83c8d88 25 AliTestShuttle::SetMainRefStorage("local://$ALICE_ROOT/SHUTTLE/TestShuttle/TestReference");
438a188c 26
da2b6160 27 printf("\n Test OCDB storage Uri: %s\n", AliShuttleInterface::GetMainCDB().Data());
28 printf(" Test Reference storage Uri: %s\n\n", AliShuttleInterface::GetMainRefStorage().Data());
438a188c 29
438a188c 30 // TODO(1)
31 //
32 // The shuttle can read DCS data, if the preprocessor should be tested to process DCS data,
33 // some fake data has to be created.
34 //
35 // The "fake" input data can be taken using either (a) or (b):
36 // (a) data from a file: Use ReadDCSAliasMap()
37 // the format of the file is explained in ReadDCSAliasMap()
38 // To use it uncomment the following line:
39 //
1a1c78d4 40 //TMap* dcsAliasMap = ReadDCSAliasMap();
438a188c 41 //
42 // (b) generated in this macro: Use CreateDCSAliasMap() and its documentation
43 // To use it uncomment the following line:
44 //
1a1c78d4 45 TMap* dcsAliasMap = CreateDCSAliasMap();
5594ae91 46 dcsAliasMap->Print("");
da2b6160 47 //WriteDCSAliasMap();
438a188c 48
49 // now give the alias map to the shuttle
50 shuttle->SetDCSInput(dcsAliasMap);
51
52 // TODO(2)
53 //
54 // The shuttle can also process files that originate from DCS, DAQ and HLT.
55 // To test it, we provide some local files and locations where these would be found when
56 // the online machinery would be there.
57 // In real life this functions would be produces by the sub-detectors
58 // calibration programs in DCS, DAQ or HLT. These files can then be retrieved using the Shuttle.
59 //
60 // Files are added with the function AliTestShuttle::AddInputFile. The syntax is:
61 // AddInputFile(<system>, <detector>, <id>, <source>, <local-file>)
62 // In this example we add a file originating from the GDC with the id PEDESTALS
63 // Three files originating from different LDCs but with the same id are also added
64 // Note that the test preprocessor name is TPC. The name of the detector's preprocessor must follow
65 // the "online" naming convention ALICE-INT-2003-039.
1ee299a5 66 //
da2b6160 67 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "PEDESTALDATA", "LDC", "ZDCPedestal.dat");
68 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "PEDESTALHISTOS", "LDC", "ZDCPedHisto.root");
321fe3b8 69 //
da2b6160 70 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "LASERDATA", "LDC", "ZDCLaserCalib.dat");
71 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "LASERHISTOS", "LDC", "ZDCLaserHisto.root");
321fe3b8 72 //
da2b6160 73 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "EMDENERGYCALIB", "LDC", "ZDCEnergyCalib.dat");
74 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "EMDTOWERCALIB", "LDC", "ZDCTowerCalib.dat");
321fe3b8 75 //
5594ae91 76 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "MBCALIB", "LDC", "ZDCMBCalib.root");
77 //
da2b6160 78 shuttle->AddInputFile(AliTestShuttle::kDAQ, "ZDC", "MAPPING", "MON", "ZDCChMapping.dat");
438a188c 79
80 // TODO(3)
81 //
82 // The shuttle can read run type stored in the DAQ logbook.
83 // To test it, we must provide the run type manually. They will be retrieved in the preprocessor
84 // using GetRunType function.
da2b6160 85 shuttle->SetInputRunType(runType);
86
438a188c 87 // TODO(4)
88 //
89 // The shuttle can read run parameters stored in the DAQ run logbook.
90 // To test it, we must provide the run parameters manually. They will be retrieved in the preprocessor
91 // using GetRunParameter function.
92 // In real life the parameters will be retrieved automatically from the run logbook;
c1ae51a3 93 shuttle->AddInputRunParameter("beamType", "A-A");
94 shuttle->AddInputRunParameter("beamEnergy", "1380");
95 //shuttle->AddInputRunParameter("beamType", "p-p");
438a188c 96 shuttle->AddInputRunParameter("totalEvents", "1000");
97 shuttle->AddInputRunParameter("NumberOfGDCs", "1");
98
952af594 99 // TODO(5) NEW!
100 //
101 // This is for preprocessor that require data from HLT.
102 // Since HLT may be switched off, the preprocessor should first query the Run logbook where
103 // the HLT status is stored. SHUTTLE implements a shortcut function (GetHLTStatus) that returns
104 // a bool directly. 1 = HLT ON, 0 = HLT OFF
105 //
106 Bool_t hltStatus=kFALSE;
d2317a93 107 shuttle->SetInputHLTStatus(hltStatus);
952af594 108
109 // TODO(6)
438a188c 110 //
111 // The shuttle can query condition parameters valid from the current run from the OCDB
112 // To test it, we must first store the object into the OCDB. It will be retrieved in the preprocessor
113 // using GetFromOCDB function.
1ee299a5 114
115 /*TObjString obj("This is a condition parameter stored in OCDB");
438a188c 116 AliCDBId id("ZDC/Calib/Data", 0, AliCDBRunRange::Infinity());
117 AliCDBMetaData md;
118 AliCDBEntry entry(&obj, id, &md);
119
120 shuttle->AddInputCDBEntry(&entry);
1ee299a5 121 */
122
438a188c 123 // TODO(6)
124 // Create the preprocessor that should be tested, it registers itself automatically to the shuttle
952af594 125 AliPreprocessor* test = new AliZDCPreprocessor(shuttle);
da2b6160 126 shuttle->Print();
127
438a188c 128 // Test the preprocessor
129 shuttle->Process();
1a1c78d4 130 //printf(" Back to test macro: final checks! \n");
438a188c 131
132 // TODO(7)
133 // In the preprocessor AliShuttleInterface::Store should be called to put the final
134 // data to the CDB. To check if all went fine have a look at the files produced in
135 // $ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB/<detector>/SHUTTLE/Data
136 //
137 // Check the file which should have been created
1ee299a5 138 AliCDBEntry* chkEntry0 = AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
dfa968c2 139 ->Get("ZDC/Calib/ChMap", 0);
da2b6160 140 TString str(runType);
86dc0caa 141 AliCDBEntry* chkEntry1;
da2b6160 142 if((str.CompareTo("STANDALONE_PEDESTAL")) == 0) chkEntry1 = AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
dfa968c2 143 ->Get("ZDC/Calib/Pedestals", 0);
da2b6160 144 else if((str.CompareTo("STANDALONE_LASER")) == 0) chkEntry1 = AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
dfa968c2 145 ->Get("ZDC/Calib/LaserCalib", 0);
da2b6160 146 else if((str.CompareTo("CALIBRATION_EMD")) == 0) chkEntry1 = AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
dfa968c2 147 ->Get("ZDC/Calib/EnergyCalib", 0);
1ee299a5 148
6024ec85 149
86dc0caa 150 if(!chkEntry0){
73bc3a3f 151 printf("No file in ZDC/Calib/ChMap\n");
438a188c 152 return;
153 }
86dc0caa 154 if(!chkEntry1){
da2b6160 155 if((str.CompareTo("STANDALONE_PEDESTAL")) == 0) printf("No file in ZDC/Calib/Pedestal\n");
156 else if((str.CompareTo("STANDALONE_LASER")) == 0) printf("No file in ZDC/Calib/LaserCalib\n");
157 else if((str.CompareTo("CALIBRATION_EMD")) == 0) printf("No file in ZDC/Calib/EnergyCalib\n");
86dc0caa 158 return;
159 }
160
438a188c 161
86dc0caa 162 /*AliTestDataDCS* output = dynamic_cast<AliTestDataDCS*> (chkEntry1->GetObject());
438a188c 163 // If everything went fine, draw the result
164 if (output)
86dc0caa 165 output->Draw();*/
438a188c 166}
167
168TMap* CreateDCSAliasMap()
169{
170 // Creates a DCS structure
171 // The structure is the following:
172 // TMap (key --> value)
173 // <DCSAlias> --> <valueList>
174 // <DCSAlias> is a string
175 // <valueList> is a TObjArray of AliDCSValue
176 // An AliDCSValue consists of timestamp and a value in form of a AliSimpleValue
177
178 // In this example 6 aliases exists: DCSAlias1 ... DCSAlias6
179 // Each contains 1000 values randomly generated by TRandom::Gaus + 5*nAlias
180
181 TMap* aliasMap = new TMap;
182 aliasMap->SetOwner(1);
183
184 TRandom random;
952af594 185 TString aliasNames[28];
438a188c 186
187 // ******************************** alignment values
188 aliasNames[0] = "ZDC_ZNA_POS.actual.position";
189 aliasNames[1] = "ZDC_ZPA_POS.actual.position";
190 aliasNames[2] = "ZDC_ZNC_POS.actual.position";
191 aliasNames[3] = "ZDC_ZPC_POS.actual.position";
192 //
57c85e6e 193 for(int nAlias=0; nAlias<4; nAlias++){
438a188c 194 TObjArray* valueSet = new TObjArray;
195 valueSet->SetOwner(1);
196
197 TString aliasName = aliasNames[nAlias];
73bc3a3f 198 //printf("\n\n alias: %s\n\n",aliasName.Data());
438a188c 199
200 Float_t simVal = (Float_t) (random.Rndm()*0.025+random.Rndm()*0.1);
57c85e6e 201 int timeStamp1[5] = {0,500,1000,1500,2000};
202 for(int i=0;i<5;i++){
438a188c 203 AliDCSValue* dcsVal = new AliDCSValue(simVal, timeStamp1[i]);
73bc3a3f 204 //printf("%s\n",dcsVal->ToString());
438a188c 205 valueSet->Add(dcsVal);
206 }
207 aliasMap->Add(new TObjString(aliasName), valueSet);
208 }
209 // ******************************** HV values
438a188c 210 aliasNames[4] = "ZDC_ZNA_HV0.actual.vMon";
211 aliasNames[5] = "ZDC_ZNA_HV1.actual.vMon";
212 aliasNames[6] = "ZDC_ZNA_HV2.actual.vMon";
213 aliasNames[7] = "ZDC_ZNA_HV3.actual.vMon";
214 aliasNames[8] = "ZDC_ZNA_HV4.actual.vMon";
215 //
216 aliasNames[9] = "ZDC_ZPA_HV0.actual.vMon";
217 aliasNames[10] = "ZDC_ZPA_HV1.actual.vMon";
218 aliasNames[11] = "ZDC_ZPA_HV2.actual.vMon";
219 aliasNames[12] = "ZDC_ZPA_HV3.actual.vMon";
220 aliasNames[13] = "ZDC_ZPA_HV4.actual.vMon";
221 //
222 aliasNames[14] = "ZDC_ZNC_HV0.actual.vMon";
223 aliasNames[15] = "ZDC_ZNC_HV1.actual.vMon";
224 aliasNames[16] = "ZDC_ZNC_HV2.actual.vMon";
225 aliasNames[17] = "ZDC_ZNC_HV3.actual.vMon";
226 aliasNames[18] = "ZDC_ZNC_HV4.actual.vMon";
227 //
228 aliasNames[19] = "ZDC_ZPC_HV0.actual.vMon";
229 aliasNames[20] = "ZDC_ZPC_HV1.actual.vMon";
230 aliasNames[21] = "ZDC_ZPC_HV2.actual.vMon";
231 aliasNames[22] = "ZDC_ZPC_HV3.actual.vMon";
232 aliasNames[23] = "ZDC_ZPC_HV4.actual.vMon";
233 //
952af594 234 aliasNames[24] = "ZDC_ZEM_HV0.actual.vMon";
235 aliasNames[25] = "ZDC_ZEM_HV1.actual.vMon";
236 //
20a3b76c 237 aliasNames[26] = "ZDC_REFA_HV.actual.vMon";
238 aliasNames[27] = "ZDC_REFC_HV.actual.vMon";
438a188c 239 //
57c85e6e 240 for(int nAlias=4;nAlias<28;nAlias++){
5594ae91 241// if(nAlias<14 || nAlias>18){
57c85e6e 242 TObjArray* valueSet = new TObjArray;
243 valueSet->SetOwner(1);
438a188c 244
57c85e6e 245 TString aliasName = aliasNames[nAlias];
1a1c78d4 246 //printf("\n\n alias: %s\n\n",aliasName.Data());
57c85e6e 247
248 for(int timeStamp=0;timeStamp<=2000;timeStamp+=500){
249 Float_t simVal = (Float_t) (random.Gaus()*600.+1800.);
250 AliDCSValue* dcsVal = new AliDCSValue(simVal, timeStamp);
251 //printf("%s\n",dcsVal->ToString());
252 valueSet->Add(dcsVal);
253 }
254 aliasMap->Add(new TObjString(aliasName), valueSet);
5594ae91 255// }
438a188c 256 }
257
258 return aliasMap;
259}
260
261TMap* ReadDCSAliasMap()
262{
263 // Open a file that contains DCS input data
264 // The CDB framework is used to open the file, this means the file is located
265 // in $ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB/<detector>/DCS/Data
266 // The file contains an AliCDBEntry that contains a TMap with the DCS structure.
267 // An explanation of the structure can be found in CreateDCSAliasMap()
268
bf867209 269 AliCDBManager *manager = AliCDBManager::Instance();
270 AliCDBStorage *sto = manager->GetStorage("local://$ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB/");
271 AliCDBId id("ZDC/DCS/Data",0,999999999);
57c85e6e 272 AliCDBEntry *entry = sto->Get("ZDC/DCS/Data", 0);
bf867209 273 if(!entry) printf("TestZDCPreprocessor.C -> ERROR! No entry found as DCS Map! \n");
438a188c 274 return dynamic_cast<TMap*> (entry->GetObject());
275}
276
277void WriteDCSAliasMap()
278{
279 // This writes the output from CreateDCSAliasMap to a CDB file
280
281 TMap* dcsAliasMap = CreateDCSAliasMap();
282
283 AliCDBMetaData metaData;
284 metaData.SetBeamPeriod(0);
d2317a93 285 metaData.SetResponsible("Chiara Oppedisano");
286 metaData.SetComment("Test object for TestZDCPreprocessor.C");
438a188c 287
bf867209 288 AliCDBId id("ZDC/DCS/Data", 0, 999999999);
438a188c 289
290 // initialize location of CDB
162637e4 291 AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB/SHUTTLE/TestShuttle/TestCDB");
438a188c 292
293 AliCDBManager::Instance()->Put(dcsAliasMap, id, &metaData);
294}