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