ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/nowaf/StudyXClean/rootScripts/MakeMetStudy.py
(Generate patch)

Comparing UserCode/nowaf/StudyXClean/rootScripts/MakeMetStudy.py (file contents):
Revision 1.1 by nowak, Thu Feb 25 06:29:24 2010 UTC vs.
Revision 1.2 by nowak, Fri Apr 9 11:22:35 2010 UTC

# Line 376 | Line 376 | if __name__ == "__main__":
376      gene_vs_ie_same.SetTitle( "newMET = oldMET" )
377      gene_vs_ie_same.Draw( "COLZ" )
378      saves["gene_vs_ie_same"] = gene_vs_ie_same
379 <    
379 >
380 >    c_ie = ROOT.TCanvas( "c_ie", "isoe vs nonie" ,1,71,600,600 )
381 >    c_ie.Divide(2,2)
382 >    c_ie.cd(1)
383 >    ie_vs_nonie_new = fin.Get( "h_ie_vs_nonie_newMETg_" + top )
384 >    ie_vs_nonie_new.SetXTitle( "number of iso elecs" )
385 >    ie_vs_nonie_new.SetYTitle( "number of non iso elecs" )
386 >    ie_vs_nonie_new.SetTitle( "newMET > oldMET" )
387 >    ie_vs_nonie_new.Draw( "COLZ" )
388 >    c_ie.cd(2)
389 >    ie_vs_nonie_old = fin.Get( "h_ie_vs_nonie_oldMETg_" + top )
390 >    ie_vs_nonie_old.SetXTitle( "number of iso elecs" )
391 >    ie_vs_nonie_old.SetYTitle( "number of non iso elecs" )
392 >    ie_vs_nonie_old.SetTitle( "newMET < oldMET" )
393 >    ie_vs_nonie_old.Draw( "COLZ" )
394 >    saves["ie_vs_nonie_new"] = ie_vs_nonie_new
395 >    saves["ie_vs_nonie_old"] = ie_vs_nonie_old
396 >    c_ie.cd(3)
397 >    ie_vs_nonie_same = fin.Get( "h_ie_vs_nonie_sameMET_" + top )
398 >    ie_vs_nonie_same.SetXTitle( "number of iso elecs" )
399 >    ie_vs_nonie_same.SetYTitle( "number of non iso elecs" )
400 >    ie_vs_nonie_same.SetTitle( "newMET = oldMET" )
401 >    ie_vs_nonie_same.Draw( "COLZ" )
402 >    saves["ie_vs_nonie_same"] = ie_vs_nonie_same
403 >
404 >    c_ie2 = ROOT.TCanvas( "c_ie2", "isoe vs dropped elecs" ,1,71,600,600 )
405 >    c_ie2.Divide(2,2)
406 >    c_ie2.cd(1)
407 >    ie_vs_de_new = fin.Get( "h_ie_vs_de_newMETg_" + top )
408 >    ie_vs_de_new.SetXTitle( "number of iso elecs" )
409 >    ie_vs_de_new.SetYTitle( "number of dropped elecs" )
410 >    ie_vs_de_new.SetTitle( "newMET > oldMET" )
411 >    ie_vs_de_new.Draw( "COLZ" )
412 >    c_ie2.cd(2)
413 >    ie_vs_de_old = fin.Get( "h_ie_vs_de_oldMETg_" + top )
414 >    ie_vs_de_old.SetXTitle( "number of iso elecs" )
415 >    ie_vs_de_old.SetYTitle( "number of dropped elecs" )
416 >    ie_vs_de_old.SetTitle( "newMET < oldMET" )
417 >    ie_vs_de_old.Draw( "COLZ" )
418 >    saves["ie_vs_de_new"] = ie_vs_de_new
419 >    saves["ie_vs_de_old"] = ie_vs_de_old
420 >    c_ie2.cd(3)
421 >    ie_vs_de_same = fin.Get( "h_ie_vs_de_sameMET_" + top )
422 >    ie_vs_de_same.SetXTitle( "number of iso elecs" )
423 >    ie_vs_de_same.SetYTitle( "number of dropped elecs" )
424 >    ie_vs_de_same.SetTitle( "newMET = oldMET" )
425 >    ie_vs_de_same.Draw( "COLZ" )
426 >    saves["ie_vs_de_same"] = ie_vs_de_same
427 >
428 >    c_ie3 = ROOT.TCanvas( "c_ie3", "ge vs re" ,1,71,600,600 )
429 >    c_ie3.Divide(2,2)
430 >    c_ie3.cd(1)
431 >    ge_vs_re_new = fin.Get( "h_ge_vs_re_newMETg_" + top )
432 >    ge_vs_re_new.SetXTitle( "number of gen elecs" )
433 >    ge_vs_re_new.SetYTitle( "number of reco elecs" )
434 >    ge_vs_re_new.SetTitle( "newMET > oldMET" )
435 >    ge_vs_re_new.Draw( "COLZ" )
436 >    c_ie3.cd(2)
437 >    ge_vs_re_old = fin.Get( "h_ge_vs_re_oldMETg_" + top )
438 >    ge_vs_re_old.SetXTitle( "number of gen elecs" )
439 >    ge_vs_re_old.SetYTitle( "number of reco elecs" )
440 >    ge_vs_re_old.SetTitle( "newMET < oldMET" )
441 >    ge_vs_re_old.Draw( "COLZ" )
442 >    saves["ge_vs_re_new"] = ge_vs_re_new
443 >    saves["ge_vs_re_old"] = ge_vs_re_old
444 >    c_ie3.cd(3)
445 >    ge_vs_re_same = fin.Get( "h_ge_vs_re_sameMET_" + top )
446 >    ge_vs_re_same.SetXTitle( "number of gen elecs" )
447 >    ge_vs_re_same.SetYTitle( "number of reco elecs" )
448 >    ge_vs_re_same.SetTitle( "newMET = oldMET" )
449 >    ge_vs_re_same.Draw( "COLZ" )
450 >    saves["ge_vs_re_same"] = ge_vs_re_same
451 >
452 >    c_met1 = ROOT.TCanvas( "c_met1", "dj vs MET" ,1,71,600,600 )
453 >    c_met1.Divide(2,2)
454 >    c_met1.cd(1)
455 >    dj_vs_MET_new = fin.Get( "h_dj_vs_MET_newMETg_" + top )
456 >    dj_vs_MET_new.SetXTitle( "number of dropped jets" )
457 >    dj_vs_MET_new.SetYTitle( "MET" )
458 >    dj_vs_MET_new.SetTitle( "newMET > oldMET" )
459 >    dj_vs_MET_new.Draw( "COLZ" )
460 >    c_met1.cd(2)
461 >    dj_vs_MET_old = fin.Get( "h_dj_vs_MET_oldMETg_" + top )
462 >    dj_vs_MET_old.SetXTitle( "number of dropped jets" )
463 >    dj_vs_MET_old.SetYTitle( "MET" )
464 >    dj_vs_MET_old.SetTitle( "newMET < oldMET" )
465 >    dj_vs_MET_old.Draw( "COLZ" )
466 >    saves["dj_vs_MET_new"] = dj_vs_MET_new
467 >    saves["dj_vs_MET_old"] = dj_vs_MET_old
468 >    c_met1.cd(3)
469 >    dj_vs_MET_same = fin.Get( "h_dj_vs_MET_sameMET_" + top )
470 >    dj_vs_MET_same.SetXTitle( "number of dropped jets" )
471 >    dj_vs_MET_same.SetYTitle( "MET" )
472 >    dj_vs_MET_same.SetTitle( "newMET = oldMET" )
473 >    dj_vs_MET_same.Draw( "COLZ" )
474 >    saves["dj_vs_MET_same"] = dj_vs_MET_same
475 >
476 >    c_met2 = ROOT.TCanvas( "c_met2", "de vs MET" ,1,71,600,600 )
477 >    c_met2.Divide(2,2)
478 >    c_met2.cd(1)
479 >    de_vs_MET_new = fin.Get( "h_de_vs_MET_newMETg_" + top )
480 >    de_vs_MET_new.SetXTitle( "number of dropped elecs" )
481 >    de_vs_MET_new.SetYTitle( "MET" )
482 >    de_vs_MET_new.SetTitle( "newMET > oldMET" )
483 >    de_vs_MET_new.Draw( "COLZ" )
484 >    c_met2.cd(2)
485 >    de_vs_MET_old = fin.Get( "h_de_vs_MET_oldMETg_" + top )
486 >    de_vs_MET_old.SetXTitle( "number of dropped elecs" )
487 >    de_vs_MET_old.SetYTitle( "MET" )
488 >    de_vs_MET_old.SetTitle( "newMET < oldMET" )
489 >    de_vs_MET_old.Draw( "COLZ" )
490 >    saves["de_vs_MET_new"] = de_vs_MET_new
491 >    saves["de_vs_MET_old"] = de_vs_MET_old
492 >    c_met2.cd(3)
493 >    de_vs_MET_same = fin.Get( "h_de_vs_MET_sameMET_" + top )
494 >    de_vs_MET_same.SetXTitle( "number of dropped elecs" )
495 >    de_vs_MET_same.SetYTitle( "MET" )
496 >    de_vs_MET_same.SetTitle( "newMET = oldMET" )
497 >    de_vs_MET_same.Draw( "COLZ" )
498 >    saves["de_vs_MET_same"] = de_vs_MET_same
499 >
500 >
501 >
502 >    c_idde = ROOT.TCanvas( "c_nde", "number of id dropped elecs" ,1,71,600,600 )
503 >    #c_idde.Divide(2,2)
504 >
505 >    hist_list_idde = [ fin.Get( "h_de_id_newMetg_" + top ), fin.Get( "h_de_id_oldMetg_" + top ), fin.Get( "h_de_id_sameMet_" + top ) ]
506 >    color_list_idde = [50,8,12]
507 >    case_list_idde = ["newMETg", "oldMETg", "sameMET"]
508 >
509 >    c_idde.cd(1)
510 >    Def.DrawSmall( c_idde,\
511 >                   hist_list_idde, color_list_idde, case_list_idde, \
512 >                   "number of id dropped elecs", "", "",\
513 >                   False, 0, False,\
514 >                   "ur", "idde", saves)
515 >
516 >
517 >    c_control =  ROOT.TCanvas( "c_control", "gen e vs sum reco e" ,1,71,600,600 )
518 >    c_control.Divide(2,2)
519 >
520 >    c_control.cd(1)
521 >    ge_vs_sre_newMETg = fin.Get( "h_gene_vs_sumrecoe_newMETg" )
522 >    ge_vs_sre_newMETg.SetXTitle( "number of gen e" )
523 >    ge_vs_sre_newMETg.SetYTitle( "sum of all reco e" )
524 >    ge_vs_sre_newMETg.SetTitle( "newMET > oldMET" )
525 >    ge_vs_sre_newMETg.Draw( "COLZ" )
526 >    saves["ge_vs_sre_newMETg"] = ge_vs_sre_newMETg
527 >
528 >    c_control.cd(2)
529 >    ge_vs_sre_oldMETg = fin.Get( "h_gene_vs_sumrecoe_oldMETg" )
530 >    ge_vs_sre_oldMETg.SetXTitle( "number of gen e" )
531 >    ge_vs_sre_oldMETg.SetYTitle( "sum of all reco e" )
532 >    ge_vs_sre_oldMETg.SetTitle( "newMET < oldMET" )
533 >    ge_vs_sre_oldMETg.Draw( "COLZ" )
534 >    saves["ge_vs_sre_oldMETg"] = ge_vs_sre_oldMETg
535 >
536 >    c_control.cd(3)
537 >    ge_vs_sre_sameMET = fin.Get( "h_gene_vs_sumrecoe_sameMET" )
538 >    ge_vs_sre_sameMET.SetXTitle( "number of gen e" )
539 >    ge_vs_sre_sameMET.SetYTitle( "sum of all reco e" )
540 >    ge_vs_sre_sameMET.SetTitle( "newMET = oldMET" )
541 >    ge_vs_sre_sameMET.Draw( "COLZ" )
542 >    saves["ge_vs_sre_sameMET"] = ge_vs_sre_sameMET
543 >
544 >    c_iso =  ROOT.TCanvas( "c_iso", "relative isolations from dropped e" ,1,71,600,600 )
545 >    c_iso.Divide( 2,2 )
546 >    hist_list_iso_track = [ fin.Get( "h_de_trackIso_newMetg"), fin.Get( "h_de_trackIso_oldMetg"), fin.Get( "h_de_trackIso_sameMet" ) ]
547 >    c_iso.cd(1)
548 >    Def.DrawSmall( c_iso,\
549 >                   hist_list_iso_track, color_list_idde, case_list_idde, \
550 >                   "relative track iso from dropped elecs", "", "",\
551 >                   True, 0, True,\
552 >                   "ur", "trackiso", saves)
553 >
554 >    hist_list_iso_ecal = [ fin.Get( "h_de_ecalIso_newMetg"), fin.Get( "h_de_ecalIso_oldMetg"), fin.Get( "h_de_ecalIso_sameMet" ) ]
555 >    c_iso.cd(2)
556 >    Def.DrawSmall( c_iso,\
557 >                   hist_list_iso_ecal, color_list_idde, case_list_idde, \
558 >                   "relative ecal iso from dropped elecs", "", "",\
559 >                   True, 0, True,\
560 >                   "ur", "ecaliso", saves)
561 >
562 >    hist_list_iso_hcal = [ fin.Get( "h_de_hcalIso_newMetg"), fin.Get( "h_de_hcalIso_oldMetg"), fin.Get( "h_de_hcalIso_sameMet" ) ]
563 >    c_iso.cd(3)
564 >    Def.DrawSmall( c_iso,\
565 >                   hist_list_iso_hcal, color_list_idde, case_list_idde, \
566 >                   "relative hcal iso from dropped elecs", "", "",\
567 >                   True, 0, True,\
568 >                   "ur", "hcaliso", saves)
569                  
570      Def.DontQuit()

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines