Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
TFile *f;
int ipad;
void plotQA(TString fName="./JpsiGamKsKlKK.root"){
f = new TFile(fName, "READ");
if(f==0){
cout << "Could not open file " << fName << endl;
return;
}
ipad=1;
TCanvas* c1 = new TCanvas("c1","c1",450,600);
c1->Divide(2,3);
TCanvas* c2 = new TCanvas("c2","c2",450,600);
c2->Divide(2,3);
ipad=1;
c1->cd(ipad);
plotHist("_PhiPhiMass");
c1->cd(ipad);
plotHist("_costGamCm");
c1->cd(ipad);
plotHist("_costKs_KsKlHeli");
c1->cd(ipad);
plotHist("_costKp_KpKmHeli");
c1->cd(ipad);
plotHist("_phiKs_KsKlHeli");
c1->cd(ipad);
plotHist("_phiKp_KpKmHeli");
ipad=1;
c2->cd(ipad);
plotHist("_costPhi_KpKm");
c2->cd(ipad);
plotHist("_phiPhi_KpKm");
c2->cd(ipad);
plotHist("_chi");
c1->Print("qaPlots.ps(");
c2->Print("qaPlots.ps)");
}
void plotHist(TString base){
TH1F* hd = (TH1F*) f->Get(base+"DataHist");
TH1F* hf = (TH1F*) f->Get(base+"FittedHist");
hd->SetMinimum(0);
hf->SetMinimum(0);
hf->SetLineWidth(2);
hf->SetLineColor(kRed);
hd->Draw("E");
hf->Draw("same");
ipad++;
return;
}