Comparison between observational and MC data
Contents
Introduction
Here, I compared the observational data (IC-22) with several MCs
(the Elbert model and corsika (Sibyll, proton and iron)) with several cuts.
NPE, ZA and CoGZ distributions without cuts
The Elbert model express the observational data well. There is a bit difference in ZA between the observational data and corsika (iron) data.
|
NPE distribution (without any cut)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (without any cut)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
|
CoGZ distribution (without any cut)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution (without any cut)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution (without any cut)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution with several cuts
CoGZ cut
There is a difference between observational data and both the Elbert model and corsika data
with CoGZ < -200m (deep ice).
|
NPE distribution (CoGZ > 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ < 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
|
NPE distribution (CoGZ > 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (0m < CoGZ < 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (-200m < CoGZ < 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ < -200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ > -200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA cut
The Elbert model express the ZA dependence of observational data relatively well.
The agreement between the observational data and corsika data becomes worse as
the ZA becomes horizontal.
|
NPE distribution (cos(ZA) > 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (0.6 < cos(ZA) < 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (0.4 < cos(ZA) < 0.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (0.2 < cos(ZA) < 0.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (0. < cos(ZA) < 0.2)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA cut for CoGZ < -200m
The Elbert model has a small disagreement with observational data for CoGZ < -200m and cos(ZA)>0.8 (vertical)
The agreement between data and corsika (proton) becomes worse as the ZA increase (horizontal).
|
NPE distribution (CoGZ < -200m, cos(ZA) > 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ < -200m, 0.6 < cos(ZA) < 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ < -200m, 0.4 < cos(ZA) < 0.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ < -200m, 0. < cos(ZA) < 0.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA cut for CoGZ > -200m
The agreement is relatively good for events with CoGZ > -200m.
|
NPE distribution (CoGZ > -200m, cos(ZA) > 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ > -200m, 0.6 < cos(ZA) < 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ > -200m, 0.4 < cos(ZA) < 0.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE distribution (CoGZ > -200m, 0. < cos(ZA) < 0.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution with several cuts
CoGZ cut
You see again the discrepancy between observational data and corsika for horizontal at deep ice.
Actually, there is not so much discrepancy for events whose CoGZ is above center of the detector.
|
ZA distribution (CoGZ > 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (CoGZ < 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
|
ZA distribution (CoGZ > 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (0m < CoGZ < 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (-200m < CoGZ < 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (CoGZ < -200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (CoGZ > -200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE cut
There is no strong NPE dependence.
|
ZA distribution (4. < log10(NPE) < 4.2)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (4.2 < log10(NPE) < 4.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (4.4 < log10(NPE) < 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (log10(NPE) > 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution with several cuts (with MC truth)
CoGZ cut
|
ZA distribution (CoGZ > 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (CoGZ < 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
|
ZA distribution (CoGZ > 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (0m < CoGZ < 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (-200m < CoGZ < 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (CoGZ < -200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA distribution (CoGZ > -200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
The zenith angle resolution (first guess)
|
ZA resolution (log(NPE) > 4.)
red: juliet
blue: corsika
|
|
CoGZ Vs ZA resolution (log(NPE) > 4.) (juliet)
|
CoGZ Vs ZA resolution (log(NPE) > 4.) (coriska)
|
|
ZA resolution Vs NPE (juliet)
|
ZA resolution Vs NPE (coriska)
|
|
ZA resolution Vs Energy (juliet)
|
ZA resolution Vs Energy (coriska)
|
|
ZA resolution (log(NPE) > 4., CoGZ > 200m)
red: juliet
blue: corsika
|
ZA resolution (log(NPE) > 4., 0 < CoGZ < 200m)
red: juliet
blue: corsika
|
ZA resolution (log(NPE) > 4., -200 < CoGZ < 0m)
red: juliet
blue: corsika
|
ZA resolution (log(NPE) > 4., CoGZ < -200m)
red: juliet
blue: corsika
|
The zenith angle resolution (first guess) (weighted)
|
ZA resolution (log(NPE) > 4.)
red: juliet
blue: corsika
|
|
CoGZ Vs ZA resolution (log(NPE) > 4.) (juliet)
|
CoGZ Vs ZA resolution (log(NPE) > 4.) (coriska)
|
|
ZA Vs ZA resolution (log(NPE) > 4.) (juliet)
|
ZA Vs ZA resolution (log(NPE) > 4.) (coriska)
|
|
ZA resolution Vs NPE (juliet)
|
ZA resolution Vs NPE (coriska)
|
|
ZA resolution Vs Energy (juliet)
|
ZA resolution Vs Energy (coriska)
|
|
ZA resolution (log(NPE) > 4., CoGZ > 250m)
red: juliet
blue: corsika
|
ZA resolution (log(NPE) > 4., 0 < CoGZ < 250m)
red: juliet
blue: corsika
|
ZA resolution (log(NPE) > 4., -250 < CoGZ < 0m)
red: juliet
blue: corsika
|
ZA resolution (log(NPE) > 4., CoGZ < -250m)
red: juliet
blue: corsika
|
CoGZ distribution with several cuts
ZA cut (4. < log10(NPE) < 5.)
Generally speaking, the Elbert model reproduce the observed CoGZ distribution nicely,
though again there is slight disagreement at deep ice for vertical events.
The shape of CoGZ for Corsika and obsevational data seems to agree.
|
CoGZ distribution (cos(ZA) > 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ distribution (0.6 < cos(ZA) < 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ distribution (0.4 < cos(ZA) < 0.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ distribution (0. < cos(ZA) < 0.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE cut
|
CoGZ distribution (4. < log10(NPE) < 4.2)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ distribution (4.2 < log10(NPE) < 4.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ distribution (4.4 < log10(NPE) < 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ distribution (log10(NPE) > 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution with several cuts
ZA cut (4. < log10(NPE) < 5.)
|
CoGX distribution (cos(ZA) > 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution (0.6 < cos(ZA) < 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution (0.4 < cos(ZA) < 0.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution (0. < cos(ZA) < 0.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE cut
|
CoGX distribution (4. < log10(NPE) < 4.2)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution (4.2 < log10(NPE) < 4.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution (4.4 < log10(NPE) < 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGX distribution (log10(NPE) > 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution with several cuts
ZA cut (4. < log10(NPE) < 5.)
|
CoGY distribution (cos(ZA) > 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution (0.6 < cos(ZA) < 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution (0.4 < cos(ZA) < 0.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution (0. < cos(ZA) < 0.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE cut
|
CoGY distribution (4. < log10(NPE) < 4.2)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution (4.2 < log10(NPE) < 4.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution (4.4 < log10(NPE) < 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGY distribution (log10(NPE) > 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution with several cuts
ZA cut (4. < log10(NPE) < 5.)
|
NDOM distribution (cos(ZA) > 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution (0.6 < cos(ZA) < 0.8)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution (0.4 < cos(ZA) < 0.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution (0. < cos(ZA) < 0.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE cut
|
NDOM distribution (4. < log10(NPE) < 4.2)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution (4.2 < log10(NPE) < 4.4)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution (4.4 < log10(NPE) < 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution (log10(NPE) > 4.6)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ cut
|
NDOM distribution
(4. < log10(NPE) < 5., CoGZ > 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution
(4. < log10(NPE) < 5., 0 < CoGZ < 200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution
(4. < log10(NPE) < 5., -200 < CoGZ < 0m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NDOM distribution
(4. < log10(NPE) < 5., CoGZ < -200m)
black: obs. data
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution with several cuts
|
Occupancy distribution (ATWD)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA cut (4. < log10(NPE) < 5.)
|
Occupancy distribution (cos(ZA) > 0.8)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (0.6 < cos(ZA) < 0.8)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (0.4 < cos(ZA) < 0.6)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (0.2 < cos(ZA) < 0.4)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (0. < cos(ZA) < 0.2)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE cut
|
Occupancy distribution (4. < log10(NPE) < 4.2)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (4.2 < log10(NPE) < 4.4)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (4.4 < log10(NPE) < 4.6)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (4.6 < log10(NPE) < 4.8)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution (log10(NPE) > 5.)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ cut
|
Occupancy distribution
(4. < log10(NPE) < 5., CoGZ > 250m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution
(4. < log10(NPE) < 5., 0 < CoGZ < 250m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution
(4. < log10(NPE) < 5., -250 < CoGZ < 0m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Occupancy distribution
(4. < log10(NPE) < 5., CoGZ < -250m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution with several cuts
|
Nearest distance distribution
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
ZA cut (4. < log10(NPE) < 5.)
|
Nearest distance distribution (cos(ZA) > 0.8)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution (0.6 < cos(ZA) < 0.8)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution (0.4 < cos(ZA) < 0.6)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution (0. < cos(ZA) < 0.4)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
NPE cut
|
Nearest distance distribution (4. < log10(NPE) < 4.2)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution (4.2 < log10(NPE) < 4.4)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution (4.4 < log10(NPE) < 4.6)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution (log10(NPE) > 4.6)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
CoGZ cut
|
Nearest distance distribution
(4. < log10(NPE) < 5., CoGZ > 200m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution
(4. < log10(NPE) < 5., 0 < CoGZ < 200m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution
(4. < log10(NPE) < 5., -200 < CoGZ < 0m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Nearest distance distribution
(4. < log10(NPE) < 5., CoGZ < -200m)
blue: Elbert model
red: corsika (proton)
purple: corsika (iron)
|
Channel-wise NPE distributions
|
Channel-wise ATWD NPE distribution (level1 (NDOM>=80))
|
Channel-wise FADC NPE distribution (level1 (NDOM>=80))
|
Summary
- The Elbert model express our data relatively well, except slight disagreement for vertical events whose CoGZ is placed at deep ice.
- There is a disagreement of horizontal events whose CoGZ is at deep ice between observational data and corsika (low rate). But, the shape (e.g. CoGZ) of corsika agrees relatively well with observational data. The problem seems that we don't have so much horizontal events at deep ice in Corsika...
Keiichi Mase
Last modified: 2008-12-22 23:22:55