Estimate of the ANTARES sensitivity for showering
Download
Report
Transcript Estimate of the ANTARES sensitivity for showering
ANTARES sensitivity for Diffuse
Fluxes and for 1000 d livetime
I. Single flavor (nm)
II. All flavors (showering events)
M. Spurio
19/06/2013
1
Sensitivity definition
90%
m90
ns
test
m90
• All analyses optimized the sensitivity using the MRF
methods from Hill-Rawlins.
•m90 (nb) depends on the
number of background
events. See plot.
• ns represent the number
of signal events. Is the sum
of all events in the final step
from all neutrino and
antineutrino channels
•The test flux test
represents the sum of the
input flux for each n and n
channel
nb
2
I- nm Analyses
1) Paper PLB 696 (2011) 16–22
2) Jutta Analysis
• Estimate of the sensitivity (nm) assuming a
reference livetime of 1000 d
3
1) Sensitivity from the PLB paper
It includes the possible contribution from prompt
nb
m90
ns
MRF= test
90%(nm)
m90/ns GeVcm-2s-1sr-1 GeVcm-2s-1sr-1
Original
(334d)
12.7
7.4
10.8
0.68
(1/2)x2x10-7
0.68x10-7
1000 d
38.1
12
32.4
0.37
10-7
0.37x10-7
1000 d
32.1
11
32.4
0.34
10-7
0.34x10-7
(no prompt)
4
2) Sensitivity from the Jutta
It does not include the contribution from prompt
m90
nb
ns
MRF= test
90%(nm)
m90/ns GeVcm-2s-1sr-1 GeVcm-2s-1sr-1
Original
(838 d)
6.68
5.7
3.22
1.77
2x10-8
0.36x10-7
1000 d
8
6.2
3.8
1.63
2x10-8
0.33x10-7
(no prompt)
5
II- Shower analyses (all flavors)
• 4 different documents analyzed
–
–
–
–
Giada Carminati PhD thesis
Qader PhD thesis
Lucia Ambrogi Master Thesis
Thomas Folger, in progress (as at the Marseille meeting)
• Estimate of the sensitivity (ne+nm+nt) assuming a
reference livetime of 1000 d
• Common problem: simulation of the nt int. missing
• Two scenarios for the nt contribution
– Pessimistic = 0 events
– Optimistic = 70% of the ne (for ne, no Glashow resonance)
6
1. Giada Carminati thesis
• Discovered a lot of bugs in the MC
• Some numbers from wrong MC (Glashow
events for upgoing events, NC bugs, etc…)
corrected “by hand”
• The analysis was supposed to be valid for 4p
sr (no dedicated cuts to remove down going
events)
• Sensitivity computed for 365 days
7
Sensitivity from Giada
nb
m90
nsh
nt
n s=
MRF= test
90%(all flavors)
nsh+nt m90/ns GeVcm-2s-1sr-1 GeVcm-2s-1sr-1
Original
(365d)
0.27
2.67
4.4
0
4.4
0.61
10-7
0.61x10-7
(single fl.)
Pessimis
0.74
3.2
12.0
0
12.0
0.26
6x10-7
1.6x10-7
Optimis
0.74
3.2
12.0
4.8*
16.8
0.19
6x10-7
1.2x10-7
* Removed the events in the Glashow resonance
8
2. Qader PhD (unblinded)
•
•
•
•
656 d livetime up to 2011
Upward going events only (<2p sr)
Unblinded: no survived events (signal =0)
Upper limit can be set (better that sensitivity)
9
Sensitivity from Qader (KNN)
nb
m90
nsh
nt
n s=
MRF= test
90%(all flavors)
nsh+nt m90/ns GeVcm-2s-1sr-1 GeVcm-2s-1sr-1
Original
(656d)
1.5
3.4
3.42
0
3.5
1
4x10-7
4x10-7
(ne+nm)
Pessimis
2.3
4.0
5.2
0
5.2
0.77
6x10-7
4.6x10-7
Optimis
2.3
4.0
5.2
3.6
8.8
0.45
6x10-7
2.7x10-7
10
3. Lucia Ambrogi Master
• Based on the reconstruction of Qader
• Different cuts to include also the downward
going events (however mostly rejected at level
of MVA)
• Sensitivity computed for 1000 days
• Slightly better than Qader
11
Sensitivity from Lucia
nb
Original
(1000 d)
1.63
m90
nsh
nt
n s=
MRF= test
90%(all flavors)
nsh+nt m90/ns GeVcm-2s-1sr-1 GeVcm-2s-1sr-1
3.5
5.18
0
5.2
0.67
4x10-7
2.7x10-7
(ne+nm)
Pessimis
1.63
3.5
5.18
0
5.2
0.67
6x10-7
4.0x10-7
Optimis
1.63
3.5
5.18
3.5
8.7
0.40
6x10-7
2.4x10-7
12
Florian Folger (Dusj in progress)
• See presentation at the last coll. meeting
13
Sensitivity from Florian
nb
m90
nsh
nt
n s=
MRF= test
90%(all flavors)
nsh+nt m90/ns GeVcm-2s-1sr-1 GeVcm-2s-1sr-1
Original
(88 days)
3.68
4.5
0.92+
0.07
0
0.99
4.5
4x4.5x10-8
8x10-7
(ne+nm)
Pessimis
42
13
11.3
0
11.3
1.15
6x4.5x10-8
3.1x10-7
Optimis
42
13
11.3
7.2*
18.5
0.70
6x4.5x10-8
1.9x10-7
* I cannot evaluate the contribution of Glashow resonance, if any
14
Summary for all flavor
Analysis
Pessimistic
Optimistic
Note
E290%(all flavors) GeVcm-2s-1sr-1
Giada
1.6x10-7
1.2x10-7
4p. Some extrapolations. Glashow
resonance contribute 40% to the signal.
Florian
3.1x10-7
1.9x10-7
Lucia
4.0x10-7
2.4x10-7
Qader
4.6x10-7
2.7x10-7
A reasonable value to be quoted seems to be
E290%(all flavors) = 3x10-7 GeVcm-2s-1sr-1
This corresponds to a single flavor: 1/3xE290%(all flavors) = 10-7 GeVcm-2s-1sr-1
15
Summary for all flavor
Analysis
Pessimistic
Optimistic
Note
E290%(all flavors) GeVcm-2s-1sr-1
1.6x10-7
Giada
1.2x10-7
4p. Some extrapolations. Glashow
resonance contribute 40% to the signal.
Florian
3.1x10-7
1.9x10-7
Lucia
4.0x10-7
2.4x10-7
Qader
4.6x10-7
2.7x10-7
and single flavor from nm
Analysis
Sensitivity
Note
E290%(nm)
GeVcm-2s-1sr-1
PRL paper+Biagi
0.37x10-7
Jutta
0.33x10-7
Prompt contribution included
16