20121009_SoftQcdAlfa_RapidityGaps
Download
Report
Transcript 20121009_SoftQcdAlfa_RapidityGaps
Diffractive analyses with gaps
Hardeep Bansil, Oldrich Kepka, Vlastimil Kus, Paul
Newman, Marek Tasevsky
Workshop on Diffractive Analyses with ALFA
09/10/2012
Contents
Measuring rapidity gaps in ATLAS
Soft diffraction
Diffractive dijets
How ALFA can help
2
pp Cross Section & Inelastic Interactions
Non Diffractive Events
Coloured exchange, Soft PT spectrum
High multiplicity final states peaking at
central rapidity
Largest cross section at LHC
Diffractive Events
Colour singlet exchange (pomeron) results in
a rapidity gap devoid of soft QCD radiation
Can involve Single or Double proton
dissociation
Size of the rapidity gap is related to the
invariant mass of the dissociated system(s)
IP
GAP
Related to pz loss of intact proton in Single
diffractive case
25-30% of the total inelastic cross section
(ξX > 5×10-6) is measured to be inelastic
diffractive
IP
GAP
3
Measuring rapidity gaps in ATLAS
Use the full tracking (|η|<2.5) and calorimetric range
(|η|<4.9) of detector
In the calorimeters electronic noise is the primary
concern
The standard ATLAS energy deposits are from Topological
clustering of cells
Seed cell required to have an energy significance σ = E/σNoise > 4
Statistically, expect 6 topological clusters per event from
noise fluctuations alone
187,616 cells multiplied by P(σ≥4) ≈ 6
Just one noise cluster can kill a gap
Additional noise suppression is employed but set the
thresholds as low as the detector will allow
Apply a statistical noise cut to the leading cell in the
cluster which comes from the LAr systems (noise from the
hadronic Tile calorimeter follows a double Gaussian)
Set Pnoise within a 0.1 η slice to be 1.4x10-4
N is the number of cells in the slice
The threshold Sth(η) varies from 5.8 σ at η = 0 to 4.8 σ
at η = 4.9
4
Gap Finding Algorithm
Detector split into bins of η
Detector Level Bin full if it
contains
Example Single
Diffractive
ΔηF:3.4 |ηStart|:4.9
one or more noise suppressed
calorimeter clusters above ET cut
of 200 MeV
- AND/OR one or more tracks reconstructed
above pT cut of 200 MeV
Generator Level Bin full if it
contains
Example NonDiffractive
ΔηF:0.4 |ηStart|:4.9
Minimum Bias Trigger Scintillators
(Physics Trigger)
one or more stable (cτ > 10 mm)
generator particles above pT cut
of 200 MeV
ΔηF = Largest region of
pseudo-rapidity from
detector edge containing no
particles within bins
Forward Rapidity Gap
Devoid of particles pT > 200 MeV
η = -4.9 to η = 0.5
ΔηF = 5.4, ξ = 1x10-4, MX = 75 GeV
5
Soft Diffractive Analysis
EPJC 72 (2012) 1926
Utilising the first stable beam physics run at √s=7 TeV
ATLAS accumulated 422,776 minimum bias events
7.1 μb-1 (at peak instantaneous luminosity 1.1x1027 cm-2s-1)
Trigger requirement as loose as possible. Require one hit in the
MBTS online, offline we required two hits with MC thresholds matched
to the efficiency observed in data
Use unfolded data up to a
forward gap size of ΔηF = 8
Exponential
Fall
Raw ΔηF plot for data and
MC at the detector level,
including trigger requirement
on MC and data
Poor Trigger
Efficiency
Diffractive Plateau
Event normalised
6
Soft Diffractive Analysis
The Raw gap size distribution is unfolded to remove
detector effects
Tune the ratios in the MCs from Tevatron
data
Data is corrected for trigger inefficiency at
large gap size
Single application of D’Agostini’s Bayesian
unfolding method
MC normalised to Default ND, DD and SD Cross
section up to ΔηF = 8
Integrated cross section in diffractive plateau:
5 < ΔηF < 8 (Approx: -5.1 < log10(ξX) < -3.1)
= 3.05 ± 0.23 mb
~4% of σInelas (From TOTEM)
Corrected ΔηF Distribution
ξ=10-5.1
ξ=10-2.5
ΔηF vs. Pythia 8
Compare to Pythia 8 4C split into sub-components
Non-Diffractive contribution dominant up to
gap size of 2, negligible for gaps larger than 3
Shape OK, overestimation of cross section in diffractive
plateau
Large Double Diffraction contribution across
entire forward gap range (large uncertainty)
7
Diffractive dijets
Work in progress
Search for single diffractive events diffraction with a hard scale set by 2 jets
Described by diffractive PDFs + pQCD cross-sections
Previous measurements of hard diffractive processes at HERA and Tevatron
Now also studied at CMS
Measure the ratio of the single diffractive to inclusive dijet events
Understand the structure of the diffractive exchange by comparison with
predictions from electron-proton data and be able to get a measure of FDjj
Gap Survival Probability – the chance of the gap between the intact proton and
diffractive system being lost due to scattering (affects measured structure function)
Tevatron have Gap Survival Probability of 0.1 relative to H1 predictions
Predict LHC to have GSP of ~ 0.03 – 0.07
Rescatter with p?
Comparison of Tevatron
diffractive PDF to H1
expectations in terms of
momentum fraction of
parton in Pomeron
Gap destruction by
secondary scattering
ξ
8
Analysis
2 medium Anti-kt jets with R=0.4 or R=0.6:
ET Jet1,2 |η| < 4.4, ET Jet1 > 30 GeV, ET Jet2 > 20 GeV
Cut values based on 2010 SM dijet analysis / JES systematic
Ask for a forward gap: |ηstart| = 4.9, ΔηF ≥ 3.0
Currently employing two separate strategies for analysis
2010 period A-B (∫L dt ≈ 7 nb-1) triggering on L1_J5 or L1_FJ5
2010 period A-F (∫L dt ≈ 3 pb-1) with pT-dependent L1 jet & forward jet triggers
Trigger below 100% efficiency plateau to collect more events than 2010 Standard
Model inclusive dijet analysis
Use POMWIG/HERWIG++ for Single Diffractive Dijets
No direct DD samples, DPE samples contribute little
Use PYTHIA 8/HERWIG++ samples as inclusive (ND) Dijets
Special MC request in preparation, filtered on gap size
Reconstruct ξ and zIP using E±pz method based on fwd gap side IP (ξ)
( E i pzi ) X
s
~
( E i p zi ) jets
~
z IP
( E i p zi ) X
9
Hard Diffraction at Generator Level
Demonstration of difference in ND and SD models for dijet events
Latest gap survival probability estimate 6% included in SD model
pTjet > 20 GeV
Hard dijets → bigger MX → smaller gaps
Like soft diffraction, have to go to bigger gaps in order to separate SD
from ND
Contrary to soft diffraction, we no longer observe diffractive plateau
10
Forward Gap Size Distributions
Different data ranges will allow for cross checks
Both provide significant statistics for forward gap sizes > 3.0
Differential cross section as a function of forward gap size for data vs.
MC models (scaled to first bin of data)
Without this scaling, difficult to distinguish Single Diffractive signal from Non-diffractive
Currently available ND samples statistics insufficient at larger gap sizes so new samples with
gap-based filter necessary
Diff. csx in ΔηF, (no noise supp.), Period B v F
Diff. csx in ΔηF (no noise supp.) – Comp. to MC
11
ξ and zIP
Aim to measure cross section as function of both of these variables as well as forward gap size
With analysis cuts applied, get good correlations between truth and reconstruction levels
Relies on picking out the correct value of E+pz or E-pz in calculations based on side gap is on
MC- Majority of time we identify side with gap correctly at truth and reconstruction levels (that
intact proton would have been on)
Data – no comparison to truth so ALFA can help by tagging proton indicating which
size gap should be on
Truth v recon ξ using Pomwig
Truth v recon zIP using Herwig++
12
What ALFA can do for us
Rapidity gaps in ATLAS data are a sensitive probe to the
dynamics of soft and hard diffractive proton dissociation
The data can be used to investigate and tune the current MC
models
ALFA can help:
Constrain relative amounts of non, single and
double diffraction by removing ND background
More precise measurements of cross sections
Study properties of soft diffractive events (multiplicities, UE)
Study properties of hard diffractive dijet events (survival
probability, jet shapes, ratios of SD to inclusive dijets)
13
Inclusive dijet s : results – full
2010
Official 2010 results
Our reproduction
We reproduced official inclusive SM2010 dijet cross section measurement
→
Excellent agreement!
However, this trigger strategy not efficient in collecting events on tail of gap spectrum
→ designing of new trigger scheme
Gap-size distribution
Different trigger strategies
ATL-COM-PHYS-2011-738
2010 inclusive dijet x-sec measurment
Big OR of jet triggers:
J20 || J30 || J35 || J50 || J75 || FJ30 || FJ50
No triggers asked for
7700 events (no triggers) vs. 40 (SM2010 triggers)
=> a lot of room for trigger strategy improvement
Large gaps require collecting events with small pT
→ we have to move below 99% trigger efficiency plateau