Signal region MC contributions

Kak W., 2018-09-18

  • Compare to Yongbin's result
  • assume integrated luminosity of 35.9 $fb^{-1}$

Sample Pre-selection

  • METFilter applied
  • Trigger: HLT_Ele27_eta2p1_WPTight_Gsf
  • Overlap removal with dR<0.4, in the preference of muon over electron over photon over jets
  • at least 1 photon object and only 1 electron object required
  • no muons in sample
  • cut-based working points
  • photon: medium, pT>20
  • electron: medium, pT>10 (will move to tight)
  • muon: tight, pT> 25

Signal region cuts (e channel)

Copied from Yongbin's slide

  • only one barrel photon $p_{T}$ > 50
  • inverse Z mass cut (center = 91GeV, width = 15)
  • photon passing electron veto (Pixel Seed/ CSEV)
  • $p^{miss}_{T}$ > 25GeV

⭕️ PU and NLO weights applied

MC samples:

  • ZGTo2LG_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/
  • DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/
  • WGamma
    • WGToLNuG_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8PhCutMax/
    • WGToLNuG_PtG-130_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8PhCutMaxPhCutMin/
    • WGToLNuG_PtG-500_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8PhCutMin/
  • WJets
    • WJetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8TrueHTOlapPhOlap/
    • WJetsToLNu_HT-100To200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • WJetsToLNu_HT-200To400_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • WJetsToLNu_HT-400To600_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • WJetsToLNu_HT-600To800_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • WJetsToLNu_HT-800To1200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • WJetsToLNu_HT-1200To2500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • WJetsToLNu_HT-2500ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
  • GJets
    • GJets_HT-40To100_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/
    • GJets_HT-100To200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/
    • GJets_HT-200To400_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/
    • GJets_HT-400To600_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/
    • GJets_HT-600ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/
  • DiPhotonJets_MGG-80toInf_13TeV_amcatnloFXFX_pythia8/
  • TTJets
    • TTJets_DiLept_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • TTJets_SingleLeptFromTbar_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/
    • TTJets_SingleLeptFromT_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap/

Import Samples

In [1]:
cd /home/kakw/efake/WG_Analysis/Plotting

/home/kakw/efake/WG_Analysis/Plotting
In [2]:
# old sample without MET filter /data/users/jkunkle/Resonances/LepGamma_elg_2018_03_28/ 
%run interactiveStackTree.py --baseDir  /data/users/fengyb/WGammaNtuple/LepGamma_elg_2018_08_19 \
--samplesConf Modules/Resonance.py --xsFile cross_sections/photon15.py \
--lumi 35900. --treeName UMDNTuple/EventTree --fileName tree.root --jupyt --quiet

Welcome to JupyROOT 6.10/05
-------------------------------------
 LOAD CROSS SECTION INFO
-------------------------------------
DYJetsToLL_M-50                xsec : 5.77e+03 pb Nevent : 49144274 lumi : 8.52e+03 scale :     4.21
ZGTo2LG                        xsec :      118 pb Nevent : 9321494  lumi : 7.91e+04 scale :    0.454
TTJets_DiLept                  xsec :     85.7 pb Nevent : 6094476  lumi : 7.11e+04 scale :    0.505
TTJets_SingleLeptFromT         xsec :      179 pb Nevent : 11957043 lumi : 6.69e+04 scale :    0.537
TTJets_SingleLeptFromTbar      xsec :      179 pb Nevent : 46455353 lumi :  2.6e+05 scale :    0.138
TTGJets                        xsec :      3.7 pb Nevent : 1577833  lumi : 4.27e+05 scale :   0.0841
WWTo2L2Nu                      xsec :     12.2 pb Nevent : 1        lumi :   0.0821 scale : 4.37e+05
WGToLNuG-amcatnloFXFX          xsec :      822 pb Nevent : 3235156  lumi : 3.94e+03 scale :     9.12
WGToLNuG-madgraphMLM           xsec :      405 pb Nevent : 6103817  lumi : 1.51e+04 scale :     2.38
WGToLNuG_PtG-130-amcatnloFXFX  xsec :     2.02 pb Nevent : 841701   lumi : 4.16e+05 scale :   0.0864
WGToLNuG_PtG-130-madgraphMLM   xsec :    0.795 pb Nevent : 1645059  lumi : 2.07e+06 scale :   0.0174
WGToLNuG_PtG-500-amcatnloFXFX  xsec :   0.0177 pb Nevent : 827560   lumi : 4.68e+07 scale : 0.000768
WGToLNuG_PtG-500-madgraphMLM   xsec :  0.00849 pb Nevent : 1393505  lumi : 1.64e+08 scale : 0.000219
WWG                            xsec :    0.215 pb Nevent : 827630   lumi : 3.85e+06 scale :  0.00931
WJetsToLNu-madgraphMLM         xsec : 6.15e+04 pb Nevent : 29705748 lumi :      483 scale :     74.4
WJetsToLNu-amcatnloFXFX        xsec : 6.15e+04 pb Nevent : 12965572 lumi :      211 scale :      170
WJetsToLNu_HT-100To200         xsec : 1.33e+03 pb Nevent : 10235198 lumi : 7.72e+03 scale :     4.65
WJetsToLNu_HT-200To400         xsec :      428 pb Nevent : 4950373  lumi : 1.16e+04 scale :      3.1
WJetsToLNu_HT-400To600         xsec :     63.9 pb Nevent : 1963464  lumi : 3.07e+04 scale :     1.17
WJetsToLNu_HT-600To800         xsec :     16.3 pb Nevent : 3779141  lumi : 2.32e+05 scale :    0.155
WJetsToLNu_HT-800To1200        xsec :     7.48 pb Nevent : 1544513  lumi : 2.06e+05 scale :    0.174
WJetsToLNu_HT-1200To2500       xsec :     1.79 pb Nevent : 244532   lumi : 1.36e+05 scale :    0.263
WJetsToLNu_HT-2500ToInf        xsec :   0.0404 pb Nevent : 253561   lumi : 6.28e+06 scale :  0.00572
WJetsToLNu_Pt-100To250         xsec :      676 pb Nevent : 3644567  lumi : 5.39e+03 scale :     6.66
GJets_HT-40To100               xsec : 2.07e+04 pb Nevent : 4467985  lumi :      216 scale :      167
GJets_HT-100To200              xsec : 9.23e+03 pb Nevent : 5131873  lumi :      556 scale :     64.5
GJets_HT-200To400              xsec :  2.3e+03 pb Nevent : 10036487 lumi : 4.36e+03 scale :     8.23
GJets_HT-400To600              xsec :      277 pb Nevent : 2529729  lumi : 9.12e+03 scale :     3.94
GJets_HT-600ToInf              xsec :     93.4 pb Nevent : 2463946  lumi : 2.64e+04 scale :     1.36
DiPhoton                       xsec :       84 pb Nevent : 19410854 lumi : 2.31e+05 scale :    0.155
ResonanceMass200               xsec :   0.0219 pb Nevent : 50000    lumi : 2.28e+05 scale :    0.157
ResonanceMass250               xsec :  0.00773 pb Nevent : 50000    lumi : 6.47e+05 scale :   0.0555
ResonanceMass300               xsec :  0.00307 pb Nevent : 50000    lumi : 1.63e+06 scale :   0.0221
ResonanceMass350               xsec :  0.00147 pb Nevent : 50000    lumi :  3.4e+06 scale :   0.0106
ResonanceMass400               xsec : 0.000679 pb Nevent : 50000    lumi : 7.36e+06 scale :  0.00488
ResonanceMass450               xsec :  0.00429 pb Nevent : 50000    lumi : 1.17e+06 scale :   0.0308
ResonanceMass500               xsec :  0.00242 pb Nevent : 50000    lumi : 2.07e+06 scale :   0.0174
ResonanceMass600               xsec :  0.00095 pb Nevent : 50000    lumi : 5.26e+06 scale :  0.00682
ResonanceMass700               xsec : 0.000472 pb Nevent : 50000    lumi : 1.06e+07 scale :  0.00339
ResonanceMass800               xsec :  0.00373 pb Nevent : 50000    lumi : 1.34e+06 scale :   0.0268
ResonanceMass900               xsec :  0.00113 pb Nevent : 50000    lumi : 4.42e+06 scale :  0.00812
ResonanceMass1000              xsec : 0.000748 pb Nevent : 50000    lumi : 6.69e+06 scale :  0.00537
ResonanceMass1200              xsec : 0.000262 pb Nevent : 50000    lumi : 1.91e+07 scale :  0.00188
ResonanceMass1400              xsec : 0.000117 pb Nevent : 50000    lumi : 4.27e+07 scale :  0.00084
ResonanceMass1600              xsec : 5.42e-05 pb Nevent : 50000    lumi : 9.22e+07 scale : 0.000389
ResonanceMass1800              xsec : 0.000252 pb Nevent : 50000    lumi : 1.98e+07 scale :  0.00181
ResonanceMass2000              xsec : 0.000172 pb Nevent : 50000    lumi :  2.9e+07 scale :  0.00124
ResonanceMass2200              xsec : 5.99e-05 pb Nevent : 50000    lumi : 8.34e+07 scale :  0.00043
ResonanceMass2400              xsec : 4.68e-05 pb Nevent : 50000    lumi : 1.07e+08 scale : 0.000336
ResonanceMass2600              xsec :  0.00023 pb Nevent : 50000    lumi : 2.17e+07 scale :  0.00165
ResonanceMass2800              xsec : 0.000136 pb Nevent : 50000    lumi : 3.69e+07 scale : 0.000974
ResonanceMass3000              xsec : 7.12e-05 pb Nevent : 50000    lumi : 7.02e+07 scale : 0.000511
ResonanceMass3500              xsec :  2.3e-05 pb Nevent : 50000    lumi : 2.17e+08 scale : 0.000165
ResonanceMass4000              xsec :  2.3e-05 pb Nevent : 50000    lumi : 2.17e+08 scale : 0.000165
-------------------------------------
BEGIN READING SAMPLES
-------------------------------------
AddSample: Reading SingleMuon                                [ Failed  ]
AddSample: Reading SingleElectron                            [ Failed  ]
AddSample: Reading DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading ZGTo2LG_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8  [ Success ]
AddSample: Reading TTJets_DiLept_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading TTJets_SingleLeptFromTbar_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading TTJets_SingleLeptFromT_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading TTGJets_TuneCUETP8M1_13TeV-amcatnloFXFX-madspin-pythia8  [ Success ]
AddSample: Reading WGToLNuG_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8  [ Success ]
AddSample: Reading WGToLNuG_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading WGToLNuG_PtG-130_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8  [ Success ]
AddSample: Reading WGToLNuG_PtG-130_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading WGToLNuG_PtG-500_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8  [ Success ]
AddSample: Reading WGToLNuG_PtG-500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading WGToLNuG_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8PhCutMax  [ Success ]
AddSample: Reading WGToLNuG_PtG-130_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8PhCutMaxPhCutMin  [ Failed  ]
AddSample: Reading WGToLNuG_PtG-500_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8PhCutMin  [ Success ]
AddSample: Reading WGToLNuG_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhCutMax  [ Success ]
AddSample: Reading WGToLNuG_PtG-130_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhCutMaxPhCutMin  [ Failed  ]
AddSample: Reading WGToLNuG_PtG-500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhCutMin  [ Success ]
AddSample: Reading WGToLNuG_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading WGToLNuG_PtG-130_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading WGToLNuG_PtG-500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading WJetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading WWG_TuneCUETP8M1_13TeV-amcatnlo-pythia8   [ Success ]
AddSample: Reading DiPhotonJets_MGG-80toInf_13TeV_amcatnloFXFX_pythia8  [ Success ]
AddSample: Reading WJetsToLNu_TuneCUETP8M1_13TeV-madgraphMLM-pythia8TrueHTOlapPhOlap  [ Failed  ]
AddSample: Reading WJetsToLNu_HT-100To200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading WJetsToLNu_HT-200To400_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading WJetsToLNu_HT-400To600_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading WJetsToLNu_HT-600To800_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading WJetsToLNu_HT-800To1200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading WJetsToLNu_HT-1200To2500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading WJetsToLNu_HT-2500ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8PhOlap  [ Success ]
AddSample: Reading GJets_HT-40To100_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading GJets_HT-100To200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading GJets_HT-200To400_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading GJets_HT-400To600_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading GJets_HT-600ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M200_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M200_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M250_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M250_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M300_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M300_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M350_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M350_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M400_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M400_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M450_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M450_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M500_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M500_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M600_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M600_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M700_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M700_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M800_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M800_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M900_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M900_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1000_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1000_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1200_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1200_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1400_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1400_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1600_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1600_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1800_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M1800_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2000_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2000_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2200_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2200_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2400_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2400_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2600_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2800_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M2800_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M3000_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M3500_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M3500_width5  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M4000_width0p01  [ Success ]
AddSample: Reading MadGraphChargedResonance_WGToLNu_M4000_width5  [ Success ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M200_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M200_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M250_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M250_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M300_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M300_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M350_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M350_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M400_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M400_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M450_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M450_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M500_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M500_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M600_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M600_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M700_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M700_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M800_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M800_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M900_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M900_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1000_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1000_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1200_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1200_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1400_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1400_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1600_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1600_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1800_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M1800_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2000_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2000_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2200_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2200_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2400_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2400_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2600_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2800_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M2800_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M3000_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M3500_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M3500_width5  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M4000_width0p01  [ Failed  ]
AddSample: Reading PythiaChargedResonance_WGToLNu_M4000_width5  [ Failed  ]
WARNING - Child sample, SingleMuon, does not exist!
WARNING - Child sample, SingleElectron, does not exist!
WARNING - Child sample, WGToLNuG_PtG-130-amcatnloFXFXPhCut, does not exist!
Grouping Wgamma
Grouping Z+jets
Grouping Zgamma
Grouping WgammaLO
WARNING - Child sample, WJetsToLNuTrueHTOlap, does not exist!
Grouping Wjets
Grouping TTG
Grouping GJets
Grouping GammaGamma
Grouping TTbar_DiLep
Grouping TTbar_SingleLep
Grouping AllTop
Grouping Others
Grouping MCBackground
Grouping MCBackgroundLO
Grouping JetBackground
Grouping EleFakeBackground
In [3]:
base = 'ph_n>=1 && el_n==1'
baseeta = base + ' && abs(ph_eta[0])<1.4'
base1 = 'ph_n==1 && el_n==1 && el_pt30_n==1 && mu_n==0'
base1eta = base1 + ' && ph_IsEB[0]'
passpix = '&& ph_hasPixSeed[0]==0'  #Pixel seed
failpix = '&& ph_hasPixSeed[0]==1'
passcsev = '&& ph_passEleVeto[0]==1' #CSEV
failcsev = '&& ph_passEleVeto[0]==0' 
ltmet = '&&met_pt<25'
gtmet = '&&met_pt>25'
phpt50 = "&&ph_pt[0]>50"
invZ = '&& abs(m_lep_ph-91)>15'
In [4]:
ROOT.gStyle.SetOptStat(0)
ROOT.gStyle.SetOptFit(1)
#samples.change_stack_order("GammaGamma","Wjets")
err = array('d',[0])

Fill histogram

Pixel Veto

In [5]:
selection = "PUWeight*NLOWeight*("+ base1eta + gtmet + invZ + passpix + phpt50+")"
samples.Draw("mt_lep_met_ph", selection, (100,0,1000),
             {"xlabel":"m_{T}(e,#gamma,p^{miss}_{T})","logy":1,'ymax':1e4,"blind":True},
             {'legendLoc':"Double","legendTranslateX":0.3});
samples.curr_canvases["top"]

mt_lep_met_ph : PUWeight*NLOWeight*(ph_n==1 && el_n==1 && el_pt30_n==1 && mu_n==0 && ph_IsEB[0]&&met_pt>25&& abs(m_lep_ph-91)>15&& ph_hasPixSeed[0]==0&&ph_pt[0]>50)
stack_sum 30762.6035345
True
Draw Signal hist  MadGraphResonanceMass300_width0p01
True
Draw Signal hist  MadGraphResonanceMass450_width5
True
Draw Signal hist  MadGraphResonanceMass800_width0p01
Out[5]:
Info in <TCanvas::Print>: png file /tmp/tmpdFElIz.png has been created

Event count by sample (pixel Veto)

  • Yongbin's number
    • Wgamma: 2.13e+4
    • Z+jets: 2.18e+3
    • GammaGamma: 6.37e+3 (incorr. xsec, should be ~ 2.15e+3)
    • Zgamma: 1.64+03
In [6]:
samples.print_stack_count()

         Wgamma 2.12e+04 +/-   586
         Z+jets 2.17e+03 +/-  97.1
     GammaGamma 2.15e+03 +/-  25.6
         Zgamma 1.64e+03 +/-  36.8
          GJets 1.06e+03 +/-   245
            TTG      990 +/-  17.1
         AllTop      879 +/-  20.7
          Wjets      662 +/-  52.5
MadGraphResonanceMass800_width0p01      166 +/-  2.18
MadGraphResonanceMass450_width5      146 +/-   2.2
MadGraphResonanceMass300_width0p01     80.4 +/-  1.38
===================================
          TOTAL 3.08e+04 +/-   646

Conversion Safe eVeto

In [7]:
selection = "PUWeight*NLOWeight*("+  base1eta + gtmet + invZ + passcsev + phpt50 +")"
samples.Draw("mt_lep_met_ph", selection,(70,0,3500),
             {"xlabel":"m_{T}(e,#gamma,p^{miss}_{T})","ymin":1e-2,"logy":1,"blind":True});
samples.curr_canvases["top"]

mt_lep_met_ph : PUWeight*NLOWeight*(ph_n==1 && el_n==1 && el_pt30_n==1 && mu_n==0 && ph_IsEB[0]&&met_pt>25&& abs(m_lep_ph-91)>15&& ph_passEleVeto[0]==1&&ph_pt[0]>50)
stack_sum 37363.9479397
True
Draw Signal hist  MadGraphResonanceMass300_width0p01
True
Draw Signal hist  MadGraphResonanceMass450_width5
True
Draw Signal hist  MadGraphResonanceMass800_width0p01
Out[7]:
Info in <TCanvas::Print>: png file /tmp/tmpn33Uaz.png has been created

Event count by sample (conversion safe eVeto)

In [8]:
samples.print_stack_count()

         Wgamma 2.25e+04 +/-   602
         Z+jets 5.68e+03 +/-   158
     GammaGamma 2.32e+03 +/-  26.6
         AllTop 2.06e+03 +/-  32.5
         Zgamma 1.78e+03 +/-  38.3
          GJets 1.17e+03 +/-   255
            TTG 1.05e+03 +/-  17.6
          Wjets      823 +/-  58.6
MadGraphResonanceMass800_width0p01      174 +/-  2.24
MadGraphResonanceMass450_width5      156 +/-  2.27
MadGraphResonanceMass300_width0p01     85.1 +/-  1.42
===================================
          TOTAL 3.74e+04 +/-   678

Conclusion

  • SM Wgamma MC passing signal selection:
    • Without weights, $3.36$ (pixel) $, 3.55$(CSEV) $\times10^{4}$
    • With NLO weight, $2.02$ (pixel) $, 2.14$(CSEV) $\times10^{4}$
    • With NLO and PU, $2.02$ (pixel) $, 2.13$(CSEV) $\times10^{4}$
    • other variables *, $2.12$ (pixel) $, 2.25$(CSEV) $\times10^{4}$
    • Yongbin's number is , $2.13$ (pixel) $, 2.25$(CSEV) $\times10^{4}$
  • The remaining difference may be due to re-reco

* mt_lep_met_ph -> mt_res, ph_eta cut -> ph_isEB

Full documentation: http://hepcms-hn.umd.edu/~kakw/2018_09_18_signal_region_check_echannel.html

BACKUP

In [9]:
selection = "PUWeight*NLOWeight*("+ base1eta  + invZ + passpix + phpt50+")"
samples.Draw("met_pt", selection, (60,0,300),
             {"xlabel":"E^{miss}_{T}","logy":1,"ymax":1e5,"blind":True},
             {"legendLoc":"Double",'legendTranslateX':0.25});
samples.curr_canvases["top"]

met_pt : PUWeight*NLOWeight*(ph_n==1 && el_n==1 && el_pt30_n==1 && mu_n==0 && ph_IsEB[0]&& abs(m_lep_ph-91)>15&& ph_hasPixSeed[0]==0&&ph_pt[0]>50)
stack_sum 43871.6556113
True
Draw Signal hist  MadGraphResonanceMass300_width0p01
True
Draw Signal hist  MadGraphResonanceMass450_width5
True
Draw Signal hist  MadGraphResonanceMass800_width0p01
Out[9]:
Info in <TCanvas::Print>: png file /tmp/tmpuqgq8s.png has been created
In [10]:
selection = "PUWeight*NLOWeight*("+ base1eta  + passpix + phpt50+")"
samples.Draw("m_lep_ph", selection, (80,0,400),
             {"xlabel":"M(e,#gamma)","logy":1,"ymax":1e6,"blind":True},
             {"legendLoc":"Double",'legendTranslateX':0.25});
samples.curr_canvases["top"]

m_lep_ph : PUWeight*NLOWeight*(ph_n==1 && el_n==1 && el_pt30_n==1 && mu_n==0 && ph_IsEB[0]&& ph_hasPixSeed[0]==0&&ph_pt[0]>50)
stack_sum 73648.3757711
True
Draw Signal hist  MadGraphResonanceMass300_width0p01
True
Draw Signal hist  MadGraphResonanceMass450_width5
True
Draw Signal hist  MadGraphResonanceMass800_width0p01
Out[10]:
Info in <TCanvas::Print>: png file /tmp/tmpOZm_sf.png has been created

Pixel Veto (WITHOUT weights)

In [11]:
samples.Draw("mt_lep_met_ph","("+ base1eta + gtmet + invZ + passpix + phpt50+")",
             (70,0,3500),{"xlabel":"m_{T}(e,#gamma,p^{miss}_{T})","ymin":1e-2,"logy":1,"blind":True});

mt_lep_met_ph : (ph_n==1 && el_n==1 && el_pt30_n==1 && mu_n==0 && ph_IsEB[0]&&met_pt>25&& abs(m_lep_ph-91)>15&& ph_hasPixSeed[0]==0&&ph_pt[0]>50)
stack_sum 50974.6809024
True
Draw Signal hist  MadGraphResonanceMass300_width0p01
True
Draw Signal hist  MadGraphResonanceMass450_width5
True
Draw Signal hist  MadGraphResonanceMass800_width0p01

Tranverse Mass

In [12]:
samples.curr_canvases["top"]
Out[12]:
Info in <TCanvas::Print>: png file /tmp/tmpKnJLwX.png has been created

Event count by sample (pixel Veto)

In [13]:
samples.print_stack_count()

         Wgamma 3.51e+04 +/-   565
     GammaGamma 4.11e+03 +/-  25.3
            TTG 3.24e+03 +/-  16.5
         Zgamma  2.9e+03 +/-  36.3
         Z+jets 2.53e+03 +/-   103
          GJets 1.42e+03 +/-   334
         AllTop      973 +/-  21.4
          Wjets      695 +/-  52.3
MadGraphResonanceMass800_width0p01      165 +/-   2.1
MadGraphResonanceMass450_width5      146 +/-  2.12
MadGraphResonanceMass300_width0p01     80.8 +/-  1.34
===================================
          TOTAL  5.1e+04 +/-   669

Conversion Safe eVeto (WITHOUT weights)

Tranverse Mass

In [14]:
histp = {"xlabel":"m_{T}(e,#gamma,p^{miss}_{T})","ymin":1e-2,"logy":1,"blind":True}
select = "("+ base1eta + gtmet + invZ + passcsev + phpt50 +")"
samples.Draw("mt_lep_met_ph",select,(70,0,3500),histp);
samples.curr_canvases["top"]

mt_lep_met_ph : (ph_n==1 && el_n==1 && el_pt30_n==1 && mu_n==0 && ph_IsEB[0]&&met_pt>25&& abs(m_lep_ph-91)>15&& ph_passEleVeto[0]==1&&ph_pt[0]>50)
stack_sum 59048.9188016
True
Draw Signal hist  MadGraphResonanceMass300_width0p01
True
Draw Signal hist  MadGraphResonanceMass450_width5
True
Draw Signal hist  MadGraphResonanceMass800_width0p01
Out[14]:
Info in <TCanvas::Print>: png file /tmp/tmp0a2_Uh.png has been created

Event count by sample (conversion safe eVeto)

In [15]:
samples.print_stack_count()

         Wgamma 3.72e+04 +/-   582
         Z+jets 6.33e+03 +/-   163
     GammaGamma 4.45e+03 +/-  26.3
            TTG 3.45e+03 +/-    17
         Zgamma 3.17e+03 +/-  37.9
         AllTop 2.21e+03 +/-  32.8
          GJets 1.38e+03 +/-   298
          Wjets      861 +/-  58.2
MadGraphResonanceMass800_width0p01      174 +/-  2.16
MadGraphResonanceMass450_width5      156 +/-  2.19
MadGraphResonanceMass300_width0p01     85.6 +/-  1.37
===================================
          TOTAL  5.9e+04 +/-   679
In [16]:
!bash /home/kakw/forjupyter/jupyterconvert.sh /home/kakw/public_html/ 2018_09_18_signal_region_check_echannel.ipynb

**** Following environment variables are going to be unset.
       DBS_CLIENT_CONFIG
       DBSCMD_HOME
       PYTHON_VALGRIND_SUPP
[NbConvertApp] Converting notebook 2018_09_18_signal_region_check_echannel.ipynb to slides
[NbConvertApp] Writing 488237 bytes to 2018_09_18_signal_region_check_echannel.slides.html
[NbConvertApp] Converting notebook 2018_09_18_signal_region_check_echannel.ipynb to html
[NbConvertApp] Writing 487256 bytes to 2018_09_18_signal_region_check_echannel.html
Wed Oct 17 08:59:33 EDT 2018