ECS run type is PEDESTAL and not PEDESTAL_RUN

[u/mrichter/AliRoot.git] / ITS / ITSSDDBASda.cxx

diff --git a/ITS/ITSSDDBASda.cxx b/ITS/ITSSDDBASda.cxx
index 6993725f28b8613f0ff77d6a5f322b7a03873784..67016265d58c01eb1516efb5bc9ddce3611ec563 100644 (file)
--- a/ITS/ITSSDDBASda.cxx
+++ b/ITS/ITSSDDBASda.cxx
@@ -1,3 +1,15 @@
+/*
+- Contact: - prino@to.infn.it
+- Link: - http://www.to.infn.it/~prino/alice/RawData/run11161.date
+- Run Type: - PEDESTAL_RUN
+- DA Type: - LDC
+- Number of events needed: 100
+- Input Files: - 
+- Output Files: - SDDbase_step1_ddl*c*_sid*.data SDDbase_step2_ddl*c*_sid*.data
+- Trigger types used: 
+*/
+
+
 //////////////////////////////////////////////////////////////////////////////
 // Detector Algorithm for analysis of SDD baseline runs.                    //
 //                                                                          //
@@ -30,6 +42,8 @@ extern "C" {
 #include <TFile.h>
 #include <TH1F.h>
 #include <TH2F.h>
+#include <TROOT.h>
+#include <TPluginManager.h>
 
 // AliRoot includes
 #include "AliRawReaderDate.h"
@@ -43,9 +57,15 @@ int main(int argc, char **argv) {
 
   int status = 0;
 
+  // line added to solve IO problems
+  gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
+                                        "*",
+                                        "TStreamerInfo",
+                                        "RIO",
+                                        "TStreamerInfo()");
 
   /* log start of process */
-  printf("ITS SDD BASE algorithm program started\n");  
+  printf("ITS SDD BASELINE+NOISE algorithm program started\n");  
 
 
   /* check that we got some arguments = list of files */
@@ -55,20 +75,25 @@ int main(int argc, char **argv) {
   }
 
 
-  Int_t maxNEvents=10; // maximum number of events to be analyzed
-  const Int_t nSDDmodules=12;  // temp for test raw data
-  AliITSOnlineSDDBase **base=new AliITSOnlineSDDBase*[2*nSDDmodules];
-  AliITSOnlineSDDCMN **corr=new AliITSOnlineSDDCMN*[2*nSDDmodules];
-  TH2F **histo=new TH2F*[2*nSDDmodules];
+  Int_t eqOffset = 256;
+  Int_t DDLrange = 24;
+  Int_t maxNEvents=18; // maximum number of events to be analyzed
+  const Int_t kTotDDL=24;
+  const Int_t kModPerDDL=12;
+  const Int_t kSides=2;
+  AliITSOnlineSDDBase **base=new AliITSOnlineSDDBase*[kTotDDL*kModPerDDL*kSides];
+  AliITSOnlineSDDCMN **corr=new AliITSOnlineSDDCMN*[kTotDDL*kModPerDDL*kSides];
+  TH2F **histo=new TH2F*[kTotDDL*kModPerDDL*kSides];
 
   Char_t hisnam[20];
-  for(Int_t imod=0; imod<nSDDmodules;imod++){
-    for(Int_t isid=0;isid<2;isid++){
-      Int_t index=2*imod+isid;
-      base[index]=new AliITSOnlineSDDBase(imod,isid);
-      corr[index]=new AliITSOnlineSDDCMN(imod,isid);
-      sprintf(hisnam,"his%03ds%d",imod,isid);
-      histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
+  for(Int_t iddl=0; iddl<kTotDDL;iddl++){
+    for(Int_t imod=0; imod<kModPerDDL;imod++){
+      for(Int_t isid=0;isid<kSides;isid++){
+       Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
+       base[index]=new AliITSOnlineSDDBase(iddl,imod,isid);
+       sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid);
+       histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
+      }
     }
   }
   
@@ -79,10 +104,12 @@ int main(int argc, char **argv) {
     printf("Start Analysis Step %d\n",iStep);
     iev=0;
     if(iStep==1){
-      for(Int_t imod=0; imod<nSDDmodules;imod++){
-       for(Int_t isid=0;isid<2;isid++){
-         Int_t index=2*imod+isid;
-         corr[index]->Reset();
+      for(Int_t iddl=0; iddl<kTotDDL;iddl++){
+       for(Int_t imod=0; imod<kModPerDDL;imod++){
+         for(Int_t isid=0;isid<kSides;isid++){
+           Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
+           corr[index]=new AliITSOnlineSDDCMN(iddl,imod,isid);
+         }
        }
       }
     }
@@ -144,8 +171,7 @@ int main(int argc, char **argv) {
        case PHYSICS_EVENT: // uncomment this line for test raw data
          printf(" event number = %i \n",iev);
          AliRawReader *rawReader = new AliRawReaderDate((void*)event);
-         rawReader->RequireHeader(kFALSE);
-         rawReader->SelectEquipment(17,101,101);  // temp for test raw data
+         rawReader->SelectEquipment(17,eqOffset+1,eqOffset+DDLrange);
 
 
          Int_t evtyp=0;
@@ -155,32 +181,32 @@ int main(int argc, char **argv) {
          }
 
          rawReader->Reset();
-         for(Int_t imod=0; imod<nSDDmodules;imod++){
-           for(Int_t isid=0; isid<2;isid++){
-             Int_t index=2*imod+isid;
-             histo[index]->Reset();
+         for(Int_t iddl=0; iddl<kTotDDL;iddl++){
+           for(Int_t imod=0; imod<kModPerDDL;imod++){
+             for(Int_t isid=0;isid<kSides;isid++){
+               Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
+               histo[index]->Reset();
+             }
            }
          }
          AliITSRawStreamSDD s(rawReader);
-         rawReader->SelectEquipment(17,101,101);  // temp for test raw data
 
          while(s.Next()){
-           Int_t iddl=rawReader->GetDDLID();
-           iddl=0; // temporary for test raw data
-           Int_t isddmod=s.GetModuleNumber(iddl,s.GetCarlosId()); 
-           isddmod-=240;  // to have SDD modules from 0 to 259
-           isddmod=s.GetCarlosId(); // temporary for test raw data
-           if(isddmod<nSDDmodules){ 
-             Int_t index=2*isddmod+s.GetChannel(); 
+           Int_t iDDL=rawReader->GetDDLID();
+           Int_t iCarlos=s.GetCarlosId();
+           if(iDDL>=0 && iDDL<kTotDDL && s.IsCompletedModule()==kFALSE){ 
+             Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s.GetChannel(); 
              histo[index]->Fill(s.GetCoord2(),s.GetCoord1(),s.GetSignal());
            }
          }
          delete rawReader;
-         for(Int_t imod=0; imod<nSDDmodules;imod++){
-           for(Int_t isid=0; isid<2;isid++){
-             Int_t index=2*imod+isid;
-             if(iStep==0) base[index]->AddEvent(histo[index]);    
-             if(iStep==1) corr[index]->AddEvent(histo[index]);    
+         for(Int_t iddl=0; iddl<kTotDDL;iddl++){
+           for(Int_t imod=0; imod<kModPerDDL;imod++){
+             for(Int_t isid=0;isid<kSides;isid++){
+               Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
+               if(iStep==0) base[index]->AddEvent(histo[index]);
+               if(iStep==1) corr[index]->AddEvent(histo[index]);
+             }
            }
          }
 
@@ -188,22 +214,19 @@ int main(int argc, char **argv) {
          free(event);
        }
       }
-    
     }
-    
-    for(Int_t imod=0; imod<nSDDmodules;imod++){
-      for(Int_t isid=0; isid<2;isid++){
-       Int_t index=2*imod+isid;
-       if(iStep==0){
-         base[index]->ValidateAnodes();
-         base[index]->WriteToASCII();
-       }
-       if(iStep==1){
-         corr[index]->ValidateAnodes();
-         corr[index]->WriteToASCII();
+
+    if(iStep==0){
+      for(Int_t iddl=0; iddl<kTotDDL;iddl++){
+       for(Int_t imod=0; imod<kModPerDDL;imod++){
+         for(Int_t isid=0;isid<kSides;isid++){
+           Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
+           base[index]->ValidateAnodes();
+           base[index]->WriteToASCII();
+         }
        }
-      }
-    }  
+      }  
+    }
   }
 
   /* write report */
@@ -215,10 +238,14 @@ int main(int argc, char **argv) {
 
 
   TFile *fh=new TFile("SDDbaseHistos.root","RECREATE");
-  for(Int_t imod=0; imod<nSDDmodules;imod++){
-    for(Int_t isid=0; isid<2;isid++){
-      Int_t index=2*imod+isid;
-      corr[index]->WriteToROOT(fh);
+  for(Int_t iddl=0; iddl<kTotDDL;iddl++){
+    for(Int_t imod=0; imod<kModPerDDL;imod++){
+      for(Int_t isid=0;isid<kSides;isid++){
+       Int_t index=kSides*(kModPerDDL*iddl+imod)+isid;
+       corr[index]->ValidateAnodes();
+       corr[index]->WriteToASCII();
+       corr[index]->WriteToROOT(fh);
+      }
     }
   }
   fh->Close();