]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TPC/TPCPULSERda.cxx
Patch for the tracker
[u/mrichter/AliRoot.git] / TPC / TPCPULSERda.cxx
index 122b17a3839b7161a88fd251dbb485dcfb05ff22..73f146f9796de5f7a61b066a53790c2165f1af14 100644 (file)
@@ -3,7 +3,7 @@ TPC DA for online calibration
 
 Contact: Haavard.Helstrup@cern.ch
 Link:
-Run Type: PULSER_RUN
+Run Type: CALIBRATION_PULSER
 DA Type: LDC
 Number of events needed: 100
 Input Files: 
@@ -19,7 +19,8 @@ TPCda_pulser.cxx - calibration algorithm for TPC pulser events
 
 10/06/2007  sylvain.chapeland@cern.ch :  first version - clean skeleton based on DAQ DA case1
 30/09/2007  haavard.helstrup@cern.ch  :  created pulser DA based on pedestal code
-
+19/09/2008  J.Wiechula@gsi.de:            Added export of the calibration data to the AMORE data base.
+                                          Added support for configuration files.
 contact: marian.ivanov@cern.ch
 
 
@@ -31,6 +32,7 @@ and save results in a file (named from RESULT_FILE define - see below).
 #define RESULT_FILE "tpcPulser.root"
 #define FILE_ID "pulser"
 #define MAPPING_FILE "tpcMapping.root"
+#define CONFIG_FILE "TPCPULSERda.conf"
 #define AliDebugLevel() -1
 
 
@@ -46,6 +48,9 @@ and save results in a file (named from RESULT_FILE define - see below).
 #include <TFile.h>
 #include "TROOT.h"
 #include "TPluginManager.h"
+#include "TString.h"
+#include "TObjString.h"
+#include "TDatime.h"
 //
 //AliRoot includes
 //
@@ -60,42 +65,62 @@ and save results in a file (named from RESULT_FILE define - see below).
 #include "TTreeStream.h"
 #include "AliLog.h"
 #include "TSystem.h"
-
+#include "AliTPCConfigDA.h"
+//
+//AMORE
+//
+#include <AmoreDA.h>
 //
 // TPC calibration algorithm includes
 //
 #include "AliTPCCalibPulser.h"
 
-
-
-
 /* Main routine
       Arguments: list of DATE raw data files
 */
 int main(int argc, char **argv) {
+  /* log start of process */
+  printf("TPC Pulser DA started - %s\n",__FILE__);
+
+  if (argc<2) {
+    printf("Wrong number of arguments\n");
+    return -1;
+  }
   AliLog::SetClassDebugLevel("AliTPCRawStream",-5);
   AliLog::SetClassDebugLevel("AliRawReaderDate",-5);
   AliLog::SetClassDebugLevel("AliTPCAltroMapping",-5);
   AliLog::SetModuleDebugLevel("RAW",-5);
 
- gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
-                                         "*",
-                                         "TStreamerInfo",
-                                         "RIO",
-                                         "TStreamerInfo()");
 gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
+                                        "*",
+                                        "TStreamerInfo",
+                                        "RIO",
+                                        "TStreamerInfo()");
 
 
+  /* declare monitoring program */
   int i,status;
+  status=monitorDeclareMp( __FILE__ );
+  if (status!=0) {
+    printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
+    return -1;
+  }
+
+  // variables 
   AliTPCmapper *mapping = 0;   // The TPC mapping
+  char localfile[255];
+  unsigned long32 runNb=0;      //run number
+  // configuration options 
+  Bool_t fastDecoding = kFALSE;
   // if  test setup get parameters from $DAQDA_TEST_DIR 
    
   if (!mapping){
     /* copy locally the mapping file from daq detector config db */
-    status = daqDA_DB_getFile(MAPPING_FILE,"./tpcMapping.root");
+    sprintf(localfile,"./%s",MAPPING_FILE);
+    status = daqDA_DB_getFile(MAPPING_FILE,localfile);
     if (status) {
       printf("Failed to get mapping file (%s) from DAQdetDB, status=%d\n", MAPPING_FILE, status);
-      printf("Continue anyway ... maybe it works?\n");              // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-      return -1;   // temporarily uncommented for testing on pcald47 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+      return -1;
     }
 
     /* open the mapping file and retrieve mapping object */
@@ -106,35 +131,35 @@ int main(int argc, char **argv) {
 
   if (mapping == 0) {
     printf("Failed to get mapping object from %s.  ...\n", MAPPING_FILE);
-    //return -1;
+    return -1;
   } else {
     printf("Got mapping object from %s\n", MAPPING_FILE);
   }
 
-
-  AliTPCCalibPulser calibPulser;   // pedestal and noise calibration
-
-  if (argc<2) {
-    printf("Wrong number of arguments\n");
+  //
+  // DA configuration from configuration file
+  //
+ //retrieve configuration file
+  sprintf(localfile,"./%s",CONFIG_FILE);
+  status = daqDA_DB_getFile(CONFIG_FILE,localfile);
+  if (status) {
+    printf("Failed to get configuration file (%s) from DAQdetDB, status=%d\n", CONFIG_FILE, status);
     return -1;
   }
-
-
-  /* log start of process */
-  printf("TPC Pulser DA started - %s\n",__FILE__);
-
-  /* set time bin range */
-  calibPulser.SetRangeTime(400,500);
-  calibPulser.SetAltroMapping(mapping->GetAltroMapping()); // Use altro mapping we got from daqDetDb
-  /* declare monitoring program */
-  status=monitorDeclareMp( __FILE__ );
-  if (status!=0) {
-    printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
-    return -1;
+  AliTPCConfigDA config(CONFIG_FILE);
+  // check configuration options
+  if ( (Int_t)config.GetValue("UseFastDecoder") == 1 ) {
+    printf("Info: The fast decoder will be used for the processing.\n");
+    fastDecoding=kTRUE;
   }
 
+  // create calibration object
+  AliTPCCalibPulser calibPulser(config.GetConfigurationMap());   // pulser calibration algorithm
+  calibPulser.SetAltroMapping(mapping->GetAltroMapping()); // Use altro mapping we got from daqDetDb
 
-  /* loop over RAW data files */
+  //===========================//
+  // loop over RAW data files //
+  //==========================//
   int nevents=0;
   for(i=1;i<argc;i++) {
 
@@ -170,24 +195,24 @@ int main(int argc, char **argv) {
         continue;
       }
       nevents++;
-
+      // get the run number
+      runNb = event->eventRunNb;
       //  Pulser calibration
-
-      AliRawReader *rawReader = new AliRawReaderDate((void*)event);
-//      calibPulser.ProcessEvent(rawReader);
-      calibPulser.ProcessEventFast(rawReader);
-      delete rawReader;
+      calibPulser.ProcessEvent(event);
 
       /* free resources */
       free(event);
     }
   }
 
-  calibPulser.Analyse(); 
+  //
+  // Analyse pulser data and write them to rootfile
+  //
+  calibPulser.Analyse();
   printf ("%d events processed\n",nevents);
 
   TFile * fileTPC = new TFile (RESULT_FILE,"recreate");
-  calibPulser.Write("calibPulser");
+  calibPulser.Write("tpcCalibPulser");
   delete fileTPC;
   printf("Wrote %s\n",RESULT_FILE);
 
@@ -197,6 +222,50 @@ int main(int argc, char **argv) {
   if (status) {
     status = -2;
   }
+  //
+  //Send objects to the AMORE DB
+  //
+  printf ("AMORE part\n");
+  const char *amoreDANameorig=gSystem->Getenv("AMORE_DA_NAME");
+  //cheet a little -- temporary solution (hopefully)
+  // 
+  //currently amoreDA uses the environment variable AMORE_DA_NAME to create the mysql
+  //table in which the calib objects are stored. This table is dropped each time AmoreDA
+  //is initialised. This of course makes a problem if we would like to store different
+  //calibration entries in the AMORE DB. Therefore in each DA which writes to the AMORE DB
+  //the AMORE_DA_NAME env variable is overwritten.  
+  //find processed sector
+  Char_t sideName='A';
+  Int_t sector = -1;
+  for ( Int_t roc = 0; roc < 72; roc++ ) {
+    if ( !calibPulser.GetCalRocT0(roc) ) continue;
+    if (mapping->GetSideFromRoc(roc)==1) sideName='C';
+    sector = mapping->GetSectorFromRoc(roc);
+  }
+//   gSystem->Setenv("AMORE_DA_NAME",Form("TPC-%c%02d-%s",sideName,sector,FILE_ID));
+  gSystem->Setenv("AMORE_DA_NAME",Form("%s-%s",gSystem->Getenv("DATE_ROLE_NAME"),FILE_ID));
+  
+  // 
+  // end cheet
+  if (sector>-1){
+    TDatime time;
+    TObjString info(Form("Run: %u; Date: %s",runNb,time.AsSQLString()));
+    
+    amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender);
+    Int_t statusDA=0;
+    statusDA+=amoreDA.Send("PulserT0",calibPulser.GetCalPadT0());
+    statusDA+=amoreDA.Send("PulserQ",calibPulser.GetCalPadQ());
+    statusDA+=amoreDA.Send("PulserRMS",calibPulser.GetCalPadRMS());
+    statusDA+=amoreDA.Send("arrayTmean",calibPulser.GetMeanTimeSectorArray());
+    statusDA+=amoreDA.Send("Info",&info);
+    if ( statusDA!=0 )
+      printf("Waring: Failed to write one of the calib objects to the AMORE database\n");
+  } else {
+    printf("Waring: No data found!\n");
+  }
+  // reset env var
+  if (amoreDANameorig) gSystem->Setenv("AMORE_DA_NAME",amoreDANameorig);
 
   return status;
 }