Fork me on GitHub

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

Last change on this file since 20 was 17, checked in by Xavier Rouby, 16 years ago

first commit
File size: 4.6 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(),"",50,0.5,1.5);
82 Analyze->Draw(temp.c_str(),mintemp.c_str(),maxtemp.c_str());
83 h->SetMarkerSize(0.6);
84 TF1 *Gauss = new TF1("Gauss","gausn",0.8,1.3);
85 h->Fit("Gauss","R");
86 h->Fit("Gauss","RI");
87 h->Fit("Gauss","RI");
88 Double_t* params = Gauss->GetParameters();
89 rms=params[2];
90 mean=params[1];
91 h->Draw("P");
92 h->GetXaxis()->SetTitle("E_{T}^{rec}/E_{T}^{MC} [GeV]");
93}
94
95
96TH1F * 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)
97{
98 string temp = histo;
99 string nom = histo.erase(0,histo.find(">>")+2);
100 TH1F *h = new TH1F(nom.c_str(),"pt_lept ",numBin,minX,maxX);
101 Analyze->Draw(temp.c_str());
102 h->SetLineColor(Lcolor);
103 h->SetFillColor(Fcolor);
104 h->SetLineStyle(Lstyle);
105 h->SetLineWidth(width);
106 if(Log)gPad->SetLogy();
107 if(h->Integral()!=0)h->Scale(1*Scale);
108 return h;
109}
110
111
112TH1F * MakeNormTH1FSum(Int_t numBin,Float_t minX,Float_t maxX,TTree * analyze1,TTree * analyze2,TTree * analyze3,
113string * histo,float * weight, Int_t Lcolor, Int_t Fcolor, Int_t Lstyle,Int_t width=2,bool Log=false)
114{
115
116 TTree *analyse[3]={analyze1,analyze2,analyze3};
117 TH1F *hq = new TH1F("hq","pt_lept ",numBin,minX,maxX);
118
119 for(unsigned int i=0; i<3; i++)
120 {
121 string temp = histo[i];
122 string nom = histo[i].erase(0,histo[i].find(">>")+2);
123 TH1F *h = new TH1F(nom.c_str(),"pt_lept ",numBin,minX,maxX);
124 analyse[i]->Draw(temp.c_str());
125 h->Scale(weight[i]);
126 hq->Add(h);
127 }
128
129 hq->SetLineColor(Lcolor);
130 hq->SetFillColor(Fcolor);
131 hq->SetLineStyle(Lstyle);
132 hq->SetLineWidth(width);
133 if(Log)gPad->SetLogy();
134 if(hq->Integral()!=0)hq->Scale(1./hq->Integral());
135
136 return hq;
137}
138
139TH1F * HiggsMakeNormTH1FSum(Int_t numBin,Float_t minX,Float_t maxX,TTree * analyze1,TTree * analyze2,TTree * analyze3,
140string * histo,float * weight, Int_t Lcolor, Int_t Fcolor, Int_t Lstyle,Int_t width=2,bool Log=false)
141{
142
143 TTree *analyse[3]={analyze1,analyze2,analyze3};
144 TH1F *hq = new TH1F("hq","pt_lept ",numBin,minX,maxX);
145
146 for(unsigned int i=0; i<3; i++)
147 {
148 string temp = histo[i];
149 string nom = histo[i].erase(0,histo[i].find(">>")+2);
150 TH1F *h = new TH1F(nom.c_str(),"pt_lept ",numBin,minX,maxX);
151 analyse[i]->Draw(temp.c_str());
152 h->Scale(weight[i]);
153 hq->Add(h);
154 }
155
156 hq->SetLineColor(Lcolor);
157 hq->SetFillColor(Fcolor);
158 hq->SetLineStyle(Lstyle);
159 hq->SetLineWidth(width);
160 if(Log)gPad->SetLogy();
161 if(hq->Integral()!=0)hq->Scale(1./hq->Integral());
162
163 return hq;
164}
165
166
167#endif
Note: See TracBrowser for help on using the repository browser.