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.7 by cgoetzma, Wed Dec 7 08:34:02 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
12	–	pGen_histo    = new TH1D("p (GeV) at generator","p (GeV) at generator",200,0.,500.);
13	–	ptGen_histo    = new TH1D("pt (GeV) at generator","pt (GeV) at generator",200,0.,160.);
14	–	etaGen_histo    = new TH1D("pseudo-rapidity at generator","pseudo-rapidity at generator",200,-3.,+3.);
15	–	phiGen_histo    = new TH1D("azimutal angle (rad) at generator","azimutal angle (rad) at generator",200,-3.6,+3.6);
16	–	mzGen_histo    = new TH1D("#mu^{#pm} invariant mass (GeV) at generator","#mu^{#pm} invariant mass (GeV) at generator" ,200,65.,115.);
17	–
18	–	p_histo    = new TH1D("p (GeV)","p (GeV)",200,0.,500.);
19	–	pt_histo    = new TH1D("pt (GeV)","pt (GeV)",200,0.,160.);
20	–	eta_histo    = new TH1D("pseudo-rapidity","pseudo-rapidity",200,-3.,+3.);
21	–	phi_histo    = new TH1D("azimutal angle (rad)","azimutal angle (rad)",200,-3.6,+3.6);
22	–	mz_histo    = new TH1D("#mu^{#pm} invariant mass (GeV)","#mu^{#pm} invariant mass (GeV)" ,200,65.,115.);
23	–
24	–	rp_histo    = new TH1D("p resolution","p resolution",200,-0.2,+0.2);
25	–	rpt_histo   = new TH1D("pT resolution","pT resolution",200,-0.2,+0.2);
26	–	reta_histo   = new TH1D("pseudo-rapidity resolution","pseudo-rapidity resolution",200,-0.2,+0.2);
27	–	rphi_histo   = new TH1D("azimutal angle resolution","azimutal angle resolution",200,-0.2,+0.2);
28	–	rmz_histo    = new TH1D("#mu^{#pm} invariant mass (GeV)","#mu^{#pm} invariant mass (GeV)" ,200,-0.2,+0.2);
51
52	<	}
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	>	}

Diff Legend

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