def fit_numerator():
'''
Fits the R9 scaling and make the plots for the situation when the nearest
muon crosses the 3x3 crystals centered on the photon seed.
'''
fitter.name = base_name + '_muInNumerator'
fitter.cuts = base_cuts + [
## Require photons to be in squre of side N-1 = 0, 1, ..., centered
## on the central crystal.
('TMath::Max(abs(muNearIEtaX - phoIEtaX), '
'abs(muNearIPhiY - phoIPhiY)) < 2'),
]
fitter.labels = [
'#mu hits 3x3',
'(R_{9} numer. and denom.)',
'',
'E_{T}^{#gamma} > 25 GeV',
'Endcaps',
]
scalings.append(fitter.main())
names.append(fitter.name)
def fit_outside():
'''
Fits the R9 scaling for a number situation when the muon does not cross
the SC and is further and further away from it.
'''
for index in range(2,10):
side = 2 * index + 1
fitter.name = base_name + '_muOutside%dx%d' % (side, side)
fitter.cuts = base_cuts + [
## Require photons to be in squre of side N-1 = 0, 1, ..., centered
## on the central crystal.
('TMath::Max(abs(muNearIEtaX - phoIEtaX), '
'abs(muNearIPhiY - phoIPhiY)) > %d' % index),
]
fitter.labels = [
'#mu outside %dx%d' % (side, side),
'',
'',
'E_{T}^{#gamma} > 25 GeV',
'Endcaps',
]
scalings.append(fitter.main())
names.append(fitter.name)
