ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/jbabb/ScriptsAndMacros/histoPlotterIsolation.C
Revision: 1.1
Committed: Sat Sep 3 21:22:35 2011 UTC (13 years, 8 months ago) by jbabb
Content type: text/plain
Branch: MAIN
CVS Tags: HEAD
Log Message:
*** empty log message ***

File Contents

# User Rev Content
1 jbabb 1.1 void patPlotter(){
2     /* Macro to plot histograms produced by PatAnalyzer
3     *
4     * Author: John Babb - UCR (stolen from Dominique Fortin)
5     */
6    
7     float nsigmas = 1.5; // Number of sigmas around mean to fit gaussian. It uses 2 iterations
8     // i.e. range is = [mu - nsigmas * sigma, mu + nsigmas * sigma]
9    
10     // Files for histogram output --> set suffixps to desired file type: e.g. .eps, .jpg, ...
11    
12     TString simType = "_full";
13     TString suffixps = ".gif";
14    
15     TString tfile = "";
16    
17     TFile *file = TFile::Open(tfile);
18    
19     TString plot2 = "IsolatedMuonPt"+simType+suffixps;
20     TString plot6 = "METPt"+simType+suffixps;
21     TString plot7 = "JetPt"+simType+suffixps;
22     TString plot8 = "Jet1Pt"+simType+suffixps;
23     TString plot9 = "Jet2Pt"+simType+suffixps;
24     TString plot10 = "Jet3Pt"+simType+suffixps;
25     TString plot11 = "Jet4Pt"+simType+suffixps;
26     TString plot37 = "IsolatedMuonTrackIsolation"+simType+suffixps;
27     TString plot40 = "IsolatedMuonCaloIsolation"+simType+suffixps;
28     TString plot45 = "MinDeltaR_IsolatedMuonJet"+simType+suffixps;
29     // TString plot82 = "h_Jet_emFraction"+simType+suffixps;
30     // TString plot83 = "h_Jet_hadFraction"+simType+suffixps;
31     // TString plot84 = "conePtDiff"+simType+suffixps;
32     // TString plot85 = "coneEtDiff"+simType+suffixps;
33     // TString plot86 = "sumPt03"+simType+suffixps;
34     // TString plot87 = "sumEt03"+simType+suffixps;
35     // TString plot88 = "sumPt05"+simType+suffixps;
36     // TString plot89 = "sumEt05"+simType+suffixps;
37     // TString plot90 = "sumEtVsPt03"+simType+suffixps;
38     // TString plot91 = "sumEtVsPt05"+simType+suffixps;
39     TString plot92 = "ptRelVsMinDr"+simType+suffixps;
40     TString plot93 = "ptRel"+simType+suffixps;
41     // TString plot94 = "genJetType_inclusive"+simType+suffixps;
42     // TString plot95 = "genJetType_2Reco"+simType+suffixps;
43     // TString plot96 = "genJetType_3Reco"+simType+suffixps;
44     // TString plot97 = "genJetType_4Reco"+simType+suffixps;
45     TString plot98 = "muonPP"+simType+suffixps;
46     TString plot99 = "nearestJetMass"+simType+suffixps;
47     TString plot100 = "mjMass"+simType+suffixps;
48     TString plot101 = "ptRelVsNearestJetPt"+simType+suffixps;
49     TString plot102 = "ptRelVsPt"+simType+suffixps;
50     TString plot103 = "ptRelVsMjMass"+simType+suffixps;
51     TString plot104 = "nearestJetPtVsMuPt"+simType+suffixps;
52     TString plot105 = "mjmetMass"+simType+suffixps;
53     TString plot106 = "MinDeltaR_IsolatedMuonJet1"+simType+suffixps;
54     TString plot107 = "MinDeltaR_IsolatedMuonJet2"+simType+suffixps;
55     // TString plot108 = "Jet1VsJet2"+simType+suffixps;
56     TString plot109 = "MinDeltaEta_IsolatedMuonJet"+simType+suffixps;
57     TString plot110 = "MinDeltaPhi_IsolatedMuonJet"+simType+suffixps;
58     TString plot111 = "nearestJetPtVsDr"+simType+suffixps;
59     TString plot112 = "nearestJetId"+simType+suffixps;
60     TString plot113 = "nearestJetEtDbMuPt"+simType+suffixps;
61     TString plot114 = "nearestJetEtDbMuPtVsPt"+simType+suffixps;
62     TString plot115 = "nearestJetEtDbMuPtVsEt"+simType+suffixps;
63     // TString plot116 = "mcMatchedPtRes"+simType+suffixps;
64     // TString plot117 = "mcMatchedDr"+simType+suffixps;
65     TString plot118 = "pTratio"+simType+suffixps;
66     TString plot119 = "IsolatedMuonRelativeIsolation"+simType+suffixps;
67     /*
68     TString plot120 = "UEvDR_noID"+simType+suffixps;
69     TString plot121 = "UEvDR_bTags"+simType+suffixps;
70     TString plot122 = "EvDR_noID"+simType+suffixps;
71     TString plot123 = "UEvE_noID"+simType+suffixps;
72     TString plot124 = "EvDR_bTags"+simType+suffixps;
73     TString plot125 = "UEvE_bTags"+simType+suffixps;
74     */
75     /*
76     TString plot126 = "MuonDr12"+simType+suffixps;
77     TString plot127 = "MuonDeta12"+simType+suffixps;
78     TString plot128 = "MuonDphi12"+simType+suffixps;
79     TString plot129 = "MuonRelPt1v2"+simType+suffixps;
80     TString plot130 = "MuonMinDr1v2"+simType+suffixps;
81    
82     TString plot131 = "MuonPt1"+simType+suffixps;
83     TString plot132 = "MuonPt2"+simType+suffixps;
84     TString plot133 = "MuonPt3"+simType+suffixps;
85     TString plot134 = "MuonEta1"+simType+suffixps;
86     TString plot135 = "MuonEta2"+simType+suffixps;
87     TString plot136 = "MuonEta3"+simType+suffixps;
88     TString plot137 = "MuonPhi1"+simType+suffixps;
89     TString plot138 = "MuonPhi2"+simType+suffixps;
90     TString plot139 = "MuonPhi3"+simType+suffixps;
91     TString plot140 = "MuonMinDr1"+simType+suffixps;
92     TString plot141 = "MuonMinDr2"+simType+suffixps;
93     TString plot142 = "MuonMinDr3"+simType+suffixps;
94     TString plot143 = "MuonRelPt1"+simType+suffixps;
95     TString plot144 = "MuonRelPt2"+simType+suffixps;
96     TString plot145 = "MuonRelPt3"+simType+suffixps;
97     TString plot146 = "MuonDeltaPhi1"+simType+suffixps;
98     TString plot147 = "MuonDeltaPhi2"+simType+suffixps;
99     TString plot148 = "MuonDeltaPhi3"+simType+suffixps;
100     TString plot149 = "MuonDeltaEta1"+simType+suffixps;
101     TString plot150 = "MuonDeltaEta2"+simType+suffixps;
102     TString plot151 = "MuonDeltaEta3"+simType+suffixps;
103     TString plot152 = "MuonPtVsMinDr1"+simType+suffixps;
104     TString plot153 = "MuonPtVsMinDr2"+simType+suffixps;
105     TString plot154 = "MuonPtVsMinDr3"+simType+suffixps;
106     */
107     TString plot155 = "ptRel_W"+simType+suffixps;
108     TString plot156 = "ptRel_B"+simType+suffixps;
109     TString plot157 = "ptRel_C"+simType+suffixps;
110     TString plot158 = "ptRel_DIF"+simType+suffixps;
111     TString plot159 = "MinDeltaR_IsolatedMuonJet_W"+simType+suffixps;
112     TString plot160 = "MinDeltaR_IsolatedMuonJet_B"+simType+suffixps;
113     TString plot161 = "MinDeltaR_IsolatedMuonJet_C"+simType+suffixps;
114     TString plot162 = "MinDeltaR_IsolatedMuonJet_DIF"+simType+suffixps;
115     TString plot163 = "ptRelVsMinDr_W"+simType+suffixps;
116     TString plot164 = "ptRelVsMinDr_B"+simType+suffixps;
117     TString plot165 = "ptRelVsMinDr_C"+simType+suffixps;
118     TString plot166 = "ptRelVsMinDr_DIF"+simType+suffixps;
119     TString plot167 = "IsolatedMuonChi2"+simType+suffixps;
120     TString plot168 = "IsolatedMuonHits"+simType+suffixps;
121     //TString plot169 = "MuonRelPtVsMinDr1"+simType+suffixps;
122     // TString plot170 = "MuonRelPtVsMinDr2"+simType+suffixps;
123     // TString plot171 = "MuonRelPtVsMinDr3"+simType+suffixps;
124     TString plot172 = "MuonD0"+simType+suffixps;
125     TString plot173 = "MuonD0Error"+simType+suffixps;
126     TString plot174 = "RelIsoEff"+simType+suffixps;
127     TString plot175 = "drMinEff"+simType+suffixps;
128     TString plot176 = "ptRelEff"+simType+suffixps;
129     TString plot177 = "TwoDCutEff_box"+simType+suffixps;
130     TString plot178 = "TwoDCutEff_lego"+simType+suffixps;
131    
132    
133    
134     // ********************************************************************
135     // Pointers to histograms
136     // ********************************************************************
137    
138    
139     r_isoMuon_Pt = (TH1F *) file->Get("r_isoMuon_Pt_");
140     r_isoMuon_chi2 = (TH1F *) file->Get("r_isoMuon_chi2");
141     r_isoMuon_nHits = (TH1F *) file->Get("r_isoMuon_nHits");
142     r_MET_Pt = (TH1F *) file->Get("r_MET_Pt_");
143     r_Jet_Pt = (TH1F *) file->Get("r_Jet_Pt_");
144     r_Jet1_Pt = (TH1F *) file->Get("r_Jet1_Pt_");
145     r_Jet2_Pt = (TH1F *) file->Get("r_Jet2_Pt_");
146     r_Jet3_Pt = (TH1F *) file->Get("r_Jet3_Pt_");
147     r_Jet4_Pt = (TH1F *) file->Get("r_Jet4_Pt_");
148     r_isoMuon_trkIso = (TH1F *) file->Get("r_isoMuon_trkIso_");
149     r_isoMuon_caloIso = (TH1F *) file->Get("r_isoMuon_caloIso_");
150     r_minDeltaR_mu_j_iso = (TH1F *) file->Get("r_minDeltaR_mu_j_iso_");
151     r_minDeltaR_mu_j_W = (TH1F *) file->Get("r_minDeltaR_mu_j_W");
152     r_minDeltaR_mu_j_b = (TH1F *) file->Get("r_minDeltaR_mu_j_b");
153     r_minDeltaR_mu_j_c = (TH1F *) file->Get("r_minDeltaR_mu_j_c");
154     r_minDeltaR_mu_j_dif = (TH1F *) file->Get("r_minDeltaR_mu_j_dif");
155     // r_Jet_emFraction = (TH1F *) file->Get("r_Jet_emFraction");
156     // r_Jet_hadFraction = (TH1F *) file->Get("r_Jet_hadFraction");
157     // r_conePtDiff = (TH1F *) file->Get("r_conePtDiff");
158     // r_coneEtDiff = (TH1F *) file->Get("r_coneEtDiff");
159     // r_sumPt03 = (TH1F *) file->Get("r_sumPt03");
160     // r_sumEt03 = (TH1F *) file->Get("r_sumEt03");
161     // r_sumPt05 = (TH1F *) file->Get("r_sumPt05");
162     // r_sumEt05 = (TH1F *) file->Get("r_sumEt05");
163     // r_sumEtVsPt03 = (TH2F *) file->Get("r_sumEtVsPt03");
164     // r_sumEtVsPt05 = (TH2F *) file->Get("r_sumEtVsPt05");
165     r_ptRelVsMinDr = (TH2F *) file->Get("r_ptRelVsMinDr");
166     r_ptRelVsMinDr_W = (TH2F *) file->Get("r_ptRelVsMinDr_W");
167     r_ptRelVsMinDr_c= (TH2F *) file->Get("r_ptRelVsMinDr_c");
168     r_ptRelVsMinDr_b = (TH2F *) file->Get("r_ptRelVsMinDr_b");
169     r_ptRelVsMinDr_dif = (TH2F *) file->Get("r_ptRelVsMinDr_dif");
170     r_ptRel = (TH1F *) file->Get("r_ptRel");
171     r_ptRel_W = (TH1F *) file->Get("r_ptRel_W");
172     r_ptRel_c = (TH1F *) file->Get("r_ptRel_c");
173     r_ptRel_b = (TH1F *) file->Get("r_ptRel_b");
174     r_ptRel_dif = (TH1F *) file->Get("r_ptRel_dif");
175     // genJetType_inclusive = (TH2F *) file->Get("genJetType_inclusive");
176     // genJetType_2Reco = (TH2F *) file->Get("genJetType_2Reco");
177     // genJetType_3Reco = (TH2F *) file->Get("genJetType_3Reco");
178     // genJetType_4Reco = (TH2F *) file->Get("genJetType_4Reco");
179     muonPP = (TH2F *) file->Get("muonPP");
180     r_mjMass = (TH1F *) file->Get("r_mjMass");
181     r_nearestJetMass = (TH1F *) file->Get("r_nearestJetMass");
182     r_ptRelVsNearestJetPt = (TH2F *) file->Get("r_ptRelVsNearestJetPt");
183     r_ptRelVsPt = (TH2F *) file->Get("r_ptRelVsPt");
184     r_ptRelVsMjMass = (TH2F *) file->Get("r_ptRelVsMjMass");
185     r_nearestJetPtVsMuPt = (TH2F *) file->Get("r_nearestJetPtVsMuPt");
186     r_mjmetMass = (TH1F *) file->Get("r_mjmetMass");
187     r_minDeltaR_mu_j1 = (TH1F *) file->Get("r_minDeltaR_mu_j1");
188     r_minDeltaR_mu_j2 = (TH1F *) file->Get("r_minDeltaR_mu_j2");
189     // r_Jet1VsJet2 = (TH1F *) file->Get("r_Jet1VsJet2");
190     r_nearestJetPtVsDr = (TH2F *) file->Get("r_nearestJetPtVsDr");
191     r_nearestJetId = (TH1F *) file->Get("r_nearestJetId");
192     eTnJbPt = (TH1F *) file->Get("eTnJbPt");
193     eTnJbPtVsEt = (TH1F *) file->Get("eTnJbPtVsEt");
194     eTnJbPtVsPt = (TH1F *) file->Get("eTnJbPtVsPt");
195     // mcPtRes = (TH1F *) file->Get("mcPtRes");
196     // mcDr = (TH1F *) file->Get("mcDr");
197     pTratio = (TH1F *) file->Get("pTratio");
198     r_isoMuon_relIso = (TH1F *) file->Get("r_isoMuon_relIso_");
199     /*
200     UEvDR_noID = (TH2F *) file->Get("UEvDR_noID");
201     EvDR_noID = (TH2F *) file->Get("EvDR_noID");
202     UEvE_noID = (TH2F *) file->Get("UEvE_noID");
203     UEvDR_bTags = (TH2F *) file->Get("UEvDR_bTags");
204     EvDR_bTags = (TH2F *) file->Get("EvDR_bTags");
205     UEvE_bTags = (TH2F *) file->Get("UEvE_bTags");
206     */
207     /*
208     muonDr1v2 = (TH1F *) file->Get("muonDr1v2");
209     muonDeta1v2 = (TH1F *) file->Get("muonDeta1v2");
210     muonDphi1v2 = (TH1F *) file->Get("muonDphi1v2");
211     muonRelPt1v2 = (TH2F *) file->Get("muonRelPt1v2");
212     muonDrMin1v2 = (TH2F *) file->Get("muonDrMin1v2");
213    
214     r_Muon1_Pt = (TH1F *) file->Get("r_Muon1_Pt");
215     r_Muon2_Pt = (TH1F *) file->Get("r_Muon2_Pt");
216     r_Muon3_Pt = (TH1F *) file->Get("r_Muon3_Pt");
217    
218     r_Muon1_Phi = (TH1F *) file->Get("r_Muon1_Phi");
219     r_Muon2_Phi = (TH1F *) file->Get("r_Muon2_Phi");
220     r_Muon3_Phi = (TH1F *) file->Get("r_Muon3_Phi");
221    
222     r_Muon1_Eta = (TH1F *) file->Get("r_Muon1_Eta");
223     r_Muon2_Eta = (TH1F *) file->Get("r_Muon2_Eta");
224     r_Muon3_Eta = (TH1F *) file->Get("r_Muon3_Eta");
225    
226     r_Muon1_dEta = (TH1F *) file->Get("r_Muon1_dEta");
227     r_Muon2_dEta = (TH1F *) file->Get("r_Muon2_dEta");
228     r_Muon3_dEta = (TH1F *) file->Get("r_Muon3_dEta");
229    
230     r_Muon1_dPhi = (TH1F *) file->Get("r_Muon1_dPhi");
231     r_Muon2_dPhi = (TH1F *) file->Get("r_Muon2_dPhi");
232     r_Muon3_dPhi = (TH1F *) file->Get("r_Muon3_dPhi");
233    
234     r_Muon1_RelPt = (TH1F *) file->Get("r_Muon1_RelPt");
235     r_Muon2_RelPt = (TH1F *) file->Get("r_Muon2_RelPt");
236     r_Muon3_RelPt = (TH1F *) file->Get("r_Muon3_RelPt");
237    
238     r_Muon1_minDr = (TH1F *) file->Get("r_Muon1_minDr");
239     r_Muon2_minDr = (TH1F *) file->Get("r_Muon2_minDr");
240     r_Muon3_minDr = (TH1F *) file->Get("r_Muon3_minDr");
241    
242     r_Muon1_ptVsminDr = (TH2F *) file->Get("r_Muon1_ptVsminDr");
243     r_Muon2_ptVsminDr = (TH2F *) file->Get("r_Muon2_ptVsminDr");
244     r_Muon3_ptVsminDr = (TH2F *) file->Get("r_Muon3_ptVsminDr");
245    
246     r_Muon1_relPtVsminDr = (TH2F *) file->Get("r_Muon1_relPtVsminDr");
247     r_Muon2_relPtVsminDr = (TH2F *) file->Get("r_Muon2_relPtVsminDr");
248     r_Muon3_relPtVsminDr = (TH2F *) file->Get("r_Muon3_relPtVsminDr");
249     */
250     r_isoMuon_d0 = (TH1F *) file->Get("r_isoMuon_d0");
251     r_isoMuon_d0Error = (TH1F *) file->Get("r_isoMuon_d0Error");
252    
253     relIsoEff_sig = (TH1F *) file->Get("relIsoEff");
254     drMinEff_sig = (TH1F *) file->Get("drMinEff");
255     ptRelEff_sig = (TH1F *) file->Get("ptRelEff");
256     ptRelVdrMinEff_sig = (TH2F *) file->Get("ptRelVdrMinEff");
257    
258    
259     // *****************************************************************
260     // Do the drawing
261     // *****************************************************************
262    
263    
264     //gStyle->SetOptStat("");
265     TCanvas *c1 = new TCanvas("c1","");
266     c1->SetFillColor(10);
267     relIsoEff_sig->SetLineColor(kRed);
268     relIsoEff_sig->SetTitle("RelIso Efficiency"); //change here
269     relIsoEff_sig->Draw();
270     relIsoEff_sig->GetYaxis()->SetTitle("Eff"); //and here
271     relIsoEff_sig->GetXaxis()->SetTitle("relIso"); //and here
272     c1->Print(plot174); //and finally here
273    
274     //gStyle->SetOptStat("");
275     TCanvas *c1 = new TCanvas("c1","");
276     c1->SetFillColor(10);
277     drMinEff_sig->SetLineColor(kRed);
278     drMinEff_sig->SetTitle("drMin Efficiency"); //change here
279     drMinEff_sig->Draw();
280     drMinEff_sig->GetYaxis()->SetTitle("Eff"); //and here
281     drMinEff_sig->GetXaxis()->SetTitle("drMin"); //and here
282     c1->Print(plot175); //and finally here
283    
284     //gStyle->SetOptStat("");
285     TCanvas *c1 = new TCanvas("c1","");
286     c1->SetFillColor(10);
287     ptRelEff_sig->SetLineColor(kRed);
288     ptRelEff_sig->SetTitle("ptRel Efficiency"); //change here
289     ptRelEff_sig->Draw();
290     ptRelEff_sig->GetYaxis()->SetTitle("Eff"); //and here
291     ptRelEff_sig->GetXaxis()->SetTitle("ptRel"); //and here
292     c1->Print(plot176); //and finally here
293    
294     //gStyle->SetOptStat("");
295     TCanvas *c1 = new TCanvas("c1","");
296     c1->SetFillColor(10);
297     ptRelVdrMinEff_sig->SetFillColor(kBlue);
298     ptRelVdrMinEff_sig->SetTitle("2D Cut Efficiency"); //change here
299     ptRelVdrMinEff_sig->Draw("BOX");
300     ptRelVdrMinEff_sig->GetYaxis()->SetTitle("drMin"); //and here
301     ptRelVdrMinEff_sig->GetXaxis()->SetTitle("ptRel"); //and here
302     c1->Print(plot177); //and finally here
303    
304     //gStyle->SetOptStat("");
305     TCanvas *c1 = new TCanvas("c1","");
306     c1->SetFillColor(10);
307     ptRelVdrMinEff_sig->SetFillColor(kBlue);
308     ptRelVdrMinEff_sig->SetTitle("2D Cut Efficiency"); //change here
309     ptRelVdrMinEff_sig->Draw("LEGO");
310     ptRelVdrMinEff_sig->GetYaxis()->SetTitle("drMin"); //and here
311     ptRelVdrMinEff_sig->GetXaxis()->SetTitle("ptRel"); //and here
312     c1->Print(plot178); //and finally here
313    
314     gStyle->SetOptStat("nemruoi");
315     TCanvas *c1 = new TCanvas("c1","");
316     c1->SetFillColor(10);
317     c1->SetFillColor(10);
318     r_isoMuon_d0->SetTitle("Muon d0"); //change here
319     r_isoMuon_d0->Draw();
320     r_isoMuon_d0->GetYaxis()->SetTitle(""); //and here
321     r_isoMuon_d0->GetXaxis()->SetTitle("d0"); //and here
322     c1->Print(plot172); //and finally here
323    
324     gStyle->SetOptStat("nemruoi");
325     TCanvas *c1 = new TCanvas("c1","");
326     c1->SetFillColor(10);
327     c1->SetFillColor(10);
328     r_isoMuon_d0Error->SetTitle("Muon d0 Error"); //change here
329     r_isoMuon_d0Error->Draw();
330     r_isoMuon_d0Error->GetYaxis()->SetTitle(""); //and here
331     r_isoMuon_d0Error->GetXaxis()->SetTitle("d0 Error"); //and here
332     c1->Print(plot173); //and finally here
333    
334     gStyle->SetOptStat("nemruoi");
335     TCanvas *c1 = new TCanvas("c1","");
336     c1->SetFillColor(10);
337     c1->SetFillColor(10);
338     muonPP->SetTitle("Muon Parent pdgId"); //change here
339     muonPP->Draw();
340     muonPP->GetYaxis()->SetTitle(""); //and here
341     muonPP->GetXaxis()->SetTitle("Particle Id"); //and here
342     c1->Print(plot98); //and finally here
343    
344    
345     gStyle->SetOptStat("nemruoi");
346     TCanvas *c1 = new TCanvas("c1","");
347     c1->SetFillColor(10);
348     c1->SetFillColor(10);
349     r_isoMuon_Pt->SetTitle("Isolated Muon Pt"); //change here
350     r_isoMuon_Pt->Draw();
351     r_isoMuon_Pt->GetYaxis()->SetTitle(""); //and here
352     r_isoMuon_Pt->GetXaxis()->SetTitle("Pt(GeV)"); //and here
353     c1->Print(plot2); //and finally here
354    
355     gStyle->SetOptStat("nemruoi");
356     TCanvas *c1 = new TCanvas("c1","");
357     c1->SetFillColor(10);
358     c1->SetFillColor(10);
359     r_MET_Pt->SetTitle("MET Pt"); //change here
360     r_MET_Pt->Draw();
361     r_MET_Pt->GetYaxis()->SetTitle(""); //and here
362     r_MET_Pt->GetXaxis()->SetTitle("Pt(GeV)"); //and here
363     c1->Print(plot6); //and finally here
364    
365     gStyle->SetOptStat("nemruoi");
366     TCanvas *c1 = new TCanvas("c1","");
367     c1->SetFillColor(10);
368     c1->SetFillColor(10);
369     r_Jet_Pt->SetTitle("Jet Pt"); //change here
370     r_Jet_Pt->Draw();
371     r_Jet_Pt->GetYaxis()->SetTitle(""); //and here
372     r_Jet_Pt->GetXaxis()->SetTitle("Pt(GeV)"); //and here
373     c1->Print(plot7); //and finally here
374    
375     gStyle->SetOptStat("nemruoi");
376     TCanvas *c1 = new TCanvas("c1","");
377     c1->SetFillColor(10);
378     c1->SetFillColor(10);
379     r_Jet1_Pt->SetTitle("Jet 1 Pt"); //change here
380     r_Jet1_Pt->Draw();
381     r_Jet1_Pt->GetYaxis()->SetTitle(""); //and here
382     r_Jet1_Pt->GetXaxis()->SetTitle("Pt(GeV)"); //and here
383     c1->Print(plot8); //and finally here
384    
385     gStyle->SetOptStat("nemruoi");
386     TCanvas *c1 = new TCanvas("c1","");
387     c1->SetFillColor(10);
388     c1->SetFillColor(10);
389     r_Jet2_Pt->SetTitle("Jet 2 Pt"); //change here
390     r_Jet2_Pt->Draw();
391     r_Jet2_Pt->GetYaxis()->SetTitle(""); //and here
392     r_Jet2_Pt->GetXaxis()->SetTitle("Pt(GeV)"); //and here
393     c1->Print(plot9); //and finally here
394    
395     gStyle->SetOptStat("nemruoi");
396     TCanvas *c1 = new TCanvas("c1","");
397     c1->SetFillColor(10);
398     c1->SetFillColor(10);
399     r_Jet3_Pt->SetTitle("Jet 3 Pt"); //change here
400     r_Jet3_Pt->Draw();
401     r_Jet3_Pt->GetYaxis()->SetTitle(""); //and here
402     r_Jet3_Pt->GetXaxis()->SetTitle("Pt(GeV)"); //and here
403     c1->Print(plot10); //and finally here
404    
405     gStyle->SetOptStat("nemruoi");
406     TCanvas *c1 = new TCanvas("c1","");
407     c1->SetFillColor(10);
408     c1->SetFillColor(10);
409     r_Jet4_Pt->SetTitle("Jet 4 Pt"); //change here
410     r_Jet4_Pt->Draw();
411     r_Jet4_Pt->GetYaxis()->SetTitle(""); //and here
412     r_Jet4_Pt->GetXaxis()->SetTitle("Pt(GeV)"); //and here
413     c1->Print(plot11); //and finally here
414    
415     gStyle->SetOptStat("nemruoi");
416     TCanvas *c1 = new TCanvas("c1","");
417     c1->SetFillColor(10);
418     c1->SetFillColor(10);
419     r_isoMuon_trkIso->SetTitle("Isolated Muon Track Isolation"); //change here
420     r_isoMuon_trkIso->Draw();
421     r_isoMuon_trkIso->GetYaxis()->SetTitle(""); //and here
422     r_isoMuon_trkIso->GetXaxis()->SetTitle("Isolation"); //and here
423     c1->Print(plot37); //and finally here
424    
425     gStyle->SetOptStat("nemruoi");
426     TCanvas *c1 = new TCanvas("c1","");
427     c1->SetFillColor(10);
428     c1->SetFillColor(10);
429     r_isoMuon_caloIso->SetTitle("Isolated Muon Calorimeter Isolation"); //change here
430     r_isoMuon_caloIso->Draw();
431     r_isoMuon_caloIso->GetYaxis()->SetTitle(""); //and here
432     r_isoMuon_caloIso->GetXaxis()->SetTitle("Isolation"); //and here
433     c1->Print(plot40); //and finally here
434    
435     gStyle->SetOptStat("nemruoi");
436     TCanvas *c1 = new TCanvas("c1","");
437     c1->SetFillColor(10);
438     c1->SetFillColor(10);
439     r_minDeltaR_mu_j_iso->SetTitle("Min Distance Between Isolated Muon and Nearest Jet"); //change here
440     r_minDeltaR_mu_j_iso->Draw();
441     r_minDeltaR_mu_j_iso->GetYaxis()->SetTitle(""); //and here
442     r_minDeltaR_mu_j_iso->GetXaxis()->SetTitle("dR"); //and here
443     c1->Print(plot45); //and finally here
444     /*
445     gStyle->SetOptStat("nemruoi");
446     TCanvas *c1 = new TCanvas("c1","");
447     c1->SetFillColor(10);
448     c1->SetFillColor(10);
449     r_Jet_emFraction->SetTitle("EM Energy Fraction"); //change here
450     r_Jet_emFraction->Draw();
451     r_Jet_emFraction->GetYaxis()->SetTitle(""); //and here
452     r_Jet_emFraction->GetXaxis()->SetTitle("em fraction"); //and here
453     c1->Print(plot82); //and finally here
454    
455     gStyle->SetOptStat("nemruoi");
456     TCanvas *c1 = new TCanvas("c1","");
457     c1->SetFillColor(10);
458     c1->SetFillColor(10);
459     r_Jet_hadFraction->SetTitle("Hadronic Energy Fraction"); //change here
460     r_Jet_hadFraction->Draw();
461     r_Jet_hadFraction->GetYaxis()->SetTitle(""); //and here
462     r_Jet_hadFraction->GetXaxis()->SetTitle("had fraction"); //and here
463     c1->Print(plot83); //and finally here
464    
465     gStyle->SetOptStat("nemruoi");
466     TCanvas *c1 = new TCanvas("c1","");
467     c1->SetFillColor(10);
468     c1->SetFillColor(10);
469     r_conePtDiff->SetTitle("Cone Pt Diff: 0.5 - 0.3"); //change here
470     r_conePtDiff->Draw();
471     r_conePtDiff->GetYaxis()->SetTitle(""); //and here
472     r_conePtDiff->GetXaxis()->SetTitle("GeV"); //and here
473     c1->Print(plot84); //and finally here
474    
475     gStyle->SetOptStat("nemruoi");
476     TCanvas *c1 = new TCanvas("c1","");
477     c1->SetFillColor(10);
478     c1->SetFillColor(10);
479     r_coneEtDiff->SetTitle("Cone Et Diff: 0.5 - 0.3"); //change here
480     r_coneEtDiff->Draw();
481     r_coneEtDiff->GetYaxis()->SetTitle(""); //and here
482     r_coneEtDiff->GetXaxis()->SetTitle("GeV"); //and here
483     c1->Print(plot85); //and finally here
484    
485     gStyle->SetOptStat("nemruoi");
486     TCanvas *c1 = new TCanvas("c1","");
487     c1->SetFillColor(10);
488     c1->SetFillColor(10);
489     r_sumPt03->SetTitle("Cone Pt: 0.3"); //change here
490     r_sumPt03->Draw();
491     r_sumPt03->GetYaxis()->SetTitle(""); //and here
492     r_sumPt03->GetXaxis()->SetTitle("GeV"); //and here
493     c1->Print(plot86); //and finally here
494    
495     gStyle->SetOptStat("nemruoi");
496     TCanvas *c1 = new TCanvas("c1","");
497     c1->SetFillColor(10);
498     c1->SetFillColor(10);
499     r_sumEt03->SetTitle("Cone Et: 0.3"); //change here
500     r_sumEt03->Draw();
501     r_sumEt03->GetYaxis()->SetTitle(""); //and here
502     r_sumEt03->GetXaxis()->SetTitle("GeV"); //and here
503     c1->Print(plot87); //and finally here
504    
505     gStyle->SetOptStat("nemruoi");
506     TCanvas *c1 = new TCanvas("c1","");
507     c1->SetFillColor(10);
508     c1->SetFillColor(10);
509     r_sumPt05->SetTitle("Cone Pt: 0.5"); //change here
510     r_sumPt05->Draw();
511     r_sumPt05->GetYaxis()->SetTitle(""); //and here
512     r_sumPt05->GetXaxis()->SetTitle("GeV"); //and here
513     c1->Print(plot88); //and finally here
514    
515     gStyle->SetOptStat("nemruoi");
516     TCanvas *c1 = new TCanvas("c1","");
517     c1->SetFillColor(10);
518     c1->SetFillColor(10);
519     r_sumEt05->SetTitle("Cone Et: 0.5"); //change here
520     r_sumEt05->Draw();
521     r_sumEt05->GetYaxis()->SetTitle(""); //and here
522     r_sumEt05->GetXaxis()->SetTitle("GeV"); //and here
523     c1->Print(plot89); //and finally here
524    
525     gStyle->SetOptStat("nemruoi");
526     TCanvas *c1 = new TCanvas("c1","");
527     c1->SetFillColor(10);
528     c1->SetFillColor(10);
529     r_sumEtVsPt03->SetTitle("Cone Et: 0.3 Vs Cone Pt: 0.3"); //change here
530     r_sumEtVsPt03->Draw();
531     r_sumEtVsPt03->GetYaxis()->SetTitle("GeV"); //and here
532     r_sumEtVsPt03->GetXaxis()->SetTitle("GeV"); //and here
533     c1->Print(plot90); //and finally here
534    
535     gStyle->SetOptStat("nemruoi");
536     TCanvas *c1 = new TCanvas("c1","");
537     c1->SetFillColor(10);
538     c1->SetFillColor(10);
539     r_sumEtVsPt05->SetTitle("Cone Et: 0.5 Vs Cone Pt: 0.5"); //change here
540     r_sumEtVsPt05->Draw();
541     r_sumEtVsPt05->GetYaxis()->SetTitle("GeV"); //and here
542     r_sumEtVsPt05->GetXaxis()->SetTitle("GeV"); //and here
543     c1->Print(plot91); //and finally here
544     */
545     gStyle->SetOptStat("nemruoi");
546     TCanvas *c1 = new TCanvas("c1","");
547     c1->SetFillColor(10);
548     c1->SetFillColor(10);
549     r_ptRelVsMinDr->SetTitle("RelPt Vs min dR"); //change here
550     r_ptRelVsMinDr->Draw();
551     r_ptRelVsMinDr->GetYaxis()->SetTitle("GeV"); //and here
552     r_ptRelVsMinDr->GetXaxis()->SetTitle("dR"); //and here
553     c1->Print(plot92); //and finally here
554    
555     gStyle->SetOptStat("nemruoi");
556     TCanvas *c1 = new TCanvas("c1","");
557     c1->SetFillColor(10);
558     c1->SetFillColor(10);
559     r_ptRel->SetTitle("Pt relative to the nearest jet"); //change here
560     r_ptRel->Draw();
561     r_ptRel->GetYaxis()->SetTitle(""); //and here
562     r_ptRel->GetXaxis()->SetTitle("GeV"); //and here
563     c1->Print(plot93); //and finally here
564     /*
565     gStyle->SetOptStat(kTRUE);
566     TCanvas *c1 = new TCanvas("c1","");
567     c1->SetFillColor(10);
568     c1->SetFillColor(10);
569     genJetType_inclusive->SetTitle("GenJet type and multiplicity per event (Inclusive)"); //change here
570     genJetType_inclusive->Draw("COLZ");
571     genJetType_inclusive->GetYaxis()->SetTitle("Jets per event"); //and here
572     genJetType_inclusive->GetXaxis()->SetTitle("Jet Type"); //and here
573     genJetType_inclusive->GetXaxis()->SetBinLabel(1, "d");
574     genJetType_inclusive->GetXaxis()->SetBinLabel(2, "u");
575     genJetType_inclusive->GetXaxis()->SetBinLabel(3, "s");
576     genJetType_inclusive->GetXaxis()->SetBinLabel(4, "c");
577     genJetType_inclusive->GetXaxis()->SetBinLabel(5, "b");
578     genJetType_inclusive->GetXaxis()->SetBinLabel(6, "g");
579     c1->Print(plot94); //and finally here
580    
581     gStyle->SetOptStat(kTRUE);
582     TCanvas *c1 = new TCanvas("c1","");
583     c1->SetFillColor(10);
584     c1->SetFillColor(10);
585     genJetType_2Reco->SetTitle("GenJet type and multiplicity per event (2Reco)"); //change here
586     genJetType_2Reco->Draw("COLZ");
587     genJetType_2Reco->GetYaxis()->SetTitle("Jets per event"); //and here
588     genJetType_2Reco->GetXaxis()->SetTitle("Jet Type"); //and here
589     genJetType_2Reco->GetXaxis()->SetBinLabel(1, "d");
590     genJetType_2Reco->GetXaxis()->SetBinLabel(2, "u");
591     genJetType_2Reco->GetXaxis()->SetBinLabel(3, "s");
592     genJetType_2Reco->GetXaxis()->SetBinLabel(4, "c");
593     genJetType_2Reco->GetXaxis()->SetBinLabel(5, "b");
594     genJetType_2Reco->GetXaxis()->SetBinLabel(6, "g");
595     c1->Print(plot95); //and finally here
596    
597     gStyle->SetOptStat(kTRUE);
598     TCanvas *c1 = new TCanvas("c1","");
599     c1->SetFillColor(10);
600     c1->SetFillColor(10);
601     genJetType_3Reco->SetTitle("GenJet type and multiplicity per event (3Reco)"); //change here
602     genJetType_3Reco->Draw("COLZ");
603     genJetType_3Reco->GetYaxis()->SetTitle("Jets per event"); //and here
604     genJetType_3Reco->GetXaxis()->SetTitle("Jet Type"); //and here
605     genJetType_3Reco->GetXaxis()->SetBinLabel(1, "d");
606     genJetType_3Reco->GetXaxis()->SetBinLabel(2, "u");
607     genJetType_3Reco->GetXaxis()->SetBinLabel(3, "s");
608     genJetType_3Reco->GetXaxis()->SetBinLabel(4, "c");
609     genJetType_3Reco->GetXaxis()->SetBinLabel(5, "b");
610     genJetType_3Reco->GetXaxis()->SetBinLabel(6, "g");
611     c1->Print(plot96); //and finally here
612    
613     gStyle->SetOptStat(kTRUE);
614     TCanvas *c1 = new TCanvas("c1","");
615     c1->SetFillColor(10);
616     c1->SetFillColor(10);
617     genJetType_4Reco->SetTitle("GenJet type and multiplicity per event (4Reco)"); //change here
618     genJetType_4Reco->Draw("COLZ");
619     genJetType_4Reco->GetYaxis()->SetTitle("Jets per event"); //and here
620     genJetType_4Reco->GetXaxis()->SetTitle("Jet Type"); //and here
621     genJetType_4Reco->GetXaxis()->SetBinLabel(1, "d");
622     genJetType_4Reco->GetXaxis()->SetBinLabel(2, "u");
623     genJetType_4Reco->GetXaxis()->SetBinLabel(3, "s");
624     genJetType_4Reco->GetXaxis()->SetBinLabel(4, "c");
625     genJetType_4Reco->GetXaxis()->SetBinLabel(5, "b");
626     genJetType_4Reco->GetXaxis()->SetBinLabel(6, "g");
627     c1->Print(plot97); //and finally here
628     */
629     gStyle->SetOptStat("nemruoi");
630     TCanvas *c1 = new TCanvas("c1","");
631     c1->SetFillColor(10);
632     c1->SetFillColor(10);
633     r_nearestJetPt->SetTitle("Pt of the jet nearest the muon"); //change here
634     r_nearestJetPt->Draw();
635     r_nearestJetPt->GetYaxis()->SetTitle(""); //and here
636     r_nearestJetPt->GetXaxis()->SetTitle("GeV"); //and here
637     c1->Print(plot99); //and finally here
638    
639     gStyle->SetOptStat("nemruoi");
640     TCanvas *c1 = new TCanvas("c1","");
641     c1->SetFillColor(10);
642     c1->SetFillColor(10);
643     r_mjMass->SetTitle("Combined mass of muon + nearest jet"); //change here
644     r_mjMass->Draw();
645     r_mjMass->GetYaxis()->SetTitle(""); //and here
646     r_mjMass->GetXaxis()->SetTitle("GeV"); //and here
647     c1->Print(plot100); //and finally here
648    
649     gStyle->SetOptStat("nemruoi");
650     TCanvas *c1 = new TCanvas("c1","");
651     c1->SetFillColor(10);
652     c1->SetFillColor(10);
653     r_ptRelVsNearestJetPt->SetTitle("RelPt Vs Nearest Jet Pt"); //change here
654     r_ptRelVsNearestJetPt->Draw();
655     r_ptRelVsNearestJetPt->GetYaxis()->SetTitle("GeV"); //and here
656     r_ptRelVsNearestJetPt->GetXaxis()->SetTitle("GeV"); //and here
657     c1->Print(plot101); //and finally here
658    
659     gStyle->SetOptStat("nemruoi");
660     TCanvas *c1 = new TCanvas("c1","");
661     c1->SetFillColor(10);
662     c1->SetFillColor(10);
663     r_ptRelVsPt->SetTitle("RelPt Vs Muon Pt"); //change here
664     r_ptRelVsPt->Draw();
665     r_ptRelVsPt->GetYaxis()->SetTitle("GeV"); //and here
666     r_ptRelVsPt->GetXaxis()->SetTitle("GeV"); //and here
667     c1->Print(plot102); //and finally here
668    
669     gStyle->SetOptStat("nemruoi");
670     TCanvas *c1 = new TCanvas("c1","");
671     c1->SetFillColor(10);
672     c1->SetFillColor(10);
673     r_ptRelVsMjMass->SetTitle("RelPt Vs M + J mass"); //change here
674     r_ptRelVsMjMass->Draw();
675     r_ptRelVsMjMass->GetYaxis()->SetTitle("GeV"); //and here
676     r_ptRelVsMjMass->GetXaxis()->SetTitle("GeV"); //and here
677     c1->Print(plot103); //and finally here
678    
679     gStyle->SetOptStat("nemruoi");
680     TCanvas *c1 = new TCanvas("c1","");
681     c1->SetFillColor(10);
682     c1->SetFillColor(10);
683     r_nearestJetPtVsMuPt->SetTitle("nearestJetPtVsMuPt"); //change here
684     r_nearestJetPtVsMuPt->Draw();
685     r_nearestJetPtVsMuPt->GetYaxis()->SetTitle("GeV"); //and here
686     r_nearestJetPtVsMuPt->GetXaxis()->SetTitle("GeV"); //and here
687     c1->Print(plot104); //and finally here
688    
689     gStyle->SetOptStat("nemruoi");
690     TCanvas *c1 = new TCanvas("c1","");
691     c1->SetFillColor(10);
692     c1->SetFillColor(10);
693     r_mjmetMass->SetTitle("Combined mass of muon + nearest jet + MET"); //change here
694     r_mjmetMass->Draw();
695     r_mjmetMass->GetYaxis()->SetTitle(""); //and here
696     r_mjmetMass->GetXaxis()->SetTitle("GeV"); //and here
697     c1->Print(plot105); //and finally here
698    
699     gStyle->SetOptStat("nemruoi");
700     TCanvas *c1 = new TCanvas("c1","");
701     c1->SetFillColor(10);
702     c1->SetFillColor(10);
703     r_minDeltaR_mu_j1->SetTitle("dR between muon and leading jet"); //change here
704     r_minDeltaR_mu_j1->Draw();
705     r_minDeltaR_mu_j1->GetYaxis()->SetTitle(""); //and here
706     r_minDeltaR_mu_j1->GetXaxis()->SetTitle("dR"); //and here
707     c1->Print(plot106); //and finally here
708    
709     gStyle->SetOptStat("nemruoi");
710     TCanvas *c1 = new TCanvas("c1","");
711     c1->SetFillColor(10);
712     c1->SetFillColor(10);
713     r_minDeltaR_mu_j2->SetTitle("dR between muon and 2nd leading jet"); //change here
714     r_minDeltaR_mu_j2->Draw();
715     r_minDeltaR_mu_j2->GetYaxis()->SetTitle(""); //and here
716     r_minDeltaR_mu_j2->GetXaxis()->SetTitle("dR"); //and here
717     c1->Print(plot107); //and finally here
718     /*
719     gStyle->SetOptStat("nemruoi");
720     TCanvas *c1 = new TCanvas("c1","");
721     c1->SetFillColor(10);
722     c1->SetFillColor(10);
723     r_Jet1VsJet2->SetTitle("Jet1 pt vs Jet2 pt"); //change here
724     r_Jet1VsJet2->Draw();
725     r_Jet1VsJet2->GetYaxis()->SetTitle("GeV"); //and here
726     r_Jet1VsJet2->GetXaxis()->SetTitle("GeV"); //and here
727     c1->Print(plot108); //and finally here
728     */
729     gStyle->SetOptStat("nemruoi");
730     TCanvas *c1 = new TCanvas("c1","");
731     c1->SetFillColor(10);
732     c1->SetFillColor(10);
733     r_minDeltaEta_mj->SetTitle("dEta between muon and nearest jet"); //change here
734     r_minDeltaEta_mj->Draw();
735     r_minDeltaEta_mj->GetYaxis()->SetTitle(""); //and here
736     r_minDeltaEta_mj->GetXaxis()->SetTitle("dEta"); //and here
737     c1->Print(plot109); //and finally here
738    
739     gStyle->SetOptStat("nemruoi");
740     TCanvas *c1 = new TCanvas("c1","");
741     c1->SetFillColor(10);
742     c1->SetFillColor(10);
743     r_minDeltaPhi_mj->SetTitle("dPhi between muon and nearest jet"); //change here
744     r_minDeltaPhi_mj->Draw();
745     r_minDeltaPhi_mj->GetYaxis()->SetTitle(""); //and here
746     r_minDeltaPhi_mj->GetXaxis()->SetTitle("dPhi"); //and here
747     c1->Print(plot110); //and finally here
748    
749     gStyle->SetOptStat("nemruoi");
750     TCanvas *c1 = new TCanvas("c1","");
751     c1->SetFillColor(10);
752     c1->SetFillColor(10);
753     r_nearestJetPtVsDr->SetTitle("nearestJetPt Vs dR"); //change here
754     r_nearestJetPtVsDr->Draw();
755     r_nearestJetPtVsDr->GetYaxis()->SetTitle("GeV"); //and here
756     r_nearestJetPtVsDr->GetXaxis()->SetTitle("dR"); //and here
757     c1->Print(plot111); //and finally here
758    
759     gStyle->SetOptStat("nemruoi");
760     TCanvas *c1 = new TCanvas("c1","");
761     c1->SetFillColor(10);
762     c1->SetFillColor(10);
763     r_nearestJetId->SetTitle("Parton flavor of the jet nearest the muon"); //change here
764     r_nearestJetId->Draw();
765     r_nearestJetId->GetYaxis()->SetTitle(""); //and here
766     r_nearestJetId->GetXaxis()->SetTitle("pdgId"); //and here
767     c1->Print(plot112); //and finally here
768    
769     gStyle->SetOptStat("nemruoi");
770     TCanvas *c1 = new TCanvas("c1","");
771     c1->SetFillColor(10);
772     c1->SetFillColor(10);
773     eTnJbPt->SetTitle("Fragmentation Function"); //change here
774     eTnJbPt->Draw();
775     eTnJbPt->GetYaxis()->SetTitle(""); //and here
776     eTnJbPt->GetXaxis()->SetTitle(""); //and here
777     c1->Print(plot113); //and finally here
778    
779     gStyle->SetOptStat("nemruoi");
780     TCanvas *c1 = new TCanvas("c1","");
781     c1->SetFillColor(10);
782     c1->SetFillColor(10);
783     eTnJbPtVsPt->SetTitle("Fragmentation Function Vs muon Pt"); //change here
784     eTnJbPtVsPt->Draw();
785     eTnJbPtVsPt->GetYaxis()->SetTitle(""); //and here
786     eTnJbPtVsPt->GetXaxis()->SetTitle("GeV"); //and here
787     c1->Print(plot114); //and finally here
788    
789     gStyle->SetOptStat("nemruoi");
790     TCanvas *c1 = new TCanvas("c1","");
791     c1->SetFillColor(10);
792     c1->SetFillColor(10);
793     eTnJbPtVsEt->SetTitle("Fragmentation Function Vs nearest jet Et"); //change here
794     eTnJbPtVsEt->Draw();
795     eTnJbPtVsEt->GetYaxis()->SetTitle(""); //and here
796     eTnJbPtVsEt->GetXaxis()->SetTitle("GeV"); //and here
797     c1->Print(plot115); //and finally here
798     /*
799     gStyle->SetOptStat("nemruoi");
800     TCanvas *c1 = new TCanvas("c1","");
801     c1->SetFillColor(10);
802     c1->SetFillColor(10);
803     mcPtRes->SetTitle("MC matched Pt Resolution"); //change here
804     mcPtRes->Draw();
805     mcPtRes->GetYaxis()->SetTitle(""); //and here
806     mcPtRes->GetXaxis()->SetTitle("GeV"); //and here
807     c1->Print(plot116); //and finally here
808    
809     gStyle->SetOptStat("nemruoi");
810     TCanvas *c1 = new TCanvas("c1","");
811     c1->SetFillColor(10);
812     c1->SetFillColor(10);
813     mcDr->SetTitle("MC matched dR"); //change here
814     mcDr->Draw();
815     mcDr->GetYaxis()->SetTitle(""); //and here
816     mcDr->GetXaxis()->SetTitle("dR"); //and here
817     c1->Print(plot117); //and finally here
818     */
819    
820     gStyle->SetOptStat("nemruoi");
821     TCanvas *c1 = new TCanvas("c1","");
822     c1->SetFillColor(10);
823     c1->SetFillColor(10);
824     pTratio->SetTitle("Pt rel / P"); //change here
825     pTratio->Draw();
826     pTratio->GetYaxis()->SetTitle(""); //and here
827     pTratio->GetXaxis()->SetTitle(""); //and here
828     c1->Print(plot118); //and finally here
829    
830     gStyle->SetOptStat("nemruoi");
831     TCanvas *c1 = new TCanvas("c1","");
832     c1->SetFillColor(10);
833     c1->SetFillColor(10);
834     r_isoMuon_relIso->SetTitle("Selected Muon Relative Isolation"); //change here
835     r_isoMuon_relIso->Draw();
836     r_isoMuon_relIso->GetYaxis()->SetTitle(""); //and here
837     r_isoMuon_relIso->GetXaxis()->SetTitle("relIso"); //and here
838     c1->Print(plot119); //and finally here
839    
840     /*
841     gStyle->SetOptStat("nemruoi");
842     TCanvas *c1 = new TCanvas("c1","");
843     c1->SetFillColor(10);
844     c1->SetFillColor(10);
845     UEvDR_noID->SetTitle("Fake Jets Uncorrected Energy Vs dR "); //change here
846     UEvDR_noID->Draw();
847     UEvDR_noID->GetYaxis()->SetTitle("GeV"); //and here
848     UEvDR_noID->GetXaxis()->SetTitle("dR"); //and here
849     c1->Print(plot120); //and finally here
850    
851     gStyle->SetOptStat("nemruoi");
852     TCanvas *c1 = new TCanvas("c1","");
853     c1->SetFillColor(10);
854     c1->SetFillColor(10);
855     UEvDR_bTags->SetTitle("b-Jets Uncorrected Energy Vs dR "); //change here
856     UEvDR_bTags->Draw();
857     UEvDR_bTags->GetYaxis()->SetTitle("GeV"); //and here
858     UEvDR_bTags->GetXaxis()->SetTitle("dR"); //and here
859     c1->Print(plot121); //and finally here
860    
861     gStyle->SetOptStat("nemruoi");
862     TCanvas *c1 = new TCanvas("c1","");
863     c1->SetFillColor(10);
864     c1->SetFillColor(10);
865     EvDR_noID->SetTitle("Fake Jets Energy Vs dR "); //change here
866     EvDR_noID->Draw();
867     EvDR_noID->GetYaxis()->SetTitle("GeV"); //and here
868     EvDR_noID->GetXaxis()->SetTitle("dR"); //and here
869     c1->Print(plot122); //and finally here
870    
871     gStyle->SetOptStat("nemruoi");
872     TCanvas *c1 = new TCanvas("c1","");
873     c1->SetFillColor(10);
874     c1->SetFillColor(10);
875     UEvE_noID->SetTitle("Fake Jets Uncorrected Energy Vs Energy "); //change here
876     UEvE_noID->Draw();
877     UEvE_noID->GetYaxis()->SetTitle("GeV"); //and here
878     UEvE_noID->GetXaxis()->SetTitle("GeV"); //and here
879     c1->Print(plot123); //and finally here
880    
881     gStyle->SetOptStat("nemruoi");
882     TCanvas *c1 = new TCanvas("c1","");
883     c1->SetFillColor(10);
884     c1->SetFillColor(10);
885     EvDR_bTags->SetTitle("b-Jets Energy Vs dR "); //change here
886     EvDR_bTags->Draw();
887     EvDR_bTags->GetYaxis()->SetTitle("GeV"); //and here
888     EvDR_bTags->GetXaxis()->SetTitle("dR"); //and here
889     c1->Print(plot124); //and finally here
890    
891     gStyle->SetOptStat("nemruoi");
892     TCanvas *c1 = new TCanvas("c1","");
893     c1->SetFillColor(10);
894     c1->SetFillColor(10);
895     UEvE_bTags->SetTitle("b-Jets Uncorrected Energy Vs Energy "); //change here
896     UEvE_bTags->Draw();
897     UEvE_bTags->GetYaxis()->SetTitle("GeV"); //and here
898     UEvE_bTags->GetXaxis()->SetTitle("GeV"); //and here
899     c1->Print(plot125); //and finally here
900     */
901     /*
902     gStyle->SetOptStat("nemruoi");
903     TCanvas *c1 = new TCanvas("c1","");
904     c1->SetFillColor(10);
905     c1->SetFillColor(10);
906     muonDphi1v2->SetTitle("dPhi between two leading muons"); //change here
907     muonDphi1v2->Draw();
908     muonDphi1v2->GetYaxis()->SetTitle(""); //and here
909     muonDphi1v2->GetXaxis()->SetTitle("dPhi"); //and here
910     c1->Print(plot126); //and finally here
911    
912     gStyle->SetOptStat("nemruoi");
913     TCanvas *c1 = new TCanvas("c1","");
914     c1->SetFillColor(10);
915     c1->SetFillColor(10);
916     muonDr1v2->SetTitle("dR between two leading muons"); //change here
917     muonDr1v2->Draw();
918     muonDr1v2->GetYaxis()->SetTitle(""); //and here
919     muonDr1v2->GetXaxis()->SetTitle("dR"); //and here
920     c1->Print(plot127); //and finally here
921    
922     gStyle->SetOptStat("nemruoi");
923     TCanvas *c1 = new TCanvas("c1","");
924     c1->SetFillColor(10);
925     c1->SetFillColor(10);
926     muonDeta1v2->SetTitle("dEta between two leading muons"); //change here
927     muonDeta1v2->Draw();
928     muonDeta1v2->GetYaxis()->SetTitle(""); //and here
929     muonDeta1v2->GetXaxis()->SetTitle(""); //and here
930     c1->Print(plot128); //and finally here
931    
932     gStyle->SetOptStat("nemruoi");
933     TCanvas *c1 = new TCanvas("c1","");
934     c1->SetFillColor(10);
935     c1->SetFillColor(10);
936     muonRelPt1v2->SetTitle("RelPt1 Vs RelPt2"); //change here
937     muonRelPt1v2->Draw();
938     muonRelPt1v2->GetYaxis()->SetTitle("GeV"); //and here
939     muonRelPt1v2->GetXaxis()->SetTitle("GeV"); //and here
940     c1->Print(plot129); //and finally here
941    
942     gStyle->SetOptStat("nemruoi");
943     TCanvas *c1 = new TCanvas("c1","");
944     c1->SetFillColor(10);
945     c1->SetFillColor(10);
946     muonDrMin1v2->SetTitle("MinDr1 vs MinDr2"); //change here
947     muonDrMin1v2->Draw();
948     muonDrMin1v2->GetYaxis()->SetTitle("dR"); //and here
949     muonDrMin1v2->GetXaxis()->SetTitle("dR"); //and here
950     c1->Print(plot130); //and finally here
951    
952     gStyle->SetOptStat("nemruoi");
953     TCanvas *c1 = new TCanvas("c1","");
954     c1->SetFillColor(10);
955     c1->SetFillColor(10);
956     r_Muon1_Pt->SetTitle("Muon 1 Pt"); //change here
957     r_Muon1_Pt->Draw();
958     r_Muon1_Pt->GetYaxis()->SetTitle(""); //and here
959     r_Muon1_Pt->GetXaxis()->SetTitle("Pt"); //and here
960     c1->Print(plot131); //and finally here
961    
962     gStyle->SetOptStat("nemruoi");
963     TCanvas *c1 = new TCanvas("c1","");
964     c1->SetFillColor(10);
965     c1->SetFillColor(10);
966     r_Muon2_Pt->SetTitle("Muon 2 Pt"); //change here
967     r_Muon2_Pt->Draw();
968     r_Muon2_Pt->GetYaxis()->SetTitle(""); //and here
969     r_Muon2_Pt->GetXaxis()->SetTitle("Pt"); //and here
970     c1->Print(plot132); //and finally here
971    
972     gStyle->SetOptStat("nemruoi");
973     TCanvas *c1 = new TCanvas("c1","");
974     c1->SetFillColor(10);
975     c1->SetFillColor(10);
976     r_Muon3_Pt->SetTitle("Muon 3 Pt"); //change here
977     r_Muon3_Pt->Draw();
978     r_Muon3_Pt->GetYaxis()->SetTitle(""); //and here
979     r_Muon3_Pt->GetXaxis()->SetTitle("Pt"); //and here
980     c1->Print(plot133); //and finally here
981    
982     gStyle->SetOptStat("nemruoi");
983     TCanvas *c1 = new TCanvas("c1","");
984     c1->SetFillColor(10);
985     c1->SetFillColor(10);
986     r_Muon1_Phi->SetTitle("Muon 1 Phi"); //change here
987     r_Muon1_Phi->Draw();
988     r_Muon1_Phi->GetYaxis()->SetTitle(""); //and here
989     r_Muon1_Phi->GetXaxis()->SetTitle("Phi"); //and here
990     c1->Print(plot134); //and finally here
991    
992     gStyle->SetOptStat("nemruoi");
993     TCanvas *c1 = new TCanvas("c1","");
994     c1->SetFillColor(10);
995     c1->SetFillColor(10);
996     r_Muon2_Phi->SetTitle("Muon 2 Phi"); //change here
997     r_Muon2_Phi->Draw();
998     r_Muon2_Phi->GetYaxis()->SetTitle(""); //and here
999     r_Muon2_Phi->GetXaxis()->SetTitle("Phi"); //and here
1000     c1->Print(plot135); //and finally here
1001    
1002     gStyle->SetOptStat("nemruoi");
1003     TCanvas *c1 = new TCanvas("c1","");
1004     c1->SetFillColor(10);
1005     c1->SetFillColor(10);
1006     r_Muon3_Phi->SetTitle("Muon 3 Phi"); //change here
1007     r_Muon3_Phi->Draw();
1008     r_Muon3_Phi->GetYaxis()->SetTitle(""); //and here
1009     r_Muon3_Phi->GetXaxis()->SetTitle("Phi"); //and here
1010     c1->Print(plot136); //and finally here
1011    
1012     gStyle->SetOptStat("nemruoi");
1013     TCanvas *c1 = new TCanvas("c1","");
1014     c1->SetFillColor(10);
1015     c1->SetFillColor(10);
1016     r_Muon1_Eta->SetTitle("Muon 1 Eta"); //change here
1017     r_Muon1_Eta->Draw();
1018     r_Muon1_Eta->GetYaxis()->SetTitle(""); //and here
1019     r_Muon1_Eta->GetXaxis()->SetTitle("Eta"); //and here
1020     c1->Print(plot137); //and finally here
1021    
1022     gStyle->SetOptStat("nemruoi");
1023     TCanvas *c1 = new TCanvas("c1","");
1024     c1->SetFillColor(10);
1025     c1->SetFillColor(10);
1026     r_Muon2_Eta->SetTitle("Muon 2 Eta"); //change here
1027     r_Muon2_Eta->Draw();
1028     r_Muon2_Eta->GetYaxis()->SetTitle(""); //and here
1029     r_Muon1_Eta->GetXaxis()->SetTitle("Eta"); //and here
1030     c1->Print(plot138); //and finally here
1031    
1032     gStyle->SetOptStat("nemruoi");
1033     TCanvas *c1 = new TCanvas("c1","");
1034     c1->SetFillColor(10);
1035     c1->SetFillColor(10);
1036     r_Muon3_Eta->SetTitle("Muon 3 Eta"); //change here
1037     r_Muon3_Eta->Draw();
1038     r_Muon3_Eta->GetYaxis()->SetTitle(""); //and here
1039     r_Muon3_Eta->GetXaxis()->SetTitle("Eta"); //and here
1040     c1->Print(plot139); //and finally here
1041    
1042     gStyle->SetOptStat("nemruoi");
1043     TCanvas *c1 = new TCanvas("c1","");
1044     c1->SetFillColor(10);
1045     c1->SetFillColor(10);
1046     r_Muon1_minDr->SetTitle("Muon 1 dR min"); //change here
1047     r_Muon1_minDr->Draw();
1048     r_Muon1_minDr->GetYaxis()->SetTitle(""); //and here
1049     r_Muon1_minDr->GetXaxis()->SetTitle("dR"); //and here
1050     c1->Print(plot140); //and finally here
1051    
1052     gStyle->SetOptStat("nemruoi");
1053     TCanvas *c1 = new TCanvas("c1","");
1054     c1->SetFillColor(10);
1055     c1->SetFillColor(10);
1056     r_Muon2_minDr->SetTitle("Muon 2 dR min"); //change here
1057     r_Muon2_minDr->Draw();
1058     r_Muon2_minDr->GetYaxis()->SetTitle(""); //and here
1059     r_Muon2_minDr->GetXaxis()->SetTitle("dR"); //and here
1060     c1->Print(plot141); //and finally here
1061    
1062     gStyle->SetOptStat("nemruoi");
1063     TCanvas *c1 = new TCanvas("c1","");
1064     c1->SetFillColor(10);
1065     c1->SetFillColor(10);
1066     r_Muon3_minDr->SetTitle("Muon 3 dR min"); //change here
1067     r_Muon3_minDr->Draw();
1068     r_Muon3_minDr->GetYaxis()->SetTitle(""); //and here
1069     r_Muon3_minDr->GetXaxis()->SetTitle("dR"); //and here
1070     c1->Print(plot142); //and finally here
1071    
1072     gStyle->SetOptStat("nemruoi");
1073     TCanvas *c1 = new TCanvas("c1","");
1074     c1->SetFillColor(10);
1075     c1->SetFillColor(10);
1076     r_Muon1_RelPt->SetTitle("Muon 1 relPt"); //change here
1077     r_Muon1_RelPt->Draw();
1078     r_Muon1_RelPt->GetYaxis()->SetTitle(""); //and here
1079     r_Muon1_RelPt->GetXaxis()->SetTitle("GeV"); //and here
1080     c1->Print(plot143); //and finally here
1081    
1082     gStyle->SetOptStat("nemruoi");
1083     TCanvas *c1 = new TCanvas("c1","");
1084     c1->SetFillColor(10);
1085     c1->SetFillColor(10);
1086     r_Muon2_RelPt->SetTitle("Muon 2 relPt"); //change here
1087     r_Muon2_RelPt->Draw();
1088     r_Muon2_RelPt->GetYaxis()->SetTitle(""); //and here
1089     r_Muon2_RelPt->GetXaxis()->SetTitle("GeV"); //and here
1090     c1->Print(plot144); //and finally here
1091    
1092     gStyle->SetOptStat("nemruoi");
1093     TCanvas *c1 = new TCanvas("c1","");
1094     c1->SetFillColor(10);
1095     c1->SetFillColor(10);
1096     r_Muon3_RelPt->SetTitle("Muon 3 relPt"); //change here
1097     r_Muon3_RelPt->Draw();
1098     r_Muon3_RelPt->GetYaxis()->SetTitle(""); //and here
1099     r_Muon3_RelPt->GetXaxis()->SetTitle("GeV"); //and here
1100     c1->Print(plot145); //and finally here
1101    
1102     gStyle->SetOptStat("nemruoi");
1103     TCanvas *c1 = new TCanvas("c1","");
1104     c1->SetFillColor(10);
1105     c1->SetFillColor(10);
1106     r_Muon1_dPhi->SetTitle("Muon 1 min dPhi "); //change here
1107     r_Muon1_dPhi->Draw();
1108     r_Muon1_dPhi->GetYaxis()->SetTitle(""); //and here
1109     r_Muon1_dPhi->GetXaxis()->SetTitle("dPhi"); //and here
1110     c1->Print(plot146); //and finally here
1111    
1112     gStyle->SetOptStat("nemruoi");
1113     TCanvas *c1 = new TCanvas("c1","");
1114     c1->SetFillColor(10);
1115     c1->SetFillColor(10);
1116     r_Muon2_dPhi->SetTitle("Muon 2 min dPhi "); //change here
1117     r_Muon2_dPhi->Draw();
1118     r_Muon2_dPhi->GetYaxis()->SetTitle(""); //and here
1119     r_Muon2_dPhi->GetXaxis()->SetTitle("dPhi"); //and here
1120     c1->Print(plot147); //and finally here
1121    
1122     gStyle->SetOptStat("nemruoi");
1123     TCanvas *c1 = new TCanvas("c1","");
1124     c1->SetFillColor(10);
1125     c1->SetFillColor(10);
1126     r_Muon3_dPhi->SetTitle("Muon 3 min dPhi "); //change here
1127     r_Muon3_dPhi->Draw();
1128     r_Muon3_dPhi->GetYaxis()->SetTitle(""); //and here
1129     r_Muon3_dPhi->GetXaxis()->SetTitle("dPhi"); //and here
1130     c1->Print(plot148); //and finally here
1131    
1132     gStyle->SetOptStat("nemruoi");
1133     TCanvas *c1 = new TCanvas("c1","");
1134     c1->SetFillColor(10);
1135     c1->SetFillColor(10);
1136     r_Muon1_dEta->SetTitle("Muon 1 min dEta "); //change here
1137     r_Muon1_dEta->Draw();
1138     r_Muon1_dEta->GetYaxis()->SetTitle(""); //and here
1139     r_Muon1_dEta->GetXaxis()->SetTitle("dEta"); //and here
1140     c1->Print(plot149); //and finally here
1141    
1142     gStyle->SetOptStat("nemruoi");
1143     TCanvas *c1 = new TCanvas("c1","");
1144     c1->SetFillColor(10);
1145     c1->SetFillColor(10);
1146     r_Muon2_dEta->SetTitle("Muon 2 min dEta "); //change here
1147     r_Muon2_dEta->Draw();
1148     r_Muon2_dEta->GetYaxis()->SetTitle(""); //and here
1149     r_Muon2_dEta->GetXaxis()->SetTitle("dEta"); //and here
1150     c1->Print(plot150); //and finally here
1151    
1152     gStyle->SetOptStat("nemruoi");
1153     TCanvas *c1 = new TCanvas("c1","");
1154     c1->SetFillColor(10);
1155     c1->SetFillColor(10);
1156     r_Muon3_dEta->SetTitle("Muon 3 min dEta "); //change here
1157     r_Muon3_dEta->Draw();
1158     r_Muon3_dEta->GetYaxis()->SetTitle(""); //and here
1159     r_Muon3_dEta->GetXaxis()->SetTitle("dEta"); //and here
1160     c1->Print(plot151); //and finally here
1161    
1162     gStyle->SetOptStat("nemruoi");
1163     TCanvas *c1 = new TCanvas("c1","");
1164     c1->SetFillColor(10);
1165     c1->SetFillColor(10);
1166     r_Muon1_ptVsminDr->SetTitle("Muon 1 Pt Vs min dR "); //change here
1167     r_Muon1_ptVsminDr->Draw();
1168     r_Muon1_ptVsminDr->GetYaxis()->SetTitle("GeV"); //and here
1169     r_Muon1_ptVsminDr->GetXaxis()->SetTitle("dR"); //and here
1170     c1->Print(plot152); //and finally here
1171    
1172     gStyle->SetOptStat("nemruoi");
1173     TCanvas *c1 = new TCanvas("c1","");
1174     c1->SetFillColor(10);
1175     c1->SetFillColor(10);
1176     r_Muon2_ptVsminDr->SetTitle("Muon 2 Pt Vs min dR "); //change here
1177     r_Muon2_ptVsminDr->Draw();
1178     r_Muon2_ptVsminDr->GetYaxis()->SetTitle("GeV"); //and here
1179     r_Muon2_ptVsminDr->GetXaxis()->SetTitle("dR"); //and here
1180     c1->Print(plot153); //and finally here
1181    
1182     gStyle->SetOptStat("nemruoi");
1183     TCanvas *c1 = new TCanvas("c1","");
1184     c1->SetFillColor(10);
1185     c1->SetFillColor(10);
1186     r_Muon3_ptVsminDr->SetTitle("Muon 3 Pt Vs min dR "); //change here
1187     r_Muon3_ptVsminDr->Draw();
1188     r_Muon3_ptVsminDr->GetYaxis()->SetTitle("GeV"); //and here
1189     r_Muon3_ptVsminDr->GetXaxis()->SetTitle("dR"); //and here
1190     c1->Print(plot154); //and finally here
1191     */
1192     gStyle->SetOptStat("nemruoi");
1193     TCanvas *c1 = new TCanvas("c1","");
1194     c1->SetFillColor(10);
1195     c1->SetFillColor(10);
1196     r_ptRel_W->SetTitle("Pt relative to the nearest jet - W"); //change here
1197     r_ptRel_W->Draw();
1198     r_ptRel_W->GetYaxis()->SetTitle(""); //and here
1199     r_ptRel_W->GetXaxis()->SetTitle("GeV"); //and here
1200     c1->Print(plot155); //and finally here
1201    
1202     gStyle->SetOptStat("nemruoi");
1203     TCanvas *c1 = new TCanvas("c1","");
1204     c1->SetFillColor(10);
1205     c1->SetFillColor(10);
1206     r_ptRel_b->SetTitle("Pt relative to the nearest jet - b"); //change here
1207     r_ptRel_b->Draw();
1208     r_ptRel_b->GetYaxis()->SetTitle(""); //and here
1209     r_ptRel_b->GetXaxis()->SetTitle("GeV"); //and here
1210     c1->Print(plot156); //and finally here
1211    
1212     gStyle->SetOptStat("nemruoi");
1213     TCanvas *c1 = new TCanvas("c1","");
1214     c1->SetFillColor(10);
1215     c1->SetFillColor(10);
1216     r_ptRel_c->SetTitle("Pt relative to the nearest jet - c"); //change here
1217     r_ptRel_c->Draw();
1218     r_ptRel_c->GetYaxis()->SetTitle(""); //and here
1219     r_ptRel_c->GetXaxis()->SetTitle("GeV"); //and here
1220     c1->Print(plot157); //and finally here
1221    
1222     gStyle->SetOptStat("nemruoi");
1223     TCanvas *c1 = new TCanvas("c1","");
1224     c1->SetFillColor(10);
1225     c1->SetFillColor(10);
1226     r_ptRel_dif->SetTitle("Pt relative to the nearest jet - DIF"); //change here
1227     r_ptRel_dif->Draw();
1228     r_ptRel_dif->GetYaxis()->SetTitle(""); //and here
1229     r_ptRel_dif->GetXaxis()->SetTitle("GeV"); //and here
1230     c1->Print(plot158); //and finally here
1231    
1232     gStyle->SetOptStat("nemruoi");
1233     TCanvas *c1 = new TCanvas("c1","");
1234     c1->SetFillColor(10);
1235     c1->SetFillColor(10);
1236     r_minDeltaR_mu_j_W->SetTitle("Min Distance Between Isolated Muon and Nearest Jet - W"); //change here
1237     r_minDeltaR_mu_j_W->Draw();
1238     r_minDeltaR_mu_j_W->GetYaxis()->SetTitle(""); //and here
1239     r_minDeltaR_mu_j_W->GetXaxis()->SetTitle("dR"); //and here
1240     c1->Print(plot159); //and finally here
1241    
1242     gStyle->SetOptStat("nemruoi");
1243     TCanvas *c1 = new TCanvas("c1","");
1244     c1->SetFillColor(10);
1245     c1->SetFillColor(10);
1246     r_minDeltaR_mu_j_b->SetTitle("Min Distance Between Isolated Muon and Nearest Jet - b"); //change here
1247     r_minDeltaR_mu_j_b->Draw();
1248     r_minDeltaR_mu_j_b->GetYaxis()->SetTitle(""); //and here
1249     r_minDeltaR_mu_j_b->GetXaxis()->SetTitle("dR"); //and here
1250     c1->Print(plot160); //and finally here
1251    
1252     gStyle->SetOptStat("nemruoi");
1253     TCanvas *c1 = new TCanvas("c1","");
1254     c1->SetFillColor(10);
1255     c1->SetFillColor(10);
1256     r_minDeltaR_mu_j_c->SetTitle("Min Distance Between Isolated Muon and Nearest Jet - c"); //change here
1257     r_minDeltaR_mu_j_c->Draw();
1258     r_minDeltaR_mu_j_c->GetYaxis()->SetTitle(""); //and here
1259     r_minDeltaR_mu_j_c->GetXaxis()->SetTitle("dR"); //and here
1260     c1->Print(plot161); //and finally here
1261    
1262     gStyle->SetOptStat("nemruoi");
1263     TCanvas *c1 = new TCanvas("c1","");
1264     c1->SetFillColor(10);
1265     c1->SetFillColor(10);
1266     r_minDeltaR_mu_j_dif->SetTitle("Min Distance Between Isolated Muon and Nearest Jet - DIF"); //change here
1267     r_minDeltaR_mu_j_dif->Draw();
1268     r_minDeltaR_mu_j_dif->GetYaxis()->SetTitle(""); //and here
1269     r_minDeltaR_mu_j_dif->GetXaxis()->SetTitle("dR"); //and here
1270     c1->Print(plot162); //and finally here
1271    
1272     gStyle->SetOptStat("nemruoi");
1273     TCanvas *c1 = new TCanvas("c1","");
1274     c1->SetFillColor(10);
1275     c1->SetFillColor(10);
1276     r_ptRelVsMinDr_W->SetTitle("RelPt Vs min dR - W"); //change here
1277     r_ptRelVsMinDr_W->Draw();
1278     r_ptRelVsMinDr_W->GetYaxis()->SetTitle("GeV"); //and here
1279     r_ptRelVsMinDr_W->GetXaxis()->SetTitle("dR"); //and here
1280     c1->Print(plot163); //and finally here
1281    
1282     gStyle->SetOptStat("nemruoi");
1283     TCanvas *c1 = new TCanvas("c1","");
1284     c1->SetFillColor(10);
1285     c1->SetFillColor(10);
1286     r_ptRelVsMinDr_c->SetTitle("RelPt Vs min dR - c"); //change here
1287     r_ptRelVsMinDr_c->Draw();
1288     r_ptRelVsMinDr_c->GetYaxis()->SetTitle("GeV"); //and here
1289     r_ptRelVsMinDr_c->GetXaxis()->SetTitle("dR"); //and here
1290     c1->Print(plot164); //and finally here
1291    
1292     gStyle->SetOptStat("nemruoi");
1293     TCanvas *c1 = new TCanvas("c1","");
1294     c1->SetFillColor(10);
1295     c1->SetFillColor(10);
1296     r_ptRelVsMinDr_b->SetTitle("RelPt Vs min dR - b"); //change here
1297     r_ptRelVsMinDr_b->Draw();
1298     r_ptRelVsMinDr_b->GetYaxis()->SetTitle("GeV"); //and here
1299     r_ptRelVsMinDr_b->GetXaxis()->SetTitle("dR"); //and here
1300     c1->Print(plot165); //and finally here
1301    
1302     gStyle->SetOptStat("nemruoi");
1303     TCanvas *c1 = new TCanvas("c1","");
1304     c1->SetFillColor(10);
1305     c1->SetFillColor(10);
1306     r_ptRelVsMinDr_dif->SetTitle("RelPt Vs min dR - dif"); //change here
1307     r_ptRelVsMinDr_dif->Draw();
1308     r_ptRelVsMinDr_dif->GetYaxis()->SetTitle("GeV"); //and here
1309     r_ptRelVsMinDr_dif->GetXaxis()->SetTitle("dR"); //and here
1310     c1->Print(plot166); //and finally here
1311    
1312     gStyle->SetOptStat("nemruoi");
1313     TCanvas *c1 = new TCanvas("c1","");
1314     c1->SetFillColor(10);
1315     c1->SetFillColor(10);
1316     r_isoMuon_chi2->SetTitle("Selected Muon chi2"); //change here
1317     r_isoMuon_chi2->Draw();
1318     r_isoMuon_chi2->GetYaxis()->SetTitle(""); //and here
1319     r_isoMuon_chi2->GetXaxis()->SetTitle("Chi2"); //and here
1320     c1->Print(plot167); //and finally here
1321    
1322     gStyle->SetOptStat("nemruoi");
1323     TCanvas *c1 = new TCanvas("c1","");
1324     c1->SetFillColor(10);
1325     c1->SetFillColor(10);
1326     r_isoMuon_nHits->SetTitle("Selected Muon nHits"); //change here
1327     r_isoMuon_nHits->Draw();
1328     r_isoMuon_nHits->GetYaxis()->SetTitle(""); //and here
1329     r_isoMuon_nHits->GetXaxis()->SetTitle("nHits"); //and here
1330     c1->Print(plot168); //and finally here
1331     /*
1332     gStyle->SetOptStat("nemruoi");
1333     TCanvas *c1 = new TCanvas("c1","");
1334     c1->SetFillColor(10);
1335     c1->SetFillColor(10);
1336     r_Muon1_relPtVsminDr->SetTitle("Muon 1 relPt Vs min dR "); //change here
1337     r_Muon1_relPtVsminDr->Draw();
1338     r_Muon1_relPtVsminDr->GetYaxis()->SetTitle("GeV"); //and here
1339     r_Muon1_relPtVsminDr->GetXaxis()->SetTitle("dR"); //and here
1340     c1->Print(plot169); //and finally here
1341    
1342     gStyle->SetOptStat("nemruoi");
1343     TCanvas *c1 = new TCanvas("c1","");
1344     c1->SetFillColor(10);
1345     c1->SetFillColor(10);
1346     r_Muon2_relPtVsminDr->SetTitle("Muon 2 relPt Vs min dR "); //change here
1347     r_Muon2_relPtVsminDr->Draw();
1348     r_Muon2_relPtVsminDr->GetYaxis()->SetTitle("GeV"); //and here
1349     r_Muon2_relPtVsminDr->GetXaxis()->SetTitle("dR"); //and here
1350     c1->Print(plot170); //and finally here
1351    
1352     gStyle->SetOptStat("nemruoi");
1353     TCanvas *c1 = new TCanvas("c1","");
1354     c1->SetFillColor(10);
1355     c1->SetFillColor(10);
1356     r_Muon3_relPtVsminDr->SetTitle("Muon 3 relPt Vs min dR "); //change here
1357     r_Muon3_relPtVsminDr->Draw();
1358     r_Muon3_relPtVsminDr->GetYaxis()->SetTitle("GeV"); //and here
1359     r_Muon3_relPtVsminDr->GetXaxis()->SetTitle("dR"); //and here
1360     c1->Print(plot171); //and finally here
1361     */
1362    
1363    
1364    
1365     gROOT->ProcessLine(".q");
1366    
1367    
1368    
1369     }