PPS CEP Studies



Analysis Information

  • Analysis Contacts: Dilson Damião (dilson@cernNOSPAMPLEASE.ch) and Eliza Melo (eliza@cernNOSPAMPLEASE.ch)
  • Authors: D. Damião, D. Figueiredo, S. Fonseca, E. Melo, L. Mundim, A. Santoro and A. Vilela
  • Status: Studing Vertex/Jets association (resolutions) and pile-up.
  • Indico: here
  • Git Analysis Code: here, version CMSSW_6_2_5_PPS_V14
  • Git Macro for signal: here,version CMSSW_6_2_5_PPS_V09
  • Git Macro for signal+bckg: here, version CMSSW_6_2_5_PPS_V11
  • PPS Fast Sim: here


DataFormat and Fastsim Software

  • Dataformats: Git
  • Responsibles: L. Mundim and M. Elena Pol


Useful Information




Software Instalation and Recipes

Recipe Step by Step

cmsrel CMSSW_6_2_5
cd CMSSW_6_2_5/src
cmsenv
git clone -b CMSSW_6_2_5 git@github.com:ForwardGroupBrazil/DataFormats.git DataFormats
scram b -j 9
git clone -b CMSSW_6_2_5_PPS_V14 git@github.com:ForwardGroupBrazil/PPS_Dijets_Skimming.git PPS_Dijets_Skimming
scram b -j 9


Analysis Software

Running CEP Analyzer

cd PPS_Dijets_Skimming/ExclusiveDijetsAnalysisUsingPPS/test
cmsRun DijetsAnalyzer_cfg.py

Running with options:
Run = MC_OOT_PU (MC OOT and with PU)
Run = MC_NO_OOT_NO_PU (MC without OOT and without PU)
Run = MC_NO_OOT_PU (MC without OOT and with PU)

Example:

cmsRun DijetsAnalyzer_cfg.py Run=MC_OOT_PU

Some git commands to make new branch

- i is number of current version of branch:

git branch CMSSW_6_2_5_PPS_V0(i+1)
git checkout CMSSW_6_2_5_PPS_V0(i+1)
git push origin CMSSW_6_2_5_PPS_V0(i+1)
git commit -m "my comment" *
git push --set-upstream origin CMSSW_6_2_5_PPS_V0(i+1)

Plotter histogram from PPS analysis

- Please should go to /Dropbox/PPSBrazil/macro
- Edit loopHist_PPS.C to put the name of input file
-run root:
root -b
.L loopHist_PPS.C 
loopHist_PPS(0) # default from root file wih output
loopHist_PPS(1) # default from png file wih output
loopHist_PPS(2) # default from pdf file wih output

SRMBasic commands for analysis5

lcg-cp -b -D srmv2 "srm://srm-eoscms.cern.ch:8443/srm/v2/server?SFN=/eos/cms/store/group/phys_diffraction/polme/Pomwig_SDP_552710710.root" "file:teste.root"

Example of PPS Analyzer

In git website repository, you can find an example of PPS analyzer.


CEP Analyzer

This software is used to create a NTuple with CEP analysis collections.

NTuple Branches

Branch Name Access, [unit] Type
VertexCMSVector_x x component [cm], all vertex at CMS std::vector
VertexCMSVector_y y component [cm], all vertex at CMS std::vector
VertexCMSVector_z z component [cm], all vertex at CMS std::vector
VertexGENVector_x x component [cm], all vertex GEN Level std::vector
VertexGENVector_y y component [cm], all vertex GEN Level std::vector
VertexGENVector_z z component [cm], all vertex GEN Level std::vector
AllDiffVertexZVector z difference between all neighboring vertices per event [cm] std::vector
DistanceBetweenJets z difference between Leading Jet and all other jets per event std::vector
MinDistanceZVertex z min. difference between neighboring vertices per event [cm] double
MaxDistanceZVertex z max. difference between neighboring vertices per event [cm] double
JetsPt All jets pT [GeV] std::vector
JetsEta All jets eta] std::vector
JetsPhi All jets phi [rad] std::vector
JetVertex X(), Y() and Z() All jets vertex position. Access info JetVertex.Z()[index] std::vector<TMath::Point>
TracksPerJet Number of tracks per all Jets std::vector
JetsSameVertex_pt pT [GeV] all jets close to leading Jet std::vector
JetsSameVertex_eta eta all jets close to leading Jet std::vector
JetsSameVertex_phi phi [rad ]all jets close to leading Jet std::vector
JetsSameVertex_x x [cm] vertex leading jet double
JetsSameVertex_y y [cm] vertex leading jet double
JetsSameVertex_z z [cm] vertex leading jet double
CandidatesJets_pt pT [GeV] all jets going out CMS/PPS vertex associated std::vector
CandidatesJets_eta eta all jets going out CMS/PPS vertex associated std::vector
CandidatesJets_phi phi [rad] all jets going out CMS/PPS vertex associated std::vector
PFCandidatePt All PF candidates pT [GeV] std::vector
PFCandidateEta All PF candidates eta [a.u] std::vector
nVertex Number of vertex per event int
nTracks Number of tracks per event int
MinDistance minimum distance between PPS/CMS vertex and all tracks double
MaxDistance maximum distance between PPS/CMS vertex and all tracks double
GoldenVertexZ CMS vertex in association with PPS, z component [cm] double
VertexZPPS PPS vertex, z component [cm] double
xiPPSArmB xi, PPS, ArmB double
xiPPSArmF xi, PPS, ArmF double
tPPSArmB t [GeV/c^2], PPS, ArmB double
tPPSArmF t [GeV/c^2], PPS, ArmF double
xPPSArmBDet1 x [mm], ArmB/Detector 1 double
xPPSArmFDet1 x [mm], ArmF/Detector 1 double
yPPSArmBDet1 y [mm], ArmB/Detector 1 double
yPPSArmFDet1 y [mm], ArmF/Detector 1 double
xPPSArmBDet2 x [mm], ArmB/Detector 2 double
xPPSArmFDet2 x [mm], ArmF/Detector 2 double
yPPSArmBDet2 y [mm], ArmB/Detector 2 double
yPPSArmFDet2 y [mm], ArmF/Detector 2 double
xPPSArmFToF x [mm], ArmF/ToF double
yPPSArmFToF y [mm], ArmF/ToF double
xPPSArmBToF x [mm], ArmB/ToF double
yPPSArmBToF y [mm], ArmB/ToF double
stopPPSArmFTrkDet1 status, p stopped at ArmF/TrkDet1 bool
stopPPSArmFTrkDet2 status, p stopped at ArmF/TrkDet2 bool
stopPPSArmBTrkDet1 status, p stopped at ArmB/TrkDet1 bool
stopPPSArmBTrkDet2 status, p stopped at ArmB/TrkDet1 bool
stopPPSArmFToF status, p stopped at ArmF/ToF bool
stopPPSArmBToF status, p stopped at ArmB/ToF bool
Mjj Leading jets invariant mass over all jets double
Mpf All Pf Candidates invariant mass double
Rjj Ratio Mjj/Mpf double
CandidatesMjj Leading jets invariant mass from CMS/PPS vertex associated double
FiducialCut 0 = No PPS Hit, 1 = PPS Hit bool
xiPPSArmFInfo xi information for Arm Forward, where vector index: 0 = signal and >1 for PU std::vector
xiPPSArmBInfo xi information for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
tPPSArmFInfo t information for Arm Forward, where vector index: 0 = signal and >1 for PU std::vector
tPPSArmBInfo t information for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmFDet1Info x axis position on Det1 for Arm Forward,where vector index 0 = signal and >1 for PU std::vector
yPPSArmFDet1Info y axis position on Det1 for Arm Forward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmFDet2Info x axis position on Det2 for Arm Forward,where vector index 0 = signal and >1 for PU std::vector
yPPSArmFDet2Info y axis position on Det2 for Arm Forward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmBDet1Info x axis position on Det1 for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
yPPSArmBDet1Info y axis position on Det1 for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmBDet2Info x axis position on Det2 for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
yPPSArmBDet2Info y axis position on Det2 for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmFToFInfo x axis position on ToF for Arm Forward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmFToFInfo y axis position on ToF for Arm Forward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmBToFInfo x axis position on ToF for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
xPPSArmBToFInfo y axis position on ToF for Arm Backward,where vector index 0 = signal and >1 for PU std::vector
deltaToF_00; Tof Combination Information double
deltaToF_0i; Tof Combination Information double
deltaToF_ii; Tof Combination Information double
stopPPSArmFToFInfo status, p stopped at ArmF/ToF std::vector
stopPPSArmBToFInfo status, p stopped at ArmB/ToF std::vector


Recipe to PPS Group Machine


  • From lxplus machine:
ssh uerjpc02
  • Common Area (/storage3, /storage4 or /storage5):
/storage3/PPSProject
  • Path to Submit Crab Jobs:
/storage3/PPSProject/MCProduction/CMSSW_6_2_0/src/CreateSampleSubmit

IMPORTANT: new users should request access to dmf@cernNOSPAMPLEASE.ch

Published Sample


Sample Type XS
/exhume_GG_13TeV_GEN/dmf-exhume_GG_13TeV_GEN-6ea6285d7af7bb71df1b704c0ccd271e/USER GEN-SIM X
/exhume_GG_13TeV_GEN/dmf-exhume_GG_13TeV_DIGI2RAW_20_May-86e7340def74f99d95fde32aa640f04e/USER DIGI-RAW X
/exhume_GG_13TeV_GEN/dmf-exhume_GG_13TeV_AODSIM_30_May-d39d75e84f73fb166fbe065ff6ec3525/USER GEN-AODSIM X
Sample at DBS phys03

CEP Algorithm Studies and Tests


Here, there is a talk to sum up how our algorithm works to match CMS/PPS vertex and to choose dijet candidate vertex from all the primary vertices.



Log Book

September 05th

Private Sample


  • ExHuMe (pp → gg → dijets) + Pythia 8 (minBias)
    • σ ~ 650.0 fb
    • 10000 events (50 PU) and 17800 (25 PU)
    • Generator cuts:
      • 300 < M < 2000 GeV
      • jet pT > 100 GeV

  • Pomwig (DPE/Inclusive) + Pythia 8 (minBias)
    • σ ~291.7 pb
    • 208224 events (50 PU) and 424408 (25 PU)
    • Generator cuts:
      • 0 < |t| < 4
      • 0.001 < xi < 0.2
      • 1 jet with pT > 80 GeV

  • Pomwig (SD dijets - Plus and Minus) + Pythia 8 (minBias)
    • σ ~2570.6 pb
    • 422711 events (50 PU) and 935831 (25 PU)
    • Generator cuts:
      • 0 < |t| < 4
      • 0.001 < xi < 0.2
      • 1 jet with pT > 80 GeV

  • QCD (non-diffractive jets) + Pythia 8 (minBias)
    • σ ~24.191 ub
    • 197943 events (50 PU) and 195300 (25 PU)
    • Generator cuts:
      • 50 < pThat < 3000 GeV

Selection Criteria


  • Select jets from the same vertex of the leading jet
    • pT (j1,j2) > 100 GeV/c
    • |η| (j1,j2) < 2.0
  • One leading proton in the CT-PPS fiducial region in both arms of the detector
    • Tracking station 1: -9.0 < y < 9.0 mm -23.15 < x < -3.15 mm
    • Tracking station 2: -9.0 < y < 9.0 mm -22.03 < x < -2.03 mm
  • No overlap hits in the same cell of the ToF detector
  • The arrival time difference at the CT-PPS depends on the z-vertex position and must be consistent with the vertex position of the central dijet system. A selection cut varying according to the z-vertex position (corresponding to the solid straight lines) keeps approximately 60% (50%) of signal events while reducing the inclusive dijet background by a factor 33 (18), for a 10 (30) ps timing resolution.
  • Dijet mass fraction, Rjj - 0.7 < Mjj / Mx < 1.15
  • Rapidity difference, (yjj - yX) < 0.1

PPS CEP plots and numbers


PU 50

All the numbers and plots bellow were obtained using the private samples! The numbers are normalized to 1 fb-1 pu50.png

PU 25

All the numbers and plots bellow were obtained using the private samples! The numbers are normalized to 1 fb-1 pu25.png

* Timing correlation plots (10 ps resolution in simulation) - Require all cuts in central detector, all combinations (signal+pileup) of 2 protons in acceptance of tracking/timing stations of both arms of PPS


-- DiegoFigueiredo - 16 Apr 2014

Topic attachments
I Attachment History Action Size Date Who Comment
Unknown file formatgz combinePlots_July17.tar.gz r1 manage 10016.0 K 2014-07-17 - 19:30 DilsonDeJesusDamiao  
PDFpdf deltaToF_00vsJet0.pdf r1 manage 26.0 K 2014-06-04 - 21:26 DilsonDeJesusDamiao  
PDFpdf deltaToF_0ivsJet0.pdf r1 manage 22.0 K 2014-06-04 - 21:26 DilsonDeJesusDamiao  
Unknown file formatgz dijetsPlots30_julho.tar.gz r1 manage 3429.9 K 2014-08-01 - 12:17 DilsonDeJesusDamiao  
Unknown file formatext minutes2105 r1 manage 1.4 K 2014-05-21 - 20:39 DilsonDeJesusDamiao  
Texttxt minutes2105.txt r1 manage 1.4 K 2014-05-21 - 20:42 DilsonDeJesusDamiao  
PNGpng pu25.png r1 manage 135.9 K 2014-09-11 - 16:53 DilsonDeJesusDamiao  
PNGpng pu50.png r1 manage 137.2 K 2014-09-11 - 16:53 DilsonDeJesusDamiao  
Compressed Zip archivezip sdpomwig_histos.zip r1 manage 980.3 K 2014-07-03 - 22:45 DilsonDeJesusDamiao  
PDFpdf vzcmsvsvzpps.pdf r1 manage 16.9 K 2014-05-07 - 15:54 ElizaMelo VertexZCMS x VertexZPPS
Edit | Attach | Watch | Print version | History: r63 < r62 < r61 < r60 < r59 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r63 - 2017-02-01 - DilsonDeJesusDamiao
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

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