]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TPC/TPCCEda.cxx
Updated flags for low flux case (A. Dainese)
[u/mrichter/AliRoot.git] / TPC / TPCCEda.cxx
index d53d1791dbb36f774d342a6a2a4b035cee32e4f1..f81332cbc62a87312c8d644e467097eda05675c9 100644 (file)
@@ -19,6 +19,7 @@ TPCCEda.cxx - calibration algorithm for TPC Central Electrode events
 
 10/06/2007  sylvain.chapeland@cern.ch :  first version - clean skeleton based on DAQ DA case1
 06/12/2007  haavard.helstrup@cern.ch  :  created CE DA based on pulser code
+19/09/2008  J.Wiechula@gsi.de:           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 "tpcCE.root"
 #define FILE_ID "CE"
 #define MAPPING_FILE "tpcMapping.root"
+#define CONFIG_FILE "TPCCEda.conf"
 
 
 #include <daqDA.h>
@@ -45,6 +47,10 @@ and save results in a file (named from RESULT_FILE define - see below).
 #include <TFile.h>
 #include "TROOT.h"
 #include "TPluginManager.h"
+#include "TSystem.h"
+#include "TString.h"
+#include "TObjString.h"
+#include "TDatime.h"
 //
 //AliRoot includes
 //
@@ -58,8 +64,11 @@ and save results in a file (named from RESULT_FILE define - see below).
 #include "AliMathBase.h"
 #include "TTreeStream.h"
 #include "AliLog.h"
-#include "TSystem.h"
-
+#include "AliTPCConfigDA.h"
+//
+//AMORE
+//
+#include <AmoreDA.h>
 //
 // TPC calibration algorithm includes
 //
@@ -72,6 +81,18 @@ and save results in a file (named from RESULT_FILE define - see below).
       Arguments: list of DATE raw data files
 */
 int main(int argc, char **argv) {
+  /* log start of process */
+  printf("TPC CE 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",
                                          "*",
@@ -79,22 +100,31 @@ int main(int argc, char **argv) {
                                          "RIO",
                                          "TStreamerInfo()");
 
-  AliLog::SetClassDebugLevel("AliTPCRawStream",-5);
-  AliLog::SetClassDebugLevel("AliRawReaderDate",-5);
-  AliLog::SetClassDebugLevel("AliTPCAltroMapping",-5);
-  AliLog::SetModuleDebugLevel("RAW",-5);
-
-  int i,status;
+ /* declare monitoring program */
+ int i,status;
+ status=monitorDeclareMp( __FILE__ );
+ if (status!=0) {
+   printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
+   return -1;
+ }
+ monitorSetNowait();
+ monitorSetNoWaitNetworkTimeout(1000);
+
+  //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 */
@@ -105,39 +135,37 @@ 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);
   }
-
-
-  AliTPCCalibCE calibCE;   // 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 CE DA started - %s\n",__FILE__);
-
-
-  /* set time bin range */
-  calibCE.SetRangeTime(800,940);
-  calibCE.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
+  if ( (Int_t)config.GetValue("UseFastDecoder") == 1 ){
+    printf("Info: The fast decoder will be used for the processing.\n");
+    fastDecoding=kTRUE;
   }
 
+  //create calibration object 
+  AliTPCCalibCE calibCE(config.GetConfigurationMap());   // central electrode calibration
+  calibCE.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++) {
+  for ( i=1; i<argc; i++) {
 
     /* define data source : this is argument i */
     printf("Processing file %s\n", argv[i]);
@@ -165,17 +193,22 @@ int main(int argc, char **argv) {
         printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
         break;
       }
+      
+      /* skip start/end of run events */
+      if ( (event->eventType != physicsEvent) && (event->eventType != calibrationEvent) )
+        continue;
 
       /* retry if got no event */
-      if (event==NULL) {
+      if (event==NULL)
         continue;
-      }
+      
       nevents++;
-
-      //  Pulser calibration
-
+      // get the run number
+      runNb = event->eventRunNb;
+      //  CE calibration
       AliRawReader *rawReader = new AliRawReaderDate((void*)event);
-      calibCE.ProcessEvent(rawReader);
+      if ( fastDecoding ) calibCE.ProcessEventFast(rawReader);
+      else calibCE.ProcessEvent(rawReader);
       delete rawReader;
 
       /* free resources */
@@ -183,7 +216,10 @@ int main(int argc, char **argv) {
     }
   }
 
-  calibCE.Analyse(); 
+  //
+  // Analyse CE data and write them to rootfile
+  //
+  calibCE.Analyse();
   printf ("%d events processed\n",nevents);
 
   TFile * fileTPC = new TFile (RESULT_FILE,"recreate");
@@ -198,5 +234,37 @@ int main(int argc, char **argv) {
     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.  
+  const char *role=gSystem->Getenv("DATE_ROLE_NAME");
+  gSystem->Setenv("AMORE_DA_NAME",Form("%s-%s",role,FILE_ID));
+  //
+  // end cheet
+  if ( role ){
+    TDatime time;
+    TObjString info(Form("Run: %u; Date: %s",runNb,time.AsString()));
+    amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender);
+    Int_t status=0;
+    status+=amoreDA.Send("CET0",calibCE.GetCalPadT0());
+    status+=amoreDA.Send("CEQ",calibCE.GetCalPadQ());
+    status+=amoreDA.Send("CERMS",calibCE.GetCalPadRMS());
+    status+=amoreDA.Send("Info",&info);
+    if ( status!=0 )
+      printf("Waring: Failed to write one of the calib objects to the AMORE database\n");
+    if (amoreDANameorig) gSystem->Setenv("AMORE_DA_NAME",amoreDANameorig);
+  } else {
+    printf ("Warning: environment variable 'AMORE_DA_NAME' not set. Cannot write to the AMORE database\n");
+  }
+
   return status;
 }