Using GRP instead of local setters. Creation of AOD for from the RAW ESDs. No need...

[u/mrichter/AliRoot.git] / ITS / AnalyzeSDDGainAllMod.C

diff --git a/ITS/AnalyzeSDDGainAllMod.C b/ITS/AnalyzeSDDGainAllMod.C
index 5d3ab2e45be7f013479a5e7c13b8c5c56f706876..2ed1c2f09544647164cf1ddbd178671189c2a3ec 100644 (file)
--- a/ITS/AnalyzeSDDGainAllMod.C
+++ b/ITS/AnalyzeSDDGainAllMod.C
@@ -25,10 +25,6 @@
 
 void AnalyzeSDDGainAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t lastEv=16, Float_t pascalDAC=100){
 
-  Int_t eqOffset = 256;
-  Int_t DDLid_range = 24;
-  //  Int_t eqOffset = 100;
-  //  Int_t DDLid_range = 1;
   const Int_t kTotDDL=24;
   const Int_t kModPerDDL=12;
   const Int_t kSides=2;
@@ -71,8 +67,6 @@ void AnalyzeSDDGainAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t l
     c0->Clear();
     c0->Divide(4,6,0.001,0.001);
     printf("Event # %d\n",iev);
-    //rd->SelectEvents(7);
-    rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range); 
     rd->Reset();
     for(Int_t iddl=0; iddl<kTotDDL;iddl++){
       for(Int_t imod=0; imod<kModPerDDL;imod++){
@@ -83,11 +77,12 @@ void AnalyzeSDDGainAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t l
       }
     }
     AliITSRawStreamSDD s(rd);
-    rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range); 
     while(s.Next()){
       Int_t iDDL=rd->GetDDLID();
       Int_t iCarlos=s.GetCarlosId();
-      if(iDDL>=0 && iDDL<kTotDDL && s.IsCompletedModule()==kFALSE){ 
+      if(s.IsCompletedModule()) continue;
+      if(s.IsCompletedDDL()) continue;
+      if(iDDL>=0 && iDDL<kTotDDL){ 
        Int_t index=kSides*(kModPerDDL*iDDL+iCarlos)+s.GetChannel(); 
        histo[index]->Fill(s.GetCoord2(),s.GetCoord1(),s.GetSignal());
        isFilled[index]=1;
@@ -113,7 +108,7 @@ void AnalyzeSDDGainAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t l
     printf(" --- OK\n");
   }while(rd->NextEvent()&&iev<=lastEv);
 
-  TH1F *htotgain=new TH1F("htotgain","",100,0.,5.);
+  TH1F *htotgain=new TH1F("htotgain","",100,0.2,4.2);
   TH1F *htotpeakpos=new TH1F("htotpeakpos","",256,-0.5,255.5);
   TH1F *hstatus=new TH1F("hstatus","",2,-0.5,1.5);
 
@@ -227,10 +222,10 @@ void AnalyzeSDDGainAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t l
 
 }
 
-void AnalyzeSDDGainAllMod(Int_t nrun, Int_t n2, Int_t nDDL=0, Int_t firstEv=10, Int_t lastEv=18, Float_t pascalDAC=100){
+void AnalyzeSDDGainAllMod(Int_t nrun, Int_t n2, Char_t* dir="LHC08d_SDD", Int_t nDDL=0, Int_t firstEv=15, Int_t lastEv=20, Float_t pascalDAC=100){
   TGrid::Connect("alien:",0,0,"t");
   Char_t filnam[200];
-  sprintf(filnam,"alien:///alice/data/2008/LHC08c_SDD/%09d/raw/08%09d%03d.10.root",nrun,nrun,n2);
+  sprintf(filnam,"alien:///alice/data/2008/%s/%09d/raw/08%09d%03d.10.root",dir,nrun,nrun,n2);
   printf("Open file %s\n",filnam);
   AnalyzeSDDGainAllMod(filnam,nDDL,firstEv,lastEv,pascalDAC);
 }