add comment, fix dipole portal MB example

Harvard-Neutrino · Jun 3, 2024 · d3adae2 · d3adae2
1 parent e86995f
commit d3adae2
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/projects/injection/public/SIREN/injection/TreeWeighter.tcc b/projects/injection/public/SIREN/injection/TreeWeighter.tcc
@@ -187,6 +187,8 @@ double ProcessWeighter<ProcessType>::NormalizedPositionProbability(std::tuple<si
     double interaction_density = detector_model->GetInteractionDensity(intersections, DetectorPosition(interaction_vertex), targets, total_cross_sections, total_decay_length); //units of m^-1
 
     double prob_density;
+    // This is equivalent to equation 11 of the SIREN paper
+    // Reach out to the authors if you disagree and we can send the derivation :)
     if(total_interaction_depth < 1e-6) {
         prob_density = interaction_density / total_interaction_depth;
     } else {

diff --git a/resources/Examples/Example2/DipolePortal_MiniBooNE.py b/resources/Examples/Example2/DipolePortal_MiniBooNE.py
@@ -17,7 +17,7 @@
 }
 
 # Number of events to inject
-events_to_inject = 10#0000
+events_to_inject = 100000
 
 # Expeirment to run
 experiment = "MiniBooNE"
@@ -68,8 +68,7 @@
     primary_type, primary_injection_distributions, primary_physical_distributions
 )
 
-controller.Initialize("output/MiniBooNE_Dipole_M%2.2e_mu%2.2e_example"
-                      % (model_kwargs["m4"], model_kwargs["mu_tr_mu4"]))
+controller.Initialize()
 
 def stop(datum, i):
     secondary_type = datum.record.signature.secondary_types[i]