No Slide Title
Download
Report
Transcript No Slide Title
D
- Calibration for jets
Reminder on the DØ detector
Jet Identification and Reconstruction
Jet Energy Scale:
results from Run 1
b-jet calibration
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
1
D
Run II Detector
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
2
D
- Calorimeter
Fine segmentation:
semi-projective towers in 0.1
x0.1
4 em layers: 2, 2, 7, 10 X0
shower-max (EM3): 0.05 x 0.05
4/5 hadronic (FH + CH)
hermetic with full coverage
|| < 4.2 ( 2o)
int > 7.2 (total)
Uranium absorber (Cu (CC) or
Steel (EC) for coarse hadronic)
compensating e/ 1
to be studied with shorter
shaping
Beam Tests of the D0 Uranium Liquid Argon Calorimeter.
NIM A324, 53 (1993)
NIM A 338 185 (1994)
ATLAS Ringberg workshop - 23/7/2002
from test beam
data
e: sE/E = 15% /E + 0.3%
: sE/E = 45% /E + 4%
Ursula Bassler, LPNHE Paris
3
D
SCA non-linearity
functional form of SCA non-linearity
correction function
correction important at low energies
electronic noise translates into
higher energy
jet become more narrow
250 MeV
0.25 ADC
count/MeV
0
~Energy/GeV
for energies
> 200MeV
non-linearity
introduces an
offset of
~250 MeV for
the gain 8
4 measurements
ATLAS Ringberg workshop - 23/7/2002
1 GeV
Gain 8
Gain 1
0
10
20
30
40
50 GeV
Ursula Bassler, LPNHE Paris
4
D
Central Jet Triggers
L2 jet
Efficiency vs jet pT
CJT(1,3)
CJT(1,5)
CJT(1,7)
CJT(1,10)
Cluster 3x3 or 5x5 trigger
towers around L1 seed
towers
L3 jet
Simple cone or tower NN
algo’s 0.1x0.1 towers
3 single jet (tower) triggers:
JT_LO L1: 5 GeV,
L3:10 GeV
JT_HI
L1:10 GeV,
L3:15 GeV
CJT40: L1:40 GeV
Efficiency
standard jet selection,
offline pT > 8 GeV
very sharp turn on
L1 Trigger efficiency CJT(1,x)
L1 Trigger efficiency CJT(2,x)
L1 single jet efficiencies
ask for one or two hadronic trigger towers
(0.2x0.2) above threshold
use -trigger as unbiased reference to
measure turn-on
ask for one and only one reconstructed jet
in ||<0.7
L1 hadronic response about 40% low for
current data set
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
5
D
NADA: noise reduction
NADA = New Anomalous Deposit Algorithm
identify isolated energy deposits in the calorimeter = “Hot Cells”
Source: electronics, Ur noise, beam splash, cosmics etc
Improve object resolution and ETmiss
Run 1: AIDA
Only examine neighbors in the same tower for Ecell > 10 GeV
99% efficient, BUT 5-10% misidentification rate
examine all cells with > 1 GeV
remove cells < -1 GeV & > 500 GeV
ET < 5 GeV removed if no neighbor with E
> 100 MeV
ET < 500 GeV removed if no neighbor with
E > 2% Ecell
high efficiency (90%) and low
misidentification
ET > 1 GeV : ~0.5%
ET > 10 GeV : ~0%
on average about 0.8 cells / event
ATLAS Ringberg workshop - 23/7/2002
ETthresold
ETneighbour> 100 MeV or 0.02Ecell
Ursula Bassler, LPNHE Paris
6
D
Jet Finding
Calorimeter
jet (cone)
jet is a collection of energy deposits with a
given cone R: R Δ 2 Δη2
cone direction maximizes the total ET of the jet
various clustering algorithms
correct for finite energy resolution
subtract underlying event
add out of cone energy
Particle jet
a spread of particles running roughly in the
same direction as the parton after hadronization
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
7
D
Jet Algorithms: Cone
Run 1 Legacy Cone
draw a cone of fixed size around a seed
compute jet axis from ET-weighted mean and jet ET from ET’s
draw a new cone around the new jet axis and recalculate axis
and new ET
iterate until stable
algorithm is sensitive to soft radiation (split & merge)
Improved Run 2 cone
use 4-vectors instead of ET
add additional midpoint seeds between pairs of close jets
split/merge after stable proto-jets found
algorithm is infrared safe
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
8
D
Jet Algorithms: kT
For each object and pair of objects:
dii k 2T,i
order all dii and dij:
If dmin=dij
2
ΔR
ij
dij min(k 2T,i ,k 2T, j ) 2
D
Soft
Collinear
merge particles
(if R<<1 )
Resolution
parameter
(D=1)
theoretically favored, no split-merge
to reduce computation time, start
with 0.2 x 0.2 pre-clusters
x-section measurement differ from
cone-jet
If dmin=dii
jet
(JETRAD)
DØ
Subjet multiplicity of gluon and quark jets reconstructed using the kT algorithm in pbarp collisions
Phys. Rev. D65 052008 (2002) hep-ex/0108054
The inclusive jet cross section in pbarp collisions at sqrt(s)=1.8 TeV using the kT algorithm
Phys. Lett. B {525}, 211 (2002) hep-ex/0109041
Ursula Bassler, LPNHE Paris
ATLAS Ringberg workshop - 23/7/2002
9
D
Hadronization effects?
• particle jets are more (less) energetic
than parton jets with kT (cone)
kT collects more energy
cone looses energy
kT jets are 7 (3)% more energetic at
60 (200) GeV than cone jets:
• consistent with HERWIG at high pT,
at 2s at low pT
applying correction to
cone-jets improves
agreement between
the 2 algorithms
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
10
D
Jet Algorithms: CellNN & Flow
Cell Nearest Neighbor
layer-by-layer clustering starting with EM3
each local maximum starts a layer-cluster then add in neighbors
energy sharing according to transverse shape parameterization
angular matching of floor clusters
search for minima in longitudinal energy distribution to separate EM
and hadronic showers
Energy Flow algorithm
use tracking information to better characterize the contributions from
charged particles
in development
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
11
D
Jet Selection
DØ Run 2 Preliminary
central jets (Run 2 cone, R=0.7)
event quality cuts
number of jets 1
Etotal in the calorimeter 2 TeV
missing ET 70% of leading jet pT
Zvtx < 50 cm
CHF
EMF
leading Jet Cuts
Jet pT > 8 GeV (offline cut)
0.05 EMF 0.95
CHF 0.4 (0.25 tight)
HotF 10 (5 tight)
(HotF = ET1st cell / ET2nd cell )
n90 > 1
(number of towers that
contain 90% of jet ET)
efficiencies from MC
loose: ~100%
~Flat in
tight: ~ 98%
HotF
n90
Data
— MC
Non-linearity of SCA included in MC
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
12
D
Jet Energy Scale
jet
correct Jet Energy to the particle level
ptcl
E jet
E meas
jet E offset
shower
R calo
R
jet
jet
Eoffset energy offset from underlying event,
pile-up, noise
determined from Min. Bias Events
Rcalo calorimeter response
using -jet events: Missing ET
Projection Fraction Method
Rshower energy contained in jet
corrections from MC - energy in
cones around the jet axis
depending on jet algorithm!
Determination of the Absolute Jet Energy Scale in the D0
Calorimeters. NIM A424, 352 (1999), hep-ex/9805009
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
13
D
Run I: Offset corrections
subtract contributions not associated to the high pt interaction:
Ur noise, pile-up, multiple interaction, underlying even
measured as ET densities D, to be multiplied by the area of a jet in
E
EOFF= EUE +NZB EUE +Enoise+Epile-up
measurement of the ET density D in zero bias event
measurement of DUE from minimum bias events
DUE=DMB-DZBno HC
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
14
D
Run I: Offset corrections
Ur noise, pile-up, multiple events
underlying event contribution
ICR
measured for different luminosities
depends on s and process
dominant error from occupancy
dependence
associated to a single event
independent of luminosity
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
15
D
Run I: response correction
using -jet events
ideal calorimeter :
ET ETjet 0
jet response (with calibrated ): R jet 1
E
T nˆT
ET
Ejetmeas:
dependent on
energy response and
resolution, threshold effects
and smearing
better: E= ET cosh jet
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
16
D
Run I: jet response
comparison
of jet
response in different
cryostat regions
CC ||<0.7
ICR 0.7<||<1.8
EC 1.8 <||<2.5
effect of finite jet
resolution at E = 10GeV
lowest response in
ICR: int < 6
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
17
D
Run I: EC/CC correction
independent of E as EC/CC
similar in construction
derived from overlap region
of CC and EC response at
60<E<180 GeV
Fncry/Fscry=0.997 0.003
EC response 2% below CC
compared to the ratio of a
fit to the 2 response
functions
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
18
D
Run I: ICR correction
inhomogeneous detector material: correction as function of and ET
high ET: jet-jet events
low ET: -jet events
or leading jet required to be central (| |< 0.5)
fit of response as Rjet = + b ln ET+ b ln (cosh )
correction derived from difference between measurement and the expectation
for an ideal detector, extrapolated from fit at ||<0.5 and 2 < | |<2.5
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
19
D
Run I: low ET bias
Etjet
> 8 GeV
jet resolution ~50%
migration of low ET jets
jets fluctuating below Etmin
are not reconstructed
bias of the response
towards higher values
• as response determined
from Etmiss and , bias
correction determined from:
Rbias
ATLAS Ringberg workshop - 23/7/2002
R jet(1)
R jet(0)
Ursula Bassler, LPNHE Paris
20
D
Run I: Response function
fit of the measured
response function
Rjet(E)=a+b ln E+c (ln E)2
logarithmic terms
justified by non
compensation at low E
fit of CC and EC
measurement for ET>30
GeV
at highest energy
prediction from MC after
tuning response on data
in measured region
error band derived
taking into account
correlations
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
21
D
Run I: Showering correction
corrects for out-of-
1%
cone energy
belonging to the jet
scales reconstructed
jet to particle level:
S=Ejet/(Ejet+EshoMC)
4%
10%
parameterizations for
different cone sizes
errors at low energy:
offset subtraction; at
high energy: stat
2.5%
5%
10%
ATLAS Ringberg workshop - 23/7/2002
shower correction
depend on jet
profiles, but not on s
Ursula Bassler, LPNHE Paris
22
D
Run II: +jet / Z+jet
/Z+jet
QCD (udsg)
signal:152 M evts
bkgd: 47.1 M evts
signal: 64.8k evts
bkgd: 650 evts
QCD (cbt)
W+jet, Z+X,
• +jet: Run I method – jet calibration possible up to 250 GeV
• Z+jet: lower statistics, but clean sample, useful at low energies, x-check!
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
23
D
b-jet calibration
naïve reconstruction
of Z-mass shows a
lower mass for selected
b-jets than light quark
jets.
energy losses from
semi-leptonic b decays
(, )
wider b-jets (due to
the large b-mass)
ATLAS Ringberg workshop - 23/7/2002
Z bb
peak: 82.6
Z qq
peak: 86.8
Ursula Bassler, LPNHE Paris
24
D
Z bb vs + b-jet
+ b-jet :
Zbb:
high statistics, allows for a tight b-jet selection (btagging).
expected number of tagged events: 1.2 M
but: sensitive fractional imbalance I= (pT() - ET(jet))/ pT()
systematics closer to physics processes (H or
Top) at high pT
resonance mass independent of
multiple interactions.
but: signal/noise~10-3 requires special trigger
(Silicon Track Trigger – operational end 2002)
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
25
D
CDF Run 1: Z bb selection
about 120 000 Z bb events produced in Run 1
expected to be observed ~ 50-100
Trigger: central muon (pT> 7.5 GeV) 5.5 M evts
Offline: request 2 tagged (0.7 cone) jets 5479 evts
QCD background rejection based on event topology:
> Z is produced by a time-like q-qbar anihilation,
> QCD produced color flow between initial and final
partons
> Z is expected to have soft radiation between the jets
> background will also have strong radiation between IS
and FS partons.
http://www-cdf.fnal.gov/physics/ewk/zbb_new.html
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
26
D
CDF Run 1: 3 ET and 12
Use 2 kinematic variables to discriminate:
3 ET : sum of ET of the clusters outside the 2
leading jets
12 : azimuthal angle difference between the 2 jets
cuts derived: 3 ET < 10 GeV, 12>3 rad
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
27
D
CDF Run 1: Z bb Signal
after cuts:
S/N=1/6 at the Z
mass peak
select/antiselect
w.r.t. the 2
variables to
determine the
tagging probability
3.2 s exces
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
28
D
CDF Run 1: Likelihood fit
Results:
MZ=90.0 2.4 GeV
sZ = 9.4 3.5 GeV
NZ=91 30(stat) 19(sys.)
Pythia: expect 12414
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
29
D
First Run 2 QCD Physics
Inclusive jet pT spectrum at 1.96 TeV
Dijet mass spectrum at 1.96 TeV
Ldt = 1.9 ± 0.2 pb-1
Ldt = 1.9 ± 0.2 pb-1
Only statistical errors
Only statistical errors
Highest 3-jet event
ETjet1 : 310 GeV
Etjet2 : 240 GeV
ETjet3 : 110 GeV
Etmiss : 8 GeV
not fully corrected distributions:
preliminary correction for jet energy scale
(but no unsmearing or resolution effects)
30-50% systematic error in cross-section
no trigger selection efficiency corrections
ATLAS Ringberg workshop - 23/7/2002
Ursula Bassler, LPNHE Paris
30