Obsolete code commented out (Annalisa)
[u/mrichter/AliRoot.git] / HMPID / AliHMPIDQADataMaker.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16
17 /* $Id$ */
18
19 //---
20 //  Produces the data needed to calculate the quality assurance. 
21 //  All data must be mergeable objects.
22 //  A. Mastroserio
23 //---
24
25 // --- ROOT system ---
26 #include <TClonesArray.h>
27 #include <TFile.h> 
28 #include <TH1F.h> 
29 #include <TH2F.h>
30 #include <TH1I.h> 
31 #include <TDirectory.h>
32 #include <Riostream.h>
33 // --- Standard library ---
34
35 // --- AliRoot header files ---
36 #include "AliESDCaloCluster.h"
37 #include "AliESDEvent.h"
38 #include "AliLog.h"
39 #include "AliHMPIDDigit.h"
40 #include "AliHMPIDHit.h"
41 #include "AliHMPIDCluster.h"
42 #include "AliHMPIDQADataMaker.h"
43
44 ClassImp(AliHMPIDQADataMaker)
45            
46 //____________________________________________________________________________ 
47   AliHMPIDQADataMaker::AliHMPIDQADataMaker() : 
48   AliQADataMaker(AliQA::GetDetName(AliQA::kHMPID), "HMPID Quality Assurance Data Maker"),
49   fhHitQdc(0x0), 
50   fhSDigits(0x0),
51   fhDigPcEvt(0x0),
52   fhDigChEvt(0x0),
53   fhDigQ(0x0),
54   fhCluEvt(0x0),
55   fhCluChi2(0x0),
56   fhCluQ(0x0),
57   fhCluFlg(0x0), 
58   fhCluSize(0x0),  
59   fhMipCluSize(0x0),
60   fhCkovP(0x0),
61   fhSigP(0x0),
62   fhMipXY(0x0),
63   fhDifXY(0x0)
64 {
65   // ctor
66   for(Int_t i=0; i<7; i++) fhHitMap[i]=0x0;
67   for(Int_t j=0; j<5; j++) fhPid[j]=0x0;
68 //   fDetectorDir = fOutput->GetDirectory(GetName()) ;  
69 //   if (!fDetectorDir) 
70 //     fDetectorDir = fOutput->mkdir(GetName()) ;  
71 }
72
73 //____________________________________________________________________________ 
74 AliHMPIDQADataMaker::AliHMPIDQADataMaker(const AliHMPIDQADataMaker& qadm) :
75   AliQADataMaker(), 
76   fhHitQdc(qadm.fhHitQdc), 
77   fhSDigits(qadm.fhSDigits),
78   fhDigPcEvt(qadm.fhDigPcEvt),
79   fhDigChEvt(qadm.fhDigChEvt),
80   fhDigQ(qadm.fhDigQ),
81   fhCluEvt(qadm.fhCluEvt),
82   fhCluChi2(qadm.fhCluChi2),
83   fhCluQ(qadm.fhCluQ),
84   fhCluFlg(qadm.fhCluFlg),
85   fhCluSize(qadm.fhCluSize),
86   fhMipCluSize(qadm.fhMipCluSize),
87   fhCkovP(qadm.fhCkovP),
88   fhSigP(qadm.fhSigP),
89   fhMipXY(qadm.fhMipXY),
90   fhDifXY(qadm.fhDifXY)
91 {
92   //copy ctor 
93   for(Int_t i=0; i<7; i++) fhHitMap[i]=qadm.fhHitMap[i];
94   for(Int_t j=0; j<5; j++) fhPid[j]=qadm.fhPid[j];
95
96   SetName((const char*)qadm.GetName()) ; 
97   SetTitle((const char*)qadm.GetTitle()); 
98 }
99
100 //__________________________________________________________________
101 AliHMPIDQADataMaker& AliHMPIDQADataMaker::operator = (const AliHMPIDQADataMaker& qadm )
102 {
103   // Equal operator.
104   this->~AliHMPIDQADataMaker();
105   new(this) AliHMPIDQADataMaker(qadm);
106   return *this;
107 }
108  
109 //____________________________________________________________________________ 
110 void AliHMPIDQADataMaker::InitHits()
111 {
112   // create Hits histograms in Hits subdir
113      fhHitQdc=new TH1F("HitQdc","HMPID Hit Qdc all chamber;QDC",500,0,4000);
114      for(Int_t iCh=0;iCh<7;iCh++) fhHitMap[iCh]=new TH2F(Form("HMPID HitMap%i",iCh),Form("Ch%i;x_{Hit};y_{Hit}",iCh),162,-1,161,146,-1,145);   
115 }
116
117 //____________________________________________________________________________ 
118 void AliHMPIDQADataMaker::InitDigits()
119 {
120   // create Digits histograms in Digits subdir
121       fhDigPcEvt=new TH1F("hDigPcEvt","PC occupancy",156,-1,77);
122       fhDigChEvt=new TH1F("hDigChEvt","Chamber occupancy",32,-1,7);
123       fhDigQ  =new TH1F("Q        "," digit charge               ",3000,0,3000);
124 }
125
126 //____________________________________________________________________________ 
127 void AliHMPIDQADataMaker::InitSDigits()
128 {
129   // create SDigits histograms in SDigits subdir
130       fhSDigits     = new TH1F("hHmpidSDigits",    "SDigits Q  distribution in HMPID",  500, 0., 5000.) ; 
131 }
132
133 //____________________________________________________________________________ 
134
135 void AliHMPIDQADataMaker::InitRecPoints()
136 {
137   // create cluster histograms in RecPoint subdir
138       fhCluEvt=new TH1F("CluPerEvt","# clusters per chamber",16,-1,7);
139       fhCluChi2  =new TH1F("CluChi2"  ,"Chi2 "               ,1000,0,100);
140       fhCluQ   =new TH1F("CluQ"   ,"Cluster charge"        ,3000,0,3000);
141       fhCluFlg   =new TH1F("CluFlg"   ,"Cluster flag"        ,14,-1.5,12.5);
142       fhCluSize  =new TH1F("CluSize"  ,"Raw cluster size    ",100,0,100);
143       fhMipCluSize =new TH1F("MipCluSize"  ,"Mip cluster size    ",100,0,100);
144 }
145 //____________________________________________________________________________
146 void AliHMPIDQADataMaker::InitESDs()
147 {
148   //create ESDs histograms in ESDs subdir
149      fhCkovP  = new TH2F("CkovP" , "#theta_{c}, [rad];P, [GeV]"   , 150,   0,  7  ,100, 0, 1)   ;
150      fhSigP   = new TH2F("SigP"  ,"#sigma_{#theta_c} [mrad];[GeV]", 150,   0,  7  ,100, 0, 1)   ;
151      fhMipXY  = new TH2F("MipXY" ,"mip position"                  , 260,   0,130  ,252, 0,126)  ;
152      fhDifXY  = new TH2F("DifXY" ,"diff"                          , 200, -10, 10  ,200,-10,10)  ;
153      fhPid[0] = new TH1F("PidE" ,"PID: e yellow #mu magenta"  ,100,0,1)                         ;
154      fhPid[1] = new TH1F("PidMu","pid of #mu"                 ,100,0,1)                         ;
155      fhPid[2] = new TH1F("PidPi","PID: #pi red K green p blue",100,0,1)                         ;
156      fhPid[3] = new TH1F("PidK" ,"pid of K"                   ,100,0,1)                         ;
157      fhPid[4] = new TH1F("PidP" ,"pid of p"                   ,100,0,1)                         ;
158 }
159
160 //____________________________________________________________________________
161 void AliHMPIDQADataMaker::MakeHits(TObject * data)
162 {
163   //fills QA histos for Hits
164   TClonesArray * hits = dynamic_cast<TClonesArray *>(data) ; 
165   if (!hits){
166     AliError("Wrong type of hits container") ; 
167   } else {
168     TIter next(hits); 
169     AliHMPIDHit * hit ; 
170     while ( (hit = dynamic_cast<AliHMPIDHit *>(next())) ) {
171       if(hit->Pid()<500000) fhHitQdc->Fill(hit->Q()) ;
172       if(hit->Pid()<500000) fhHitMap[hit->Ch()]->Fill(hit->LorsX(),hit->LorsY());
173     }
174   } 
175 }
176
177 //____________________________________________________________________________
178 void AliHMPIDQADataMaker::MakeDigits( TObject * data)
179 {
180   //fills QA histos for Digits
181   TObjArray *chambers = dynamic_cast<TObjArray*>(data);
182   if ( !chambers) {
183     AliError("Wrong type of digits container") ; 
184   } else {
185     for(Int_t i =0; i< chambers->GetEntries(); i++)
186       {
187         TClonesArray * digits = dynamic_cast<TClonesArray*>(chambers->At(i)); 
188         fhDigChEvt->Fill(i,digits->GetEntriesFast()/(48.*80.*6.));
189         TIter next(digits); 
190         AliHMPIDDigit * digit; 
191         while ( (digit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
192           fhDigPcEvt->Fill(10.*i+digit->Pc(),1./(48.*80.));
193           fhDigQ->Fill(digit->Q());
194         }  
195       }
196   }
197 }
198
199 //____________________________________________________________________________
200 void AliHMPIDQADataMaker::MakeSDigits( TObject * data)
201 {
202   //fills QA histos for SDigits
203   TClonesArray * sdigits = dynamic_cast<TClonesArray *>(data) ; 
204   if (!sdigits) {
205     AliError("Wrong type of sdigits container") ; 
206   } else {
207     AliHMPIDDigit *ref = (AliHMPIDDigit *)sdigits->At(0);
208     Float_t zero = ref->GetTrack(0); 
209     TIter next(sdigits) ; 
210     AliHMPIDDigit * sdigit ; 
211     while ( (sdigit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
212       fhSDigits->Fill(sdigit->Q()) ;
213       if(zero == sdigit->GetTrack(0)) continue;
214       else zero = sdigit->GetTrack(0);
215     } 
216   }
217 }
218
219 //____________________________________________________________________________
220 void AliHMPIDQADataMaker::MakeRecPoints(TTree * clustersTree)
221 {
222   //fills QA histos for clusters
223
224   TClonesArray *clusters = new TClonesArray("AliHMPIDCluster");
225   for(int i=AliHMPIDParam::kMinCh;i<=AliHMPIDParam::kMaxCh;i++){
226     TBranch *branch = clustersTree->GetBranch(Form("HMPID%d",i));
227     branch->SetAddress(&clusters);
228     branch->GetEntry(0);
229
230     fhCluEvt->Fill(i,clusters->GetEntries());
231     TIter next(clusters);
232     AliHMPIDCluster *clu;
233     while ( (clu = dynamic_cast<AliHMPIDCluster *>(next())) ) {;
234       fhCluFlg->Fill(clu->Status());  fhCluChi2->Fill(clu->Chi2());  fhCluSize->Fill(clu->Size());
235       fhCluQ->Fill(clu->Q()); 
236       Int_t qCut=100;
237       if(clu->Q()>qCut) {
238         fhMipCluSize->SetTitle(Form("Mip cluster size at a Qcut = %i ADC",qCut));
239         fhMipCluSize->Fill(clu->Size());
240       }
241     }
242   }
243
244   clusters->Delete();
245   delete clusters;
246 }
247
248 //____________________________________________________________________________
249 void AliHMPIDQADataMaker::MakeESDs(AliESDEvent * esd)
250 {
251   //fills QA histos for ESD
252   for(Int_t iTrk = 0 ; iTrk < esd->GetNumberOfTracks() ; iTrk++){
253     AliESDtrack *pTrk = esd->GetTrack(iTrk) ;
254     fhCkovP->Fill(pTrk->GetP(),pTrk->GetHMPIDsignal());
255     fhSigP->Fill( pTrk->GetP(),TMath::Sqrt(pTrk->GetHMPIDchi2()));
256     Float_t xm,ym; Int_t q,np;  
257     pTrk->GetHMPIDmip(xm,ym,q,np);                       //mip info
258     fhMipXY->Fill(xm,ym);
259     Float_t xRad,yRad,th,ph;        
260     pTrk->GetHMPIDtrk(xRad,yRad,th,ph);              //track info at the middle of the radiator
261     Float_t xPc = xRad+9.25*TMath::Tan(th)*TMath::Cos(ph); // temporar: linear extrapol (B=0!)
262     Float_t yPc = yRad+9.25*TMath::Tan(th)*TMath::Sin(ph); // temporar:          "
263     fhDifXY->Fill(xm-xPc,ym-yPc); //track info
264     Double_t pid[5] ;      pTrk->GetHMPIDpid(pid) ;
265     for(Int_t i = 0 ; i < 5 ; i++) fhPid[i]->Fill(pid[i]) ;
266   }
267 }
268