Merge branch 'master' into LocalDev
[u/mrichter/AliRoot.git] / PHOS / PHOSGAINda.cxx
index 06598d0..803281f 100644 (file)
@@ -1,12 +1,12 @@
 /*
 contact: Boris.Polishchuk@cern.ch
 link: see comments in the $ALICE_ROOT/PHOS/AliPHOSRcuDA1.cxx
-reference run: /castor/cern.ch/alice/phos/2007/10/04/18/07000008249001.1000.root
+reference run: /alice/data/2009/LHC09b_PHOS/000075883/raw/09000075883017.20.root
 run type: PHYSICS
 DA type: MON 
 number of events needed: 1000
-input files: RCU0.data  RCU1.data  RCU2.data  RCU3.data
-Output files: PHOS_Module2_Calib.root
+input files: Mod0RCU0.data Mod0RCU1.data Mod0RCU2.data Mod0RCU3.data Mod1RCU0.data Mod1RCU1.data Mod1RCU2.data Mod1RCU3.data Mod2RCU0.data Mod2RCU1.data Mod2RCU2.data Mod2RCU3.data Mod3RCU0.data Mod3RCU1.data Mod3RCU2.data Mod3RCU3.data Mod4RCU0.data Mod4RCU1.data Mod4RCU2.data Mod4RCU3.data zs.txt
+Output files: PHOS_Calib_Total.root contains cumulative statistics for a number of runs. 
 Trigger types used: PHYSICS
 */
 
@@ -30,7 +30,7 @@ extern "C" {
 #include "AliPHOSRawFitterv0.h"
 #include "AliCaloAltroMapping.h"
 #include "AliCaloRawStreamV3.h"
-
+#include "AliLog.h"
 
 /* Main routine
       Arguments: 
@@ -44,6 +44,9 @@ int main(int argc, char **argv) {
                                        "RIO",
                                        "TStreamerInfo()");
 
+  AliLog::SetGlobalDebugLevel(0) ;
+  AliLog::SetGlobalLogLevel(AliLog::kFatal);
+  
   int status;
   
   if (argc!=2) {
@@ -63,9 +66,30 @@ int main(int argc, char **argv) {
   }
   
   /* Retrieve mapping files from DAQ DB */ 
-  const char* mapFiles[4] = {"RCU0.data","RCU1.data","RCU2.data","RCU3.data"};
+  const char* mapFiles[20] = {
+    "Mod0RCU0.data",
+    "Mod0RCU1.data",
+    "Mod0RCU2.data",
+    "Mod0RCU3.data",
+    "Mod1RCU0.data",
+    "Mod1RCU1.data",
+    "Mod1RCU2.data",
+    "Mod1RCU3.data",
+    "Mod2RCU0.data",
+    "Mod2RCU1.data",
+    "Mod2RCU2.data",
+    "Mod2RCU3.data",
+    "Mod3RCU0.data",
+    "Mod3RCU1.data",
+    "Mod3RCU2.data",
+    "Mod3RCU3.data",
+    "Mod4RCU0.data",
+    "Mod4RCU1.data",
+    "Mod4RCU2.data",
+    "Mod4RCU3.data"
+  };
   
-  for(Int_t iFile=0; iFile<4; iFile++) {
+  for(Int_t iFile=0; iFile<20; iFile++) {
     int failed = daqDA_DB_getFile(mapFiles[iFile], mapFiles[iFile]);
     if(failed) { 
       printf("Cannot retrieve file %s from DAQ DB. Exit.\n",mapFiles[iFile]);
@@ -74,15 +98,26 @@ int main(int argc, char **argv) {
   }
   
   /* Open mapping files */
-  AliAltroMapping *mapping[4];
+  AliAltroMapping *mapping[20];
   TString path = "./";
-  path += "RCU";
+  
+  path += "Mod";
   TString path2;
-  for(Int_t i = 0; i < 4; i++) {
+  TString path3;
+  Int_t iMap = 0;
+  
+  for(Int_t iMod = 0; iMod < 5; iMod++) {
     path2 = path;
-    path2 += i;
-    path2 += ".data";
-    mapping[i] = new AliCaloAltroMapping(path2.Data());
+    path2 += iMod;
+    path2 += "RCU";
+    
+    for(Int_t iRCU=0; iRCU<4; iRCU++) {
+      path3 = path2;
+      path3 += iRCU;
+      path3 += ".data";
+      mapping[iMap] = new AliCaloAltroMapping(path3.Data());
+      iMap++;
+    }
   }
   
   /* define data source : this is argument 1 */  
@@ -178,7 +213,21 @@ int main(int argc, char **argv) {
       
       while (stream.NextDDL()) {
        while (stream.NextChannel()) {
-
+         
+         /* Retrieve ZS parameters from data*/
+         if(failZS) {
+           short value = stream.GetAltroCFG1();
+           bool ZeroSuppressionEnabled = (value >> 15) & 0x1;
+           bool AutomaticBaselineSubtraction = (value >> 14) & 0x1;
+           if(ZeroSuppressionEnabled) {
+             offset = (value >> 10) & 0xf;
+             threshold = value & 0x3ff;
+             fitter.SubtractPedestals(kFALSE);
+             fitter.SetAmpOffset(offset);
+             fitter.SetAmpThreshold(threshold);
+           }
+         }
+         
          cellX    = stream.GetCellX();
          cellZ    = stream.GetCellZ();
          caloFlag = stream.GetCaloFlag();  // 0=LG, 1=HG, 2=TRU
@@ -198,7 +247,7 @@ int main(int argc, char **argv) {
            fitter.Eval(stream.GetSignals(),sigStart,sigLength);
          } // End of NextBunch()
          
-         if (nBunches>1) continue;
+         if (nBunches != 1) continue;
          
          e[cellX][cellZ][caloFlag] = fitter.GetEnergy();
          t[cellX][cellZ][caloFlag] = fitter.GetTime();
@@ -209,7 +258,7 @@ int main(int argc, char **argv) {
        if(dAs[stream.GetModule()])
          dAs[stream.GetModule()]->FillHistograms(e,t);
        else {
-         dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1);
+         dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1,0);
          dAs[stream.GetModule()]->FillHistograms(e,t);
        }
        
@@ -243,22 +292,44 @@ int main(int argc, char **argv) {
     }
   }
   
-  for(Int_t i = 0; i < 4; i++) delete mapping[i];  
+  for(Int_t i = 0; i < 20; i++) delete mapping[i];  
   
   /* Be sure that all histograms are saved */
+
+  char h2name[80];
+  char totfile[80];
   
-  char localfile[128];
+  //Write the Total file (accumulated statistics for number of runs)
+  sprintf(totfile,"PHOS_Calib_Total.root");
+  TFile * ftot = new TFile(totfile,"recreate");
   
-  for(Int_t iMod=0; iMod<5; iMod++) {
-    if(!dAs[iMod]) continue;
+//   if (!ftot->IsZombie()){
+//     printf("Updating file %s.\n",ftot->GetName());
+
+//     for(Int_t iMod=0; iMod<5; iMod++) {
+//       if(!dAs[iMod]) continue;
     
-    dAs[iMod]->UpdateHistoFile();
-    dAs[iMod]->SetWriteToFile(kFALSE);    
+//       printf("DA1 for module %d detected.\n",iMod);
     
-    /* Store output files to the File Exchange Server */
-    sprintf(localfile,"PHOS_Module%d_Calib.root",iMod);
-    daqDA_FES_storeFile(localfile,"AMPLITUDES");
-  }
+//       for(Int_t iX=0; iX<64; iX++) {
+//     for(Int_t iZ=0; iZ<56; iZ++) {
+       
+//       for(Int_t iGain=0; iGain<2; iGain++) {
+//         sprintf(h2name,"%d_%d_%d_%d",iMod,iX,iZ,iGain);
+//         TH2F* h2tot = (TH2F*)ftot->Get(h2name);
+//         const TH2F* h2run = dAs[iMod]->GetTimeEnergyHistogram(iX,iZ,iGain); // Time vs Energy
+//         if(!h2tot && h2run) h2run->Write();
+//         if(h2tot && h2run) { h2tot->Add(h2run); h2tot->Write(h2tot->GetName(),TObject::kWriteDelete); }
+//       }
+//     }
+//       }
+//     }
+//   }
+  
+  ftot->Close();
+  
+  /* Store output files to the File Exchange Server */
+  daqDA_FES_storeFile(totfile,"AMPLITUDES");
   
   return status;
 }