Previous page

Exclusive tt analysis using full run II CMS data

Overview

This page contain analysis instructions (in development)

On this page:

Usefull Links

  • Exclusive ttbar talk (LHC WG on Forward Physics and Diffraction): link

Data and MC samples

Sample Xsec
/TTJets_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIFall17MiniAODv2-PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/MINIAODSIM 832

These samples produced privately and can be found in CERN EOS

Sample Xsec [fb]
/eos/cms/store/cmst3/group/top/exclusive_top/excl_ttbar_QED/MINIAOD/excl_ttbar_WmLep_QED/ 0.064
/eos/cms/store/cmst3/group/top/exclusive_top/excl_ttbar_QED/MINIAOD/excl_ttbar_WpLep_QED/ 0.064

Configuration file can be found in the attachment . It consist of several steps:

  • Simulation with FPMC of exclusive γγ→tt events
  • Interfacing of FPMC output with MadSpin such that only one top quark decays leptonically
  • Interfacing MadSpin output with Pythia8
  • Processing events after hadronization with full CMS detector simulation following by reconstruction
Execute run.sh script with correct arguments. The output files will be of a miniAOD.root format.

TO UPDATE
Stream Lumi
/SingleMuon/Run2017H-17Nov2017-v2/MINIAOD XXX
/HighEGJet/Run2017H-17Nov2017-v1/MINIAOD XXX

Analysis

Ntuple production

Production is performed using TopLJets2015 analysis package which use PPS reconstruction data.

PPS fast simulation

The full instruction for simulating the forward protons can be found in TaggedProtonsDirectSimulation twiki. The configuration file is attached (ProtonReco.zip). The following lines should be executed (CMSSW should point to your CMSSW setup as explained here):

cd ${CMSSW}/src
eval `scram r -sh`
cd -
cd ${CMSSW_BASE}/../cms-data
export CMSSW_SEARCH_PATH="$CMSSW_SEARCH_PATH:$(pwd -P)"
export CMSSW_DATA_PATH="/afs"
cd -
cp ${CMSSW_BASE}/src/proton_simulation_validation/settings/2017_postTS2/* .
cmsRun step5_cfg.py

NOTE: to process miniAOD, one should change the collection name to: cms.InputTag("prunedGenParticles")

Making validation plots

miniAOD

To plot top - anti-top invariant mass from the samples, either directly from tops or from the protons

scram project CMSSW CMSSW_9_4_7
cd CMSSW_9_4_7/src
cmsenv

root -l  /eos/cms/store/cmst3/group/top/exclusive_top/excl_ttbar_QED/MINIAOD/excl_ttbar_WpLep_QED//excl_ttbar_WpLep_QED_663.root
Events->SetAlias("genpar","recoGenParticles_prunedGenParticles__PAT.obj")
Events->Draw("sqrt((genpar.energy()[2]+genpar.energy()[3])**2-(genpar.px()[2]+genpar.px()[3])**2-(genpar.py()[2]+genpar.py()[3])**2-(genpar.pz()[2]+genpar.pz()[3])**2)")
Events->Draw("2*sqrt( (6500.-genpar.pz()[0])*(6500.+genpar.pz()[1]))")

miniAOD with protons

To plot proton pZ vs reconstructed pZ (similarly can be made for ξ):

cd ${CMSSW}/src
cmsenv

root -l  miniAOD_withProtons.root
Events->SetAlias("genpar","recoGenParticles_prunedGenParticles__PAT.obj")
Events->SetAlias("PUproton","recoGenParticles_genPUProtons_genPUProtons_HLT.obj")
Events->SetAlias("proton","recoForwardProtons_ctppsProtons_singleRP_CTPPSTestAcceptance.obj")
Events->Draw("proton.pz() : genpar.pz()[0]","proton.pz() > 0")

For more variables see reco::ForwardProton.

Analysis

To process nanoAODs execute the following:

ln -s /eos/user/m/mpitt/PPS_ttbar/FPMC/Merged/excl_ttbar_WmLep_QED.root MiniEvents.root
python $CMSSW_BASE/src/TopLJets2015/TopAnalysis/scripts/runLocalAnalysis.py -i MiniEvents.root -o ntuple.root --njobs 1 -q local --era era2017 -m RunExclusiveTop2020 --skimtree True

PPS reconstruction (data)

Submit to grid

After the local tests are succesfull, one can process all the data on the grid:

Plots

execute this and that

-- MichaelPitt - 2019-12-07

  • GEN.zip: production of exclusive ttbar events with CMS

Topic attachments
I Attachment History Action Size Date Who Comment
Compressed Zip archivezip GEN.zip r1 manage 17.1 K 2020-01-01 - 15:40 MichaelPitt production of exclusive ttbar events with CMS
Compressed Zip archivezip ProtonReco.zip r2 r1 manage 2.4 K 2020-01-21 - 18:51 MichaelPitt cfg file for FastSim of tagged protons
Edit | Attach | Watch | Print version | History: r11 < r10 < r9 < r8 < r7 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r11 - 2020-01-21 - MichaelPitt
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback