ZDC/TestZDCPreprocessor.C

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