-
Notifications
You must be signed in to change notification settings - Fork 3
Expand file tree
/
Copy pathdraw_B_monte_carlo.c
More file actions
117 lines (101 loc) · 3.57 KB
/
draw_B_monte_carlo.c
File metadata and controls
117 lines (101 loc) · 3.57 KB
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
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
{
TFile *tfile = new TFile("B_monte_carlo_feeding.root");
Int_t rebin = 1;
Int_t rebin_realtime = 1;
h_A_realtime->SetLineColor(kBlack);
h_D_realtime->SetLineColor(kBlack);
h_D_cyctime ->SetLineColor(kBlack);
h_T1_realtime->SetLineColor(kGreen);
h_U1_realtime->SetLineColor(kGreen);
h_DT1_realtime->SetLineColor(kGreen);
h_DU1_realtime->SetLineColor(kGreen);
h_DT1_cyctime ->SetLineColor(kGreen);
h_DU1_cyctime ->SetLineColor(kGreen);
h_T2_realtime->SetLineColor(kBlue);
h_U2_realtime->SetLineColor(kBlue);
h_DT2_realtime->SetLineColor(kBlue);
h_DU2_realtime->SetLineColor(kBlue);
h_DT2_cyctime ->SetLineColor(kBlue);
h_DU2_cyctime ->SetLineColor(kBlue);
h_T3_realtime->SetLineColor(kRed);
h_U3_realtime->SetLineColor(kRed);
h_DT3_realtime->SetLineColor(kRed);
h_DU3_realtime->SetLineColor(kRed);
h_DT3_cyctime ->SetLineColor(kRed);
h_DU3_cyctime ->SetLineColor(kRed);
h_U1_realtime->SetLineStyle(7);
h_U2_realtime->SetLineStyle(7);
h_U3_realtime->SetLineStyle(7);
h_DU1_realtime->SetLineStyle(7);
h_DU2_realtime->SetLineStyle(7);
h_DU3_realtime->SetLineStyle(7);
h_DU1_cyctime->SetLineStyle(7);
h_DU2_cyctime->SetLineStyle(7);
h_DU3_cyctime->SetLineStyle(7);
h_D_realtime ->Rebin(rebin_realtime);
h_DT1_realtime ->Rebin(rebin_realtime);
h_DT2_realtime ->Rebin(rebin_realtime);
h_DT3_realtime ->Rebin(rebin_realtime);
h_DU1_realtime ->Rebin(rebin_realtime);
h_DU2_realtime ->Rebin(rebin_realtime);
h_DU3_realtime ->Rebin(rebin_realtime);
h_D_cyctime ->Rebin(rebin);
h_DT1_cyctime ->Rebin(rebin);
h_DT2_cyctime ->Rebin(rebin);
h_DT3_cyctime ->Rebin(rebin);
h_DU1_cyctime ->Rebin(rebin);
h_DU2_cyctime ->Rebin(rebin);
h_DU3_cyctime ->Rebin(rebin);
/*
TCanvas *c_pops_realtime = new TCanvas("c_pops_realtime","Ion populations versus real time",945,600);
h_A_realtime->Draw();
h_T1_realtime->Draw("SAME");
h_T2_realtime->Draw("SAME");
h_T3_realtime->Draw("SAME");
h_U1_realtime->Draw("SAME");
h_U2_realtime->Draw("SAME");
h_U3_realtime->Draw("SAME");
TCanvas *c_decays_realtime = new TCanvas("c_decays_realtime","Decays versus cycle time",945,600);
h_D_realtime->Draw();
h_DT1_realtime->Draw("SAME");
h_DT2_realtime->Draw("SAME");
h_DT3_realtime->Draw("SAME");
h_DU1_realtime->Draw("SAME");
h_DU2_realtime->Draw("SAME");
h_DU3_realtime->Draw("SAME");
*/
TCanvas *c_decays_cyctime = new TCanvas("c_decays_cyctime","Decays versus cycle time",945,600);
h_D_cyctime->Draw();
h_DT1_cyctime->Draw("SAME");
h_DT2_cyctime->Draw("SAME");
h_DT3_cyctime->Draw("SAME");
h_DU1_cyctime->Draw("SAME");
h_DU2_cyctime->Draw("SAME");
h_DU3_cyctime->Draw("SAME");
// X and Y pops
h_DX2_cyctime ->Rebin(rebin);
h_DX3_cyctime ->Rebin(rebin);
h_DY2_cyctime ->Rebin(rebin);
h_DY3_cyctime ->Rebin(rebin);
h_DX2_cyctime->SetLineStyle(8);
h_DY2_cyctime->SetLineStyle(3);
h_DX2_cyctime->SetLineColor(kBlue);
h_DY2_cyctime->SetLineColor(kBlue);
h_DX3_cyctime->SetLineStyle(8);
h_DY3_cyctime->SetLineStyle(3);
h_DX3_cyctime->SetLineColor(kRed);
h_DY3_cyctime->SetLineColor(kRed);
TH1D *h_DU2_cyctime_1 = (TH1D*)h_DU2_cyctime->Rebin(rebin,"h_DU2_cyctime_1");
TH1D *h_DU3_cyctime_1 = (TH1D*)h_DU3_cyctime->Rebin(rebin,"h_DU3_cyctime_1");
h_DU2_cyctime->SetLineStyle(7);
h_DU3_cyctime->SetLineStyle(7);
h_DU2_cyctime->SetLineColor(kBlue);
h_DU3_cyctime->SetLineColor(kRed);
TCanvas *c_feeding_cyctime = new TCanvas("c_feeding_cyctime","Decays from pops X and Y versus cycle time",945,600);
h_DU3_cyctime_1->Draw();
h_DX3_cyctime->Draw("SAME");
h_DY3_cyctime->Draw("SAME");
h_DU2_cyctime_1->Draw("SAME");
h_DX2_cyctime->Draw("SAME");
h_DY2_cyctime->Draw("SAME");
}