Merge branch 'master' into LocalDev
[u/mrichter/AliRoot.git] / PHOS / PHOSGAINda.cxx
index 74da65b..803281f 100644 (file)
@@ -1,13 +1,13 @@
 /*
 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
-run type: PHYSICS, STANDALONE
+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
-trigger types used: L0
+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
 */
 
 
@@ -27,9 +27,10 @@ extern "C" {
 #include "AliRawReader.h"
 #include "AliRawReaderDate.h"
 #include "AliPHOSRcuDA1.h"
-#include "AliPHOSRawDecoder.h"
+#include "AliPHOSRawFitterv0.h"
 #include "AliCaloAltroMapping.h"
-
+#include "AliCaloRawStreamV3.h"
+#include "AliLog.h"
 
 /* Main routine
       Arguments: 
@@ -43,17 +44,52 @@ int main(int argc, char **argv) {
                                        "RIO",
                                        "TStreamerInfo()");
 
+  AliLog::SetGlobalDebugLevel(0) ;
+  AliLog::SetGlobalLogLevel(AliLog::kFatal);
+  
   int status;
   
   if (argc!=2) {
     printf("Wrong number of arguments\n");
     return -1;
   }
-
+  
+  /* Retrieve ZS parameters from DAQ DB */
+  const char* zsfile = "zs.txt";
+  int failZS = daqDA_DB_getFile(zsfile, zsfile);
+  
+  Int_t offset,threshold;
+  
+  if(!failZS) {
+    FILE *f = fopen(zsfile,"r");
+    int scan = fscanf(f,"%d %d",&offset,&threshold);
+  }
+  
   /* Retrieve mapping files from DAQ DB */ 
-  const char* mapFiles[4] = {"RCU0.data","RCU1.data","RCU2.data","RCU3.data"};
-
-  for(Int_t iFile=0; iFile<4; iFile++) {
+  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<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]);
@@ -62,18 +98,28 @@ 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 */  
   status=monitorSetDataSource( argv[1] );
   if (status!=0) {
@@ -99,15 +145,29 @@ int main(int argc, char **argv) {
   int nevents_total=0;
 
   AliRawReader *rawReader = NULL;
+  AliPHOSRcuDA1* dAs[5];
 
-  AliPHOSRcuDA1 da1(2,-1); // DA1 (Calibration DA) for module2
+  for(Int_t iMod=0; iMod<5; iMod++) {
+    dAs[iMod] = 0;
+  }
   
   Float_t e[64][56][2];
   Float_t t[64][56][2];
 
-  Int_t gain = -1;
-  Int_t X = -1;
-  Int_t Z = -1;
+  for(Int_t iX=0; iX<64; iX++) {
+    for(Int_t iZ=0; iZ<56; iZ++) {
+      for(Int_t iGain=0; iGain<2; iGain++) {
+       e[iX][iZ][iGain] = 0.;
+       t[iX][iZ][iGain] = 0.;
+      }
+    }
+  }
+  
+  Int_t cellX    = -1;
+  Int_t cellZ    = -1;
+  Int_t nBunches =  0;
+  Int_t sigStart, sigLength;
+  Int_t caloFlag;
 
   /* main loop (infinite) */
   for(;;) {
@@ -140,39 +200,84 @@ int main(int argc, char **argv) {
     
     if (eventT==PHYSICS_EVENT) {
       
-      for(Int_t iX=0; iX<64; iX++) {
-       for(Int_t iZ=0; iZ<56; iZ++) {
-         for(Int_t iGain=0; iGain<2; iGain++) {
-           e[iX][iZ][iGain] = 0.;
-           t[iX][iZ][iGain] = 0.;
-         }
-       }
-      }
-
       rawReader = new AliRawReaderDate((void*)event);
-//       AliPHOSRawDecoderv1 dc(rawReader,mapping);
-      AliPHOSRawDecoder dc(rawReader,mapping);
-      dc.SubtractPedestals(kTRUE);
+      AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping);
+      AliPHOSRawFitterv0 fitter;
+      fitter.SubtractPedestals(kTRUE); // assume that data is non-ZS
       
-      while(dc.NextDigit()) {
-
-       X = dc.GetRow() - 1;
-       Z = dc.GetColumn() - 1;
-
-       if(dc.IsLowGain()) gain = 0;
-       else
-         gain = 1;
+      if(!failZS) {
+       fitter.SubtractPedestals(kFALSE);
+       fitter.SetAmpOffset(offset);
+       fitter.SetAmpThreshold(threshold);
+      }
+      
+      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
+         
+         if(caloFlag!=0 && caloFlag!=1) continue; //TRU data!
+         
+         // In case of oscillating signals with ZS, 
+         //a channel can have several bunches.
+         
+         nBunches = 0;
+         while (stream.NextBunch()) {
+           nBunches++;
+           if (nBunches > 1) continue;
+           sigStart  = stream.GetStartTimeBin();
+           sigLength = stream.GetBunchLength();
+           fitter.SetChannelGeo(stream.GetModule(),cellX,cellZ,caloFlag);
+           fitter.Eval(stream.GetSignals(),sigStart,sigLength);
+         } // End of NextBunch()
+         
+         if (nBunches != 1) continue;
+         
+         e[cellX][cellZ][caloFlag] = fitter.GetEnergy();
+         t[cellX][cellZ][caloFlag] = fitter.GetTime();
+       }
+       
+       if(stream.GetModule()<0 || stream.GetModule()>4) continue;
        
-       e[X][Z][gain] = dc.GetEnergy();
-       t[X][Z][gain] = dc.GetTime();
+       if(dAs[stream.GetModule()])
+         dAs[stream.GetModule()]->FillHistograms(e,t);
+       else {
+         dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1,0);
+         dAs[stream.GetModule()]->FillHistograms(e,t);
+       }
        
+       for(Int_t iX=0; iX<64; iX++) {
+         for(Int_t iZ=0; iZ<56; iZ++) {
+           for(Int_t iGain=0; iGain<2; iGain++) {
+             e[iX][iZ][iGain] = 0.;
+             t[iX][iZ][iGain] = 0.;
+           }
+         }
+       }
+
       }
 
-      da1.FillHistograms(e,t);
-      //da1.UpdateHistoFile();
+//       da1.FillHistograms(e,t);
+//     //da1.UpdateHistoFile();
       
-      delete rawReader;     
-      nevents_physics++;
+       delete rawReader;     
+       nevents_physics++;
     }
     
     nevents_total++;
@@ -187,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 */
 
-  da1.UpdateHistoFile();
-  da1.SetWriteToFile(kFALSE);
+  char h2name[80];
+  char totfile[80];
   
-  /* Store output files to the File Exchange Server */
-  char localfile[128];
-  char fesfileid[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++) {
-    sprintf(localfile,"PHOS_Module%d_Calib.root",iMod);
-    sprintf(fesfileid,"PHOS_Module%d_AMPLITUDES",iMod);
-    daqDA_FES_storeFile(localfile,fesfileid);
-  }
+//   if (!ftot->IsZombie()){
+//     printf("Updating file %s.\n",ftot->GetName());
+
+//     for(Int_t iMod=0; iMod<5; iMod++) {
+//       if(!dAs[iMod]) continue;
+    
+//       printf("DA1 for module %d detected.\n",iMod);
+    
+//       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;
 }