Fork me on GitHub

source: svn/trunk/interface/FuncDef.h@ 25

Last change on this file since 25 was 23, checked in by severine ovyn, 16 years ago

Working on the resolution
File size: 4.7 KB
Line 
1#ifndef __FuncDef__
2#define __FuncDef__
3
4#include "TStyle.h"
5#include "TCanvas.h"
6#include "TH1F.h"
7#include "TLegend.h"
8#include "TPad.h"
9#include "TTree.h"
10#include "TFile.h"
11#include "THStack.h"
12#include "TLine.h"
13#include "TArrow.h"
14#include "TMarker.h"
15#include "TStyle.h"
16#include "TROOT.h"
17
18#include <iomanip>
19#include <vector>
20#include <iostream>
21#include <string>
22#include "TPaveText.h"
23#include "TProfile.h"
24#include "TF1.h"
25#include "TGraph.h"
26#include "TGraphErrors.h"
27
28#include "interface/TdrStyle.h"
29
30using namespace std;
31
32TPaveText * MakeTPave(Float_t xmin,Float_t ymin,Float_t xmax,Float_t ymax,string text)
33{
34
35 TPaveText *leg1 = new TPaveText(xmin,ymin,xmax,ymax,"NDC");
36 leg1->AddText(text.c_str());
37 leg1->SetFillColor(10);
38 leg1->SetBorderSize(0);
39 leg1->SetTextFont(42);
40 leg1->SetTextSize(0.03791469);
41 return leg1;
42
43
44}
45
46
47TLegend * MakeLegend(Int_t FilledColor,Int_t BorderSize,Float_t TextSize,Double_t x1 = 0.60, Double_t y1 = 0.71, Double_t x2= 0.92, Double_t y2 = 0.92)
48{
49 TLegend *leg = new TLegend(x1,y1,x2,y2,NULL,"NDC");
50 leg->SetFillColor(FilledColor);
51 leg->SetBorderSize(BorderSize);
52 leg->SetTextSize(TextSize);
53 leg->SetTextFont(42);
54
55 return leg;
56}
57
58TH1F * MakeNormTH1F(Int_t numBin,Float_t minX,Float_t maxX, TTree * Analyze,string histo, Int_t Lcolor, Int_t Fcolor, Int_t Lstyle,Int_t width=2,bool Log=false)
59{
60 string temp = histo;
61 string nom = histo.erase(0,histo.find(">>")+2);
62 TH1F *h = new TH1F(nom.c_str(),"pt_lept ",numBin,minX,maxX);
63 //Analyze->Draw(temp.c_str(),"(JetPTResol.NonSmearePT >0 && JetPTResol.NonSmearePT < 20)");
64 Analyze->Draw(temp.c_str());
65 h->SetLineColor(Lcolor);
66 h->SetFillColor(Fcolor);
67 h->SetLineStyle(Lstyle);
68 h->SetLineWidth(width);
69 if(Log)gPad->SetLogy();
70 if(h->Integral()!=0)h->Scale(1./h->Integral());
71 return h;
72}
73
74void GaussValues(TTree * Analyze,string histo,double &rms, double &mean, string min, string max)
75{
76 string temp = histo;
77 string mintemp = min;
78 string maxtemp = max;
79
80 string nom = histo.erase(0,histo.find(">>")+2);
81 TH1F *h = new TH1F(nom.c_str(),"",100,0,2);
82 Analyze->Draw(temp.c_str(),mintemp.c_str(),maxtemp.c_str());
83 h->SetMarkerSize(0.6);
84 double MeanFix = h->GetMean();
85 double RangMin = h->GetMean()-h->GetRMS();
86 double RangMax = h->GetMean()+h->GetRMS();
87 TF1 *Gauss = new TF1("Gauss","gausn",RangMin,RangMax);
88 Gauss->FixParameter(1,MeanFix);
89 h->Fit("Gauss","R");
90 h->Fit("Gauss","RI");
91 h->Fit("Gauss","RI");
92 Double_t* params = Gauss->GetParameters();
93 rms=params[2];
94 mean=params[1];
95 h->Draw("P");
96 h->GetXaxis()->SetTitle("E_{T}^{rec}/E_{T}^{MC} [GeV]");
97}
98
99TH1F * HiggsMakeNormTH1F(Float_t Scale, Int_t numBin,Float_t minX,Float_t maxX, TTree * Analyze,string histo, Int_t Lcolor, Int_t Fcolor, Int_t Lstyle,Int_t width=2,bool Log=false)
100{
101 string temp = histo;
102 string nom = histo.erase(0,histo.find(">>")+2);
103 TH1F *h = new TH1F(nom.c_str(),"pt_lept ",numBin,minX,maxX);
104 Analyze->Draw(temp.c_str());
105 h->SetLineColor(Lcolor);
106 h->SetFillColor(Fcolor);
107 h->SetLineStyle(Lstyle);
108 h->SetLineWidth(width);
109 if(Log)gPad->SetLogy();
110 if(h->Integral()!=0)h->Scale(1*Scale);
111 return h;
112}
113
114
115TH1F * MakeNormTH1FSum(Int_t numBin,Float_t minX,Float_t maxX,TTree * analyze1,TTree * analyze2,TTree * analyze3,
116string * histo,float * weight, Int_t Lcolor, Int_t Fcolor, Int_t Lstyle,Int_t width=2,bool Log=false)
117{
118
119 TTree *analyse[3]={analyze1,analyze2,analyze3};
120 TH1F *hq = new TH1F("hq","pt_lept ",numBin,minX,maxX);
121
122 for(unsigned int i=0; i<3; i++)
123 {
124 string temp = histo[i];
125 string nom = histo[i].erase(0,histo[i].find(">>")+2);
126 TH1F *h = new TH1F(nom.c_str(),"pt_lept ",numBin,minX,maxX);
127 analyse[i]->Draw(temp.c_str());
128 h->Scale(weight[i]);
129 hq->Add(h);
130 }
131
132 hq->SetLineColor(Lcolor);
133 hq->SetFillColor(Fcolor);
134 hq->SetLineStyle(Lstyle);
135 hq->SetLineWidth(width);
136 if(Log)gPad->SetLogy();
137 if(hq->Integral()!=0)hq->Scale(1./hq->Integral());
138
139 return hq;
140}
141
142TH1F * HiggsMakeNormTH1FSum(Int_t numBin,Float_t minX,Float_t maxX,TTree * analyze1,TTree * analyze2,TTree * analyze3,
143string * histo,float * weight, Int_t Lcolor, Int_t Fcolor, Int_t Lstyle,Int_t width=2,bool Log=false)
144{
145
146 TTree *analyse[3]={analyze1,analyze2,analyze3};
147 TH1F *hq = new TH1F("hq","pt_lept ",numBin,minX,maxX);
148
149 for(unsigned int i=0; i<3; i++)
150 {
151 string temp = histo[i];
152 string nom = histo[i].erase(0,histo[i].find(">>")+2);
153 TH1F *h = new TH1F(nom.c_str(),"pt_lept ",numBin,minX,maxX);
154 analyse[i]->Draw(temp.c_str());
155 h->Scale(weight[i]);
156 hq->Add(h);
157 }
158
159 hq->SetLineColor(Lcolor);
160 hq->SetFillColor(Fcolor);
161 hq->SetLineStyle(Lstyle);
162 hq->SetLineWidth(width);
163 if(Log)gPad->SetLogy();
164 if(hq->Integral()!=0)hq->Scale(1./hq->Integral());
165
166 return hq;
167}
168
169
170#endif
Note: See TracBrowser for help on using the repository browser.