# Presentation: 167708

## Inclusive Jets at the Tevatron

- Sally Seidel
- University of New Mexico
- ISMD’99
*for the CDF and D0      *
*Collaborations*

## 1.   Jets at CDF and D0

- 2.   The Inclusive Jet Cross Section      	versus ET at
- 3.   The Inclusive Jet Cross Section 	versus ET at
- 4.   The Ratio of the Dimensionless 	Cross Sections versus xT

## The motivation:

- Jet distributions at colliders can:
-  signal new particles + interactions
-  test QCD predictions
-  check parton distribution functions
*The complementary Tevatron detectors D0 and CDF have studied the inclusive jet cross section at*
* * 					....

## Common features of the two detectors:

- Luminosity monitors near beamline to monitor collision products.  At each site,
- 						     .
- CDF + D0 use slightly different        cross sections  the CDF jet cross section ends up 2.7% higher.
- Both calorimeters have projective towers that point to the nominal IP.
- The *p* beam is   ,  is measured relative     to   , + azimuth  is measured from the Tevatron plane.

## D0 in Run 1:

##  Central Tracking (radius to 75 cm, 270 cm 	length, ||  2), to find vertices:

- vertex drift chamber
- transition radiation detector
- 3-barrel central + 2 forward drift chambers
-  Calorimetry, to find jets:
- thick, uniform, finely segmented, unit gain LAr with U + Cu absorbers; ~compensating (Gaussian response, E/E = 5 – 7%)
- central calorimeter: ||  1.2
- 32 EM (U; 4 samples),                          16 fine HA (U; 3 samples),                    16 coarse HA (Cu, 1 sample) modules
- 7.2 absorption lengths @  = 0
-    = 0.1  0.1 (0.05  0.05 @ shower max)

-  Endcap calorimeter: ||  4.4
- inner EM (U), inner HA (U), middle HA (U + stainless), outer HA (stainless)
- 10 0, 8–9 segments @  = 2

## D0 jet energy calibration procedure corrects for:

- calorimeter hadronic response
- electromagnetic showering out of the reconstruction cone (see below)
- offsets due to instrumental effects, multiple event pile-up, multiple interactions, underlying event
- EM response calibrated with Zee and        0  
- HA response calibrated from central -jet events by ET balancing
- Calibration of jets
- w/ ET < 150 GeV: directly from events 					w/central , jet
- w/ ET < 300 GeV: with central  + high- 				jet
- Extend to ET > 300 GeV by simulation

## CDF in Run 1:

## Central tracking (||  1):

-  4-layer Si strip vertex chamber for *r-***
-  vertex drift chamber for *r-z* + IP
-  3.2 m central drift chamber                    @ 31  *r*  132 cm
-  all inside 1.5 T *B* field
-  Calorimeters :
- Central Calorimeter
-  EM (scintillator + Pb, ||  1.1);    HA (scintillator + Fe, ||  1.3)
-     = 15°  0.1
-  4.5 absorption lengths

- Plug + forward calorimeters (2.2    4.2)
- gas proportional chambers
-     0.1  5°

## CDF jet energy calibration corrects for:

- calorimeter response
- multiple interactions
- neutrinos + muons
- underlying event
- Corrections are based on
- Monte Carlo tuned on test beam (for          10 GeV  Ebeam  227 GeV).
- CTC (drift chamber) tracks w/                  400 MeV/c  pT  10 GeV/c.
- MC fragmentation parameters tuned to match observed particle momenta + numbers.
- Calorimeter resolution:  = 0.1ET + 1GeV.
- Forward + plug calorimeter response calibrated by balancing dijets in which 1 jet of the pair is central.

## The data: 

- CDF + D0 reconstruct jets using an *iterative cone algorithm* with cone radius

- A typical jet crosses 40 calorimeter cells in CDF + 100 cells in D0.
- Both experiments use algorithms based on the Snowmass one, which says: if *i *indexes all towers within cone radius *R,*

## The D0 iterative cone algorithm:

- Examine all calorimeter towers with             ET > 1 GeV.
- Begin w/ highest ET tower, form preclusters by adding all neighboring towers within          *R* = 0.3.
- Calculate jet direction (,) for each precluster.
- Sum all energy within *R *= 0.7 about precluster jet direction; recalculate (,).
- Iterate until jet direction is stable.
- If two jets overlap –
- Merge them if  50% of the smaller jet’s ET is in the overlap region; otherwise assign the overlap energy to the nearest jet.
-  Then recalculate (,).

## The CDF iterative cone algorithm:

- Examine all towers with ET > 1 GeV.
- Form preclusters from continuous groups of towers with monotonically decreasing ET.
- If a tower is outside a window of 7 x 7 towers from the seed of its cluster, start a new precluster with it.
- For each precluster, find the ET-weighted centroid for all towers with ET > 100 MeV within *R* = 0.7.
- Define the centroid to be the new cluster axis.
- Iterate until the tower list is stable.
- If 2 jets overlap –
- Merge them if  75% of the smaller jet’s ET is in the overlap region; otherwise assign overlap energy to the nearest jet.
- Recalculate (,).

## For D0, 

- For CDF,                               , where
- ....the D0 jets are massless (pT = ET) whereas the CDF jets have mass.

## The Inclusive Jet Cross Section, E·d3/dp3

- For jet transverse energies achievable  at the Tevatron, this probes distances down to 10-17 cm.
- For massless jets + 2 acceptance in ,
- 					       
- This is what we measure.

## C

- CDF (D0) begin with similar data quality requirements:
- |zvertex | < 60 (50) cm to maintain projective geometry of calorimeter towers.
- 0.1 (0.0)  |detector|  0.7 (0.5) for full containment of energy in central barrel.
- (CDF only:) Etotal < 1800 GeV to reject accelerator loss events.
- To reject cosmic rays + misvertexed events, define     = missing *E**T*.  Require
- 						    (CDF)
-         < (30 GeV or 0.3*E**T*leading jet,
-   		  whichever is larger).	    (D0)
- Apply EM/HA + jet shape cuts to reject noise fakes.

## Both next correct for

- Pre-scaling of triggers.
- Detection efficiencies (typically 94 –100%).
- underlying event + multiple interactions.
- “Smearing”: energy mismeasurement + detector resolution.
*No correction is made for **jet energy deposited outside the cone** by the fragmentation process, as this is **included in the NLO calculation** to which the data are ultimately compared.*

## The CDF Unsmearing Procedure:

*Simultaneous correction for detector response + resolution produces a result that is indepen-dent of binning but preserves the statistical uncertainty on the measured cross section.*
- A smooth function is proposed representing a trial “pre-detector” cross section.  This is called the “physics function.”  For the inclusive cross section, the physics function has the form:

## The physics function is convoluted with measured ET-loss and resolution functions and binned in ET. The response functions have long tails due to uninstrumented regions.  The response function is parameterized by an exponential combined with a Gaussian and depends on the choice of cone opening angle.   The response function requires

- Gaussian mean M = C1ET3 + C2ET + C3ET2 + C4ET + C5ET + C6,
- Gaussian standard deviation  = C7(ET + C8)    - C10, and
- Exponential decay constant S = C11ET + C12  (where C11 and C12 depend upon the ET range).

## The Ci are determined by comparing measured calorimeter jet response to jets in Monte Carlo simulation.  The Monte Carlo is based on single particle response from test beam data and isolated tracks from minimum bias events.

- The result of the convolution is compared with the measured cross section.
- The parameters of the physics function are iterated to obtain the best match between the convoluted function and the data.
- The ratio of (data):(physics function) for each ET bin is tabulated.  The correction factors (typically 10%) are subsequently applied to data to obtain the “unsmeared physics cross section.”

## The D0 Unsmearing Procedure:

*This procedure corrects for detector response and resolution sequentially*.
- Correct the *E**T* of each measured jet for the average response of the calorimeter.
- Assume an ansatz function of the form
- Convolute this with the (ET-dependent, Gaussian) measured resolution function.  Compare the result with the measured cross section.
- Iterate to optimize *A*, *B*, and *C*.
- Tabulate ratio (initial physics function) :  (smeared physics function) for each *E**T* bin + apply it to the cross section.

## The inclusive jet cross section results, for

## The CDF Result:

**(88.8 pb****-1****)**
**(20.0 pb****-1****)**

*The CDF Result:*

## Slide 26

## Slide 27

## Slide 28

## Systematic uncertainties (all uncorrelated)

- on the CDF inclusive jet cross section:
- i. Calorimeter response to high-pT charged hadrons
- ii. Calorimeter response to low-pT charged hadrons
- iii. Energy scale stability (±1%)
- iv. Jet fragmentation model used in the simulation
- v. Energy of the underlying event in the jet cone (±30%)
- vi. Calorimeter response to electrons + photons
- vii. Modelling of the jet energy resolution function
- viii. Normalization (3.8%)

## Slide 30

## The high-ET excess present in the Run 1b data is consistent with what was previously observed in the Run 1a data.  No single experimental source of systematic uncertainty can account for the excess. The significance of the excess was analyzed after Run 1a with 4 normalization-independent, shape-dependent statistical tests:

- signed Kolmogorov-Smirnov
- unsigned Kolmogorov-Smirnov
- Smirnov-Cramer-VonMises
- Anderson-Darling
- The MRSD0' PDF was used since it describes the low-ET data best. All 8 sources of systematic error + the effect of finite binning was included.  Two ET ranges were tested independently:

## For 40 < ET < 160 GeV:

- 	*All 4 tests show > 80% agreement **	**	**between data + theory.*
- For ET > 160 GeV:
- 	*Each test indicates that the probability **	**that the excess is due to fluctuation is **	**	**1% .*
- For other PDF’s (CTEQ2M, CTEQ2ML, GRV94, MRSA', + MRSG), agreement at low ET is reduced if the agreement at high ET is improved.
- The best agreement at high ET is for CTEQ2M: 8% likelihood for high ET but only 23% for low ET.
- Quantitative comparison of the CDF Run 1b result with theory is now underway.

## The D0 result:

*The D0 result:*

## Slide 34

## Systematic uncertainties on the D0 inclusive jet cross section:

- i. The calorimeter energy scale (including response + corrections).
- ii. The jet selection procedure (including jet,      , and vertex cut efficiencies).
- iii. Relative luminosity (uncertainties on prescale values of component triggers).
- iv. Determination of the jet energy resolution function needed for unsmearing (including fitting + choice of q-g fractions).
- v. Luminosity (including hodoscope acceptance, Level 0 hardware efficiency, and                  uncertainty).

## D0 systematic uncertainties

- 		 @

## 2 Comparison of the D0 data and theory:

- i.) Define
- ii.) Construct the C*ij* by analyzing the correlation of uncertainties between each pair of ET values. (Bin-to-bin correlations for representative ET bins are ~ 40% + positive.)
- iii.) There are 24 d.o.f.

## iv.)  Compare the data to JETRAD with 5 PDF’s:

- iv.)  Compare the data to JETRAD with 5 PDF’s:
- 2/dof = 0.65 – 1.0 (Prob. 47 – 90%) for data in range ||  0.5.
- 2/dof = 0.8 – 1.0 (Prob. 24 – 72%) for data in range 0.1  ||  0.7.
- v.)  Compare the data to EKS using CTEQ3M, *R**sep*= 1.3*R, *and * *scale  = cETmax or  = cETjet     (c = 0.25, 0.5, or 1.0):
- Prob.  57% in all cases.
- 		Good agreement.

## Comparison between D0 + CDF measurements, both relative to JETRAD: 

*Comparison between D0 + CDF measurements, **both relative to JETRAD:*

## Excellent agreement between nominal values for ET  350 GeV.

- To quantify the level of agreement over the full ET range, D0 carried out a 2 comparison between D0 data + the nominal curve describing the central values of the CDF Run 1b data:
- 2 = 41.5 for 24 d.o.f.
- To approximate a  “statistical error only” comparison, find the value of the CDF curve at the D0 ET points, multiply the D0 stat errors by        to make them equivalent to CDF’s, + remove the 2.7% relative normalization difference...2 = 35.1 (prob. 5.4%).
- Add syst. error info by expanding the covariance matrix to include both D0 + CDF...2 = 13.1 (prob. 96%).

## The Ratio of Dimensionless Cross Sections

- Consider the dimensionless cross section,
- Again for massless partons, this can be written as

- 

## The Naïve Parton Model hypothesis that d is independent of      ,  when plotted versus xT, is called “scaling.”

- QCD predicts scaling violation due to the energy-scale dependence of the probability for gluon radiation from a primary parton in the collision.
-  This energy scale dependence leads to the running of s + the evolution of the PDF’s.
-  Comparison of d measured at 2 different energies by the same experiment suppresses many theoretical + experimental uncertainties.

## Suppression of theoretical uncertainties:

- was calculated by JETRAD, for the case of   = 0.5ETmax and CTEQ3M, then compared to a range of other options for , PDF, + s:

## the result of this study on the dependence of the calculation on input parameters:

*the result of this study on the dependence of the calculation on input parameters:*

- All show  10% variation for *x**T* < 0.4.

## D0 + CDF each collected 600 nb-1 of data @        	         in Dec. ‘95.

- Results for a previous data set of 8.6 nb-1 @         		   had been published by CDF in ‘93.
- These data were analyzed like the 1800 GeV sample + compared to JETRAD w/ MRSA +    = 0.5ETjet.
-  Principal difference between the treatment of the 630 GeV + 1800 GeV data sets is the underlying event correction, which increases by 50% as       increases from 630 to 1800 GeV.

## The D0 cross section result, for 

- The D0 cross section result, for

## The CDF Result:

## Slide 48

## Slide 49

## Systematic Errors on the CDF 630 GeV measurement: 

- Systematic Errors on the CDF 630 GeV measurement:

## Comparison of CDF + D0 measurements of the  cross section @ 630 GeV:

## Result + Preliminary Conclusions:

-  CDF + D0 measurements of the cross section @ 630 GeV agree with each other above ~ 80 GeV.
- Conclusions about agreement between the 2 experiments below 80 GeV cannot be drawn yet; pending further studies of systematic errors, the data may diverge in the lowest few ET bins.
-  Both measurements are consistent w/theory above ~ 80 GeV; theory somewhat higher below.
*...**additional studies are needed for energy scale determination @ low E**T**.*

## The cross section ratio,

- for D0 + CDF, is then formed + compared to
- JETRAD w/ 3 pdf’s, and
- EKS w/ 2 pdf’s:

## The D0 result for the cross section ratio:

- 2       _Prob._
- 50.8  0.02%
- 51.7  0.01%

- 2    _Prob._
- 29.9  7.19%
- 41.8  2.91%
- 50.8  0.02%
- 30.9  5.66%

## Uncertainties on the D0 cross section ratio:

*Uncertainties on the D0 cross section ratio:*

## D0 has considered the case in which the scale  is      -dependent.  Using CTEQ3M:

- D0 has considered the case in which the scale  is      -dependent.  Using CTEQ3M:

-   __2         _Prob._
- 12.1      91.3%
- 28.5        9.8%
- 10.9      94.7%

- The apparent preference for different scales may indicate that the NNLO terms, when calculated, will not be negligible.

## The CDF measurement of the cross section ratio, compared to EKS: 

- The CDF measurement of the cross section ratio, compared to EKS:

## CDF Systematics:

## Comparison of CDF + D0 Results:

## The ratio measurements are consistent for      xT  0.1.

- Discrepancy below xT = 0.1 may be traced to the d(630 GeV) results + was also seen in the 546 GeV data.
- Variation among theoretical predictions is small unless  is      -dependent.
- The theoretical uncertainty on the ratio is 1/2 what it is on the inclusive cross section.