-
Notifications
You must be signed in to change notification settings - Fork 0
/
plothelper.cxx
67 lines (48 loc) · 2.1 KB
/
plothelper.cxx
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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
{
// // ====================================== Pythia spectrum
// TFile * in=new TFile("Results/Groom_Rhic.Zg.root","READ");
// TH1D* h= (TH1D*) in->Get("LeadPt");
// h->SetAxisRange(0, 60,"x");
// TCanvas* c = new TCanvas( "c","",400, 500 );
// gPad->SetLogy();
// h->Draw();
// gPad->SaveAs("tmp/PythiaSpec.png");
// // ====================================== rho
// // TFile * in=new TFile("Results/Groom_AuAu.root","READ");
// // TFile * in=new TFile("Results/Groom_Rhic.root","READ");
// TFile * in=new TFile("../PaperAj/BasicAj/AjResults/Fresh_NicksList_HC100_AuAu.root","READ");
// TCanvas* c = new TCanvas( "c","",500, 500 );
// TH1D* h= new TH1D ( "h","", 60, 0, 120 );
// ResultTree->Draw("rho>>h","refmult>=269","");
// h->SetTitle(";#rho [GeV/c]");
// h->Draw();
// gPad->SaveAs("tmp/Rho.png");
// // ====================================== rhoerr
// TFile * in=new TFile("../PaperAj/BasicAj/AjResults/Fresh_NicksList_HC100_AuAu.root","READ");
// TCanvas* c = new TCanvas( "c","",500, 500 );
// TH1D* h= new TH1D ( "h","", 40, 0, 20);
// ResultTree->Draw("rhoerr>>h","refmult>=269","");
// h->SetTitle(";#sigma ( #rho ) [GeV/c]");
// h->Draw();
// gPad->SaveAs("tmp/RhoErr.png");
// // ====================================== delta pt
// TFile * in=new TFile("Results/Groom_Rhic.root","READ");
// TCanvas* c = new TCanvas( "c","",500, 500 );
// gPad->SetLogy();
// TH1D* h= new TH1D ( "h","", 60, -30, 30 );
// ResultTree->Draw("-Jets[0].Pt() + EmbJets[0].Pt()>>h");
// h->SetTitle(";#Delta p_{T} [GeV/c]");
// h->Draw();
// gPad->SaveAs("tmp/DeltaPt.png");
// // ====================================== delta pt
// TFile * in=new TFile("Results/GroomUnfoldingTest.root","READ");
// TCanvas* c = new TCanvas( "c","",550, 500 );
// gPad->SetLogz();
// TH2D* h = (TH2D*) LeadPtResponse->Hresponse();
// h->SetTitle("; p_{T}^{meas} [GeV/c];p_{T}^{meas} [GeV/c]");
// h->SetAxisRange(10, 40, "x");
// h->SetAxisRange( 0, 40, "y");
// gStyle->SetOptStat(0);
// h->Draw("colz");
// gPad->SaveAs("tmp/response.png");
}