Removed unnecessary calls to SetEquipmentId (F. Prino)
[u/mrichter/AliRoot.git] / ITS / ITSSDDQAMaker.C
1 void ITSSDDQAMaker(char *iFile, Int_t MaxEvts=1000000, Int_t FirstEvt=0) {
2
3   //To have Baseline Histos uncomment parts with " // BL!!! " comment
4
5 cout << "SDD Quality Assurance Prototype Macro" << endl; 
6
7 const Int_t nSDDmodules= 260;
8 const Int_t imodoffset = 240;
9 const Int_t modtotSDD  = nSDDmodules*2;
10 const Int_t anode = 256;
11
12 Float_t xi = -0.5;
13 Float_t xf = xi + nSDDmodules;
14 TH1F *hModulePattern = new TH1F("hModulePattern","Modules pattern",nSDDmodules,xi,xf);
15 xf = xi + modtotSDD;
16 TH1F *hModuleSidePattern = new TH1F("hModuleSidePattern","Modules/Side pattern",modtotSDD,xi,xf);
17
18 TH2F *hModuleChargeMap[modtotSDD];  //260 dx e 260 sx  with A, T, Q
19 TH2F *hModuleCountsMap[modtotSDD];  //260 dx e 260 sx  with A, T, Ncounts
20 TH2I *hModuleCarlos = new TH2I("hModuleCarlos","hModuleCarlos",modtotSDD,xi,xf,101,-0.5,100.5);
21 /*
22   TH1F *hModuleSideBL[modtotSDD][anode];                                // BL!!!
23 */ 
24
25
26 //-------histograms definition 
27 Char_t *hisnam = new Char_t[50];
28 Char_t *histit = new Char_t[50];  
29 Char_t *hisnam2 = new Char_t[50];
30 Char_t *histit2 = new Char_t[50];    
31 Char_t *hisnam3 = new Char_t[50];
32  for(Int_t imod=0; imod<nSDDmodules;imod++){
33    for(Int_t isid=0;isid<2;isid++){
34      Int_t index=2*imod+isid;       //260*2 position
35
36      sprintf(hisnam,"chargeMap%d",index);
37      sprintf(histit,"Total Charge, module number %d",index);
38      hModuleChargeMap[index]=new TH2F(hisnam,histit,256,-0.5,255.5,256,-0.5,255.5);  
39
40      sprintf(hisnam2,"countsMap%d",index);
41      sprintf(histit2,"Number of Counts, module number %d",index);
42      hModuleCountsMap[index] = new TH2F(hisnam2,histit2,256,-0.5,255.5,256,-0.5,255.5);
43      /*
44        for(Int_t ianode=0; ianode<anode; ianode++){                         // BL!!!
45        sprintf(hisnam3,"BL_module_%d_%d",index,ianode); 
46        //cout<<hisnam3 <<endl;
47        hModuleSideBL[index][ianode] = new TH1F(hisnam3,hisnam3,256,0.,1024.);
48      }
49      */
50    }
51  }
52
53   TString strFile = iFile;
54   strFile += "?EventType=7";
55   AliRawReader *rd = new AliRawReaderDate(strFile.Data(),FirstEvt);  // open run
56   Int_t evCounter = 0;
57   do{                       // start loop on events
58     if(++evCounter > MaxEvts) { cout << MaxEvts << " events read, stop" << endl; evCounter--; break; }  
59     cout << "Read Event: " << evCounter+FirstEvt-1 << endl;
60
61     rd->RequireHeader(kFALSE);             
62     rd->Reset();                           // reset the current position to the beginning of the event
63  
64     Int_t nSkip = 0;                     // number of skipped signals
65     AliITSRawStreamSDD s(rd);            //This class provides access to ITS SDD digits in raw data.
66     Int_t iddl;
67     Int_t isddmod;
68     Int_t moduleSDD;
69     gStyle->SetPalette(1);
70     while(s.Next()){                     //read the next raw digit; returns kFALSE if there is no digit left
71       if(s.IsCompletedModule()) continue;
72       iddl=rd->GetDDLID()-2; // -2 is temporary for test raw data
73
74         isddmod=s.GetModuleNumber(iddl,s.GetCarlosId());        //this is the FEE Carlos
75         //cout<<"DDLID= "<<iddl <<"; Module number= " <<isddmod  <<endl;
76         if(isddmod >= imodoffset) { 
77           hModulePattern->Fill(isddmod-imodoffset);             // 0 to 259    so  240 to 499
78           moduleSDD=2*(isddmod-imodoffset)+s.GetChannel();
79           hModuleSidePattern->Fill(moduleSDD);                  // 0 to 519
80           hModuleCarlos->Fill(isddmod-imodoffset,s.GetCarlosId());
81           //cout << "anode " << s.GetCoord1() << ", time bin: " << s.GetCoord2() << ", charge: " << s.GetSignal() << endl;        
82           Int_t coord1 = s.GetCoord1();
83           Int_t coord2 = s.GetCoord2();
84           Int_t signal = s.GetSignal();
85           hModuleChargeMap[moduleSDD]->Fill(coord2, coord1,signal);
86           hModuleCountsMap[moduleSDD]->Fill(coord2, coord1 );   
87           //hModuleSideBL[moduleSDD][coord1]->Fill(signal);             // BL  !!!
88         } else {
89           nSkip++;
90         }
91     }    
92     cout << "End of Event " << evCounter+FirstEvt-1 << ", " << nSkip << " wrong module numbers" << endl;
93   } while(rd->NextEvent()); // end loop on events
94   delete rd;
95
96   cout << "end after " << evCounter << " events" << endl;
97   /*
98   TNtuple *Baseline = new TNtuple("Baseline","Baseline","HalfModule:Anode:Mean:RMS");      // BL!!!
99   Float_t meanBL;
100   Float_t rmsBL;  
101   */
102   for(Int_t i=0; i<modtotSDD; i++){   
103     if(hModuleSidePattern->GetBinContent(i+1)){              //check if they're not empty
104       hModuleChargeMap[i]->GetXaxis()->SetTitle("Time Bin");
105       hModuleChargeMap[i]->GetYaxis()->SetTitle("Anode"); 
106       hModuleCountsMap[i]->GetXaxis()->SetTitle("Time Bin");
107       hModuleCountsMap[i]->GetYaxis()->SetTitle("Anode");  
108       /*
109       for(Int_t ianode=0; ianode<anode; ianode++ ){                                      // BL!!!
110         hModuleSideBL[i][ianode]->GetXaxis()->SetTitle("ADC counts");
111         hModuleSideBL[i][ianode]->GetYaxis()->SetTitle("#"); 
112         meanBL = hModuleSideBL[i][ianode]->GetMean();
113         rmsBL = hModuleSideBL[i][ianode]->GetRMS();
114         gaussfitBL = hModuleSideBL[i][ianode]->Fit("gaus");     
115         Baseline->Fill(i,ianode,meanBL,rmsBL);
116       }
117       */
118     }
119   }  
120   
121   hModuleSidePattern->GetXaxis()->SetTitle("2*(Module Number-1)+Side");
122   hModuleSidePattern->GetYaxis()->SetTitle("Counts");  
123   hModulePattern->GetXaxis()->SetTitle("Module Number");  
124   hModulePattern->GetYaxis()->SetTitle("Counts");  
125
126
127   //-------store Histograms
128   cout << "Store Histograms" << endl;
129   TString oFileName(iFile);
130   oFileName.Append(".root");
131   TFile *oFile = TFile::Open(oFileName,"recreate");
132   hModulePattern->Write();
133   hModuleSidePattern->Write();
134   hModuleCarlos->Write();
135   for(Int_t i=0; i<modtotSDD; i++){ 
136       if(hModuleSidePattern->GetBinContent(i+1)){            //check if they're not empty
137         hModuleChargeMap[i]->Write();
138         hModuleCountsMap[i]->Write();     
139         /* 
140         for(Int_t ianode=0; ianode<anode; ianode++ ){                           // BL!!!
141           hModuleSideBL[i][ianode]->Write();
142           Baseline->Write();
143         }
144         */
145       }
146   }
147   
148   oFile->Close();
149   cout << "Clear memory" << endl;
150   for(Int_t imod=0; imod<nSDDmodules;imod++){
151     for(Int_t isid=0;isid<2;isid++){
152       Int_t index=2*imod+isid;       //260*2 position
153       delete hModuleChargeMap[index];
154       delete hModuleCountsMap[index];       
155       /*
156       for(Int_t ianode=0; ianode<anode; ianode++ ){                              // BL!!!
157         delete hModuleSideBL[index][ianode]; 
158         delete Baseline;
159       }
160       */
161     }
162   }
163   delete hModulePattern;
164   delete hModuleSidePattern;
165   delete hModuleCarlos;
166
167 }