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.1 by econte, Fri Nov 25 16:33:02 2011 UTC vs.
Revision 1.8 by cgoetzma, Wed Dec 7 18:31:52 2011 UTC

#	Line 6 | Line 6
6
7		#define FITMODE 1
8
9	+
10		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(),100,0.,500.);
15	+	name = "pt (GeV) at generator";
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(),100,-3.,+3.);
19	+	name = "azimutal angle (rad) at generator";
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(),100,0.,500.);
25	+	name = "pt (GeV)";
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(),100,-3.,+3.);
29	+	name = "azimutal angle (rad)";
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(),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	<	pGen_histo    = new TH1D("p (GeV) at generator","p (GeV) at generator",200,0.,500.);
53	<	ptGen_histo    = new TH1D("pt (GeV) at generator","pt (GeV) at generator",200,0.,160.);
54	<	etaGen_histo    = new TH1D("pseudo-rapidity at generator","pseudo-rapidity at generator",200,-3.,+3.);
55	<	phiGen_histo    = new TH1D("azimutal angle (rad) at generator","azimutal angle (rad) at generator",200,-3.6,+3.6);
56	<	mzGen_histo    = new TH1D("#mu^{#pm} invariant mass (GeV) at generator","#mu^{#pm} invariant mass (GeV) at generator" ,200,65.,115.);
57	<
58	<	p_histo    = new TH1D("p (GeV)","p (GeV)",200,0.,500.);
59	<	pt_histo    = new TH1D("pt (GeV)","pt (GeV)",200,0.,160.);
60	<	eta_histo    = new TH1D("pseudo-rapidity","pseudo-rapidity",200,-3.,+3.);
61	<	phi_histo    = new TH1D("azimutal angle (rad)","azimutal angle (rad)",200,-3.6,+3.6);
62	<	mz_histo    = new TH1D("#mu^{#pm} invariant mass (GeV)","#mu^{#pm} invariant mass (GeV)" ,200,65.,115.);
63	<
64	<	rp_histo    = new TH1D("p resolution","p resolution",200,-0.2,+0.2);
65	<	rpt_histo   = new TH1D("pT resolution","pT resolution",200,-0.2,+0.2);
66	<	reta_histo   = new TH1D("pseudo-rapidity resolution","pseudo-rapidity resolution",200,-0.2,+0.2);
67	<	rphi_histo   = new TH1D("azimutal angle resolution","azimutal angle resolution",200,-0.2,+0.2);
68	<	rmz_histo    = new TH1D("#mu^{#pm} invariant mass (GeV)","#mu^{#pm} invariant mass (GeV)" ,200,-0.2,+0.2);
52	>	name = "pT resolution";
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(),100,-0.003,+0.003);
57	>	name = "azimutal angle resolution";
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	<	}
71	>
72	>	}

Diff Legend

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