]>
Commit | Line | Data |
---|---|---|
08cb00f2 | 1 | void ITSlego(Float_t nchtheta, Float_t nchphi) { |
2 | // macro to visualize the lego plots generated by gAlive->RunLego | |
3 | ||
4 | gROOT->Reset(); | |
5 | TFile *file = new TFile("galice.root"); | |
6 | ||
7 | ||
8 | /* | |
9 | Float_t theta = 10; | |
10 | Float_t phi = 170; | |
11 | Int_t ncont = 50; | |
12 | ||
13 | TCanvas *cradl = new TCanvas("cradl","radl",10,10,800,600); | |
14 | cradl->SetTheta(theta); | |
15 | cradl->SetPhi(phi); | |
16 | TH2F *hradl = (TH2F*)file->Get("hradl"); | |
17 | hradl->SetStats(kFALSE); | |
18 | hradl->GetXaxis()->SetTitle("Phi (degrees)"); | |
19 | hradl->GetYaxis()->SetTitle("Theta (degrees)"); | |
20 | hradl->SetFillColor(2); | |
21 | hradl->SetContour(ncont); | |
22 | // hradl->SetMaximum(1); | |
23 | hradl->Draw("colz"); | |
24 | // hradl->Draw("lego2sphe"); | |
25 | */ | |
26 | ||
27 | TCanvas *cradlx = new TCanvas("cradl1","radl",50,50,800,600); | |
28 | hradl->ProjectionX(); | |
29 | hradl_px->SetStats(kFALSE); | |
30 | // hradl_px->SetOptLogY(); | |
31 | // hradl_px->SetMinimum(0.001); | |
32 | // hradl_px->SetMaximum(1); | |
33 | hradl_px->Scale(1./nchtheta); | |
34 | hradl_px->GetXaxis()->SetTitle("Phi (degrees)"); | |
35 | hradl_px->GetYaxis()->SetTitle("X/X0"); | |
36 | hradl_px->Draw(); | |
37 | // cout << "Average over Phi: " << hradl_px->GetSumOfWeights()/nchphi << " X/X0" << endl; | |
38 | ||
39 | ||
40 | TCanvas *cradly = new TCanvas("cradl2","radl",100,100,800,600); | |
41 | hradl->ProjectionY(); | |
42 | hradl_py->SetStats(kFALSE); | |
43 | // hradl_py->SetOptLogY(); | |
44 | // hradl_py->SetMinimum(0.001); | |
45 | // hradl_py->SetMaximum(1); | |
46 | hradl_py->Scale(1./nchphi); | |
47 | hradl_py->GetXaxis()->SetTitle("Theta (degrees)"); | |
48 | hradl_py->GetYaxis()->SetTitle("X/X0"); | |
49 | hradl_py->Draw(); | |
50 | // cout << "Average over Theta: " << hradl_py->GetSumOfWeights()/nchtheta << " X/X0" << endl; | |
51 | cout << "Average: " << hradl_py->GetSumOfWeights()/nchtheta << " X/X0" << endl; | |
52 | } |