GRP/TestGRPPreprocessor.C

   1
   2 // This macro runs the GRP test preprocessor
   3 // It uses AliTestShuttle to simulate a full Shuttle process
   4
   5 // The input data is created in the functions
   6 //   CreateDCSAliasMap() creates input that would in the same way come from DCS
   7 //   ReadDCSAliasMap() reads from a file
   8 //   CreateInputFilesMap() creates a list of local files, that can be accessed by the shuttle
   9
  10 #include <iostream>
  11 #include <fstream>
  12 using namespace std;
  13
  14 void TestGRPPreprocessor()
  15 {
  16   // load library
  17 //  gSystem->Load("libSTEER.so");        // needed for AliGRPPreprocessor
  18   gSystem->Load("$ALICE_ROOT/SHUTTLE/TestShuttle/libTestShuttle.so");
  19
  20   // create AliTestShuttle instance
  21   // The parameters are run, startTime, endTime
  22   Int_t kRun = 7;
  23   AliTestShuttle* shuttle = new AliTestShuttle(kRun, 1, 10);
  24
  25   // by default they are set to $ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB and TestReference
  26   AliTestShuttle::SetMainCDB("local://$ALICE_ROOT/OCDB/SHUTTLE/TestShuttle/TestCDB");
  27   AliTestShuttle::SetMainRefStorage("local://$ALICE_ROOT/OCDB/SHUTTLE/TestShuttle/TestReference");
  28
  29   printf("Test OCDB storage Uri: %s\n", AliShuttleInterface::GetMainCDB().Data());
  30   printf("Test Reference storage Uri: %s\n", AliShuttleInterface::GetMainRefStorage().Data());
  31
  32
  33   // The shuttle can read DCS data, if the preprocessor should be tested to process DCS data,
  34   // some fake data has to be created.
  35   //
  36   // The "fake" input data can be taken using either (a) or (b):
  37   // (a) data from a file: Use ReadDCSAliasMap()
  38   //     the format of the file is explained in ReadDCSAliasMap()
  39   //     To use it uncomment the following line:
  40   //
  41   //TMap* dcsAliasMap = ReadDCSAliasMap();
  42   //
  43   // (b) generated in this macro: Use CreateDCSAliasMap() and its documentation
  44   //     To use it uncomment the following line:
  45   //
  46   TMap* dcsAliasMap = CreateDCSAliasMap();
  47
  48   // now give the alias map to the shuttle
  49   shuttle->SetDCSInput(dcsAliasMap);
  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   //shuttle->AddInputFile(AliShuttleInterface::kDAQ, "GRP", "Period_test000.tag.root", "GDC0", "runTags01.root");
  64   //shuttle->AddInputFile(AliShuttleInterface::kDAQ, "GRP", "Period_test001.tag.root", "GDC1", "runTags02.root");
  65   //shuttle->AddInputFile(AliShuttleInterface::kDAQ, "GRP", "Period_test002.tag.root", "GDC2", "runTags03.root");
  66
  67   shuttle->AddInputFile(AliShuttleInterface::kDAQ, "GRP", "run000021390_GRP_gdc0_Period_TDSMtest.Seq_0.tag.root", "GDC0", "$ALICE_ROOT/GRP/ShuttleInput/run000021390_GRP_gdc0_Period_TDSMtest.Seq_0.tag.root");
  68   shuttle->AddInputFile(AliShuttleInterface::kDAQ, "GRP", "run000021390_GRP_gdc1_Period_TDSMtest.Seq_0.tag.root", "GDC1", "$ALICE_ROOT/GRP/ShuttleInput/run000021390_GRP_gdc0_Period_TDSMtest.Seq_0.tag.root");
  69   shuttle->AddInputFile(AliShuttleInterface::kDAQ, "GRP", "run000021391_GRP_gdc0_Period_TDSMtest.Seq_0.tag.root", "GDC0", "$ALICE_ROOT/GRP/ShuttleInput/run000021391_GRP_gdc0_Period_TDSMtest.Seq_0.tag.root");
  70   shuttle->AddInputFile(AliShuttleInterface::kDAQ, "GRP", "run000021391_GRP_gdc1_Period_TDSMtest.Seq_0.tag.root", "GDC1", "$ALICE_ROOT/GRP/ShuttleInput/run000021391_GRP_gdc1_Period_TDSMtest.Seq_0.tag.root");
  71
  72   shuttle->AddInputFile(AliShuttleInterface::kDCS, "GRP", "CTP_runconfig", "", gSystem->ExpandPathName("$ALICE_ROOT/GRP/CTP/p-p.cfg"));
  73   shuttle->AddInputFile(AliShuttleInterface::kDCS, "GRP", "CTP_xcounters", "DCS FXS", gSystem->ExpandPathName("$ALICE_ROOT/GRP/CTP/xcounters.txt"));
  74
  75   Char_t * filename = gSystem->ExpandPathName("$ALICE_ROOT/GRP/CTP/p-p.cfg");
  76   ifstream is;
  77   is.open(filename);
  78   is.seekg(0,ios::end);
  79   int length = is.tellg();
  80   const char *buffer = new char[length];
  81   is.seekg(0,ios::beg);
  82   is.read(buffer,length);
  83   is.close();
  84
  85   shuttle->SetInputTriggerConfiguration(buffer);
  86
  87   // The shuttle can read run type stored in the DAQ logbook.
  88   // To test it, we must provide the run type manually. They will be retrieved in the preprocessor
  89   // using GetRunType function.
  90   shuttle->SetInputRunType("PHYSICS");
  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   shuttle->AddInputRunParameter("time_start", "1233213.22");
  96   shuttle->AddInputRunParameter("time_end",   "1345645.22");
  97   shuttle->AddInputRunParameter("beamEnergy", "1400.");
  98   shuttle->AddInputRunParameter("beamType",    "p-p");
  99   shuttle->AddInputRunParameter("numberOfDetectors", "5");
 100   shuttle->AddInputRunParameter("detectorMask", "34555");
 101   shuttle->AddInputRunParameter("LHCperiod",    "LHC08b");
 102
 103 //  shuttle->AddInputRunParameter("totalEvents", "30000");
 104 //  shuttle->AddInputRunParameter("NumberOfGDCs", "15");
 105
 106   // TODO(5) NOT NEEDED
 107   //
 108   // This is for preprocessor that require data from HLT.
 109   // Since HLT may be switched off, the preprocessor should first query the Run logbook where
 110   // the HLT status is stored. SHUTTLE implements a shortcut function (GetHLTStatus) that returns
 111   // a bool directly. 1 = HLT ON, 0 = HLT OFF
 112   //
 113
 114   Bool_t hltStatus = kTRUE;
 115   shuttle->SetInputHLTStatus(hltStatus);
 116
 117
 118   // Create the preprocessor that should be tested, it registers itself automatically to the shuttle
 119   AliPreprocessor* test = new AliGRPPreprocessor(shuttle);
 120
 121   // Test the preprocessor
 122   shuttle->Process();
 123
 124   printf("\n\n");
 125
 126   // In the preprocessor AliShuttleInterface::Store should be called to put the final
 127   // data to the CDB. To check if all went fine have a look at the files produced in
 128   // $ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB/<detector>/SHUTTLE/Data
 129   //
 130   // Check the file which should have been created
 131   AliCDBEntry* chkEntry = AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
 132                         ->Get("GRP/GRP/Data", kRun);
 133   if (!chkEntry) {
 134     printf("The file is not there. Something went wrong.\n");
 135     return;
 136   }
 137   chkEntry->PrintId();
 138   chkEntry->GetObject()->Print();
 139   printf("\n\n");
 140
 141   AliDCSSensor* sen = (AliDCSSensor*)(((TMap*)chkEntry->GetObject())->GetValue("fP2Pressure"));
 142   if(sen)
 143     sen->GetFit()->MakeGraph(1.5, 9.5, 15)->Draw();
 144
 145 //  AliTestDataDCS* output = dynamic_cast<AliTestDataDCS*> (chkEntry->GetObject());
 146   // If everything went fine, draw the result
 147 //  if (output) {
 148 //    output->Print();
 149 //  }
 150
 151
 152 }
 153
 154 TMap* CreateDCSAliasMap()
 155 {
 156   // Creates a DCS structure
 157   // The structure is the following:
 158   //   TMap (key --> value)
 159   //     <DCSAlias> --> <valueList>
 160   //     <DCSAlias> is a string
 161   //     <valueList> is a TObjArray of AliDCSValue
 162   //     An AliDCSValue consists of timestamp and a value in form of a AliSimpleValue
 163
 164   const Int_t fgknDCSDP = 10;
 165   const char* fgkDCSDataPoints[AliGRPPreprocessor::fgknDCSDP] = {
 166                    "LHCState",
 167                    "L3Polarity",
 168                    "DipolePolarity",
 169                    "LHCLuminosity",
 170                    "BeamIntensity",
 171                    "L3Current",
 172                    "DipoleCurrent",
 173                    "CavernTemperature",
 174                    "CavernAtmosPressure",
 175                    "SurfaceAtmosPressure"
 176                  };
 177
 178   TMap* aliasMap;
 179   TObjArray* valueSet;
 180   AliDCSValue* dcsVal;
 181
 182   aliasMap = new TMap;
 183   aliasMap->SetOwner(1);
 184
 185   // LHCState
 186   valueSet = new TObjArray;
 187   valueSet->SetOwner(1);
 188   dcsVal = new AliDCSValue( 'F', 0 );
 189   valueSet->Add(dcsVal);
 190 //  aliasMap->Add( new TObjString(fgkDCSDataPoints[0]), valueSet );
 191
 192   // L3Polarity
 193   valueSet = new TObjArray;
 194   valueSet->SetOwner(1);
 195   dcsVal = new AliDCSValue( kTRUE, 0 );
 196   valueSet->Add(dcsVal);
 197   aliasMap->Add( new TObjString(fgkDCSDataPoints[1]), valueSet );
 198
 199   // DipolePolarity
 200   valueSet = new TObjArray;
 201   valueSet->SetOwner(1);
 202   dcsVal = new AliDCSValue( kTRUE, 0 );
 203   valueSet->Add(dcsVal);
 204   aliasMap->Add( new TObjString(fgkDCSDataPoints[2]), valueSet );
 205
 206   TRandom random;
 207
 208   //  for( int nAlias=3; nAlias<fgknDCSDP-1; nAlias++)  {
 209   for( int nAlias=3; nAlias<fgknDCSDP; nAlias++)  {
 210     valueSet = new TObjArray;
 211     valueSet->SetOwner(1);
 212
 213     for (int timeStamp=0; timeStamp<100; timeStamp++) {
 214       dcsVal = new AliDCSValue((Float_t) (random.Gaus()+5*nAlias), timeStamp);
 215       //printf("%s\n",dcsVal->ToString().Data());
 216       valueSet->Add(dcsVal);
 217     }
 218     aliasMap->Add( new TObjString( fgkDCSDataPoints[nAlias]), valueSet );
 219   }
 220
 221   return aliasMap;
 222 }
 223
 224 /*
 225 TMap* ReadDCSAliasMap()
 226 {
 227   // Open a file that contains DCS input data
 228   // The CDB framework is used to open the file, this means the file is located
 229   // in $ALICE_ROOT/SHUTTLE/TestShuttle/TestCDB/<detector>/DCS/Data
 230   // The file contains an AliCDBEntry that contains a TMap with the DCS structure.
 231   // An explanation of the structure can be found in CreateDCSAliasMap()
 232
 233   AliCDBEntry *entry = AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
 234                         ->Get("DET/DCS/Data", 0);
 235   return dynamic_cast<TMap*> (entry->GetObject());
 236 }
 237 */
 238
 239
 240 void WriteDCSAliasMap()
 241 {
 242   // This writes the output from CreateDCSAliasMap to a CDB file
 243
 244   TMap* dcsAliasMap = CreateDCSAliasMap();
 245
 246   AliCDBMetaData metaData;
 247   metaData.SetBeamPeriod(0);
 248   metaData.SetResponsible("Ernesto Lopez Torres");
 249   metaData.SetComment("Test object for TestGRPPreprocessor.C");
 250
 251   AliCDBId id("GRP/Data", 0, 0);
 252
 253   // look into AliTestShuttle's CDB main folder
 254
 255   AliCDBManager::Instance()->GetStorage(AliShuttleInterface::GetMainCDB())
 256                         ->Put(dcsAliasMap, id, &metaData);
 257 }