]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/TRDPEDESTALda.cxx
- fixing bug in cluster analyser
[u/mrichter/AliRoot.git] / TRD / TRDPEDESTALda.cxx
index 15534635c18fb074f3a84e0c81e0a03c3c0c3dec..23cd6403bb83816ee64775f7ca6e577926dd658c 100644 (file)
@@ -1,17 +1,19 @@
 /*
 
 Contact: r.bailhache@gsi.de
-Link: http://www-linux.gsi.de/~bailhach/PEDESTAL/raw.root.date
-Run Type: dedicated pedestals run
+Link:
+Reference run: 12170
+Run Type: PEDESTAL
 DA Type: LDC
 Number of events needed: 100
-Input Files: no config files, no previous result files, RAW DATA files from DDL = 1024 to DDL = 1041 included
-Output Files: trdCalibration.root, trdCalibration.root, no persitent file over runs
-Trigger types used: PEDESTAL_RUN
+Input Files: TRD raw files
+Output Files: trdCalibration.root
+Trigger types used: PHYSICS_EVENT
 
 */
 
-#define RESULT_FILE "trdCalibration.root"
+#define RESULT_FILE  "trdPedestal.root"
+#define FILE_ID "PADSTATUS"
 
 extern "C" {
 #include <daqDA.h>
@@ -27,26 +29,42 @@ extern "C" {
 //
 #include <TFile.h>
 #include <TStopwatch.h>
-
+#include "TROOT.h"
+#include "TPluginManager.h"
 //
 // AliRoot includes
 //
 #include "AliRawReader.h"
 #include "AliRawReaderDate.h"
-#include "AliTRDrawStreamTB.h"
+#include "AliTRDrawFastStream.h"
 #include "AliTRDrawStreamBase.h"
 #include "AliCDBManager.h"
+
+//
+//AMORE
+//
+#include <AmoreDA.h>
+
 //
 // AliRoot TRD calib classes
 //
 #include "AliTRDCalibPadStatus.h"
 
+//functios, implementation below
+//void SendToAmoreDB(TObject *o, unsigned long32 runNb);
 
 /* Main routine
       Arguments: list of DATE raw data files
 */
 int main(int argc, char **argv) {
 
+  /* magic line from Rene */
+  gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
+                                       "*",
+                                       "TStreamerInfo",
+                                       "RIO",
+                                       "TStreamerInfo()");
+  
   int status;
 
 
@@ -76,27 +94,16 @@ int main(int argc, char **argv) {
  
   //Instance of AliCDBManager: needed by AliTRDRawStream
   //AliCDBManager *man = AliCDBManager::Instance();
-  //man->SetDefaultStorage("local://$ALICE_ROOT");
+  //man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
   //man->SetRun(0);
   // AliTRDCalibPadStatus object
   AliTRDCalibPadStatus calipad = AliTRDCalibPadStatus();
   Bool_t passpadstatus = kTRUE;
+  //unsigned long32 runNb=0;      //run number
 
-
-  /*see the time*/
-  TStopwatch timer;
-  timer.Start();
-  
   // setting
-  // AliTRDrawStreamTB::SetNoDebug();
-  AliTRDrawStreamTB::SetNoErrorWarning();
-  AliTRDrawStreamTB::SetForceCleanDataOnly();
-  AliTRDrawStreamTB::AllowCorruptedData();
-  //AliTRDrawStreamTB::SetSkipCDH();
-  //AliTRDrawStreamTB::SetExtraWordsFix();
-  //AliTRDrawStreamTB::EnableDebugStream();
-  //AliTRDrawStreamTB::SetDumpHead(320);
-  //AliTRDrawStreamTB::SetDumpHead(80);
+  AliTRDrawFastStream::DisableSkipData();
 
   /* read the data files */
   int n;
@@ -113,6 +120,7 @@ int main(int argc, char **argv) {
     /* read the file  until EOF */
     for(;;) {
       struct eventHeaderStruct *event;
+      eventTypeType eventT;
       
       /* get next event */
       status=monitorGetEventDynamic((void **)&event);
@@ -130,21 +138,27 @@ int main(int argc, char **argv) {
         break;
       }
 
-      if(passpadstatus){
+      /* use event - here, just write event id to result file */
+      eventT=event->eventType;
+
+      if((eventT==PHYSICS_EVENT) && (passpadstatus)){
 
        AliRawReader *rawReader = new AliRawReaderDate((void*)event);
        rawReader->Select("TRD");
-       AliTRDrawStreamTB *trdRawStream = new AliTRDrawStreamTB((AliRawReader *)rawReader);
-       //trdRawStream->Init();
-       if(!calipad.ProcessEvent((AliTRDrawStreamBase *)trdRawStream,(Bool_t)nevents_total)) passpadstatus = kFALSE;
-       nevents++;
-       delete trdRawStream;
+       
+       Int_t result = calipad.ProcessEvent2((AliRawReader *) rawReader);
+       // 0 error, 1 no input, 2 output
+       if(result == 2) nevents++;
+       if(result == 0) passpadstatus = kFALSE;
        delete rawReader;
-     
+      
       }
 
       nevents_total++;
 
+      // get the run number
+      //runNb = event->eventRunNb;
+
       /* free resources */
       free(event);
     }
@@ -154,34 +168,96 @@ int main(int argc, char **argv) {
   /* report progress */
   printf("%d events processed and %d used\n",nevents_total,nevents);
 
-  /*see the time*/
-  timer.Stop();
-  timer.Print();
-
   /* write file in any case to see what happens in case of problems*/
-  /*see the time*/
-  TStopwatch timer1;
-  timer1.Start();
   TFile *fileTRD = new TFile(RESULT_FILE,"recreate");
+  if(nevents < 30) calipad.Destroy();
   calipad.AnalyseHisto();
   calipad.Write("calibpadstatus");
   fileTRD->Close();   
-  printf("Wrote local file %s\n",RESULT_FILE);
-  /*see the time*/
-  timer1.Stop();
-  timer1.Print();
-  
+     
   /* store the result file on FES */
-  status=daqDA_FES_storeFile(RESULT_FILE,RESULT_FILE);
+  status=daqDA_FES_storeFile(RESULT_FILE,FILE_ID);
   if (status) {
     printf("Failed to export file : %d\n",status);
     return -1;
   }
-  
-   
-  /* report progress */
-  daqDA_progressReport(100);
 
   
   return status;
+
+  /*
+
+  // send objects to AMORE
+  //SendToAmoreDB(calibpad,runNb);
+  
+  void SendToAmoreDB(TObject *calipad, unsigned long32 runNb)
+  {
+  //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.
+  
+  
+  //
+  // The reference data are stored in:
+  // PadStatus1 for sm-00-01-02
+  // PadStatus2 for sm-03-04-05
+  // PadStatus3 for sm-06-07-08
+  // PadStatus4 for sm-09-10-11
+  // PadStatus5 for sm-12-13-14
+  // PadStatus6 for sm-15-16-17
+  // PadStatus0 if nothing found..means problems
+  //
+  
+  ///////////////////
+  // Find wich LDC
+  ///////////////////
+  Int_t ldcnumber = -1;
+  Int_t sm = -1;
+  for (Int_t idet=0; idet<540; idet++) {
+  AliTRDCalROC *rocMean  = calipad->GetCalRocMean(idet, kFALSE);
+    if ( rocMean )  {
+      sm  = AliTRDgeometry::GetSector(idet);  
+      if((sm==0) || (sm==1) || (sm==2)) ldcnumber = 1;
+      if((sm==3) || (sm==4) || (sm==5)) ldcnumber = 2;
+      if((sm==6) || (sm==7) || (sm==8)) ldcnumber = 3;
+      if((sm==9) || (sm==10) || (sm==11)) ldcnumber = 1;
+      if((sm==12) || (sm==13) || (sm==14)) ldcnumber = 2;
+      if((sm==15) || (sm==16) || (sm==17)) ldcnumber = 3;
+    }
+  }
+  const char *amoreDANameorig=gSystem->Getenv("AMORE_DA_NAME");
+  
+  gSystem->Setenv("AMORE_DA_NAME",Form("TRD-dataQA-%02d-%s",ldcnumber,FILE_ID));
+  /////////////////////
+  // Send the stuff
+  ////////////////////
+  if (ldcnumber>-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("Pedestals",calibpad);
+  statusDA+=amoreDA.Send("Info",&info);
+  if ( statusDA )
+  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);
+  
+  }
+  
+  */
+
+
+
 }
+