ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/IPHCalignment2/analysis/histoAlgo.cpp

Comparing UserCode/IPHCalignment2/analysis/histoAlgo.cpp (file contents):
Revision 1.3 by econte, Sun Nov 27 11:11:49 2011 UTC vs.
Revision 1.9 by cgoetzma, Thu Dec 8 10:23:04 2011 UTC

#	Line 11 | Line 11 | void histoAlgo::initialize()
11		{
12		std::string name;
13		name = "p (GeV) at generator";
14	<	pGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,0.,500.);
14	>	pGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,500.);
15		name = "pt (GeV) at generator";
16	<	ptGen_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,0.,160.);
16	>	ptGen_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,160.);
17		name = "pseudo-rapidity at generator";
18	<	etaGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-3.,+3.);
18	>	etaGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.,+3.);
19		name = "azimutal angle (rad) at generator";
20	<	phiGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-3.6,+3.6);
21	<	name = "#mu^{#pm} invariant mass (GeV) at generator";
22	<	mzGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str() ,200,65.,115.);
20	>	phiGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.6,+3.6);
21	>	name = "#mu^{+}/#mu^{-} invariant mass (GeV) at generator";
22	>	mzGen_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str() ,100,65.,115.);
23		name = "p (GeV)";
24	<	p_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,0.,500.);
24	>	p_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,500.);
25		name = "pt (GeV)";
26	<	pt_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,0.,160.);
26	>	pt_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,160.);
27		name = "pseudo-rapidity";
28	<	eta_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-3.,+3.);
28	>	eta_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.,+3.);
29		name = "azimutal angle (rad)";
30	<	phi_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-3.6,+3.6);
31	<	name = "#mu^{#pm} invariant mass (GeV)";
32	<	mz_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,65.,115.);
30	>	phi_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.6,+3.6);
31	>	name = "#mu^{+}/#mu^{-} invariant mass (GeV)";
32	>	mz_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,65.,115.);
33
34		name = "p resolution";
35	<	rp_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-0.2,+0.2);
35	>	rp_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
36	>	name = "px resolution";
37	>	rpx_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
38	>	name = "py resolution";
39	>	rpy_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
40	>	name = "pz resolution";
41	>	rpz_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
42	>
43	>	name = "px resolution over px";
44	>	rpxOverpx_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.,+1);
45	>	name = "py resolution over py";
46	>	rpyOverpy_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.,+1);
47	>	name = "pz resolution over pz";
48	>	rpzOverpz_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.,+1);
49	>
50	>
51
52		name = "pT resolution";
53	<	rpt_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-0.2,+0.2);
53	>	rpt_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
54
55		name = "pseudo-rapidity resolution";
56	<	reta_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-0.2,+0.2);
56	>	reta_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.003,+0.003);
57		name = "azimutal angle resolution";
58	<	rphi_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),200,-0.2,+0.2);
59	<	name = "Z mass resolution";
60	<	rmz_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str() ,200,-0.2,+0.2);
58	>	rphi_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.002,+0.002);
59	>	name = "Z mass resolution";
60	>	rmz_histo    = new TH1D((name+"_"+algoName).c_str(),name.c_str() ,100,-0.2,+0.2);
61	>
62	>	name = "p resolution for Z";
63	>	Z_rp_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
64	>	name = "pT resolution for Z";
65	>	Z_rpt_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
66	>	name = "pseudo-rapidity resolution for Z";
67	>	Z_reta_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.05,+0.05);
68	>	name = "azimutal angle resolution for Z";
69	>	Z_rphi_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.05,+0.05);
70	>
71	>	// Alessio plots
72	>
73	>	name = "q/p";
74	>	QoverP_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.1,+0.1);
75	>
76	>	name = "lambda";
77	>	lambda_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-1.5,+1.5);
78	>
79	>	name = "dxy";
80	>	dxy_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.02,+0.02);
81	>
82	>	name = "dsz";
83	>	dsz_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-10.,+10.);
84	>
85	>	name = "dxy0";
86	>	dxy0_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.6,+0.6);
87	>
88	>	name = "dsz0";
89	>	dsz0_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-10.,+10.);
90	>
91	>
92	>	name = "q/p resolution";
93	>	rQoverP_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-1,+1);
94	>
95	>	name = "lambda resolution";
96	>	rlambda_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.005,+0.005);
97	>
98	>	name = "dxy resolution";
99	>	rdxy_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-100,+100);
100	>
101	>	name = "dsz resolution";
102	>	rdsz_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.04,+0.04);
103	>
104	>	name = "dxy0 resolution";
105	>	rdxy0_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-100,+100);
106	>
107	>	name = "dsz0 resolution";
108	>	rdsz0_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.04,+0.04);
109	>
110	>
111	>
112	>	name = "q/p";
113	>	QoverP_error_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.1,+0.1);
114	>
115	>	name = "lambda";
116	>	lambda_error_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-1.5,+1.5);
117	>
118	>	name = "dxy";
119	>	dxy_error_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.02,+0.02);
120	>
121	>	name = "dsz";
122	>	dsz_error_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-10.,+10.);
123	>
124	>	name = "dxy0";
125	>	dxy0_error_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.6,+0.6);
126	>
127	>	name = "dsz0";
128	>	dsz0_error_histo   = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-10.,+10.);
129	>
130
131		}

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines