ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/Development/Plotting/Modules/FSRStudy.C
(Generate patch)

Comparing UserCode/cbrown/Development/Plotting/Modules/FSRStudy.C (file contents):
Revision 1.1 by buchmann, Thu Oct 11 09:11:50 2012 UTC vs.
Revision 1.6 by buchmann, Tue Nov 6 17:00:04 2012 UTC

# Line 4 | Line 4
4   samplecollection PhotonSamples("PhotonSamples");
5  
6   namespace GenPhotons {
7 <  void CompareZmassToFinalStateLeptons();
8 <  void GenLevelStudies();
9 <  
10 <  
7 >    void CompareZmassToFinalStateLeptons();
8 >    void GenLevelStudies();
9 >    
10 >    void NumberOfPhotonsPerEvent();
11 >    void CharacterizingFSRPhotons();
12 >    void MakeCharacterizationPlot(string variable, int nbins, float min, float max, string label, string saveas);
13 >    void dR2DPlots();
14 >    void CompareGenJZBDistributions();
15 >    void EfficiencyPurity();
16 >    
17 >    TCut FSRcut;
18 >    TCut NoFSRcut;
19 >    TCut BaseCut;
20   }
21  
22   namespace RecoPhotons {
23 <  void RecoLevelStudies();
24 <  
25 <  
26 <  
23 >    void RecoLevelStudies();
24 >    
25 >    void NumberOfPhotonsPerEvent();
26 >    void CharacterizingFSRPhotons();
27 >    void MakeCharacterizationPlot(string variable, int nbins, float min, float max, string label, string saveas);
28 >    void dR2DPlots();
29 >    void CompareRecoJZBDistributions();
30 >    void EfficiencyPurity();
31 >    void FSRPhotonVsUnmatchedVsMatchedPhoton();
32 >    void FSRPhotonVsUnmatchedVsMatchedPhoton2D();
33 >    void JZBvsmllClosure();
34 >    
35 >    TCut FSRcut;
36 >    TCut NoFSRcut;
37 >    TCut BaseCut;
38 > }
39 >
40 > void RecoPhotons::JZBvsmllClosure() {
41 >    TCanvas *can = new TCanvas("can","can",1000,1000);
42 >    can->Divide(2,2);
43 >    can->cd(1);
44 >    
45 >    write_warning(__FUNCTION__,"JZB distribution at gen level seems to be screwed up");
46 >    
47 >    //we only want Z->ee and Z->mm for this (not Z->tautau)
48 >    
49 >    TCut em("!EventZToTaus");
50 >    
51 >    //no FSR recovery
52 >    TH1F *negmll_noFSR = PhotonSamples.Draw("negmll_noFSR", "genFSRmll", 36, 20, 200, "m_{ll}^{gen}", "events",BaseCut && TCut("genFSRjzb<0.1") && em,   mc,PlottingSetup::luminosity);
53 >    TH1F *posmll_noFSR = PhotonSamples.Draw("posmll_noFSR", "genFSRmll", 36, 20, 200, "m_{ll}^{gen}", "events",BaseCut && TCut("genFSRjzb>0.1") && em,   mc,PlottingSetup::luminosity);
54 >    
55 >    
56 >    //with FSR recovery
57 >    TH1F *negmll_FSR = PhotonSamples.Draw("negmll_FSR", "genFSRmllg", 36, 20, 200, "m_{ll}^{gen}", "events",BaseCut && TCut("genFSRjzbG<0.1") && em,   mc,PlottingSetup::luminosity);
58 >    TH1F *posmll_FSR = PhotonSamples.Draw("posmll_FSR", "genFSRmllg", 36, 20, 200, "m_{ll}^{gen}", "events",BaseCut && TCut("genFSRjzbG>0.1") && em,   mc,PlottingSetup::luminosity);
59 >    
60 >    cout << "Integrals: " << endl;
61 >    cout << "     no FSR:   JZB<0 : " << negmll_noFSR->Integral() << endl;
62 >    cout << "     no FSR:   JZB>0 : " << posmll_noFSR->Integral() << endl;
63 >    cout << "     with FSR: JZB<0 : " << negmll_FSR->Integral() << endl;
64 >    cout << "     with FSR: JZB>0 : " << posmll_FSR->Integral() << endl;
65 >    
66 >    negmll_FSR->Scale(1.0/negmll_FSR->Integral());
67 >    posmll_FSR->Scale(1.0/posmll_FSR->Integral());
68 >    
69 >    negmll_noFSR->Scale(1.0/negmll_noFSR->Integral());
70 >    posmll_noFSR->Scale(1.0/posmll_noFSR->Integral());
71 >    
72 >    
73 >    can->cd(1)->SetLogy(1);
74 >    can->cd(2)->SetLogy(1);
75 > //    can->cd(3)->SetLogy(1);
76 > //    can->cd(4)->SetLogy(1);
77 >    
78 >    negmll_noFSR->SetLineColor(kRed);
79 >    posmll_noFSR->SetLineColor(kBlack);
80 >    
81 >    negmll_FSR->SetLineColor(kRed);
82 >    posmll_FSR->SetLineColor(kBlack);
83 >    
84 >    TLegend *leg = make_legend();
85 >    leg->AddEntry(negmll_noFSR,"prediction","L");
86 >    leg->AddEntry(posmll_noFSR,"observation","P");
87 >    leg->SetHeader("#splitline{DY MC, gen level}{Only Z#rightarrow ee,mm}");
88 >    
89 >    can->cd(1);
90 >    negmll_noFSR->Draw("histo");
91 >    posmll_noFSR->Draw("e1,same");
92 >    TText *noFSR = write_title("No FSR recovery");
93 >    noFSR->Draw();
94 >    leg->Draw();
95 >    
96 >    can->cd(2);
97 >    negmll_FSR->Draw("histo");
98 >    posmll_FSR->Draw("e1,same");
99 >    TText *FSR = write_title("With FSR recovery");
100 >    FSR->Draw();
101 >    leg->Draw();
102 >    
103 >    can->cd(3);
104 >    TH1F *ratioNoFSR = (TH1F*)posmll_noFSR->Clone("ratioNoFSR");
105 >    ratioNoFSR->Divide(negmll_noFSR);
106 >    ratioNoFSR->GetYaxis()->SetRangeUser(0,2);
107 >    ratioNoFSR->GetYaxis()->SetTitle("obs/pred");
108 >    ratioNoFSR->Draw("e1");
109 >    TLine *oneline = new TLine(ratioNoFSR->GetBinLowEdge(1),1.0,ratioNoFSR->GetBinLowEdge(ratioNoFSR->GetNbinsX())+ratioNoFSR->GetBinWidth(ratioNoFSR->GetNbinsX()),1.0);
110 >    oneline->SetLineColor(kBlue);
111 >    oneline->Draw();
112 >    noFSR->Draw();
113 >    
114 >    can->cd(4);
115 >    TH1F *ratioFSR = (TH1F*)posmll_FSR->Clone("ratioFSR");
116 >    ratioFSR->Divide(negmll_FSR);
117 >    ratioFSR->GetYaxis()->SetTitle("obs/pred");
118 >    ratioFSR->GetYaxis()->SetRangeUser(0,2);
119 >    ratioFSR->Draw("e1");
120 >    oneline->Draw();
121 >    FSR->Draw();
122 >    
123 >    
124 >    CompleteSave(can,"PhotonStudies/Recophotons/ComparePredObsZmass");
125 >    
126 >    
127 >    cout << __FUNCTION__ << " : Still need to delete objects" << endl;
128 >    delete can;
129 > }
130 >
131 > void RecoPhotons::FSRPhotonVsUnmatchedVsMatchedPhoton() {
132 >    
133 >    TCanvas *can = new TCanvas("can","can");
134 >    
135 >    can->SetLogy(1);
136 >    TH1F *MatchedFSR  = PhotonSamples.Draw("MatchedFSR" ,"min(recoPhotonsDR1,recoPhotonsDR2)",20,0,2,"min(#Delta R (#gamma,l_{1}), #Delta R (#gamma,l_{2}))" ,  "events", "recoPhotonGenIndex>=0 && genPhotonsIsFSR[recoPhotonGenIndex]",mc,PlottingSetup::luminosity);
137 >    TH1F *MatchednoFSR  = PhotonSamples.Draw("MatchednoFSR" ,"min(recoPhotonsDR1,recoPhotonsDR2)",20,0,2,"min(#Delta R (#gamma,l_{1}), #Delta R (#gamma,l_{2}))" ,  "events", "recoPhotonGenIndex>=0 && !genPhotonsIsFSR[recoPhotonGenIndex]",mc,PlottingSetup::luminosity);
138 >    TH1F *UnMatched  = PhotonSamples.Draw("UnMatched" ,"min(recoPhotonsDR1,recoPhotonsDR2)",20,0,2,"min(#Delta R (#gamma,l_{1}), #Delta R (#gamma,l_{2}))" ,  "events", "recoPhotonGenIndex<0",mc,PlottingSetup::luminosity);
139 >        
140 >    MatchedFSR->SetLineColor(TColor::GetColor("#01DF01"));//green
141 >    MatchednoFSR->SetLineColor(TColor::GetColor("#FF4000"));//red
142 >    UnMatched->SetLineColor(TColor::GetColor("#194fff"));//blue
143 >    
144 >    TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
145 >    leg->SetFillColor(kWhite);
146 >    leg->AddEntry(MatchedFSR,"Matched: FSR","l");
147 >    leg->AddEntry(MatchednoFSR,"Matched: no FSR","l");
148 >    leg->AddEntry(UnMatched,"Unmatched #gamma","l");
149 >    
150 >    MatchedFSR->DrawNormalized("histo");
151 >    MatchednoFSR->DrawNormalized("same,histo");
152 >    UnMatched->DrawNormalized("same,histo");
153 >    leg->SetHeader("DY MC:");
154 >    leg->Draw("same,histo");
155 >    DrawMCPrelim();
156 >    
157 >    CompleteSave(can,"PhotonStudies/Recophotons/DeltaR_Comparison_Matched_Unmatched_FSR");
158 >    
159 >    delete can;
160 >    
161 >    
162 > }
163 >
164 > void RecoPhotons::FSRPhotonVsUnmatchedVsMatchedPhoton2D() {
165 >    
166 >    TCanvas *can = new TCanvas("can","can");
167 >    can->SetLogy(1);
168 >    TH1F *MatchedFSR  = PhotonSamples.Draw("MatchedFSR" ,"min(recoPhotonsDR1,recoPhotonsDR2)",100,0,2,"min(#Delta R (#gamma,l_{1}), #Delta R (#gamma,l_{2}))" ,  "events", "recoPhotonGenIndex>=0 && genPhotonsIsFSR[recoPhotonGenIndex]&&recoPhotonsDR1>0&&recoPhotonsDR2>0",mc,PlottingSetup::luminosity);
169 >    TH1F *MatchednoFSR  = PhotonSamples.Draw("MatchednoFSR" ,"min(recoPhotonsDR1,recoPhotonsDR2)",100,0,2,"min(#Delta R (#gamma,l_{1}), #Delta R (#gamma,l_{2}))" ,  "events", "recoPhotonGenIndex>=0 && !genPhotonsIsFSR[recoPhotonGenIndex]&&recoPhotonsDR1>0&&recoPhotonsDR2>0",mc,PlottingSetup::luminosity);
170 >    TH1F *UnMatched  = PhotonSamples.Draw("UnMatched" ,"min(recoPhotonsDR1,recoPhotonsDR2)",100,0,2,"min(#Delta R (#gamma,l_{1}), #Delta R (#gamma,l_{2}))" ,  "events", "recoPhotonGenIndex<0&&recoPhotonsDR1>0&&recoPhotonsDR2>0",mc,PlottingSetup::luminosity);
171 >    
172 >    MatchedFSR->SetLineColor(TColor::GetColor("#01DF01"));//green
173 >    MatchednoFSR->SetLineColor(TColor::GetColor("#FF4000"));//red
174 >    UnMatched->SetLineColor(TColor::GetColor("#194fff"));//blue
175 >    
176 >    TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
177 >    leg->SetFillColor(kWhite);
178 >    leg->AddEntry(MatchedFSR,"Matched: FSR","l");
179 >    leg->AddEntry(MatchednoFSR,"Matched: no FSR","l");
180 >    leg->AddEntry(UnMatched,"Unmatched #gamma","l");
181 >    
182 >    MatchedFSR->DrawNormalized("histo");
183 >    MatchednoFSR->DrawNormalized("same,histo");
184 >    UnMatched->DrawNormalized("same,histo");
185 >    leg->SetHeader("DY MC:");
186 >    leg->Draw("same,histo");
187 >    DrawMCPrelim();
188 >    
189 >    CompleteSave(can,"PhotonStudies/Recophotons/DeltaR_Comparison_Matched_Unmatched_FSR");
190 >    
191 >    
192 >    delete can;
193 >    
194 >    TCanvas *can2 = new TCanvas("can2","can2",1000,1000);
195 >    can2->Divide(2,2);
196 >    can2->cd(1);
197 >    
198 >    
199 >    vector<float> binX, binY;
200 >    int nBins=40;
201 >    for(int i=0;i<nBins;i++) {
202 >        binX.push_back(i*(3.0/nBins));
203 >        binY.push_back(i*(3.0/nBins));
204 >    }
205 >    
206 >    cout << __LINE__ << endl;
207 >    TH2F *tMatchedFSR   = PhotonSamples.Draw("tMatchedFSR",  "recoPhotonsDR1:recoPhotonsDR2", binX,binY, "#Delta R (#gamma^{FSR},l_{1})","#Delta R (#gamma^{FSR},l_{2})",BaseCut&&TCut("recoPhotonGenIndex>=0 && genPhotonsIsFSR[recoPhotonGenIndex]&&recoPhotonsDR1>0&&recoPhotonsDR2>0"),mc,PlottingSetup::luminosity, PhotonSamples.FindSample("o"),false);
208 >    cout << __LINE__ << endl;
209 >    TH2F *tMatchednoFSR   = PhotonSamples.Draw("tMatchednoFSR",  "recoPhotonsDR1:recoPhotonsDR2", binX,binY, "#Delta R (#gamma^{FSR},l_{1})","#Delta R (#gamma^{FSR},l_{2})",BaseCut&&TCut("recoPhotonGenIndex>=0 && !genPhotonsIsFSR[recoPhotonGenIndex]&&recoPhotonsDR1>0&&recoPhotonsDR2>0"),mc,PlottingSetup::luminosity, PhotonSamples.FindSample("o"),false);
210 >    cout << __LINE__ << endl;
211 >    TH2F *tUnMatched   = PhotonSamples.Draw("tUnMatched",  "recoPhotonsDR1:recoPhotonsDR2", binX,binY, "#Delta R (#gamma^{FSR},l_{1})","#Delta R (#gamma^{FSR},l_{2})",BaseCut&&TCut("recoPhotonGenIndex<0&&recoPhotonsDR1>0&&recoPhotonsDR2>0"),mc,PlottingSetup::luminosity, PhotonSamples.FindSample("o"),false);
212 >    cout << __LINE__ << endl;
213 >    
214 >    
215 >    
216 >    
217 >    can2->cd(1);
218 >    can2->cd(1)->SetRightMargin(0.14);
219 >    tMatchedFSR->Draw("COLZ");
220 >    TText *ma = write_title("Matched FSR photons");
221 >    ma->Draw();
222 >    can2->cd(2);
223 >    can2->cd(2)->SetRightMargin(0.14);
224 >    cout << __LINE__ << endl;
225 >    tMatchednoFSR->Draw("COLZ");
226 >    TText *man = write_title("Matched non-FSR phtons");
227 >    man->Draw();
228 >    can2->cd(3);
229 >    can2->cd(3)->SetRightMargin(0.14);
230 >    cout << __LINE__ << endl;
231 >    TText *unm = write_title("Unmatched photons");
232 >    tUnMatched->Draw("COLZ");
233 >    unm->Draw();
234 >    can2->cd(4);
235 >    can2->cd(4)->SetLogy(1);
236 >    MatchedFSR->DrawNormalized("histo");
237 >    MatchednoFSR->DrawNormalized("same,histo");
238 >    UnMatched->DrawNormalized("same,histo");
239 >    leg->SetHeader("DY MC:");
240 >    leg->Draw("same,histo");
241 >    DrawMCPrelim();
242 >    
243 >    cout << __LINE__ << endl;
244 >    CompleteSave(can2,"PhotonStudies/Recophotons/DeltaR_Comparison_Matched_Unmatched_FSR__2D");
245 >    
246 >    delete can2;
247 > }
248 >
249 > void RecoPhotons::NumberOfPhotonsPerEvent() {
250 >    cout << __FUNCTION__ << " not implemented." << endl;
251 > }
252 >
253 > void RecoPhotons::CharacterizingFSRPhotons() {
254 >    cout << __FUNCTION__ << " not implemented." << endl;
255 > }
256 >
257 > void RecoPhotons::dR2DPlots() {
258 >    cout << __FUNCTION__ << " not implemented." << endl;
259 > }
260 >
261 > void RecoPhotons::MakeCharacterizationPlot(string variable, int nbins, float min, float max, string label, string saveas) {
262 >    cout << __FUNCTION__ << " not implemented." << endl;
263 > }
264 >
265 > void RecoPhotons::CompareRecoJZBDistributions() {
266 >    cout << __FUNCTION__ << " not implemented." << endl;
267 > }
268 >
269 > void RecoPhotons::EfficiencyPurity() {
270 >    cout << __FUNCTION__ << " not implemented." << endl;
271 > }
272 >
273 > void CompareZmassToFinalStateLeptons() {
274 >    
275 >    TCanvas *can = new TCanvas("can","can");
276 >    can->SetLogy(1);
277 >    
278 >    TH1F *FullZ  = PhotonSamples.Draw("FullZ" ,"genFSRmllg",36,20,200,"m_{ll} [GeV]" ,  "events", "genFSRmllg>20&&genFSRmllg<200",mc,PlottingSetup::luminosity);
279 >    TH1F *RecoLL = PhotonSamples.Draw("RecoLL","mll",       36,20,200,"m_{ll} [GeV]"      , "events","mll>20&&mll<200",mc,PlottingSetup::luminosity);
280 >    TH1F *FSRRll = PhotonSamples.Draw("FSRRll","mllg",      36,20,200,"m_{ll#gamma} [GeV]", "events","mllg>20&&mllg<200",mc,PlottingSetup::luminosity);
281 >    
282 >    FullZ->SetLineColor(TColor::GetColor("#01DF01"));//green
283 >    RecoLL->SetLineColor(TColor::GetColor("#FF4000"));//red
284 >    FSRRll->SetLineColor(TColor::GetColor("#194fff"));//blue
285 >    
286 >    TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
287 >    leg->SetFillColor(kWhite);
288 >    leg->AddEntry(FullZ,"gen Z mass","l");
289 >    leg->AddEntry(RecoLL,"reco m_{ll}","l");
290 >    leg->AddEntry(FSRRll,"reco m_{ll} w/ FSR recovery","l");
291 >    
292 >    FullZ->DrawNormalized("histo");
293 >    RecoLL->DrawNormalized("same,histo");
294 >    FSRRll->DrawNormalized("same,histo");
295 >    leg->SetHeader("DY MC:");
296 >    leg->Draw("same,histo");
297 >    DrawMCPrelim();
298 >    
299 >    CompleteSave(can,"PhotonStudies/Recophotons/CompareZmassToFinalStateLeptons");
300 >    
301 >    
302 >    TH1F *RecoLL2 = PhotonSamples.Draw("RecoLL2","mll-genFSRmllg",       50,-100,100,"m_{ll}-m_{Z}^{gen} [GeV]", "events","mll>20&&genFSRmllg>20", mc,PlottingSetup::luminosity);
303 >    TH1F *FSRRll2 = PhotonSamples.Draw("FSRRll2","mllg-genFSRmllg",      50,-100,100,"m_{ll}-m_{Z}^{gen} [GeV]", "events","mllg>20&&genFSRmllg>20",mc,PlottingSetup::luminosity);
304 >    
305 >    RecoLL2->SetLineColor(TColor::GetColor("#FF4000"));//red
306 >    FSRRll2->SetLineColor(TColor::GetColor("#194fff"));//blue
307 >    
308 >    TH1F *phony = new TH1F("phony","",1,0,1);
309 >    phony->SetLineColor(kWhite);
310 >    
311 >    stringstream summary;
312 >    summary << "Before: Mean @ " << std::setprecision(3) << RecoLL2->GetMean() << ", RMS: " << std::setprecision(3) << RecoLL2->GetRMS();
313 >    stringstream summary2;
314 >    summary2 << "After:  Mean @ " << std::setprecision(3) << FSRRll2->GetMean() << ", RMS: " << std::setprecision(3) << FSRRll2->GetRMS() << "";
315 >    
316 >    TLegend *leg2 = new TLegend(0.6,0.7,0.89,0.89);
317 >    leg2->SetFillColor(kWhite);
318 >    leg2->AddEntry(FullZ,"gen Z mass","l");
319 >    leg2->AddEntry(RecoLL,"m_{ll}-m{Z}^{gen}: no FSRR","l");
320 >    leg2->AddEntry(FSRRll,"m_{ll}-m{Z}^{gen}: FSRR","l");
321 >    leg2->AddEntry(phony, summary.str().c_str(),"l");
322 >    leg2->AddEntry(phony, summary2.str().c_str(),"l");
323 >    
324 >    RecoLL2->DrawNormalized("histo");
325 >    FSRRll2->DrawNormalized("same,histo");
326 >    leg2->SetHeader("DY MC:");
327 >    leg2->Draw("same,histo");
328 >    DrawMCPrelim();
329 >    
330 >    cout << "Mean (before recovery) : " << RecoLL2->GetMean() << " , RMS: " << RecoLL2->GetRMS() << endl;
331 >    cout << "Mean (after recovery) : " <<  FSRRll2->GetMean() << " , RMS: " << FSRRll2->GetRMS() << endl;
332 >    
333 >    
334 >    CompleteSave(can,"PhotonStudies/Recophotons/CompareZmassToFinalStateLeptons_MassDiff");
335 >    
336 >    
337 >    
338 >    delete RecoLL;
339 >    delete RecoLL2;
340 >    delete FullZ;
341 >    delete FSRRll;
342 >    delete FSRRll2;
343 >    delete phony;
344 >    delete can;
345   }
346  
347   void RecoPhotons::RecoLevelStudies() {
348 <  cout << "RecoLevelStudies have not been implemented yet" << endl;
348 >    TCut essential_bkp = essentialcut;
349 >    essentialcut=TCut("");
350 >    
351 >    JZBvsmllClosure();
352 > /*    CompareZmassToFinalStateLeptons();
353 >    FSRPhotonVsUnmatchedVsMatchedPhoton();
354 >    FSRPhotonVsUnmatchedVsMatchedPhoton2D();
355 >    NumberOfPhotonsPerEvent();
356 >    CharacterizingFSRPhotons();
357 >    CompareRecoJZBDistributions();
358 >    
359 >    dR2DPlots();
360 >    
361 >    EfficiencyPurity();*/
362 >    
363 >    essentialcut=essential_bkp;
364   }
365  
366  
367   void GenPhotons::CompareZmassToFinalStateLeptons() {
368 <  
369 <  TCanvas *can = new TCanvas("can","can");
370 <  can->SetLogy(1);
371 <  
372 <  TH1F *FullZ = PhotonSamples.Draw("FullZ","genFSRmllg",36,20,200,"m_{Z}^{gen} (GeV)",  "events","genFSRmllg>20&&genFSRmllg<200",mc,PlottingSetup::luminosity);
373 <  TH1F *GenLL = PhotonSamples.Draw("GenLL","genFSRmll", 36,20,200,"m_{ll}^{gen} (GeV)", "events","genFSRmllg>20&&genFSRmllg<200",mc,PlottingSetup::luminosity);
374 <  
375 <  FullZ->SetLineColor(TColor::GetColor("#01DF01"));
376 <  GenLL->SetLineColor(TColor::GetColor("#FF4000"));
377 <  
378 <  TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
379 <  leg->AddEntry(FullZ,"gen Z mass","l");
380 <  leg->AddEntry(GenLL,"Mass from final state gen leptons","l");
381 <  
382 <  FullZ->Draw();
383 <  GenLL->Draw("same");
384 <  leg->Draw("same");
385 <  
386 <  CompleteSave(can,"PhotonStudies/CompareZmassToFinalStateLeptons");
387 <  
388 < }  
368 >    
369 >    TCanvas *can = new TCanvas("can","can");
370 >    can->SetLogy(1);
371 >    
372 >    TH1F *FullZ = PhotonSamples.Draw("FullZ","genFSRmllg",36,20,200,"m_{Z}^{gen} [GeV]",  "events", "genFSRmllg>20&&genFSRmllg<200",mc,PlottingSetup::luminosity);
373 >    TH1F *GenLL = PhotonSamples.Draw("GenLL","genFSRmll", 36,20,200,"m_{ll}^{gen} [GeV]", "events","genFSRmllg>20&&genFSRmllg<200",mc,PlottingSetup::luminosity);
374 >    
375 >    FullZ->SetLineColor(TColor::GetColor("#01DF01"));
376 >    GenLL->SetLineColor(TColor::GetColor("#FF4000"));
377 >    
378 >    TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
379 >    leg->SetFillColor(kWhite);
380 >    leg->AddEntry(FullZ,"gen Z mass","l");
381 >    leg->AddEntry(GenLL,"gen m_{ll}","l");
382 >    
383 >    FullZ->Draw("histo");
384 >    GenLL->Draw("same,histo");
385 >    leg->SetHeader("DY MC:");
386 >    leg->Draw("same,histo");
387 >    DrawMCPrelim();
388 >    
389 >    CompleteSave(can,"PhotonStudies/GenPhotons/CompareZmassToFinalStateLeptons");
390 >    delete can;
391 > }
392 >
393 >
394 > void GenPhotons::NumberOfPhotonsPerEvent() {
395 >    TCanvas *can = new TCanvas("can","can");
396 >    can->SetLogy(1);
397 >    
398 >    TH1F *FSR = PhotonSamples.Draw("FSR","genFSRNparticles"  ,11,-0.5,10.5,"N(#gamma)", "events", BaseCut,mc,PlottingSetup::luminosity);
399 >    TH1F *all = PhotonSamples.Draw("all","genPhotonsNPhotons",11,-0.5,10.5,"N(#gamma)", "events", BaseCut,mc,PlottingSetup::luminosity);
400 >    
401 >    FSR->SetLineColor(TColor::GetColor("#01DF01"));
402 >    all->SetLineColor(TColor::GetColor("#FF4000"));
403 >    
404 >    TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
405 >    leg->SetFillColor(kWhite);
406 >    leg->AddEntry(FSR,"FSR photons","l");
407 >    leg->AddEntry(all,"all photons","l");
408 >    
409 >    FSR->Draw("histo");
410 >    all->Draw("same,histo");
411 >    leg->SetHeader("DY MC:");
412 >    leg->Draw("same,histo");
413 >    DrawMCPrelim();
414 >    
415 >    CompleteSave(can,"PhotonStudies/GenPhotons/NPhotons");
416 >    
417 >    
418 >    TH1F *PhotonPt[10];
419 >    for(int i=0;i<10;i++) {
420 >        PhotonPt[i] = PhotonSamples.Draw("FSR","genPhotonsPt",25,0,100,"N(#gamma)",  "events", FSRcut && TCut(("genFSRNparticles=="+any2string(i)).c_str()),mc,PlottingSetup::luminosity);
421 >    }
422 >    
423 >    int DrawCounter=0;
424 >    
425 >    TLegend *leg2 = make_legend();
426 >    for(int i=0;i<10;i++) {
427 >        if(PhotonPt[i]->Integral()>50) {
428 >            PhotonPt[i]->SetLineColor(DrawCounter+1);
429 >            leg2->AddEntry(PhotonPt[i],(any2string(i)+" FSR Photons").c_str(),"l");
430 >            if(DrawCounter) PhotonPt[i]->DrawNormalized("same,hist");
431 >            else PhotonPt[i]->DrawNormalized("hist");
432 >            DrawCounter++;
433 >        }
434 >    }
435 >    
436 >    leg2->Draw();
437 >    
438 >    leg2->SetHeader("DY MC (N_{jets} #geq 0):");
439 >    DrawMCPrelim();
440 >    
441 >    CompleteSave(can,"PhotonStudies/GenPhotons/PhotonPtvsNPhotons");
442 >    
443 >    delete can;
444 > }
445 >
446 > void GenPhotons::MakeCharacterizationPlot(string variable, int nbins, float min, float max, string label, string saveas) {
447 >    TCanvas *can = new TCanvas("can","can");
448 >    can->SetLogy(1);
449 >    
450 >    TH1F *dR_FSR  = PhotonSamples.Draw("dR_FSR", variable, nbins, min, max, label, "events",BaseCut && FSRcut,  mc,PlottingSetup::luminosity);
451 >    TH1F *dR_nFSR = PhotonSamples.Draw("dR_nFSR", variable, nbins, min, max, label, "events",BaseCut && NoFSRcut, mc,PlottingSetup::luminosity);
452 >    
453 >    dR_FSR->SetLineColor(TColor::GetColor("#01DF01"));
454 >    dR_nFSR->SetLineColor(TColor::GetColor("#FF4000"));
455 >    
456 >    TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
457 >    leg->SetFillColor(kWhite);
458 >    leg->AddEntry(dR_FSR,"FSR #gamma's","l");
459 >    leg->AddEntry(dR_nFSR,"non-FSR #gamma's","l");
460 >    
461 >    if(dR_nFSR->GetMinimum()<dR_FSR->GetMinimum()) dR_FSR->SetMinimum(dR_nFSR->GetMinimum());
462 >    if(dR_nFSR->GetMaximum()>dR_FSR->GetMaximum()) dR_FSR->SetMaximum(dR_nFSR->GetMaximum());
463 >    
464 >    
465 >    dR_FSR->Draw("histo");
466 >    dR_nFSR->Draw("same,histo");
467 >    leg->SetHeader("DY MC:");
468 >    leg->Draw("same,histo");
469 >    DrawMCPrelim();
470 >    
471 >    CompleteSave(can,"PhotonStudies/GenPhotons/CharacterizingFSRPhotons/"+saveas);
472 >    delete dR_FSR;
473 >    delete dR_nFSR;
474 >    delete can;
475 > }
476 >
477 >
478 > void GenPhotons::CharacterizingFSRPhotons() {
479 >    MakeCharacterizationPlot("min(genPhotonsdR1,genPhotonsdR2)",100,0,5,"min(#Delta R (#gamma,l_{1}), #Delta R (#gamma,l_{2}))", "MinDR");
480 >    MakeCharacterizationPlot("genPhotonsPt",100,0,100,"p_{T}^{#gamma}", "PhotonPt");
481 >    MakeCharacterizationPlot("genPhotonsEta",100,-5,5,"#eta^{#gamma}", "PhotonEta");
482 >    MakeCharacterizationPlot("genPhotonsPhi",40,-4,4,"#phi^{#gamma}", "PhotonPhi");
483 > }
484 >
485 > void GenPhotons::CompareGenJZBDistributions() {
486 >    TCanvas *can = new TCanvas("can","can");
487 >    can->SetLogy(1);
488 >    
489 >    TH1F *dR_FSR  = PhotonSamples.Draw("dR_FSR", "genFSRjzbG", 30, -30, 30, "JZB^{gen}", "events",BaseCut && TCut("abs(genFSRjzbG)>0.1"),   mc,PlottingSetup::luminosity);
490 >    TH1F *dR_nFSR = PhotonSamples.Draw("dR_nFSR", "genFSRjzb", 30, -30, 30, "JZB^{gen}", "events",BaseCut && TCut("abs(genFSRjzb)>0.1") , mc,PlottingSetup::luminosity);
491 >    
492 >    dR_FSR->SetLineColor(TColor::GetColor("#01DF01"));
493 >    dR_nFSR->SetLineColor(TColor::GetColor("#FF4000"));
494 >    
495 >    TLegend *leg = new TLegend(0.7,0.7,0.89,0.89);
496 >    leg->SetFillColor(kWhite);
497 >    leg->AddEntry(dR_FSR,"With FSR recovery","l");
498 >    leg->AddEntry(dR_nFSR,"Without FSR recovery","l");
499 >    
500 >    if(dR_nFSR->GetMinimum()<dR_FSR->GetMinimum()) dR_FSR->SetMinimum(dR_nFSR->GetMinimum());
501 >    if(dR_nFSR->GetMaximum()>dR_FSR->GetMaximum()) dR_FSR->SetMaximum(dR_nFSR->GetMaximum());
502 >    
503 >    
504 >    dR_FSR->Draw("histo");
505 >    dR_nFSR->Draw("same,histo");
506 >    leg->SetHeader("DY MC:");
507 >    leg->Draw("same,histo");
508 >    DrawMCPrelim();
509 >    
510 >    CompleteSave(can,"PhotonStudies/GenPhotons/GenJZBDistributions");
511 >    delete dR_FSR;
512 >    delete dR_nFSR;
513 >    delete can;
514 > }
515 >
516 > void GenPhotons::dR2DPlots() {
517 >    TCanvas *can = new TCanvas("can","can",1000,500);
518 >    can->Divide(2,1);
519 >    can->cd(1);
520 >    
521 >    vector<float> binX, binY;
522 >    for(int i=0;i<=31;i++) {
523 >        binX.push_back(i*0.1);
524 >        binY.push_back(i*0.1);
525 >    }
526 >    
527 >    TH2F *dR_FSR   = PhotonSamples.Draw("dR_FSR",  "genPhotonsdR1:genPhotonsdR2", binX,binY, "#Delta R (#gamma^{FSR},l_{1})","#Delta R (#gamma^{FSR},l_{2})",BaseCut && FSRcut   ,mc,PlottingSetup::luminosity, PhotonSamples.FindSample("o"),false);
528 >    TH2F *dR_nFSR  = PhotonSamples.Draw("dR_nFSR", "genPhotonsdR1:genPhotonsdR2", binX,binY, "#Delta R (#gamma^{non-FSR},l_{1})","#Delta R (#gamma^{non-FSR},l_{2})",BaseCut && NoFSRcut ,mc,PlottingSetup::luminosity, PhotonSamples.FindSample("o"),false);
529 >    
530 >    can->cd(1);
531 >    dR_FSR->Draw("COLZ");
532 >    DrawMCPrelim();
533 >    
534 >    can->cd(2);
535 >    dR_nFSR->Draw("COLZ");
536 >    DrawMCPrelim();
537 >    
538 >    CompleteSave(can,"PhotonStudies/GenPhotons/dR_in_2D");
539 >    delete dR_FSR;
540 >    delete dR_nFSR;
541 >    delete can;
542 > }
543 >
544 >
545 > void GenPhotons::EfficiencyPurity() {
546 >    TCanvas *can = new TCanvas("can","can");
547 >    vector<float> bindR, binPt;
548 >    for(int i=0;i<=100;i++) bindR.push_back(i*0.05);
549 >    binPt.push_back(0);binPt.push_back(10);binPt.push_back(50);binPt.push_back(100);binPt.push_back(1000);binPt.push_back(1010);
550 >    //    binPt.push_back(0);binPt.push_back(1000);binPt.push_back(1001);
551 >    
552 >    
553 >    TGraph *purity[binPt.size()];
554 >    TGraph *efficiency[binPt.size()];
555 >    TGraph *purity_n_efficiency[binPt.size()];
556 >    TH1F *FSR[binPt.size()];
557 >    TH1F *Any[binPt.size()];
558 >    for(int i=1;i<binPt.size();i++) {
559 >        //do this for each pt range
560 >        float ptlow=binPt[i-1];
561 >        float pthi =binPt[i];
562 >        stringstream ptcut;
563 >        ptcut << "genPhotonsPt>" << ptlow << " && genPhotonsPt<" << pthi << " && min(genPhotonsdR1,genPhotonsdR2)<5";
564 >        
565 >        FSR[i-1]  = PhotonSamples.Draw("FSR", "min(genPhotonsdR1,genPhotonsdR2)", bindR,"min DR","events",TCut(ptcut.str().c_str()) && BaseCut && FSRcut ,mc,PlottingSetup::luminosity);
566 >        Any[i-1]  = PhotonSamples.Draw("Any", "min(genPhotonsdR1,genPhotonsdR2)", bindR,"min DR","events",TCut(ptcut.str().c_str()) && BaseCut  ,mc,PlottingSetup::luminosity);
567 >        
568 >        float NPhotons=0;
569 >        float NFSRPhotons=0;
570 >        float totalNFSRphotons=0;
571 >        
572 >        for(int k=1;k<FSR[i-1]->GetNbinsX();k++) totalNFSRphotons+=FSR[i-1]->GetBinContent(k);
573 >        purity[i-1]=new TGraph(FSR[i-1]->GetNbinsX());
574 >        efficiency[i-1]=new TGraph(FSR[i-1]->GetNbinsX());
575 >        purity_n_efficiency[i-1]=new TGraph(FSR[i-1]->GetNbinsX());
576 >        for(int j=1;j<=FSR[i-1]->GetNbinsX();j++) {
577 >            NPhotons+=(int)Any[i-1]->GetBinContent(j);
578 >            NFSRPhotons+=FSR[i-1]->GetBinContent(j);
579 >            //            cout << "       ::::::::" << j << " : " << NPhotons << " : " << NFSRPhotons << " : " << totalNFSRphotons << endl;
580 >            if(NPhotons>0) purity[i-1]->SetPoint(j-1,FSR[i-1]->GetXaxis()->GetBinCenter(j)+FSR[i-1]->GetBinWidth(j),(float)NFSRPhotons/NPhotons);
581 >            else purity[i-1]->SetPoint(j-1,FSR[i-1]->GetXaxis()->GetBinCenter(j)+FSR[i-1]->GetBinWidth(j),-1);
582 >            if(totalNFSRphotons>0) efficiency[i-1]->SetPoint(j,FSR[i-1]->GetXaxis()->GetBinCenter(j)+FSR[i-1]->GetBinWidth(j),(float)NFSRPhotons/totalNFSRphotons);
583 >            else efficiency[i-1]->SetPoint(j-1,FSR[i-1]->GetXaxis()->GetBinCenter(j)+FSR[i-1]->GetBinWidth(j),-1);
584 >            if(totalNFSRphotons>0&&NPhotons>0) purity_n_efficiency[i-1]->SetPoint(j-1,(float)NFSRPhotons/totalNFSRphotons,(float)NFSRPhotons/NPhotons);
585 >        }
586 >        purity[i-1]->GetXaxis()->SetTitle("min(#Delta R(#gamma,l_{1}), #Delta R(#gamma,_{2}))");
587 >        purity[i-1]->GetXaxis()->CenterTitle();
588 >        purity[i-1]->GetYaxis()->SetTitle("purity");
589 >        purity[i-1]->GetYaxis()->CenterTitle();
590 >        
591 >        efficiency[i-1]->GetXaxis()->SetTitle("min(#Delta R(#gamma,l_{1}), #Delta R(#gamma,_{2}))");
592 >        efficiency[i-1]->GetXaxis()->CenterTitle();
593 >        efficiency[i-1]->GetYaxis()->SetTitle("efficiency");
594 >        efficiency[i-1]->GetYaxis()->CenterTitle();
595 >        
596 >        purity_n_efficiency[i-1]->GetXaxis()->SetTitle("efficiency");
597 >        purity_n_efficiency[i-1]->GetXaxis()->CenterTitle();
598 >        purity_n_efficiency[i-1]->GetYaxis()->SetTitle("purity");
599 >        purity_n_efficiency[i-1]->GetYaxis()->CenterTitle();
600 >    }
601 >    
602 >    TLegend *leg = make_legend();
603 >    
604 >    
605 >    for(int i=0;i<binPt.size()-2;i++) {
606 >        stringstream description;
607 >        description << binPt[i] << " < p_{T}^{#gamma} < " << binPt[i+1];
608 >        purity[i]->SetLineColor(i+1);
609 >        purity_n_efficiency[i]->SetLineColor(i+1);
610 >        efficiency[i]->SetLineColor(i+1);
611 >        leg->AddEntry(purity[i],description.str().c_str(),"l");
612 >    }
613 >    
614 >    for(int i=0;i<binPt.size()-2;i++) {
615 >        if(i==0) purity[i]->Draw("AC");
616 >        else purity[i]->Draw("C");
617 >    }
618 >    leg->Draw();
619 >    CompleteSave(can,"PhotonStudies/GenPhotons/Purity");
620 >    
621 >    for(int i=0;i<binPt.size()-2;i++) {
622 >        if(i==0) efficiency[i]->Draw("AC");
623 >        else efficiency[i]->Draw("C");
624 >    }
625 >    leg->Draw();
626 >    CompleteSave(can,"PhotonStudies/GenPhotons/Efficiency");
627 >    
628 >    for(int i=0;i<binPt.size()-1;i++) {
629 >        if(i==0) purity_n_efficiency[i]->Draw("AC");
630 >        else purity_n_efficiency[i]->Draw("C");
631 >    }
632 >    leg->SetX2(0.8);
633 >    leg->SetX1(0.65);
634 >    leg->SetY1(0.1);
635 >    leg->SetY2(0.6);
636 >    leg->Draw();
637 >    CompleteSave(can,"PhotonStudies/GenPhotons/Purity_n_efficiency");
638 > }
639 >
640  
641   void GenPhotons::GenLevelStudies() {
642 <  CompareZmassToFinalStateLeptons();
643 <  
644 <  
642 >    TCut essential_bkp = essentialcut;
643 >    essentialcut=TCut("");
644 >    
645 >    CompareZmassToFinalStateLeptons();
646 >    NumberOfPhotonsPerEvent();
647 >    CharacterizingFSRPhotons();
648 >    CompareGenJZBDistributions();
649 >    
650 >    dR2DPlots();
651 >    
652 >    EfficiencyPurity();
653 >    
654 >    essentialcut=essential_bkp;
655   }
656  
657   void InitializePhotonSamples() {
658 <  float ZJetsCrossSection         = 3532.8;
659 <  float LowZJetsCrossSection      = 11050*0.069*1.15; // LO xs * filter eff * k-factor (same as ZJets)
660 <
661 <  PhotonSamples.AddSample("/shome/buchmann/ntuples/MC8tev/PhotonStudies/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12_DR53X-PU_S10_START53_V7A-v1.root","Z+Jets",-30346766,ZJetsCrossSection,false,false,1,kYellow);
662 <  PhotonSamples.AddSample("//shome/buchmann/ntuples/MC8tev/PhotonStudies/DYJetsToLL_M-10To50filter_8TeV-madgraph-Summer12_DR53X-PU_S10_START53_V7A-v1.root","Z+Jets",-6955271,LowZJetsCrossSection,false,false,1,kYellow);
663 <  
664 <  write_info(__FUNCTION__,"Photon samples have been initialized!");
665 <  PhotonSamples.ListSamples();
666 <  
658 >    float ZJetsCrossSection         = 3532.8;
659 >    float LowZJetsCrossSection      = 11050*0.069*1.15; // LO xs * filter eff * k-factor (same as ZJets)
660 >    
661 >    PhotonSamples.AddSample("/shome/buchmann/ntuples/MC8tev/PhotonStudies/Skim3genJets_DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12_DR53X-PU_S10_START53_V7A-v1.root","Z+Jets",-30346766,ZJetsCrossSection,false,false,1,kYellow);
662 >    PhotonSamples.AddSample("/shome/buchmann/ntuples/MC8tev/PhotonStudies/Skim3genJets_DYJetsToLL_M-10To50filter_8TeV-madgraph-Summer12_DR53X-PU_S10_START53_V7A-v1.root","Z+Jets",-6955271,LowZJetsCrossSection,false,false,1,kYellow);
663 >    
664 >    write_info(__FUNCTION__,"Photon samples have been initialized!");
665 >    PhotonSamples.ListSamples();
666 >    
667 >    GenPhotons::FSRcut   = TCut("genPhotonsIsFSR");
668 >    GenPhotons::NoFSRcut = TCut("!genPhotonsIsFSR");
669 >    GenPhotons::BaseCut = TCut("genNjets>=3 && genPhotonsPt > 5 && NgenLeps>0");
670   }
671  
672   void FSRstudy() {
673 <  InitializePhotonSamples();
674 <  
675 <  GenPhotons::GenLevelStudies();
676 <  RecoPhotons::RecoLevelStudies();
673 >    InitializePhotonSamples();
674 >    
675 >    //  GenPhotons::GenLevelStudies();
676 >    RecoPhotons::RecoLevelStudies();
677   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines