Download - Photon energy smearing from events
Glen CowanRHUL Physics
1
Photon energy smearing from events
Glen Cowan, Sudan Paramesvaran, David Hopkins, Henning Flaecher
Royal Holloway, University of London
EMC Calibration Meeting 6 Sept 2006
Glen CowanRHUL Physics
2
Context
Follows on from resolution studies by David, Henning, Sudan.
Same n-tuples and selection as before, see e.g. David Hopkins talk 5 April 06 at EMC software meeting:
2 ‘good’ measured tracks 1 ‘good’ measured photon 1 identified (loose) Kinematic fit 2 prob > 0.05
184 fb→ 1.2 million events
Glen CowanRHUL Physics
3
Goal of study
For the photon in the events we have
Emeas measured by the calorimeter
Efit from the kinematic fit
Goal: (Roodman, Kocian, ...) try smearing MC photon energies to give better data/MC agreement.
Histograms of x = Emeas/Efit
found to have high-x tail, especially for higher E
Glen CowanRHUL Physics
4
Method (1)
Find Try to smear MC so that it looks like data, i.e.,
find a pdf s(z) such that y + z ~ f(x)
Scale to samearea as data histogram
For x = y + z
Glen CowanRHUL Physics
5
Method (2)
Find
For x = y + z we have
In terms of the histograms this is where
For a parameterized pdf s(z;) we therefore have
original MC
smearing matrix
smeared MC
Glen CowanRHUL Physics
6
Method (3)
Try for smearing pdf: Gaussian
Student’s t
Use binned ML but for now ignore MC statistical errors,equivalent to minimizing
For Student’s t, controls extent of tails
= ∞ is Gaussian, = 1 is Cauchy
Glen CowanRHUL Physics
7
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
8
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
9
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
10
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
11
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
12
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
13
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
14
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
15
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
16
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
17
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
18
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
19
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
20
Gaussian fit Student’s t fit
Glen CowanRHUL Physics
21
Gaussian fit parameters
Glen CowanRHUL Physics
22
Student’s t fit parameters
For 1.0 < E < 1.5 GeV, →∞(consistent with Gaussian).
Glen CowanRHUL Physics
23
Andrew Wagner, 29 Aug 06 Neutrals meeting, showedsimilar study with L/R asymmetric Gaussian smearing:
Glen CowanRHUL Physics
24
Andrew Wagner, 29 Aug 06 Neutrals meeting:
Glen CowanRHUL Physics
25
From here,
Parameterize the fitted parameters vs energy, (e.g. polynomial),then for each photon
generate z ~ s(z;(E)) and replace E → E (1 + z)
Resulting distribution will not be exactly equal to smeared MCfrom fit due to several approximations made, but should be close.
To what extent are the tails due solely to EMC response? Are theyin part caused by, e.g., modeling of tracking, backgrounds,...?
Some further steps:
investigate angular dependence, refine E binning,
need more MC data (also use more real data?),
check effect on other quantities, e.g., 0 peak