New CORSIKA extended history simulations
New CORSIKA datasets are simulated and stored at:
- /data/sim/IceCube/2023/generated/CORSIKA_EHISTORY/
(test datasets: 30010, 30011, 30012, 30013)
- /data/sim/IceCube/2023/generated/CORSIKA_EHISTORY/generated/
(official datasets: 22774, 22775, 22776, 22777, 22778)
At first, in 2023 test simulations were performed using the software simulation_scripts. This software is a framework similar to IceProd, but it is used for local simulations on our Madison cluster. It imports all the icetray functions that are used in any other IceCube simulation. The main difference is that this framework is maintained locally in Dortmund and thus it’s easier to adapt to our needs.
The test datasets are named 30010, 30011, 30012 and 30013. Caution, these numbers do not correspond to the official dataset numbers.
We used these simulations to develop and test this analysis of the prompt component of the atmospheric muon flux.
First dataset explorations are stored in Appendix/First extended history simulations.
For the official analysis, we provide an official dataset simulated with IceProd to make it re-producible for the entire IceCube community.
Official CORSIKA Ehist IceProd simulation
A preliminary large-scale simulation is currently performed with Iceprod. The datasets are:
22774: 1e4 GeV - 1e6 GeV
22775: 1e6 GeV - 1e8 GeV
22776: 1e8 GeV - 1e9 GeV
22777: 1e9 GeV - 1e10 GeV
22778: 1e10 GeV - 1e11 GeV
The following settings are used:
CORSIKA version 77500
SIBYLL 2.3d
Icetray 1.11.0-rc1
5 components (p, He, N, Al, Fe)
Zenith angle: 0 - 90 degrees
Polyplopia: False (no coincident events are simulated)
Ecuts1: 273 GeV (hadron min energy)
Ecuts2: 273 GeV (muon min energy)
Ecuts3: 273 GeV (electron min energy)
Ecuts4: 273 GeV (photon min energy)
TrimShower: False
Atmosphere: all 12 seasons (defined here)
The detailed settings can be found in the config files at IceProd
In the following two figures, the primary energy distribution of the simulated CORSIKA Ehist IceProd datasets is shown on level 2 for the different primary particles. No cuts or filters are applied. This is not the entire statistics, but a subset of the data. The relative ratios of the primaries are different for each dataset to match the contributions of the primaries to the weighted, physical primary flux.
Fig. 2 : Primary energy distribution of the simulated CORSIKA Ehist IceProd datasets is shown on level 2. Not cuts or filters are applied. This is not the entire statistics, but a subset of the data.
Fig. 3 : Primary energy distribution of the simulated CORSIKA Ehist IceProd datasets is shown on level 2. The weights are applied using the GaisserH3a flux model.
v1.11.0-rc1 code fix
Note
In the icetray version 1.11.0-rc1, a bug was introduced that we fixed locally without committing to move on with our simulation. This fix is mentioned below.
The following code modifications were made:
diff --git a/KalmanFilter/python/KalmanFilter.py b/KalmanFilter/python/KalmanFilter.py
index 8e6c04d..1616a90 100644
--- a/KalmanFilter/python/KalmanFilter.py
+++ b/KalmanFilter/python/KalmanFilter.py
@@ -30,8 +30,12 @@ class SlowMPHit:
def hitfilter(frame, inputMapName, ignoreDC=False):
- inputMap = dataclasses.I3RecoPulseSeriesMap.from_frame(frame, inputMapName)
-
+ #inputMap = dataclasses.I3RecoPulseSeriesMap.from_frame(frame, inputMapName)
+ if type(frame[inputMapName]) == dataclasses.I3RecoPulseSeriesMapMask:
+ inputMap = frame[inputMapName].apply(frame)
+ else:
+ inputMap = frame[inputMapName]
+
hitlist = [
SlowMPHit(omkey, pulse.time, pulse)
for omkey, recoPulseVector in inputMap
diff --git a/sim-services/private/sim-services/MCPEMCPulseTools.cxx b/sim-services/private/sim-services/MCPEMCPulseTools.cxx
index b29234f..1d02365 100644
--- a/sim-services/private/sim-services/MCPEMCPulseTools.cxx
+++ b/sim-services/private/sim-services/MCPEMCPulseTools.cxx
@@ -118,11 +118,11 @@ public:
PushFrame(frame);
return;
}
- if(frame->Has(inputName+"ParticleIDMap")){
- log_warn_stream("Frame already contains " << (inputName+"ParticleIDMap")
- << " suggesting that " << inputName << " is already compressed!"
- << " Make sure that merging again is what you intent do.");
- }
+ //if(frame->Has(inputName+"ParticleIDMap")){
+ // log_warn_stream("Frame already contains " << (inputName+"ParticleIDMap")
+ // << " suggesting that " << inputName << " is already compressed!"
+ // << " Make sure that merging again is what you intent do.");
+ //}
//We have to copy all of the input data. This is gonna hurt.
boost::shared_ptr<I3Map<OMKey,std::vector<I3MCPE>>> hits=
diff --git a/topeventcleaning/private/topeventcleaning/I3TankPulseMerger.cxx b/topeventcleaning/private/topeventcleaning/I3TankPulseMerger.cxx
index e57d647..6977da4 100644
--- a/topeventcleaning/private/topeventcleaning/I3TankPulseMerger.cxx
+++ b/topeventcleaning/private/topeventcleaning/I3TankPulseMerger.cxx
@@ -144,10 +144,10 @@ void I3TankPulseMerger::GrabBadTanksAndDOMs(I3FramePtr frame)
if (badDOMList) {
badDOMs_.clear();
badDOMs_.insert(badDOMList->begin(), badDOMList->end());
- } else {
- log_warn("Bad DOM list '%s' not in the frame. Not updating.",
- badDOMListName_.c_str());
- }
+ }// else {
+ // log_warn("Bad DOM list '%s' not in the frame. Not updating.",
+ // badDOMListName_.c_str());
+ //}
}
if (!badTankListName_.empty()) {
@@ -156,10 +156,10 @@ void I3TankPulseMerger::GrabBadTanksAndDOMs(I3FramePtr frame)
if (badTankList) {
badTanks_.clear();
badTanks_.insert(badTankList->begin(), badTankList->end());
- } else {
- log_warn("Bad tank list '%s' not in the frame. Not updating.",
- badTankListName_.c_str());
- }
+ }// else {
+ // log_warn("Bad tank list '%s' not in the frame. Not updating.",
+ // badTankListName_.c_str());
+ //}
}