adding neutron threshold and cleaning up configurations

Author: pgreen135

sbncode/SBNEventWeight/Calculators/Geant4/Geant4WeightCalc.cxx

@@ -186,6 +186,8 @@ void Geant4WeightCalc::Configure(fhicl::ParameterSet const& p,
     FitParNominals.push_back(theNominal);
     FitParSigmas.push_back(theSigma);
 
+    if (fDebug) std::cout << "Name: " << theName << ", Nominal: " << theNominal << ", Sigma: " << theSigma << std::endl;
+
     theNominals[theName] = theNominal;
 
     fParameterSet.AddParameter(theName, theSigma);
@@ -514,8 +516,7 @@ std::vector<float> Geant4WeightCalc::GetWeight(art::Event& e, size_t itruth ) {
         }
         std::cout << std::endl;
       }
-
-    } // if ( ( TMath::Abs(p_PDG) == 211 || p_PDG == 2212 ) )
+    }
     if (fMakeOutputTrees) fOutTree_Particle->Fill();
   } // loop over mcparticles (i)
   if (fMakeOutputTrees) fOutTree_MCTruth->Fill();

sbncode/SBNEventWeight/jobs/geant4/eventweight_geant4_sbn.fcl

@@ -138,7 +138,6 @@ sbn_eventweight_geant4: {
 # - Sigma defines the 1-sigma range for multisims and unisims. Currently the
 #   code can only accept a single value for sigma, giving symmetric
 #   uncertainty bounds: nominal-sigma and nominal+sigma 
-#   [Only in use for fitting purposes, instead range defined through input files - Jul 2025]
 # - Cut must be one of the following:
 #      "total" <- total scattering cross section, pulled from Geant4
 #      "abs" <- absorption cross section (pi+ and pi- only at the moment -- Oct 2019)

sbncode/SBNEventWeight/jobs/geant4/neutron_reweight_parameters.fcl

@@ -6,7 +6,7 @@ NeutronParameters: [
   {
     Name:       "fNeutronTotal"
     Cut:        "total"
-    Range:      [10.0, 2005.0]
+    Range:      [100.0, 2005.0]
     Nominal:    1.0
     Sigma:      0.2
   }

sbncode/SBNEventWeight/jobs/geant4/piminus_reweight_parameters.fcl

@@ -4,15 +4,15 @@ BEGIN_PROLOG
 
 PiMinusParameters: [ 
   {
-    Name:       "fPiMinusReacLow"
-    Cut:        "reac"
+    Name:       "fPiMinusInelLow"
+    Cut:        "inel"
     Range:      [10., 200.]
     Nominal:    1
     Sigma:      0.2
   } ,
   {
-    Name:       "fPiMinusReacHigh"
-    Cut:        "reac"
+    Name:       "fPiMinusInelHigh"
+    Cut:        "inel"
     Range:      [700., 2005.]
     Nominal:    1.0
     Sigma:      0.2
@@ -21,7 +21,6 @@ PiMinusParameters: [
   {
     Name:       "fPiMinusAbs"
     Cut:        "abs"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
@@ -30,7 +29,6 @@ PiMinusParameters: [
   {
     Name:       "fPiMinusCex"
     Cut:        "cex"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
@@ -39,7 +37,6 @@ PiMinusParameters: [
   {
     Name:       "fPiMinusDCex"
     Cut:        "dcex"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
@@ -48,7 +45,6 @@ PiMinusParameters: [
   {
     Name:       "fPiMinusPiProd"
     Cut:        "prod"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2

sbncode/SBNEventWeight/jobs/geant4/piplus_reweight_parameters.fcl

@@ -4,15 +4,15 @@ BEGIN_PROLOG
 
 PiPlusParameters: [
   {
-    Name:       "fPiPlusReacLow"
-    Cut:        "reac"
+    Name:       "fPiPlusInelLow"
+    Cut:        "inel"
     Range:      [10., 200.]
     Nominal:    1.0
     Sigma:      0.2
   } ,
   {
-    Name:       "fPiPlusReacHigh"
-    Cut:        "reac"
+    Name:       "fPiPlusInelHigh"
+    Cut:        "inel"
     Range:      [700., 2005.]
     Nominal:    1.0
     Sigma:      0.2
@@ -21,7 +21,6 @@ PiPlusParameters: [
   {
     Name:       "fPiPlusAbs"
     Cut:        "abs"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
@@ -30,7 +29,6 @@ PiPlusParameters: [
   {
     Name:       "fPiPlusCex"
     Cut:        "cex"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
@@ -39,7 +37,6 @@ PiPlusParameters: [
   {
     Name:       "fPiPlusDCex"
     Cut:        "dcex"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
@@ -48,7 +45,6 @@ PiPlusParameters: [
   {
     Name:       "fPiPlusPiProd"
     Cut:        "prod"
-    # Range:      [200.0, 700.00]
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2