ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/IPHCalignment2/analysis/histoAlgo.cpp
Revision: 1.6
Committed: Wed Nov 30 08:31:34 2011 UTC (13 years, 5 months ago) by cgoetzma
Branch: MAIN
Changes since 1.5: +15 -15 lines
Log Message: 
chris' mods

File Contents

# User Rev Content
1 econte 1.1 #include "histoAlgo.h"
2     #include <math.h>
3     #include <TColor.h>
4    
5    
6    
7     #define FITMODE 1
8    
9 econte 1.2
10 econte 1.1 void histoAlgo::initialize()
11     {
12 econte 1.2 std::string name;
13     name = "p (GeV) at generator";
14 cgoetzma 1.6 pGen_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,500.);
15 econte 1.2 name = "pt (GeV) at generator";
16 cgoetzma 1.6 ptGen_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,160.);
17 econte 1.2 name = "pseudo-rapidity at generator";
18 cgoetzma 1.6 etaGen_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.,+3.);
19 econte 1.2 name = "azimutal angle (rad) at generator";
20 cgoetzma 1.6 phiGen_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.6,+3.6);
21 cgoetzma 1.5 name = "#mu^{+}/#mu^{-} invariant mass (GeV) at generator";
22 cgoetzma 1.6 mzGen_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str() ,100,65.,115.);
23 econte 1.2 name = "p (GeV)";
24 cgoetzma 1.6 p_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,500.);
25 econte 1.2 name = "pt (GeV)";
26 cgoetzma 1.6 pt_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,0.,160.);
27 econte 1.2 name = "pseudo-rapidity";
28 cgoetzma 1.6 eta_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.,+3.);
29 econte 1.2 name = "azimutal angle (rad)";
30 cgoetzma 1.6 phi_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-3.6,+3.6);
31 cgoetzma 1.5 name = "#mu^{+}/#mu^{-} invariant mass (GeV)";
32 cgoetzma 1.6 mz_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,65.,115.);
33 econte 1.1
34 econte 1.2 name = "p resolution";
35 cgoetzma 1.6 rp_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
36 econte 1.1
37 econte 1.2 name = "pT resolution";
38 cgoetzma 1.6 rpt_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
39 econte 1.1
40 econte 1.2 name = "pseudo-rapidity resolution";
41 cgoetzma 1.6 reta_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
42 econte 1.2 name = "azimutal angle resolution";
43 cgoetzma 1.6 rphi_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str(),100,-0.2,+0.2);
44 econte 1.4 name = "Z mass resolution";
45 cgoetzma 1.6 rmz_histo = new TH1D((name+"_"+algoName).c_str(),name.c_str() ,100,-0.2,+0.2);
46 econte 1.2
47     }