Too wide range for charge.

[u/mrichter/AliRoot.git] / HMPID / HMPIDda.cxx

diff --git a/HMPID/HMPIDda.cxx b/HMPID/HMPIDda.cxx
index 724b14605fcd78828a37aa9302cecf41e1f04772..b7012824456b23068314c2d432eb421f47865fb7 100644 (file)
--- a/HMPID/HMPIDda.cxx
+++ b/HMPID/HMPIDda.cxx
@@ -17,12 +17,14 @@ extern "C" {
 
 //AliRoot
 #include "AliHMPIDRawStream.h"
+#include "AliHMPIDCalib.h"
 #include "AliRawReaderDate.h"
 #include "AliBitPacking.h"
 #include "TMath.h"
 
 //ROOT
 #include "TFile.h"
+#include "TSystem.h"
 #include "TKey.h"
 #include "TH2S.h"
 #include "TObject.h"
@@ -30,6 +32,7 @@ extern "C" {
 #include "TMath.h"
 #include "TRandom.h"
 
+
 int main(int argc, char **argv){ 
 
   int status;
@@ -55,25 +58,14 @@ int main(int argc, char **argv){
   daqDA_progressReport(10);
 
 
-  /* init some counters */
-  Double_t SummQ[14][48][11][25], SummQ2[14][48][11][25];
-  Int_t   nEntries[14][48][11][25];
-  Bool_t  isDDLOn[14];
-  for(Int_t ddl=0;ddl<=13;ddl++) {
-    isDDLOn[ddl] = kFALSE;
-    for(Int_t row=1;row<=24;row++)
-      for(Int_t dil=1;dil<=10;dil++)
-       for(Int_t adr=0;adr<=47;adr++)
-         {
-           SummQ[ddl][adr][dil][row]=0;
-           SummQ2[ddl][adr][dil][row]=0;
-           nEntries[ddl][adr][dil][row]=0;
-         }
-  }
-
+  /* init the pedestal calculation */
+  AliHMPIDCalib *pCal=new AliHMPIDCalib();
+  
   /* init event counter */
   Int_t iEvtNcal=0;
 
+  Int_t cnt=0;
+  
   int n;
   for (n=1;n<argc;n++) {
 
@@ -93,7 +85,13 @@ int main(int argc, char **argv){
 
       /* get next event */
       status=monitorGetEventDynamic((void **)&event);
-      if (status==MON_ERR_EOF) break; /* end of monitoring file has been reached */
+      if (status==MON_ERR_EOF)                                              /* end of monitoring file has been reached */
+      {
+        printf("End of monitoring file has been reached! \n");
+        break;
+        }
+      
+       
       if (status!=0) {
         printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
        return -1;
@@ -101,16 +99,17 @@ int main(int argc, char **argv){
 
       /* retry if got no event */
       if (event==NULL) {
-        break;
+        //break;
+        continue;
       }
 
       /* use event - here, just write event id to result file */
       eventT=event->eventType;
 
-      if (eventT==PHYSICS_EVENT) {
+      if (eventT==PHYSICS_EVENT) {                                                //we use PHYSICS_EVENT for pedestal not CALIBRATION_EVENT
        
        iEvtNcal++;
-
+        
        AliRawReader *reader = new AliRawReaderDate((void*)event);
 
        // Temporary there. Shall be removed as soon as the equipment ID is set correctly
@@ -121,20 +120,14 @@ int main(int argc, char **argv){
        AliHMPIDRawStream stream(reader);
 
        while(stream.Next()) {
-
-         Int_t ddl = stream.GetDDLNumber();
-
-         for(Int_t row=1; row < 25; row++)
-           for(Int_t dil=1; dil < 11; dil++)
-             for(Int_t adr=0; adr < 48; adr++) {
-               Short_t q = stream.GetCharge(row,dil,adr);
-               if (q<0) continue;
-
-               SummQ[ddl][adr][dil][row] += q;
-               SummQ2[ddl][adr][dil][row] += (q*q);
-               nEntries[ddl][adr][dil][row]++;
-               isDDLOn[ddl] = kTRUE;
-             }
+          Int_t nDDL=stream.GetDDLNumber();
+            for(Int_t row = 1; row <=AliHMPIDRawStream::kNRows; row++){
+              for(Int_t dil = 1; dil <=AliHMPIDRawStream::kNDILOGICAdd; dil++){
+                for(Int_t pad = 0; pad < AliHMPIDRawStream::kNPadAdd; pad++){
+                  pCal->FillPedestal(nDDL,row,dil,pad,stream.GetCharge(nDDL,row,dil,pad));
+                }//pad
+              }//dil
+            }//row
        } //raw data loop
        
        delete reader;
@@ -151,7 +144,7 @@ int main(int argc, char **argv){
       free(event);
     }
 
-  }
+  }//arg
 
   /* write report */
   printf("Run #%s, received %d calibration events\n",getenv("DATE_RUN_NUMBER"),iEvtNcal);
@@ -164,41 +157,21 @@ int main(int argc, char **argv){
   /* report progress */
   daqDA_progressReport(90);
 
-  for(Int_t ddl=0; ddl < 14; ddl++) {
-    if (!isDDLOn[ddl]) continue;
-
-    ofstream out;
-    out.open(Form("./HmpidPedDdl%02i.txt",ddl));
-
-    for(Int_t row=1; row < 25; row++)
-      for(Int_t dil=1; dil < 11; dil++)
-       for(Int_t adr=0; adr < 48; adr++) {
-
-         Int_t n = nEntries[ddl][adr][dil][row];
-         if (!n) {
-           printf("No data for channel: %d %d %d %d\n",ddl,adr,dil,row);
-           continue;
-         }
-         Double_t mean = SummQ[ddl][adr][dil][row]/n;
-         if ((SummQ2[ddl][adr][dil][row]/n - (SummQ[ddl][adr][dil][row]/n)*(SummQ[ddl][adr][dil][row]/n)) < 0) {
-           printf("Invalid sums: %d %d %d   %d  %f %f\n",row,dil,adr,n,SummQ[ddl][adr][dil][row],SummQ2[ddl][adr][dil][row]);
-           return -1;
-         }
-          Float_t sigma = TMath::Sqrt(SummQ2[ddl][adr][dil][row]/n - (mean*mean));
-         Int_t inhard=((Int_t(mean))<<9)+Int_t(mean+3*sigma);
-
-         out << Form("%2i %2i %2i %5.2f %5.2f %x\n",row,dil,adr,mean,sigma,inhard);
-       }
-
+  for(Int_t nDDL=0; nDDL < AliHMPIDCalib::kNDDL; nDDL++) {
+    
+    /* Calculate pedestal for the given ddl, if there is no ddl go t next */
+    if(!pCal->CalcPedestal(nDDL,Form("./HmpidPedDdl%02i.txt",nDDL),iEvtNcal)) continue;
+    
     /* store the result file on FES */
-    status=daqDA_FES_storeFile(Form("./HmpidPedDdl%02i.txt",ddl),Form("HMPID_DA_Pedestals_ddl=%02i",ddl));
+    
+    status=daqDA_FES_storeFile(Form("./HmpidPedDdl%02i.txt",nDDL),Form("HMPID_DA_Pedestals_ddl=%02i",nDDL));
     if (status) {
       printf("Failed to export file : %d\n",status);
       return -1;
     }
-  }
-
-
+    
+  }//nDDL
+  
   /* report progress */
   daqDA_progressReport(100);