Corrected identifier for tpcCalibRaw
[u/mrichter/AliRoot.git] / TPC / AliTPCPreprocessor.cxx
index 179d9f8..11b2b02 100644 (file)
@@ -27,7 +27,9 @@
 #include "AliTPCCalibPedestal.h"
 #include "AliTPCCalibPulser.h"
 #include "AliTPCCalibCE.h"
+#include "AliTPCCalibRaw.h"
 #include "AliTPCdataQA.h"
+#include "ARVersion.h"
 #include "TFile.h"
 #include "TTree.h"
 #include "TGraph.h" 
@@ -325,6 +327,35 @@ UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap)
     }
   }
 
+
+// raw calibration processing
+
+  if(runType == kPhysicsRunType) {
+    Int_t numSources = 1;
+    Int_t rawSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ;
+    TString source = fConfEnv->GetValue("Raw","DAQ");
+    source.ToUpper();
+    if ( source != "OFF") { 
+     if ( source == "HLT") rawSource[0] = AliShuttleInterface::kHLT;
+     if (!GetHLTStatus()) rawSource[0] = AliShuttleInterface::kDAQ;
+     if (source == "HLTDAQ" ) {
+         numSources=2;
+        rawSource[0] = AliShuttleInterface::kHLT;
+        rawSource[1] = AliShuttleInterface::kDAQ;
+     }
+     if (source == "DAQHLT" ) numSources=2;
+     UInt_t rawResult=0;
+     for (Int_t i=0; i<numSources; i++ ) {     
+       rawResult = ExtractRaw(rawSource[i]);
+       if ( rawResult == 0 ) break;
+     }
+     result += rawResult;
+     status = new TParameter<int>("rawResult",rawResult);
+     resultArray->Add(status);
+    }
+  }
+
+
   // Altro configuration
 
 
@@ -341,9 +372,8 @@ UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap)
   // Central Electrode processing
 
   if( runType == kPhysicsRunType || 
-      runType == kDaqRunType ) {    
+      runType == kLaserRunType ) {    
 
-//   if (true) {                 // do CE processing for all run types
     Int_t numSources = 1;
     Int_t ceSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ;
     TString source = fConfEnv->GetValue("CE","DAQ");
@@ -362,7 +392,11 @@ UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap)
        ceResult = ExtractCE(ceSource[i]);
        if ( ceResult == 0 ) break;
      }
-     result += ceResult;
+
+   // only flag error if CE result is missing from LASER runs
+   //    -- for PHYSICS run do CE processing if data available
+   
+     if ( runType == kLaserRunType ) result += ceResult;
      status = new TParameter<int>("ceResult",ceResult);
      resultArray->Add(status);
 
@@ -396,6 +430,7 @@ UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap)
     AliCDBMetaData metaData;
     metaData.SetBeamPeriod(0);
     metaData.SetResponsible("Haavard Helstrup");
+    metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
     metaData.SetComment("Preprocessor AliTPC status.");
     Store("Calib", "PreprocStatus", resultArray, &metaData, 0, kFALSE);
     resultArray->Delete();
@@ -432,6 +467,7 @@ UInt_t AliTPCPreprocessor::MapTemperature(TMap* dcsAliasMap)
         AliCDBMetaData metaData;
        metaData.SetBeamPeriod(0);
        metaData.SetResponsible("Haavard Helstrup");
+       metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
        metaData.SetComment("Preprocessor AliTPC data base entries.");
 
        Bool_t storeOK = Store("Calib", "Temperature", fTemp, &metaData, 0, kFALSE);
@@ -489,6 +525,7 @@ UInt_t AliTPCPreprocessor::MapHighVoltage(TMap* dcsAliasMap)
         AliCDBMetaData metaData;
        metaData.SetBeamPeriod(0);
        metaData.SetResponsible("Haavard Helstrup");
+       metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
        metaData.SetComment("Preprocessor AliTPC data base entries.");
 
        Bool_t storeOK = Store("Calib", "HighVoltage", fHighVoltage, &metaData, 0, kFALSE);
@@ -524,6 +561,7 @@ UInt_t AliTPCPreprocessor::MapGoofie(TMap* dcsAliasMap)
         AliCDBMetaData metaData;
        metaData.SetBeamPeriod(0);
        metaData.SetResponsible("Haavard Helstrup");
+        metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
        metaData.SetComment("Preprocessor AliTPC data base entries.");
 
        Bool_t storeOK = Store("Calib", "Goofie", fGoofie, &metaData, 0, kFALSE);
@@ -612,6 +650,7 @@ UInt_t AliTPCPreprocessor::ExtractPedestals(Int_t sourceFXS)
      AliCDBMetaData metaData;
      metaData.SetBeamPeriod(0);
      metaData.SetResponsible("Haavard Helstrup");
+     metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
      metaData.SetComment("Preprocessor AliTPC data base entries."); 
  
      Bool_t storeOK = Store("Calib", "Pedestals", calPadPed, &metaData, 0, kTRUE);
@@ -717,6 +756,7 @@ UInt_t AliTPCPreprocessor::ExtractPulser(Int_t sourceFXS)
      AliCDBMetaData metaData;
      metaData.SetBeamPeriod(0);
      metaData.SetResponsible("Haavard Helstrup");
+     metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
      metaData.SetComment("Preprocessor AliTPC data base entries.");
 
      Bool_t storeOK = Store("Calib", "Pulser", pulserObjects, &metaData, 0, kTRUE);
@@ -730,58 +770,103 @@ UInt_t AliTPCPreprocessor::ExtractPulser(Int_t sourceFXS)
   return result;
 }
 
+//______________________________________________________________________________________________
+
+
+UInt_t AliTPCPreprocessor::ExtractRaw(Int_t sourceFXS)
+{
+ //
+ //  Read Raw calibration file from file exchage server
+ //
+ UInt_t result=0;
+ TObjArray* rawArray = new TObjArray;
+
+ TList* list = GetFileSources(sourceFXS,"tpcCalibRaw");
+ if (list && list->GetEntries()>0) {
+
+//  loop through all files
+
+    UInt_t index = 0;
+    while (list->At(index)!=NULL) {
+     TObjString* fileNameEntry = (TObjString*) list->At(index);
+     if (fileNameEntry!=NULL) {
+        TString fileName = GetFile(sourceFXS, "tpcCalibRaw",
+                                        fileNameEntry->GetString().Data());
+        TFile *f = TFile::Open(fileName);
+        if (!f) {
+         Log ("Error opening raw file.");
+         result =2;
+         break;
+       }
+        AliTPCCalibRaw *calRaw;
+       f->GetObject("tpcCalibRaw",calRaw);
+        if ( !calRaw ) {
+         Log ("No raw calibration object in file.");
+         result = 2;
+         break;
+       }
+       rawArray->Add(calRaw);
+       delete calRaw;
+       f->Close();
+      }
+     ++index;
+    }  // while(list)
+//
+//  Store updated pedestal entry to OCDB
+//
+     AliCDBMetaData metaData;
+     metaData.SetBeamPeriod(0);
+     metaData.SetResponsible("Haavard Helstrup");
+     metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
+     metaData.SetComment("Preprocessor AliTPC data base entries.");
+
+     Bool_t storeOK = Store("Calib", "Raw", rawArray, &metaData, 0, kTRUE);
+     if ( !storeOK ) ++result;
+  } else {
+    Log ("Error: no entries in input file list!");
+    result = 1;
+  }
+
+  return result;
+}
+//______________________________________________________________________________________________
+
 UInt_t AliTPCPreprocessor::ExtractCE(Int_t sourceFXS)
 {
  //
  //  Read Central Electrode file from file exchage server
- //  Keep original entry from OCDB in case no new CE calibration is available
+ //  
  //
- TObjArray    *ceObjects=0;
  AliTPCCalPad *ceTmean=0;
  AliTPCCalPad *ceTrms=0;
  AliTPCCalPad *ceQmean=0;
  TObjArray    *rocTtime=0;  
  TObjArray    *rocQtime=0;  
 
- AliCDBEntry* entry = GetFromOCDB("Calib", "CE");
- if (entry) ceObjects = (TObjArray*)entry->GetObject();
- if ( ceObjects==NULL ) {
-     Log("AliTPCPreprocsessor: No previous TPC central electrode entry available.\n");
-     ceObjects = new TObjArray;    
- }
+ TObjArray    *ceObjects= new TObjArray;
+  
 
  Int_t nSectors = fROC->GetNSectors();
 
- ceTmean = (AliTPCCalPad*)ceObjects->FindObject("CETmean");
- if ( !ceTmean ) {
-    ceTmean = new AliTPCCalPad("CETmean","CETmean");
-    ceObjects->Add(ceTmean);
- }
- ceTrms = (AliTPCCalPad*)ceObjects->FindObject("CETrms");
- if ( !ceTrms )  { 
-    ceTrms = new AliTPCCalPad("CETrms","CETrms");
-    ceObjects->Add(ceTrms);
- }
- ceQmean = (AliTPCCalPad*)ceObjects->FindObject("CEQmean");
- if ( !ceQmean )  { 
-    ceQmean = new AliTPCCalPad("CEQmean","CEQmean");
-    ceObjects->Add(ceQmean);
- }
- //!new from here please have a look!!!
- rocTtime = (TObjArray*)ceObjects->FindObject("rocTtime");
- if ( !rocTtime ) {
-     rocTtime = new TObjArray(nSectors);
-     rocTtime->SetName("rocTtime");
-     ceObjects->Add(rocTtime);
- }
- rocQtime = (TObjArray*)ceObjects->FindObject("rocQtime");
- if ( !rocQtime ) {
-     rocQtime = new TObjArray(nSectors);
-     rocQtime->SetName("rocQtime");
-     ceObjects->Add(rocQtime);
- }
+ ceTmean = new AliTPCCalPad("CETmean","CETmean");
+ ceObjects->Add(ceTmean);
 
+ ceTrms = new AliTPCCalPad("CETrms","CETrms");
+ ceObjects->Add(ceTrms);
+
+ ceQmean = new AliTPCCalPad("CEQmean","CEQmean");
+ ceObjects->Add(ceQmean);
+ rocTtime = new TObjArray(nSectors+2);   // also make room for A and C side average
+ rocTtime->SetName("rocTtime");
+ ceObjects->Add(rocTtime);
+ rocQtime = new TObjArray(nSectors);
+ rocQtime->SetName("rocQtime");
+ ceObjects->Add(rocQtime);
 
  UInt_t result=0;
 
@@ -825,6 +910,13 @@ UInt_t AliTPCPreprocessor::ExtractCE(Int_t sourceFXS)
           TGraph *grQ=calCE->MakeGraphTimeCE(sector,0,3); // Q time graph
            if ( grQ ) rocTtime->AddAt(grQ,sector);         
         }
+
+       TGraph *grT=calCE->MakeGraphTimeCE(-1,0,2); // A side average
+       if ( grT ) rocTtime->AddAt(grT,nSectors);         
+       grT=calCE->MakeGraphTimeCE(-2,0,2); // C side average
+       if ( grT ) rocTtime->AddAt(grT,nSectors+1);         
+
+
        delete calCE;
        f->Close();
       }
@@ -836,6 +928,7 @@ UInt_t AliTPCPreprocessor::ExtractCE(Int_t sourceFXS)
     AliCDBMetaData metaData;
     metaData.SetBeamPeriod(0);
     metaData.SetResponsible("Haavard Helstrup");
+    metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
     metaData.SetComment("Preprocessor AliTPC data base entries.");
 
     Bool_t storeOK = Store("Calib", "CE", ceObjects, &metaData, 0, kTRUE);
@@ -886,6 +979,7 @@ UInt_t AliTPCPreprocessor::ExtractQA(Int_t sourceFXS)
            AliCDBMetaData metaData;
            metaData.SetBeamPeriod(0);
            metaData.SetResponsible("Haavard Helstrup");
+          metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
            metaData.SetComment("Preprocessor AliTPC data base entries.");
 
            Bool_t storeOK = Store("Calib", "QA", calQA, &metaData, 0, kTRUE);
@@ -1038,6 +1132,7 @@ UInt_t AliTPCPreprocessor::ExtractAltro(Int_t sourceFXS)
      AliCDBMetaData metaData;
      metaData.SetBeamPeriod(0);
      metaData.SetResponsible("Haavard Helstrup");
+     metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH);
      metaData.SetComment("Preprocessor AliTPC data base entries.");
 
      Bool_t storeOK = Store("Calib", "AltroConfig", altroObjects, &metaData, 0, kTRUE);