Secondary beam analysis – PS and SPS
Started May 22, 2012
rpc/PS_pion
pion_2.f
pion_11.kumac
or
pion_11_1GeV.kumac
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1
pions: C1 = C2 = 0 and BC<8
muons: N > 20 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Not enough memory for run 109
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = C2 = 0 and BC<8 and R > 2.0
muons: N > 10 and C1 = C2 = 0 and BC>20
Hits in last 2 layers
rpc/PS_pion/PS_summary.kumac
No calibration whatsoever
Calibration with muon response
Fits

E

E

E




E

E

E
c

β

E
c
α
c
Fit 1
65.67± 0.1
Fit 2
61.3± 0.3
11.0± 0.3
61.3 ± 0.2
11.0 ± 0.2
Fit 3
0.0 ± 2.9
Calibration study
rpc/PS_pion/PS_summary_alpha.kumac
runs calib_1.f
 N
 N 
N
 
0
Calib
N





α=0.4
This is interesting. It appears
that a linear calibration does
not work as well as a calibration
attenuated by a power…
Positive beam
rpc/PS_pion
pion_2.f
pion_11_pos.kumac
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = 0 and C2 = 1 and BC<8 and
R > 2.0
protons: C1 = 0 and C2 = 0 and BC<8 and
R > 2.0
muons: N > 10 and C1 = 0 and C2 = 1 and
BC>20 and Hits in last 2 layers
and r<2.0
Run 660177
Run 660178
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = 0 and C2 = 1 and BC<8 and
R > 2.0
protons: C1 = 0 and C2 = 0 and BC<8 and
R > 2.0
muons: N > 10 and C1 = 0 and C2 = 1 and
BC>20 and Hits in last 2 layers
and r<2.0
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = 0 and C2 = 1 and BC<8 and
R > 2.0
protons: C1 = 0 and C2 = 0 and BC<8 and
R > 2.0
muons: N > 10 and C1 = 0 and C2 = 1 and
BC>20 and Hits in last 2 layers
and r<2.0
Event selection
all: 1 cluster in layer 0 with at most 3 hits
electrons: C1 = C2 =1 and BC < 8
pions: C1 = 0 and C2 = 1 and BC<8 and
R > 2.0
protons: C1 = 0 and C2 = 0 and BC<8 and
R > 2.0
muons: N > 10 and C1 = 0 and C2 = 1 and
BC>20 and Hits in last 2 layers
and r<2.0
β
α
c
Fit 1
65.67± 0.1
Fit 2
61.3± 0.3
11.0± 0.3
61.3 ± 0.2
11.0 ± 0.2
54.4 ± 0.6
20.2 ± 1.7
Fit 3
p fit
0.0 ± 2.9
SPS Run Plan
Rates
Spills of 9 seconds every 45 seconds -> 720 seconds/hour or 17,280 seconds/day
At 50 GeV assume maximum of 100 Hz -> 72,000 events/hour or 1,728,000 events/day
At 100 GeV assume maximum of 50 Hz -> 36,000 events/hour or 864,000 events/day
Data sets
Assume minimum of 500,000 events per energy point
Overall plan
Start June 4th in evening (days start at 8 am)
4,5,7,8 -> -50 GeV
6 -> Machine development
9,10, 11, 12 -> -180
13,14,15,16,17 -> +180
July 30 - August 6 -> 7 days: muons, (+180 GeV), high energy run (<300 GeV)
November 2-9 -> 7 days: muons (+180 GeV), reserve or high energy again
-120 GeV running (4.5 days)
50, 40 (e and pions), {30 (e and pions),} 60, 80, 100, 120
-180 GeV running (2 days)
120, 150, 180
-30 GeV running (1 day)
{30 (e and pions),} 20 (e and pions), {10 (e and pions)}
+180 GeV running (5 days)
Muons
60, 70, 80, 90 (with Cerenkov tagging of pions)
120, 150, 180
Event selection
all: 1 cluster in layer 0 with at most 3 hits
𝜎 95.5%
=
𝐸
√𝐸
660211
pions: C1 = 0 and C2 = 0 and IL>2 and
and no hits in last 4 layers
muons: BC>20 and Hits in last 4 layers
and r<3.0
𝜎 89.2%
=
𝐸
√𝐸
660217
Event selection
all: 1 cluster in layer 0 with at most 3 hits
𝜎 83.4%
=
𝐸
√𝐸
pions: C1 = 0 and C2 = 0 and IL>2 and
and no hits in last 4 layers
muons: BC>20 and Hits in last 4 layers
and r<3.0
Posi.701.atuple
Octagon not properly implemented
T3B – RPC not implemented
T3B – support not implemented
Resolution just a bit better
Peak at 303, instead of 255
Simulation
Event selection
all: 1 cluster in layer 0 with at most 3 hits
𝜎 98.2%
=
𝐸
√𝐸
660223
pions: No hits in last 4 layers
muons: BC>20 and hits in last 4 layers
and r<3.0
𝜎 97.3%
=
𝐸
√𝐸
660224
Event selection
all: 1 cluster in layer 0 with at most 3 hits
𝜎 106.6%
=
𝐸
√𝐸
660226
pions: No hits in last 4 layers
muons: BC>20 and hits in last 4 layers
and r<3.0
𝜎 105.6%
=
𝐸
√𝐸
660227
Event selection
all: 1 cluster in layer 0 with at most 3 hits
𝜎 106.7%
=
𝐸
√𝐸
660228
pions: No hits in last 4 layers
muons: BC>20 and hits in last 4 layers
and r<3.0
𝜎 105.9%
=
𝐸
√𝐸
660229
June 7, 2012
Study of time difference
fort.91
pion_19.kumac
- Some hits far off in both
+ and – direction
- From now on, only
take -18,-19
Event selection
all: 1 cluster in layer 0 with at most 3 hits
𝜎 117.9%
=
𝐸
√𝐸
660233
pions: No hits in last 4 layers
muons: BC>20 and hits in last 4 layers
and r<3.0
𝜎 116.8%
=
𝐸
√𝐸
660234
Momentum
# of events (18
mm Pb)
# of events
(Air)
20
41,336
438,356
30
594,132
410,731
40
325,803
50
886,201
60
497,738
80
722,268
100
526,323
120
505,465
TOTAL
Event selection
all: 1 cluster in layer 0 with at most 3 hits
𝜎 141.1%
=
𝐸
√𝐸
pions: No hits in last 4 layers
muons: BC>20 and hits in last 4 layers
and r<3.0
𝜎 116.7%
=
𝐸
√𝐸
Tail Catcher Calibration Study
DHCAL
TCMT_1
TCMT_2
8-June-2012
Run 660239
120 GeV
root_1.C
root
.L root_1.C++
.x root_1
Interaction Layer <8
χ2 fit to establish
calibration constants
DHCAL weight = 0.149
TCMT 1 weight = 0.048
TCMT 2 weight = 0.401
No apparent improvement
Run 660239
120 GeV
root_1.C
χ2 fit to establish
calibration constants
DHCAL weight = 0.151
TCMT 1 weight = 0.083
TCMT 2 weight = 0.520
No apparent improvement
Cross checks: program seems to work. Source of tail appears to be not related to
energy in tail catcher. Need to scan events in the tail.
Event selection
𝜎 128.8%
=
𝐸
√𝐸
all: 1 cluster in layer 0 with at most 3 hits
pions: No hits in last 4 layers
muons: BC>20 and hits in last 4 layers
and r<3.0
𝜎 158.4%
=
𝐸
√𝐸
Adjusted cuts
(Nhit in layer 0)-cut has minimal effect on resolution
Allowing for (Ncluster > 1 in layer 0) increases resolution
New cut: Nhit in layer 0 < 99
→ Gain a factor of ~20 in statistics
Event selection
𝜎 61.1%
=
𝐸
√𝐸
Peak at
103.0
all: 1 cluster in layer 0 with at most 99 hits
electrons: No hits in last 4 layers, r>8.
muons: BC>20 and hits in last 4 layers
and r<3.0
Peak at
106.1
𝜎 59.1%
=
𝐸
√𝐸
Event selection
all: 1 cluster in layer 0 with at most 99 hits
𝜎 56.9%
=
𝐸
√𝐸
electrons: No hits in last 4 layers, r>8.
muons: BC>20 and hits in last 4 layers
and r<3.0
Peak at
106.2
660279
𝜎 64.3%
=
𝐸
√𝐸
Peak at
137.6
660294