230 |
|
// Exclusion in M0 - M1/2 |
231 |
|
TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion", |
232 |
|
40,200,600,26,140,400); |
233 |
< |
TGraph * gexpexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 2,2); |
234 |
< |
TGraph * gobsexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL, 3,0, 2,1); |
235 |
< |
TGraph * gMCMCexpexcl = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,2); |
236 |
< |
TGraph * gMCMCobsexcl = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1); |
237 |
< |
TGraph * gCLsExpExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 1,2); |
238 |
< |
TGraph * gCLsObsExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL, 3,0, 1,1); |
239 |
< |
TGraph * gCLsExpExclHTm1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLm1sigma,3,0, 5,2); |
240 |
< |
TGraph * gCLsExpExclHTp1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLp1sigma,3,0, 5,2); |
233 |
> |
TGraph * gexpexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2); |
234 |
> |
TGraph * gexpexcl_LO = plotTools ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 2,4); |
235 |
> |
TGraph * gobsexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1); |
236 |
> |
TGraph * gMCMCexpexcl_LO = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,4); |
237 |
> |
TGraph * gMCMCobsexcl_LO = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1); |
238 |
> |
TGraph * gCLsExpExclHT_LO = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 1,4); |
239 |
> |
TGraph * gCLsExpExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2); |
240 |
> |
TGraph * gCLsObsExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1); |
241 |
> |
TGraph * gCLsExpExclHTm1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2); |
242 |
> |
TGraph * gCLsExpExclHTp1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2); |
243 |
> |
Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 ); |
244 |
> |
Smooth( gCLsObsExclHT, 33); gCLsObsExclHT->SetLineWidth( 3 ); |
245 |
> |
Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 ); |
246 |
> |
Smooth( gexpexcl, 33 ); |
247 |
> |
Smooth( gobsexcl, 33 ); |
248 |
> |
Smooth( gMCMCexpexcl_LO, 33 ); |
249 |
> |
Smooth( gCLsExpExclHTm1, 33 ); |
250 |
> |
Smooth( gCLsExpExclHTp1, 33 ); |
251 |
|
TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1); |
252 |
< |
hexcl->Draw("colz"); |
253 |
< |
if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf"); |
244 |
< |
if (gexpexcl) gexpexcl->Draw("l"); |
245 |
< |
if (gCLsExpExclHT) gCLsExpExclHT->Draw("l"); |
246 |
< |
if (gMCMCexpexcl) gMCMCexpexcl->Draw("l"); |
247 |
< |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
248 |
< |
//if (gobsexcl) gobsexcl->Draw("l"); |
249 |
< |
//if (gPLobsexcl) gPLobsexcl->Draw("l"); |
250 |
< |
//if (gFCobsexcl) gFCobsexcl->Draw("l"); |
251 |
< |
//if (gMCMCobsexcl) gMCMCobsexcl->Draw("l"); |
252 |
< |
//if (gFCexpexcl) gFCexpexcl->Draw("l"); |
253 |
< |
//if (gPLexpexcl) gPLexpexcl->Draw("l"); |
254 |
< |
//if (gFCexpexcl) gFCexpexcl->Draw("l"); |
252 |
> |
hexcl->Draw("colz"); |
253 |
> |
gCLsExp1Sigma->Draw("f"); |
254 |
|
//set old exclusion Limits |
255 |
|
TGraph* LEP_ch = set_lep_ch(50); |
256 |
|
TGraph* LEP_sl = set_lep_sl(50);//slepton curve |
270 |
|
First->GetYaxis()->SetTitle("m_{1/2} (GeV)"); |
271 |
|
TSpline3 *sFirst = new TSpline3("sFirst",First); |
272 |
|
sFirst->SetLineColor(kBlue); |
273 |
< |
sFirst->SetLineWidth(3); |
273 |
> |
//sFirst->SetLineWidth(3); |
274 |
|
TSpline3 *sSecond = new TSpline3("sSecond",Second); |
275 |
|
sSecond->SetLineColor(kBlue); |
276 |
|
sSecond->SetLineStyle(2); |
277 |
< |
sSecond->SetLineWidth(3); |
277 |
> |
//sSecond->SetLineWidth(3); |
278 |
|
TSpline3 *sThird = new TSpline3("sThird",Third); |
279 |
< |
sThird->SetLineColor(kGreen+2); |
279 |
> |
sThird->SetLineColor(kBlue); |
280 |
|
sThird->SetLineStyle(4); |
281 |
< |
sThird->SetLineWidth(3); |
281 |
> |
//sThird->SetLineWidth(3); |
282 |
|
TEV_sn_d0_1->Draw("fsame"); |
283 |
|
TEV_sn_d0_2->Draw("fsame"); |
284 |
|
TEV_sg_d0->Draw("fsame"); |
289 |
|
sFirst->Draw("same"); |
290 |
|
sSecond->Draw("same"); |
291 |
|
sThird->Draw("same"); |
292 |
+ |
if (gexpexcl) gexpexcl->Draw("l"); |
293 |
+ |
if (gCLsExpExclHT) gCLsExpExclHT->Draw("l"); |
294 |
+ |
if (gCLsExpExclHT_LO) gCLsExpExclHT_LO->Draw("l"); |
295 |
+ |
if (gMCMCexpexcl_LO) gMCMCexpexcl_LO->Draw("l"); |
296 |
+ |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
297 |
|
TLatex b; b.SetTextSize(0.02); b.SetTextColor(1); |
298 |
|
//b.DrawLatex( 10,420,"#tilde{#tau} LSP"); |
299 |
|
b.DrawLatex(230,145,"LEP2"); |
304 |
|
//b.DrawLatex( 80,180,"D0 #tilde{#nu}"); |
305 |
|
TLegend * leg = new TLegend(0.45,0.7,0.85,0.89); |
306 |
|
leg->SetBorderSize(0);leg->SetFillColor(0); |
307 |
+ |
TGraph * expLeg = (TGraph*)gCLsExpExclHT->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor()); |
308 |
|
//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l"); |
309 |
< |
if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"LO Observed (HT, CLs)","l"); |
310 |
< |
if (gCLsExpExclHT) leg->AddEntry(gCLsExpExclHT,"LO Expected (HT, CLs)","l"); |
311 |
< |
if (gexpexcl) leg->AddEntry(gexpexcl,"LO Expected (MHT, CLs)","l"); |
312 |
< |
//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l"); |
313 |
< |
//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l"); |
314 |
< |
//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"Observed (FC, RooStat)","l"); |
315 |
< |
//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, RooStat)","l"); |
316 |
< |
//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l"); |
312 |
< |
if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l"); |
313 |
< |
if (sFirst) leg->AddEntry(sFirst, "NLO RA1 Observed"); |
314 |
< |
if (sSecond) leg->AddEntry(sSecond,"NLO RA1 Expected"); |
315 |
< |
if (sThird) leg->AddEntry(sThird,"LO RA1 Observed"); |
309 |
> |
if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT, "NLO Observed (HT, CLs)","l"); |
310 |
> |
if (gCLsExpExclHT) leg->AddEntry(expLeg , "NLO Expected (HT, CLs)","l"); |
311 |
> |
if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l"); |
312 |
> |
if (gMCMCexpexcl_LO) leg->AddEntry(gMCMCexpexcl_LO, "LO Expected (HT, MCMC, RooStat)","l"); |
313 |
> |
if (gexpexcl) leg->AddEntry(gexpexcl, "NLO Expected (MHT, CLs)","l"); |
314 |
> |
if (sSecond) leg->AddEntry(sSecond, "NLO #alpha_{T} Expected"); |
315 |
> |
if (sFirst) leg->AddEntry(sFirst, "NLO #alpha_{T} Observed"); |
316 |
> |
if (sThird) leg->AddEntry(sThird, "LO #alpha_{T} Expected"); |
317 |
|
leg->Draw(); |
318 |
|
gPad->RedrawAxis(); |
319 |
|
c1->SaveAs("results_tb50/Exclusion_m0_m12_tb50.pdf"); |
340 |
|
if (gobsexcl) gobsexcl->Draw("l"); |
341 |
|
c1->SaveAs("results_tb50/SignalContamination_MHT_m0_m12_tb50.pdf"); |
342 |
|
|
343 |
+ |
// Signal JEC Uncertainty MHT |
344 |
+ |
c1->SetLogz(0); |
345 |
+ |
TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty", |
346 |
+ |
40,200,600,26,140,400); |
347 |
+ |
plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC); |
348 |
+ |
hsig_jec_mht->SetMinimum(0.0); |
349 |
+ |
hsig_jec_mht->SetMaximum(0.25); |
350 |
+ |
hsig_jec_mht->SetContour(10); |
351 |
+ |
hsig_jec_mht->Draw("colz"); |
352 |
+ |
if (gobsexcl) gobsexcl->Draw("l"); |
353 |
+ |
c1->SaveAs("results_tb50/SigJEC_MHT_m0_m12_tb50.pdf"); |
354 |
+ |
|
355 |
+ |
// Signal JEC Uncertainty HT |
356 |
+ |
c1->SetLogz(0); |
357 |
+ |
TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty", |
358 |
+ |
40,200,600,26,140,400); |
359 |
+ |
plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC); |
360 |
+ |
hsig_jec_ht->SetMinimum(0.0); |
361 |
+ |
hsig_jec_ht->SetMaximum(0.25); |
362 |
+ |
hsig_jec_ht->SetContour(10); |
363 |
+ |
hsig_jec_ht->Draw("colz"); |
364 |
+ |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
365 |
+ |
c1->SaveAs("results_tb50/SigJEC_HT_m0_m12_tb50.pdf"); |
366 |
+ |
|
367 |
+ |
|
368 |
+ |
// Signal MuIso Uncertainty MHT |
369 |
+ |
c1->SetLogz(0); |
370 |
+ |
TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
371 |
+ |
40,200,600,26,140,400); |
372 |
+ |
plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso); |
373 |
+ |
hsig_MuIso_mht->SetMinimum(0.0); |
374 |
+ |
hsig_MuIso_mht->SetMaximum(0.05); |
375 |
+ |
hsig_MuIso_mht->SetContour(5); |
376 |
+ |
hsig_MuIso_mht->Draw("colz"); |
377 |
+ |
if (gobsexcl) gobsexcl->Draw("l"); |
378 |
+ |
c1->SaveAs("results_tb50/SigMuIso_MHT_m0_m12_tb50.pdf"); |
379 |
+ |
|
380 |
+ |
// Signal MuIso Uncertainty HT |
381 |
+ |
c1->SetLogz(0); |
382 |
+ |
TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
383 |
+ |
40,200,600,26,140,400); |
384 |
+ |
plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso); |
385 |
+ |
hsig_MuIso_ht->SetMinimum(0.0); |
386 |
+ |
hsig_MuIso_ht->SetMaximum(0.05); |
387 |
+ |
hsig_MuIso_ht->SetContour(5); |
388 |
+ |
hsig_MuIso_ht->Draw("colz"); |
389 |
+ |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
390 |
+ |
c1->SaveAs("results_tb50/SigMuIso_HT_m0_m12_tb50.pdf"); |
391 |
+ |
|
392 |
+ |
|
393 |
+ |
// Signal kFactor Uncertainty MHT |
394 |
+ |
c1->SetLogz(0); |
395 |
+ |
TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
396 |
+ |
40,200,600,26,140,400); |
397 |
+ |
plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor); |
398 |
+ |
hsig_kFactor_mht->SetMinimum(0.1); |
399 |
+ |
hsig_kFactor_mht->SetMaximum(0.2); |
400 |
+ |
hsig_kFactor_mht->SetContour(10); |
401 |
+ |
hsig_kFactor_mht->Draw("colz"); |
402 |
+ |
if (gobsexcl) gobsexcl->Draw("l"); |
403 |
+ |
c1->SaveAs("results_tb50/SigkFactor_MHT_m0_m12_tb50.pdf"); |
404 |
+ |
|
405 |
+ |
// Signal kFactor Uncertainty HT |
406 |
+ |
c1->SetLogz(0); |
407 |
+ |
TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
408 |
+ |
40,200,600,26,140,400); |
409 |
+ |
plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor); |
410 |
+ |
hsig_kFactor_ht->SetMinimum(0.1); |
411 |
+ |
hsig_kFactor_ht->SetMaximum(0.2); |
412 |
+ |
hsig_kFactor_ht->SetContour(10); |
413 |
+ |
hsig_kFactor_ht->Draw("colz"); |
414 |
+ |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
415 |
+ |
c1->SaveAs("results_tb50/SigkFactor_HT_m0_m12_tb50.pdf"); |
416 |
+ |
|
417 |
|
|
418 |
|
//----------------------------------------------------------------------------------- |
419 |
|
c1->SetLogz(1); |
538 |
|
// Exclusion in squark - gluino mass ---------------------------------------- |
539 |
|
TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion", |
540 |
|
40,200,900,30,200,900); |
541 |
< |
TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 2,2); |
542 |
< |
TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL, 3,0, 2,1); |
541 |
> |
TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino, NLOExpExclCL, 3,0, 2,2); |
542 |
> |
TGraph * gexpexcl_qg_LO = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino, ExpExclCL, 3,0, 2,4); |
543 |
> |
TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino, NLOObsExclCL, 3,0, 2,1); |
544 |
|
TGraph * gMCMCexpexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2); |
545 |
|
TGraph * gMCMCobsexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1); |
546 |
< |
TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,2); |
547 |
< |
TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL, 3,0, 1,1); |
548 |
< |
TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLm1sigma,3,0, 5,2); |
549 |
< |
TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLp1sigma,3,0, 5,2); |
546 |
> |
TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 1,2); |
547 |
> |
TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,4); |
548 |
> |
TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 1,1); |
549 |
> |
TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2); |
550 |
> |
TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2); |
551 |
|
TGraph * gCLsExp1Sigma_qg = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg); |
552 |
|
hexcl_qg->Draw(""); |
553 |
|
// |
555 |
|
if (gobsexcl_qg) gobsexcl_qg->Draw("l"); |
556 |
|
if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l"); |
557 |
|
if (gexpexcl_qg) gexpexcl_qg->Draw("l"); |
558 |
+ |
if (gexpexcl_qg_LO) gexpexcl_qg_LO->Draw("l"); |
559 |
|
if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l"); |
560 |
+ |
if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l"); |
561 |
|
// |
562 |
|
TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89); |
563 |
|
leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0); |
564 |
< |
if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"LO Observed (MHT, CLs)","l"); |
565 |
< |
if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"LO Expected (MHT, CLs)","l"); |
566 |
< |
if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"LO Observed (HT, CLs)","l"); |
567 |
< |
if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"LO Expected (HT, CLs)","l"); |
564 |
> |
TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor()); |
565 |
> |
if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l"); |
566 |
> |
if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l"); |
567 |
> |
if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l"); |
568 |
> |
if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l"); |
569 |
> |
if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg_qg,"NLO Expected (HT, CLs)","lf"); |
570 |
> |
if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","l"); |
571 |
|
//if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l"); |
572 |
|
//if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l"); |
573 |
|
leg_qg->Draw(); |