Changes related to the extraction of the V0 finder into a separate class (A. Dainese...
[u/mrichter/AliRoot.git] / ITS / ITSSDDQA.C
1 void ITSSDDQA(char *iFile, Int_t MaxEvts=1000000, Int_t FirstEvt=0) {
2
3 cout << "SDD Quality Assurance Prototype Macro" << endl; 
4
5 const Int_t nSDDmodules= 260;
6 const Int_t imodoffset = 240;
7 const Int_t modtotSDD  = nSDDmodules*2;
8
9 Float_t xi = 0.5;
10 Float_t xf = xi + nSDDmodules;
11 TH1F *modulePattern = new TH1F("patternModule","Modules pattern",nSDDmodules,xi,xf);
12 xf = xi + modtotSDD;
13 TH1F *moduleSidePattern = new TH1F("patternSide","Modules/Side pattern",modtotSDD,xi,xf);
14
15 TH2F *hismap[modtotSDD];  //260 dx e 260 sx  with A, T, Q
16 TH2F *hispop[modtotSDD];  //260 dx e 260 sx  with A, T, Ncounts
17  
18 Char_t *cindex = new Char_t[5];
19 for(Int_t imod=0; imod<nSDDmodules;imod++){
20   for(Int_t isid=0;isid<2;isid++){
21     Int_t index=2*imod+isid;       //260*2 position
22
23     sprintf(cindex,"%d",index+1); // imod,isid);
24     TString sindex((const char *) cindex);
25
26     TString histnam("chargeMap");
27     TString histit("Total Charge, module number ");
28     histnam.Append(sindex);
29     histit.Append(sindex);
30     hismap[index]=new TH2F(histnam.Data(),histit.Data(),256,-0.5,255.5,256,-0.5,255.5);
31       
32     TString hisnam2("countsMap");
33     TString histit2("Number of Counts, module number ");
34     hisnam2.Append(sindex);
35     histit2.Append(sindex);
36     hispop[index]=new TH2F(hisnam2.Data(),histit2.Data(),256,-0.5,255.5,256,-0.5,255.5);
37   
38     /*
39       sprintf(hisnam,"hisprojX%03ds%d",imod,isid);
40       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
41       hisprojX[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
42
43       sprintf(hisnam,"hisprojY%03ds%d",imod,isid);
44       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
45       hisprojY[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
46
47       sprintf(hisnam,"hisprofX%03ds%d",imod,isid);
48       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
49       hisprofX[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
50
51       sprintf(hisnam,"hisprofY%03ds%d",imod,isid);
52       sprintf(histitle,"layer , ladder, module position, channel, %d, %d", imod, isid);
53       hisprofY[index]=new TH2F(hisnam,histitle,256,-0.5,255.5,256,-0.5,255.5);
54       */
55  
56   }
57 }
58 delete [] cindex;
59
60   TString strFile = iFile;
61   strFile += "?EventType=7";
62   AliRawReader *rd = new AliRawReaderDate(strFile.Data(),FirstEvt);  // open run
63   Int_t evCounter = 0;
64
65   //AliITS *itsRun = new AliITS();
66   //  TGeoManager::Import("$ALICE_ROOT/EVE/alice-data/alice_fullgeo.root");
67   /*
68     AliITSInitGeometry *initgeom = new AliITSInitGeometry("AliITSvPPRasymmFMD",2);
69     geom = initgeom->CreateAliITSgeom();
70     delete initgeom;
71   */
72
73   do{       // start loop on events
74     if(++evCounter > MaxEvts) { cout << MaxEvts << " events read, stop" << endl; evCounter--; break; }  
75      cout << "Read Event: " << evCounter+FirstEvt-1 << endl;
76
77     rd->RequireHeader(kFALSE);             
78
79
80     rd->Reset();    // reset the current position to the beginning of the event
81     AliITSRawStreamSDD s(rd);    //This class provides access to ITS SDD digits in raw data.
82     Int_t iddl;
83     Int_t isddmod;
84     Int_t moduleSDD;
85     while(s.Next()){                       //read the next raw digit; returns kFALSE if there is no digit left
86
87         iddl=rd->GetDDLID();
88         isddmod=s.GetModuleNumber(iddl,s.GetCarlosId());  //this is the FEE Carlos
89         //cout<<"DDLID= "<<iddl <<"; Module number= " <<isddmod  <<endl;
90         modulePattern->Fill(isddmod-imodoffset+1);     // 1 to 260
91         moduleSDD=2*(isddmod-imodoffset)+s.GetChannel();
92         moduleSidePattern->Fill(moduleSDD+1);          // 1 to 520
93         //cout << "anode " << s.GetCoord1() << ", time bin: " << s.GetCoord2() << ", charge: " << s.GetSignal() << endl;
94         hismap[moduleSDD]->Fill(s.GetCoord1(), s.GetCoord2(),s.GetSignal() );
95         hispop[moduleSDD]->Fill(s.GetCoord1(), s.GetCoord2() );  
96     }    
97
98   } while(rd->NextEvent()); // end loop on events
99   delete rd;
100
101   cout << "end after " << evCounter << " events" << endl;
102
103    TString oFileName(iFile);
104    oFileName.Append(".root");
105    TFile *oFile = TFile::Open(oFileName,"recreate");
106    modulePattern->Write();
107    moduleSidePattern->Write();
108    for(Int_t i=0; i<modtotSDD; i++){
109      hismap[i]->Write();
110      hispop[i]->Write();
111    }
112    oFile->Close();
113
114    for(Int_t imod=0; imod<nSDDmodules;imod++){
115      for(Int_t isid=0;isid<2;isid++){
116        Int_t index=2*imod+isid;       //260*2 position
117        delete hismap[index];
118        delete hispop[index];
119      }
120    }
121    delete modulePattern;
122    delete moduleSidePattern;
123
124 }