diff --git a/IceModel.cc b/IceModel.cc
index 52ae2e12..d966df1c 100644
--- a/IceModel.cc
+++ b/IceModel.cc
@@ -271,23 +271,22 @@ if ( file.is_open() ) {
 // read depth in positive value and return attenuation length (m) at the depth
 double IceModel::GetARAIceAttenuLength(double depth, Settings *settings1) {
 
-    double AttenL=0.0;
-
-    // check if depth is positive value
-    if ( depth < 0. ) {// whether above the ice or wrong value!
-        #ifdef VERBOSE_MODE
-        std::cerr << "depth negative! " << depth << std::endl;
-        #endif
-    }
-    else {
-
-        AttenL = Tools::SimpleLinearInterpolation_extend_Single(ARA_IceAtten_bin, ARA_IceAtten_Depth, ARA_IceAtten_Length, depth );
-        AttenL *= GetIceAttenSystematicsFactor(depth, settings1); //Systematics
-
-    }
+  double AttenL = 0.0;
 
+  // check if depth is positive value
+  if ( depth < 0. ) {// whether above the ice or wrong value!
+    #ifdef VERBOSE_MODE
+    std::cerr << "depth negative! " << depth << std::endl;
+    #endif
     return AttenL;
+  }
+  if ( depth > UNRELIABLE_DEPTH_M ) return 0.;
+  if ( depth > BEDROCK_DEPTH_M )    return 0.;
+
+  AttenL = Tools::SimpleLinearInterpolation_extend_Single(ARA_IceAtten_bin, ARA_IceAtten_Depth, ARA_IceAtten_Length, depth );
+  AttenL *= GetIceAttenSystematicsFactor(depth, settings1); //Systematics
 
+  return AttenL;
 }
 
 
@@ -299,17 +298,34 @@ double IceModel::GetARAIceAttenuLength(double depth, Settings *settings1) {
     And is described here: https://icecube.wisc.edu/~araproject/radio/#figure1a
     and https://icecube.wisc.edu/~araproject/radio/temp/
 
+    The "strict" parameter controls whether or not this function will let you call
+    it above the ice or below the bedrock.
+    It can be useful for *plotting* to have access to unphysical depths,
+    but generally the code should fail loudly if you try and access this function
+    somewhere unphysical.
+    Default is therefore "true" on strict behavior.
+
     \param z depth as a positive number in meters
+    \param strict whether or not to encorce strict rules
     \return ice temperature in Celsius
 */
-double IceModel::temperature(double z){
-  if( z < 0.){
-    #ifdef VERBOSE_MODE
-    std::cerr << "depth negative! " << z << std::endl;
-    #endif
+double IceModel::temperature(double z, bool strict){
+  if( strict ){
+    if( z < 0. ){
+      throw std::runtime_error(
+        "IceModel::temperature: depth " + std::to_string(z) +
+        " m is above the ice surface (z < 0).");
+    }
+    if( z > BEDROCK_DEPTH_M ){
+      throw std::runtime_error(
+        "IceModel::temperature: depth " + std::to_string(z) +
+        " m is below bedrock (> " + std::to_string(BEDROCK_DEPTH_M) + " m).");
+    }
   }
   double temp = (-51.0696) + (0.00267687 * z) + (-1.59061E-08 * pow(z,2.)) + (1.83415E-09 * pow(z,3.));
-  return temp;
+
+  // we also clamp the output to freezing to make sure this can't return a *positive* temperature
+  return std::min(temp, 0.);
 }
 
 //! Get the ice attenuation length at a depth and frequency
@@ -320,34 +336,51 @@ double IceModel::temperature(double z){
 
     \param depth depth as a positive number in meters
     \param freq frequency as a number in GHz
+    \param settings a settings pointer so we know whether or not to apply systematics factors
     \return ice temperature in Celsius
 */
 double IceModel::GetFreqDepIceAttenuLength(double depth, double freq, Settings *settings1) {
+  // Constants for the unreliable / out-of-ice regimes.
+  // Per D. Besson, the attenuation model becomes unphysical below ~2.5 km depth.
+  // The bedrock check is a fail-safe in case the 2.5 km cutoff is ever relaxed.
+
   double AttenL = 0.0;
-  if ( depth < 0. ) {
+  if (depth < 0.) {
     #ifdef VERBOSE_MODE
     std::cerr << "depth negative! " << depth << std::endl;
     #endif
+    return AttenL;
   }
-  else {
-    double t = temperature(depth);
-    const double f0=0.0001, f2=3.16;
-    const double w0=log(f0), w1=0.0, w2=log(f2), w=log(freq);
-    const double b0=-6.74890+t*(0.026709-t*0.000884);
-    const double b1=-6.22121-t*(0.070927+t*0.001773);
-    const double b2=-4.09468-t*(0.002213+t*0.000332);
-    double a,bb;
-    if(freq<1.){
-      a=(b1*w0-b0*w1)/(w0-w1);
-      bb=(b1-b0)/(w1-w0);
-    }
-    else{
-      a=(b2*w1-b1*w2)/(w1-w2);
-      bb=(b2-b1)/(w2-w1);
-    }
-    AttenL = 1./exp(a+bb*w);
-    AttenL *= GetIceAttenSystematicsFactor(depth, settings1); //Systematics
+  if (depth > UNRELIABLE_DEPTH_M){
+    #ifdef VERBOSE_MODE
+    std::cerr << "depth is below where we trust the ice model. Depth = " << depth << std::endl;
+    #endif
+    return 0.;
+  }
+  if (depth > BEDROCK_DEPTH_M){
+    #ifdef VERBOSE_MODE
+    std::cerr << "depth is below bedrock at SP. Depth = " << depth << std::endl;
+    #endif
+    return 0.;
+  }
+
+  double t = temperature(depth);
+  const double f0=0.0001, f2=3.16;
+  const double w0=log(f0), w1=0.0, w2=log(f2), w=log(freq);
+  const double b0=-6.74890+t*(0.026709-t*0.000884);
+  const double b1=-6.22121-t*(0.070927+t*0.001773);
+  const double b2=-4.09468-t*(0.002213+t*0.000332);
+  double a,bb;
+  if(freq<1.){
+    a=(b1*w0-b0*w1)/(w0-w1);
+    bb=(b1-b0)/(w1-w0);
+  }
+  else{
+    a=(b2*w1-b1*w2)/(w1-w2);
+    bb=(b2-b1)/(w2-w1);
   }
+  AttenL = 1./exp(a+bb*w);
+  AttenL *= GetIceAttenSystematicsFactor(depth, settings1); //Systematics
   return AttenL;
 }
 
@@ -1384,135 +1417,66 @@ double IceModel::GetEffectiveN(const Position &pos) const{
   return GetEffectiveN(n_local);
 } //GetEffectiveN(Position)
 
-double IceModel::EffectiveAttenuationLength(const Position &pos,const int &whichray) const {
+double IceModel::EffectiveAttenuationLength(Settings *settings1, const Position &pos, const int &whichray) const {
+
   double localmaxdepth = IceThickness(pos);
   double depth = Surface(pos) - pos.Mag();
 
-  int depth_index=0;
-  double attenuation_length=0.0;
-//   if (inu<10) {
-//     cout << "pos is ";pos.Print();
-//     cout << "surface is " << Surface(pos) << "\n";
-//   }
-  if(WestLand(pos) && !CONSTANTICETHICKNESS)
-    {
-      depth_index=int(depth*419.9/localmaxdepth);//use 420 m ice shelf attenuation length data as the standard, squeeze or stretch if localmaxdepth is longer or shorter than 420m.
-      if(RossIceShelf(pos) || RonneIceShelf(pos))
-	{
-	  if(whichray==0)
-	    attenuation_length=l_shelfup[depth_index];
-	  else if(whichray==1)
-	    attenuation_length=l_shelfdown[depth_index];
-	  else
-	    cerr << " wrong attenuation length " <<endl;
-
-	  //for sanity check
-	  if((depth_index+0.5)!=d_shelfup[depth_index])
-	    {
-	      cerr << "the index of the array l_iceshelfup is wrong!" << endl;
-	      exit(1);
-	    }
-	}
-      else //in ice sheet of westland
-	{
-	  if(whichray==0)
-	    attenuation_length=l_westlandup[depth_index];
-	  else if(whichray==1)
-	    attenuation_length=l_westlanddown[depth_index];
-	  else
-	    cerr << " wrong attenuation length " <<endl;
-      	}
-
-      if(mooreBayFlag)//if use Moore's Bay measured data for the west land
-	attenuation_length*=1.717557; //about 450 m (field attenuation length) for one whole way when assuming -3dB for the power loss at the bottom
+  if (depth < 0.) {
+    #ifdef VERBOSE_MODE
+    std::cerr << "depth negative! " << depth << std::endl;
+    #endif
+    return 0.;
+  }
+  if (depth > UNRELIABLE_DEPTH_M) return 0.;
+  if (depth > BEDROCK_DEPTH_M)    return 0.;
+
+  int depth_index = 0;
+  double attenuation_length = 0.0;
+
+  if (WestLand(pos) && !CONSTANTICETHICKNESS) {
+    depth_index = int(depth * 419.9 / localmaxdepth);  // use 420 m ice shelf attenuation length data as the standard, squeeze or stretch if localmaxdepth is longer or shorter than 420m.
+    if (RossIceShelf(pos) || RonneIceShelf(pos)) {
+      if (whichray == 0)
+        attenuation_length = l_shelfup[depth_index];
+      else if (whichray == 1)
+        attenuation_length = l_shelfdown[depth_index];
+      else
+        cerr << " wrong attenuation length " << endl;
+
+      // for sanity check
+      if ((depth_index + 0.5) != d_shelfup[depth_index]) {
+        cerr << "the index of the array l_iceshelfup is wrong!" << endl;
+        exit(1);
+      }
+    }
+    else {  // in ice sheet of westland
+      if (whichray == 0)
+        attenuation_length = l_westlandup[depth_index];
+      else if (whichray == 1)
+        attenuation_length = l_westlanddown[depth_index];
+      else
+        cerr << " wrong attenuation length " << endl;
     }
-  else //in east antarctica or constant ice thickness
-     {
-//        if (inu<10) {
-//        cout << "localmaxdepth is " << localmaxdepth << "\n";
-//        cout << "depth is " << depth << "\n";
-//        }
-       depth_index =int(depth*(2809.9/localmaxdepth));
-       //if (inu<10)
-	 //       cout << "depth_index is " << depth_index << "\n";
-
-       if(whichray==0)
-	 attenuation_length =l_sheetup[depth_index];
-       else if(whichray==1)
-	 attenuation_length =l_sheetdown[depth_index];
-       else
-	 cerr << " wrong attenuation length " <<endl;
-     } //else
-
-  return attenuation_length;
-} //EffectiveAttenuationLengthUp
 
+    if (settings1->MOOREBAY)  // if use Moore's Bay measured data for the west land
+      attenuation_length *= 1.717557;  // about 450 m (field attenuation length) for one whole way when assuming -3dB for the power loss at the bottom
+  }
+  else {  // in east antarctica or constant ice thickness
+    depth_index = int(depth * (2809.9 / localmaxdepth));
 
-double IceModel::EffectiveAttenuationLength(Settings *settings1, const Position &pos,const int &whichray) const {
-  double localmaxdepth = IceThickness(pos);
-  double depth = Surface(pos) - pos.Mag();
+    if (whichray == 0)
+      attenuation_length = l_sheetup[depth_index];
+    else if (whichray == 1)
+      attenuation_length = l_sheetdown[depth_index];
+    else
+      cerr << " wrong attenuation length " << endl;
+  }
 
-  int depth_index=0;
-  double attenuation_length=0.0;
-//   if (inu<10) {
-//     cout << "pos is ";pos.Print();
-//     cout << "surface is " << Surface(pos) << "\n";
-//   }
-  if(WestLand(pos) && !CONSTANTICETHICKNESS)
-    {
-      depth_index=int(depth*419.9/localmaxdepth);//use 420 m ice shelf attenuation length data as the standard, squeeze or stretch if localmaxdepth is longer or shorter than 420m.
-      if(RossIceShelf(pos) || RonneIceShelf(pos))
-	{
-	  if(whichray==0)
-	    attenuation_length=l_shelfup[depth_index];
-	  else if(whichray==1)
-	    attenuation_length=l_shelfdown[depth_index];
-	  else
-	    cerr << " wrong attenuation length " <<endl;
-
-	  //for sanity check
-	  if((depth_index+0.5)!=d_shelfup[depth_index])
-	    {
-	      cerr << "the index of the array l_iceshelfup is wrong!" << endl;
-	      exit(1);
-	    }
-	}
-      else //in ice sheet of westland
-	{
-	  if(whichray==0)
-	    attenuation_length=l_westlandup[depth_index];
-	  else if(whichray==1)
-	    attenuation_length=l_westlanddown[depth_index];
-	  else
-	    cerr << " wrong attenuation length " <<endl;
-      	}
-
-      //if(mooreBayFlag)//if use Moore's Bay measured data for the west land
-      if(settings1->MOOREBAY)//if use Moore's Bay measured data for the west land
-	attenuation_length*=1.717557; //about 450 m (field attenuation length) for one whole way when assuming -3dB for the power loss at the bottom
-    }
-  else //in east antarctica or constant ice thickness
-     {
-//        if (inu<10) {
-//        cout << "localmaxdepth is " << localmaxdepth << "\n";
-//        cout << "depth is " << depth << "\n";
-//        }
-       depth_index =int(depth*(2809.9/localmaxdepth));
-       //if (inu<10)
-	 //       cout << "depth_index is " << depth_index << "\n";
-
-       if(whichray==0)
-	 attenuation_length =l_sheetup[depth_index];
-       else if(whichray==1)
-	 attenuation_length =l_sheetdown[depth_index];
-       else
-	 cerr << " wrong attenuation length " <<endl;
-     } //else
-
-    attenuation_length *= GetIceAttenSystematicsFactor(depth, settings1);
+  attenuation_length *= GetIceAttenSystematicsFactor(depth, settings1);
 
   return attenuation_length;
-} //EffectiveAttenuationLengthUp
+}  // EffectiveAttenuationLength
 
 
 
diff --git a/IceModel.h b/IceModel.h
index 4a620053..e8638241 100644
--- a/IceModel.h
+++ b/IceModel.h
@@ -111,6 +111,11 @@ void ENtoLonLat(int e_coord,
   double ground_elevation[1068][869]; //elevation above geoid at which ice starts
   double water_depth[1200][1000]; //depth of water under ice
 
+  static constexpr double BEDROCK_DEPTH_M = 2850.;  // South Pole nominal
+  static constexpr double UNRELIABLE_DEPTH_M = 2500.; // the depth below which we don't trust our ARA ice models according to D. Besson
+  static constexpr double MIN_PHYSICAL_ATTEN_M = 1.; // the smallest physical attenuation length we think we'd have (to guard against divison by zero issues)
+
+
   void IceENtoLonLat(int e,
 		     int n,
 		     
@@ -165,7 +170,6 @@ void ENtoLonLat(int e_coord,
   double GetN(const Position &pos) const;
   double GetEffectiveN(double n_local) const;
   double GetEffectiveN(const Position &pos) const;
-  double EffectiveAttenuationLength(const Position &pos, const int &whichray) const;
   double EffectiveAttenuationLength(Settings *settings1, const Position &pos, const int &whichray) const;
   
   void IceLonLattoEN(double lon, 
@@ -200,7 +204,7 @@ void GetFresnel (
   double GetFreqDepIceAttenuLength(double depth, double freq, Settings *settings1);
   bool LoadIceAttenPercentTable(const std::string& filename);
   double GetIceAttenSystematicsFactor(double depth, Settings *settings1) const;
-  double temperature(double z);
+  double temperature(double z, bool strict = true);
 
   //  void FillArraysforTree(double lon_ground[1068][869],double lat_ground[1068][869],double lon_ice[1200][1000],double lat_ice[1200][1000],double lon_water[1200][1000],double lat_water[1200][1000]);
 
diff --git a/Report.cc b/Report.cc
index 93e78a84..a4964322 100644
--- a/Report.cc
+++ b/Report.cc
@@ -965,23 +965,45 @@ void Report::ModelRay(
 
             dx = RayStep[ray_idx][0][steps - 1] - RayStep[ray_idx][0][steps];
             dz = RayStep[ray_idx][1][steps - 1] - RayStep[ray_idx][1][steps];
-            dl = sqrt((dx *dx) + (dz *dz));
+            dl = sqrt((dx * dx) + (dz * dz));
 
-            // Skipping attenuation calculation when the distance between two RaySteps is 0. 
-            // PrevenTing adds -nan into the IceAttenFactor.
+            // Skipping attenuation calculation when the distance between two RaySteps is 0.
+            // Preventing adds -nan into the IceAttenFactor.
             if (dl > 0) {
-                // use new ice model
-                // use the midpoint of the array to calculate the attenuation length, instead of the end of the ray 
-                IceAttenFactor *= (
-                    exp(-dl / icemodel->GetARAIceAttenuLength(-RayStep[ray_idx][1][steps], settings)) + 
-                    exp(-dl / icemodel->GetARAIceAttenuLength(-RayStep[ray_idx][1][steps - 1], settings))
-                ) / 2;
+                double L1 = icemodel->GetARAIceAttenuLength(-RayStep[ray_idx][1][steps],     settings);
+                double L2 = icemodel->GetARAIceAttenuLength(-RayStep[ray_idx][1][steps - 1], settings);
+
+                // If either endpoint is in unmodelable ice, treat the whole
+                // ray as fully attenuated and stop accumulating.
+                if (L1 < IceModel::MIN_PHYSICAL_ATTEN_M ||
+                    L2 < IceModel::MIN_PHYSICAL_ATTEN_M) {
+                    // if (-RayStep[ray_idx][1][steps - 1] < 2500. && -RayStep[ray_idx][1][steps] < 2500){
+                    //     printf("Start (%.2f, %.2f), End (%.2f, %.2f)\n", 
+                    //         RayStep[ray_idx][0][steps - 1], RayStep[ray_idx][1][steps - 1],
+                    //         RayStep[ray_idx][0][steps], RayStep[ray_idx][1][steps]
+                    //     );
+                    //     printf("This is firing! Tighter check! L2 and L1 are %f, %f, and dl is %f \n", L2, L1, dl);
+                    //     printf("----\n");
+                    // }
+                    // IceAttenFactor = 0.;
+                    IceAttenFactor *= 0.;
+                    // break;
+                }
+                else{
+                    // use new ice model
+                    // use the midpoint of the array to calculate the attenuation length, instead of the end of the ray
+                    IceAttenFactor *= (exp(-dl / L1) + exp(-dl / L2)) / 2;
+                }
             }
-
         }
     }
     else if (settings->USE_ARA_ICEATTENU == 0) { // use old method
-        IceAttenFactor = exp(-ray_output[0][ray_idx] / icemodel->EffectiveAttenuationLength(settings, event->Nu_Interaction[interaction_idx].posnu, 0));
+        double L = icemodel->EffectiveAttenuationLength(settings, event->Nu_Interaction[interaction_idx].posnu, 0);
+        if (L < IceModel::MIN_PHYSICAL_ATTEN_M) {
+            IceAttenFactor = 0.;
+        } else {
+            IceAttenFactor = exp(-ray_output[0][ray_idx] / L);
+        }
     }
 
     // If the simulation is not in debug mode, calculate the field from the
@@ -1048,27 +1070,39 @@ void Report::ModelRay(
                 }
                 else if (settings->USE_ARA_ICEATTENU == 2) {
 
+
                     double IceAttenFactor = 1.;
                     double dx, dz, dl;
                     for (int steps = 1; steps < (int) RayStep[ray_idx][0].size(); steps++) {
                         dx = RayStep[ray_idx][0][steps - 1] - RayStep[ray_idx][0][steps];
                         dz = RayStep[ray_idx][1][steps - 1] - RayStep[ray_idx][1][steps];
-                        dl = sqrt((dx *dx) + (dz *dz));
+                        dl = sqrt((dx * dx) + (dz * dz));
 
-                        // Skipping attenuation calculation when the distance between two RaySteps is 0. 
-                        // Prevening adds -nan into the IceAttenFactor.
+                        // Skipping attenuation calculation when the distance between two RaySteps is 0.
+                        // Preventing adds -nan into the IceAttenFactor.
                         if (dl > 0)
                         {
                             // use ray midpoint for attenuation calculation
-                            IceAttenFactor *= (exp(-dl / icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps], detector->GetFreq(l) / 1e9, settings)) +
-                                exp(-dl / icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps - 1], detector->GetFreq(l) / 1e9,  settings))
-                            ) / 2.;  // 1e9 to convert to GHz
+                            // But, if either endpoint is in unmodelable ice (L < 1 m), the ray
+                            // segment can't be evaluated with any confidence. Treat the
+                            // whole ray as fully attenuated and stop accumulating.
+
+                            double freq_GHz = detector->GetFreq(l) / 1e9;
+                            double L1 = icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps],     freq_GHz, settings);
+                            double L2 = icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps - 1], freq_GHz, settings);
+                
+                            if (L1 < IceModel::MIN_PHYSICAL_ATTEN_M || L2 < IceModel::MIN_PHYSICAL_ATTEN_M) {
+                                IceAttenFactor = 0.;
+                                break;
+                            }
+                
+                            IceAttenFactor *= (exp(-dl / L1) + exp(-dl / L2)) / 2.;
                         }
                     }
-                    
+                
                     // assume whichray = 0, now vmmhz1m_tmp has all factors except for the detector properties (antenna gain, etc)
-                    vmmhz1m_tmp = vmmhz1m_tmp / ray_output[0][ray_idx] *IceAttenFactor *mag * fresnel;
-
+                    vmmhz1m_tmp = vmmhz1m_tmp / ray_output[0][ray_idx] * IceAttenFactor * mag * fresnel;
+                
                 }
 
                 vmmhz1m_sum += vmmhz1m_tmp;
@@ -1644,22 +1678,31 @@ void Report::PropagateSignal(
                 for (int steps = 1; steps < (int) RayStep[ray_idx][0].size(); steps++) {
                         dx = RayStep[ray_idx][0][steps - 1] - RayStep[ray_idx][0][steps];
                         dz = RayStep[ray_idx][1][steps - 1] - RayStep[ray_idx][1][steps];
-                        dl = sqrt((dx *dx) + (dz *dz));
+                        dl = sqrt((dx * dx) + (dz * dz));
 
-                        // Skipping attenuation calculation when the distance between two RaySteps is 0. 
-                        // Prevening adds -nan into the IceAttenFactor. 
+                        // Skipping attenuation calculation when the distance between two RaySteps is 0.
+                        // Preventing adds -nan into the IceAttenFactor.
                         if (dl > 0) {
+                            double freq_GHz = freq_tmp * 1.E-9;
+                            double L1 = icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps],     freq_GHz, settings);
+                            double L2 = icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps - 1], freq_GHz, settings);
+
+                            // If either endpoint is in unmodelable ice, treat the whole
+                            // ray as fully attenuated and stop accumulating.
+                            if (L1 < IceModel::MIN_PHYSICAL_ATTEN_M ||
+                                L2 < IceModel::MIN_PHYSICAL_ATTEN_M) {
+                                IceAttenFactor = 0.;
+                                break;
+                            }
+
                             // use ray midpoint for attenuation calculation
-                            IceAttenFactor *= (
-                                exp(-dl / icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps],     freq_tmp *1.E-9, settings)) +
-                                exp(-dl / icemodel->GetFreqDepIceAttenuLength(-RayStep[ray_idx][1][steps - 1], freq_tmp *1.E-9, settings))
-                            ) / 2.;  // 1e9 for conversion to GHz
+                            IceAttenFactor *= (exp(-dl / L1) + exp(-dl / L2)) / 2.;
                         }
 
                 }
-                
-                V_forfft[2 *n] *= IceAttenFactor;   // apply IceAttenFactor to the real part of fft
-                V_forfft[2 *n + 1] *= IceAttenFactor;   // apply IceAttenFactor to the imag part of fft
+
+                V_forfft[2 * n] *= IceAttenFactor;       // apply IceAttenFactor to the real part of fft
+                V_forfft[2 * n + 1] *= IceAttenFactor;   // apply IceAttenFactor to the imag part of fft
 
             }
 
diff --git a/examples/atten_length/README.md b/examples/atten_length/README.md
new file mode 100644
index 00000000..616c7905
--- /dev/null
+++ b/examples/atten_length/README.md
@@ -0,0 +1,9 @@
+Run these two functions to see plots of the attenuation length as a function of depth.
+
+Run as:
+```
+python get_atten.py
+python plot_atten.py
+```
+
+They do the obvious thing you think.
\ No newline at end of file
diff --git a/examples/atten_length/get_atten.py b/examples/atten_length/get_atten.py
new file mode 100644
index 00000000..0b80b82e
--- /dev/null
+++ b/examples/atten_length/get_atten.py
@@ -0,0 +1,35 @@
+#!/usr/bin/env python
+"""Compute AraSim ice attenuation length vs depth at a few frequencies."""
+
+import numpy as np
+import ROOT
+from ROOT import gSystem
+
+# change this to your use case!
+gSystem.Load('/home/baclark/scratch/ARA/AraSim/libAra.so')
+
+icemodel = ROOT.IceModel(10, 0, 0)
+settings = ROOT.Settings()
+
+depths_m = np.linspace(1.0, 3000.0, 301)
+frequencies_mhz = [100.0, 300.0, 700.0]
+
+attens = {f: np.empty_like(depths_m) for f in frequencies_mhz}
+for f_mhz in frequencies_mhz:
+    f_ghz = f_mhz / 1000.0
+    for i, z in enumerate(depths_m):
+        attens[f_mhz][i] = icemodel.GetFreqDepIceAttenuLength(
+            float(z), float(f_ghz), settings
+        )
+
+temperatures_C = np.array([icemodel.temperature(float(z), False) for z in depths_m])
+
+csv_path = "atten_vs_depth.csv"
+header = "depth_m,temperature_C," + ",".join(
+    f"L_{int(f)}MHz_m" for f in frequencies_mhz
+)
+stacked = np.column_stack(
+    [depths_m, temperatures_C] + [attens[f] for f in frequencies_mhz]
+)
+np.savetxt(csv_path, stacked, delimiter=",", header=header, comments="", fmt="%.4f")
+print(f"Wrote {csv_path}")
diff --git a/examples/atten_length/plot_atten.py b/examples/atten_length/plot_atten.py
new file mode 100644
index 00000000..5a96dade
--- /dev/null
+++ b/examples/atten_length/plot_atten.py
@@ -0,0 +1,53 @@
+#!/usr/bin/env python
+"""Plot ice attenuation length and temperature vs depth from CSV."""
+
+import numpy as np
+import pandas as pd
+import matplotlib.pyplot as plt
+
+df = pd.read_csv("atten_vs_depth.csv")
+depths_m = df["depth_m"].values
+temperatures_C = df["temperature_C"].values
+
+freq_cols = [c for c in df.columns if c.startswith("L_") and c.endswith("MHz_m")]
+frequencies_mhz = [float(c.split("_")[1].replace("MHz", "")) for c in freq_cols]
+
+linestyles = ["-", "--", ":", "-."]
+line_color = "tab:blue"
+temp_color = "tab:red"
+
+fig, ax = plt.subplots(figsize=(7.0, 6.5))
+
+for col, f_mhz, ls in zip(freq_cols, frequencies_mhz, linestyles):
+    ax.plot(df[col].values, depths_m,
+            color=line_color, ls=ls, lw=1.8, label=f"{int(f_mhz)} MHz")
+
+ax.set_ylabel("Depth (m)")
+ax.set_xlabel("Attenuation length (m)", color=line_color)
+ax.tick_params(axis="x", colors=line_color)
+ax.spines["bottom"].set_color(line_color)
+ax.invert_yaxis()
+ax.grid(True, alpha=0.3)
+
+south_pole_ice_depth_m = 2850.0
+ax.axhline(south_pole_ice_depth_m, color="k", lw=1.0, ls="-.", alpha=0.7)
+ax.text(ax.get_xlim()[1], south_pole_ice_depth_m,
+        f" South Pole bedrock ({south_pole_ice_depth_m:.0f} m)",
+        color="k", fontsize=8, va="bottom", ha="right")
+
+ax_T = ax.twiny()
+ax_T.plot(temperatures_C, depths_m, color=temp_color, lw=1.8, label="Temperature")
+ax_T.set_xlabel("Temperature (\N{DEGREE SIGN}C)", color=temp_color)
+ax_T.tick_params(axis="x", colors=temp_color)
+ax_T.spines["top"].set_color(temp_color)
+ax_T.spines["bottom"].set_visible(False)
+
+lines_L, labels_L = ax.get_legend_handles_labels()
+lines_T, labels_T = ax_T.get_legend_handles_labels()
+ax.legend(lines_L + lines_T, labels_L + labels_T, loc="center right")
+
+ax.set_title("AraSim ice attenuation length and temperature vs depth", pad=28)
+
+fig.tight_layout()
+fig.savefig("atten_vs_depth.png", dpi=150)
+print("Wrote atten_vs_depth.png")
diff --git a/log.txt b/log.txt
index 789cfa0d..6a756a92 100644
--- a/log.txt
+++ b/log.txt
@@ -1913,3 +1913,31 @@ Removed ReadVgain, ReadVgainTop, ReadHgain, and ReadTxgain functions and replace
 Added new option ELECTRONICS_ANTENNA_CONSISTENCY and corresponding functionality to ensure antenna gain model and electronics gain model are self-consistent within a simulation. This is done by calculating the electronics gain model amplitude on-the-fly from the in situ noise model, antenna gain model, and amplifier noise figure. The default (ELECTRONICS_ANTENNA_CONSISTENCY=1) is to ensure these are self-consistent whenever an in-situ electronics gain model is requested and a corresponding in-situ noise model is used. The electronics gain phase read-in from the in situ model file is left unchanged. 
 
 ==============================================================================
+
+2026/5/13 Brian Clark
+
+Hardened ice attenuation models against unphysical behavior at large depths.
+Per D. Besson, the attenuation model is unreliable beyond ~2.5 km depth.
+
+IceModel: Added public static constexpr BEDROCK_DEPTH_M (2850 m) 
+and MIN_PHYSICAL_ATTEN_M (1 m).
+These set the depth of bedrock at SP and the minimum attenuation length
+we think is physical. This is important to help protect against nans
+(see below in the report class).
+GetFreqDepIceAttenuLength, GetARAIceAttenuLength, and 
+EffectiveAttenuationLength now return 0 past 2.5 km and past bedrock. 
+temperature() gained an opt-in strict mode (throws on out-of-ice depths) 
+and clamps positive values to 0.
+I also removed an unused version of EffectiveAttenuationLength,
+which must have been cruft.
+
+Report: All uses of the attenuation length in e.g. exp(-dl/L)
+now guard against L < MIN_PHYSICAL_ATTEN_M, 
+treating such steps as fully attenuating to 
+avoid NaNs from division by zero.
+
+NB: I've hardcoded in the depth of bedrock, which means
+if someone later wanted to run this at a site with deeper ice,
+they'll need to adjust it.
+
+==============================================================================
\ No newline at end of file
diff --git a/test/veff_ara3/run_veff_ara3_test.sh b/test/veff_ara3/run_veff_ara3_test.sh
index 19225bbd..d95e0af1 100644
--- a/test/veff_ara3/run_veff_ara3_test.sh
+++ b/test/veff_ara3/run_veff_ara3_test.sh
@@ -5,6 +5,6 @@
 # then, do a comparison to check for consistency
 
 EXPECTED_GLOBAL_PASS=3
-EXPECTED_TOTAL_WEIGHT=1.8148
+EXPECTED_TOTAL_WEIGHT=0.9879
 EXPECTED_TOTAL_WEIGHT_SIGMA=0.0001
 python3 test/check_sim.py test/veff_ara3/veff_ara3_test_output.txt $EXPECTED_GLOBAL_PASS $EXPECTED_TOTAL_WEIGHT $EXPECTED_TOTAL_WEIGHT_SIGMA
diff --git a/test/veff_ara3/veff_ara3_test_output.txt b/test/veff_ara3/veff_ara3_test_output.txt
index ce2cc04d..2ecb0512 100644
--- a/test/veff_ara3/veff_ara3_test_output.txt
+++ b/test/veff_ara3/veff_ara3_test_output.txt
@@ -1,4 +1,4 @@
-The Git Commit Hash: db07246e2a18a11943b09f4afd46b505a92583dc
+The Git Commit Hash: b7a845875883896465a00bfbc7c7c0ee83371f9e
 
 	Default values!
 NNU : 100
@@ -30,6 +30,7 @@ DETECTOR : 4
 POSNU_RADIUS : 5000
 EVENT_GENERATION_MODE: 0
 first random : 0.553
+Loaded ice attenuation percent table from /home/baclark/scratch/ARA/AraSim_old//data/iceattenuation_systematics.txt with 3001 rows.
 
 Surface at (log:0, lat:0) : 6.35973e+06
 SurfaceAboveGeoid at (log:0, lat:0) : 2977
@@ -38,7 +39,7 @@ read stations_per_side
 	Error, station number not match !
 InstalledStations[StationID].nChannels is 20
 Opening 2013-2017 SQliteDB using AraGeomTool::getStationInfo()
-Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /home/mmuzio/AraRoot/build/share/araCalib/AntennaInfo.sqlite
+Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /cvmfs/ara.opensciencegrid.org/trunk/RHEL_8_x86_64/ara_build/share/araCalib/AntennaInfo.sqlite
 Opening default 2013-2017 SQLite tables for all stations 
 Borehole ch: 0 inserted station: 0 station: 3 string: 3 ant: 2 Type: 0
 Borehole ch: 1 inserted station: 0 station: 3 string: 0 ant: 2 Type: 0
@@ -95,7 +96,7 @@ Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 0 : 10004.8 : 9990.5
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 1 : 10004.8 : 9990.58 : -188.408 : 14140.1 : 1.58412 : 0.784686 : 6.37813e+06 : 
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 2 : 10004.8 : 9990.58 : -172.607 : 14139.9 : 1.583 : 0.784686 : 6.37813e+06 : 
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 3 : 10004.8 : 9990.58 : -169.486 : 14139.9 : 1.58278 : 0.784686 : 6.37813e+06 : 
-The number of channels: 16
+Index of refraction for antenna gain file /home/baclark/scratch/ARA/AraSim_old//data/antennas/realizedGain/PVA_RealizedGainAndPhase_Copol_Kansas2024.txt could not be determined. Index will adjust to always match local value.
 DATA_BIN_SIZE: 16384 TIMESTEP: 6.25e-10 df_fft: 97656.2
 number of f bins: 8192
 System temp ch0 : 300
@@ -116,14 +117,14 @@ System temp ch14 : 325
 System temp ch15 : 325
 System temp ch16 : 0
 Reading in situ ARA noise for this station and configuration from file: 
-/home/mmuzio/AraSimMain/AraSim/data/noise/sigmavsfreq_A_3_config_2.csv
+/home/baclark/scratch/ARA/AraSim_old//data/noise/sigmavsfreq_A_3_config_2.csv
 start read elect chain
  Reading standard ARA electronics response from file:
-/home/mmuzio/AraSimMain/AraSim/data/gain/ARA_Electronics_TotalGain_TwoFilters.csv
+/home/baclark/scratch/ARA/AraSim_old//data/gain/ARA_Electronics_TotalGain_TwoFilters.csv
  Recalculating in situ electronic response amplitude to ensure consistency with antenna model used here.
 done read elect chain
  Reading trigger formation values for this station and configuration from file:
-/home/mmuzio/AraSimMain/AraSim/data/trigger/delays_masking_A3_C2.csv
+/home/baclark/scratch/ARA/AraSim_old//data/trigger/delays_masking_A3_C2.csv
 Sharesurface: 0 : 0 : 0 : 10010.2 : 10003 : 6.35944e+06 : 
 Sharesurface: 0 : 0 : 1 : 10010.2 : 10003 : 6.35944e+06 : 
 Sharesurface: 0 : 0 : 2 : 10010.2 : 10003 : 6.35946e+06 : 
@@ -211,14 +212,12 @@ station[0].strings[3].antennas[1] no_ch:14
 station[0].strings[3].antennas[2] no_ch:15
 station[0].strings[3].antennas[3] no_ch:16
 *Thrown 0
-**********
-trigger passed at bin 5380  passed ch : 3 (1type) Direct dist btw posnu : 1823.63 noiseID : 0 ViewAngle : 55.1341 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 5413  passed ch : 5 (1type) Direct dist btw posnu : 1811.83 noiseID : 0 ViewAngle : 55.7054 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 5314  passed ch : 6 (0type) Direct dist btw posnu : 1817.96 noiseID : 0 ViewAngle : 55.5114 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 5434  passed ch : 10 (0type) Direct dist btw posnu : 1831.08 noiseID : 0 ViewAngle : 55.5173 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 5416  passed ch : 11 (1type) Direct dist btw posnu : 1832.25 noiseID : 0 ViewAngle : 55.4806 LikelyTrigSignal : interaction 0, ray 1
-Global_Pass : 5380 evt : 10 added weight : 0.81612
-0.81612 : 9
+**************************
+trigger passed at bin 4374  passed ch : 6 (0type) Direct dist btw posnu : 1945 noiseID : 0 ViewAngle : 52.6573 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4246  passed ch : 12 (0type) Direct dist btw posnu : 1930.14 noiseID : 0 ViewAngle : 52.3707 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4359  passed ch : 14 (0type) Direct dist btw posnu : 1943.87 noiseID : 0 ViewAngle : 52.7658 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4374 evt : 26 added weight : 0.987877
+0.987877 : 9
 
 Making useful event
 StationID: 3
@@ -230,44 +229,53 @@ OR
 c) Call UsefulAtriStationEvent *realAtriEvPtr = new UsefulAtriStationEvent(rawAtriEvPtr, AraCalType::kLatestCalib); ***BEFORE*** 
 AraGeomTool::getStationInfo(3,2017)
 Opening 2013-2017 SQliteDB using AraGeomTool::getStationInfo()
-Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /home/mmuzio/AraRoot/build/share/araCalib/AntennaInfo.sqlite
+Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /cvmfs/ara.opensciencegrid.org/trunk/RHEL_8_x86_64/ara_build/share/araCalib/AntennaInfo.sqlite
 Opening default 2013-2017 SQLite tables for all stations 
-******
-trigger passed at bin 4387  passed ch : 0 (0type) Direct dist btw posnu : 3201.86 noiseID : 0 ViewAngle : 56.3795 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 4468  passed ch : 6 (0type) Direct dist btw posnu : 3218.74 noiseID : 0 ViewAngle : 56.9636 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 4590  passed ch : 10 (0type) Direct dist btw posnu : 3219.54 noiseID : 0 ViewAngle : 56.9094 LikelyTrigSignal : interaction 0, ray 1
-Global_Pass : 4387 evt : 16 added weight : 4.29438e-27
-4.29438e-27 : 0
+****************************************************************
+trigger passed at bin 4215  passed ch : 5 (1type) Direct dist btw posnu : 1806.47 noiseID : 0 ViewAngle : 54.3836 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4294  passed ch : 7 (1type) Direct dist btw posnu : 1814.86 noiseID : 0 ViewAngle : 54.2689 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4168  passed ch : 9 (1type) Direct dist btw posnu : 1801.01 noiseID : 0 ViewAngle : 53.9506 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4245  passed ch : 11 (1type) Direct dist btw posnu : 1809.35 noiseID : 0 ViewAngle : 53.8391 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4273  passed ch : 13 (1type) Direct dist btw posnu : 1812.32 noiseID : 0 ViewAngle : 53.7993 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4355  passed ch : 15 (1type) Direct dist btw posnu : 1820.63 noiseID : 0 ViewAngle : 53.6875 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4215 evt : 90 added weight : 1.19837e-07
+1.19837e-07 : 0
 
 Making useful event
 StationID: 3
 StationID_AraRoot: 3
-************************************************************************************Thrown 100
-***********************************************************************************
-trigger passed at bin 4704  passed ch : 4 (0type) Direct dist btw posnu : 1911.64 noiseID : 0 ViewAngle : 51.8715 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 4651  passed ch : 6 (0type) Direct dist btw posnu : 1915.05 noiseID : 0 ViewAngle : 51.2907 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 4609  passed ch : 12 (0type) Direct dist btw posnu : 1899.4 noiseID : 0 ViewAngle : 51.7349 LikelyTrigSignal : interaction 0, ray 1
-Global_Pass : 4704 evt : 183 added weight : 0.998722
-0.998722 : 9
+**********Thrown 100
+****************************************************************************************************Thrown 200
+****************************************************************************************************Thrown 300
+**********************************************************************************************
+trigger passed at bin 4565  passed ch : 2 (0type) Direct dist btw posnu : 1797.57 noiseID : 0 ViewAngle : 57.8673 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4557  passed ch : 3 (1type) Direct dist btw posnu : 1798.2 noiseID : 0 ViewAngle : 57.9643 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4691  passed ch : 5 (1type) Direct dist btw posnu : 1803.56 noiseID : 0 ViewAngle : 57.3868 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4634  passed ch : 6 (0type) Direct dist btw posnu : 1806.8 noiseID : 0 ViewAngle : 57.8891 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4727  passed ch : 10 (0type) Direct dist btw posnu : 1818.1 noiseID : 0 ViewAngle : 58.2522 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4720  passed ch : 11 (1type) Direct dist btw posnu : 1818.72 noiseID : 0 ViewAngle : 58.3448 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4694  passed ch : 13 (1type) Direct dist btw posnu : 1805.52 noiseID : 0 ViewAngle : 57.761 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4645  passed ch : 14 (0type) Direct dist btw posnu : 1808.8 noiseID : 0 ViewAngle : 58.271 LikelyTrigSignal : interaction 0, ray 1
+Global_Pass : 4565 evt : 394 added weight : 6.95039e-29
+6.95039e-29 : 0
 
 Making useful event
 StationID: 3
 StationID_AraRoot: 3
-*****************Thrown 200
-****************************************************************************************************Thrown 300
-****************************************************************************************************Thrown 400
+******Thrown 400
 *************************************************************************************************** end loop
 Total Events Thrown: 500
 Total_Global_Pass : 3
-Total_Weight : 1.81484
-Total_Probability : -1619.74
-weight bin values : 1, 0, 0, 0, 0, 0, 0, 0, 0, 2
+Total_Weight : 0.987878
+Total_Probability : 3614.66
+weight bin values : 2, 0, 0, 0, 0, 0, 0, 0, 0, 1
 
 Radius: 5000 [m]
 IceVolume : 2.35619e+11
-test Veff(ice) : 1.07471e+10 m3sr, 10.7471 km3sr
-test Veff(water eq.) : 9.85506e+09 m3sr, 9.85506 km3sr
-And Veff(water eq.) error plus : 9.11416 km3sr and error minus : 4.84385 km3sr
+test Veff(ice) : 5.84998e+09 m3sr, 5.84998 km3sr
+test Veff(water eq.) : 5.36443e+09 m3sr, 5.36443 km3sr
+And Veff(water eq.) error plus : 7.31507 km3sr and error minus : 2.42195 km3sr
+Warning in <THistPainter::MakeChopt>: option SCAT is deprecated.
 Info in <TCanvas::Print>: pdf file sigmaCrossSection.pdf has been created
 outputdir= outputs/.
 This AraSim run is complete and will exit with code 0
diff --git a/test/veff_birefringence/run_veff_birefringence_test.sh b/test/veff_birefringence/run_veff_birefringence_test.sh
index 7d39393d..d541ce20 100644
--- a/test/veff_birefringence/run_veff_birefringence_test.sh
+++ b/test/veff_birefringence/run_veff_birefringence_test.sh
@@ -4,7 +4,7 @@
 
 # then, do a comparison to check for consistency
 
-EXPECTED_GLOBAL_PASS=6
-EXPECTED_TOTAL_WEIGHT=2.89898
+EXPECTED_GLOBAL_PASS=3
+EXPECTED_TOTAL_WEIGHT=0.999202
 EXPECTED_TOTAL_WEIGHT_SIGMA=0.0001
 python3 test/check_sim.py test/veff_birefringence/veff_birefringence_test_output.txt $EXPECTED_GLOBAL_PASS $EXPECTED_TOTAL_WEIGHT $EXPECTED_TOTAL_WEIGHT_SIGMA
diff --git a/test/veff_birefringence/veff_birefringence_test_output.txt b/test/veff_birefringence/veff_birefringence_test_output.txt
index c79e845e..2381c845 100644
--- a/test/veff_birefringence/veff_birefringence_test_output.txt
+++ b/test/veff_birefringence/veff_birefringence_test_output.txt
@@ -1,4 +1,4 @@
-The Git Commit Hash: 30328d5436ef38ee5ee2ee10ed82ae595675a984
+The Git Commit Hash: b7a845875883896465a00bfbc7c7c0ee83371f9e
 
 	Default values!
 NNU : 100
@@ -30,6 +30,7 @@ DETECTOR : 4
 POSNU_RADIUS : 5000
 EVENT_GENERATION_MODE: 0
 first random : 0.994531
+Loaded ice attenuation percent table from /home/baclark/scratch/ARA/AraSim_old//data/iceattenuation_systematics.txt with 3001 rows.
 
 Surface at (log:0, lat:0) : 6.35973e+06
 SurfaceAboveGeoid at (log:0, lat:0) : 2977
@@ -38,7 +39,7 @@ read stations_per_side
 	Error, station number not match !
 InstalledStations[StationID].nChannels is 20
 Opening 2013-2017 SQliteDB using AraGeomTool::getStationInfo()
-Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /home/abishop/analysis/AraRoot_Install//share/araCalib/AntennaInfo.sqlite
+Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /cvmfs/ara.opensciencegrid.org/trunk/RHEL_8_x86_64/ara_build/share/araCalib/AntennaInfo.sqlite
 Opening default 2013-2017 SQLite tables for all stations 
 Borehole ch: 0 inserted station: 0 station: 2 string: 3 ant: 2 Type: 0
 Borehole ch: 1 inserted station: 0 station: 2 string: 0 ant: 2 Type: 0
@@ -95,7 +96,7 @@ Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 0 : 10010.6 : 10002.
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 1 : 10010.6 : 10002.3 : -186.546 : 14152.5 : 1.58398 : 0.784986 : 6.37813e+06 : 
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 2 : 10010.6 : 10002.3 : -170.247 : 14152.3 : 1.58283 : 0.784986 : 6.37813e+06 : 
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 3 : 10010.6 : 10002.3 : -167.492 : 14152.2 : 1.58263 : 0.784985 : 6.37813e+06 : 
-The number of channels: 16
+Index of refraction for antenna gain file /home/baclark/scratch/ARA/AraSim_old//data/antennas/realizedGain/PVA_RealizedGainAndPhase_Copol_Kansas2024.txt could not be determined. Index will adjust to always match local value.
 DATA_BIN_SIZE: 16384 TIMESTEP: 6.25e-10 df_fft: 97656.2
 number of f bins: 8192
 System temp ch0 : 300
@@ -116,14 +117,14 @@ System temp ch14 : 325
 System temp ch15 : 325
 System temp ch16 : 0
 Reading in situ ARA noise for this station and configuration from file: 
-/home/abishop/AraSim-AraSoft//data/noise/sigmavsfreq_A_2_config_4.csv
+/home/baclark/scratch/ARA/AraSim_old//data/noise/sigmavsfreq_A_2_config_4.csv
 start read elect chain
  Reading standard ARA electronics response from file:
-/home/abishop/AraSim-AraSoft//data/gain/ARA_Electronics_TotalGain_TwoFilters.csv
+/home/baclark/scratch/ARA/AraSim_old//data/gain/ARA_Electronics_TotalGain_TwoFilters.csv
  Recalculating in situ electronic response amplitude to ensure consistency with antenna model used here.
 done read elect chain
  Reading trigger formation values for this station and configuration from file:
-/home/abishop/AraSim-AraSoft//data/trigger/delays_masking_A2_C4.csv
+/home/baclark/scratch/ARA/AraSim_old//data/trigger/delays_masking_A2_C4.csv
 Sharesurface: 0 : 0 : 0 : 10004.4 : 9989.17 : 6.35944e+06 : 
 Sharesurface: 0 : 0 : 1 : 10004.4 : 9989.17 : 6.35945e+06 : 
 Sharesurface: 0 : 0 : 2 : 10004.4 : 9989.2 : 6.35946e+06 : 
@@ -211,15 +212,13 @@ station[0].strings[3].antennas[1] no_ch:14
 station[0].strings[3].antennas[2] no_ch:15
 station[0].strings[3].antennas[3] no_ch:16
 *Thrown 0
-*************************************************short chord 0.166532
-***************************************************Thrown 100
-****************************************************************************************************Thrown 200
-*************************************
-trigger passed at bin 4228  passed ch : 0 (0type) Direct dist btw posnu : 1591.57 noiseID : 0 ViewAngle : 55.4854 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4257  passed ch : 2 (0type) Direct dist btw posnu : 1595.82 noiseID : 0 ViewAngle : 56.1332 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4241  passed ch : 12 (0type) Direct dist btw posnu : 1592.92 noiseID : 0 ViewAngle : 55.3397 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4257 evt : 237 added weight : 0.999199
-0.999199 : 9
+****************************************************************************************************Thrown 100
+
+trigger passed at bin 4267  passed ch : 1 (1type) Direct dist btw posnu : 1916.69 noiseID : 0 ViewAngle : 56.8507 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4301  passed ch : 3 (1type) Direct dist btw posnu : 1920.9 noiseID : 0 ViewAngle : 56.7298 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4251  passed ch : 5 (1type) Direct dist btw posnu : 1914.85 noiseID : 0 ViewAngle : 57.4716 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4301 evt : 100 added weight : 7.18397e-06
+7.18397e-06 : 0
 
 Making useful event
 StationID: 2
@@ -231,68 +230,46 @@ OR
 c) Call UsefulAtriStationEvent *realAtriEvPtr = new UsefulAtriStationEvent(rawAtriEvPtr, AraCalType::kLatestCalib); ***BEFORE*** 
 AraGeomTool::getStationInfo(3,2017)
 Opening 2013-2017 SQliteDB using AraGeomTool::getStationInfo()
-Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /home/abishop/analysis/AraRoot_Install//share/araCalib/AntennaInfo.sqlite
+Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /cvmfs/ara.opensciencegrid.org/trunk/RHEL_8_x86_64/ara_build/share/araCalib/AntennaInfo.sqlite
 Opening default 2013-2017 SQLite tables for all stations 
-******************************************************
-trigger passed at bin 4213  passed ch : 0 (0type) Direct dist btw posnu : 1675.84 noiseID : 0 ViewAngle : 54.4786 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4222  passed ch : 2 (0type) Direct dist btw posnu : 1678.91 noiseID : 0 ViewAngle : 53.8541 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4290  passed ch : 12 (0type) Direct dist btw posnu : 1685.29 noiseID : 0 ViewAngle : 54.4365 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4290 evt : 291 added weight : 1.21462e-26
-1.21462e-26 : 0
-
-Making useful event
-StationID: 2
-StationID_AraRoot: 2
-*********Thrown 300
-*****
-trigger passed at bin 4300  passed ch : 4 (0type) Direct dist btw posnu : 1500.3 noiseID : 0 ViewAngle : 59.7475 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4233  passed ch : 8 (0type) Direct dist btw posnu : 1493.54 noiseID : 0 ViewAngle : 59.2655 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4308  passed ch : 10 (0type) Direct dist btw posnu : 1501.89 noiseID : 0 ViewAngle : 59.452 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4308 evt : 305 added weight : 0.900189
-0.900189 : 9
+****************************************************************************************************Thrown 200
+***************************************
+trigger passed at bin 4214  passed ch : 4 (0type) Direct dist btw posnu : 1954.96 noiseID : 0 ViewAngle : 58.1995 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4202  passed ch : 5 (1type) Direct dist btw posnu : 1955.6 noiseID : 0 ViewAngle : 58.1527 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4231  passed ch : 7 (1type) Direct dist btw posnu : 1958.94 noiseID : 0 ViewAngle : 57.914 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4166  passed ch : 12 (0type) Direct dist btw posnu : 1949.68 noiseID : 0 ViewAngle : 58.5309 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4180  passed ch : 15 (1type) Direct dist btw posnu : 1953.56 noiseID : 0 ViewAngle : 58.2548 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4231 evt : 239 added weight : 5.31138e-08
+5.31138e-08 : 0
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
-***********************************************************************************************Thrown 400
-*******************************
-trigger passed at bin 4397  passed ch : 4 (0type) Direct dist btw posnu : 1791.8 noiseID : 0 ViewAngle : 57.8268 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 4341  passed ch : 6 (0type) Direct dist btw posnu : 1792.52 noiseID : 0 ViewAngle : 56.9911 LikelyTrigSignal : interaction 0, ray 1
-trigger passed at bin 4415  passed ch : 10 (0type) Direct dist btw posnu : 1801.08 noiseID : 0 ViewAngle : 57.0021 LikelyTrigSignal : interaction 0, ray 1
-Global_Pass : 4415 evt : 431 added weight : 0.999593
-0.999593 : 9
+************
+trigger passed at bin 4666  passed ch : 3 (1type) Direct dist btw posnu : 1734.72 noiseID : 0 ViewAngle : 56.08 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4687  passed ch : 11 (1type) Direct dist btw posnu : 1734.33 noiseID : 0 ViewAngle : 56.4274 LikelyTrigSignal : interaction 0, ray 1
+trigger passed at bin 4634  passed ch : 13 (1type) Direct dist btw posnu : 1745.22 noiseID : 0 ViewAngle : 46.2013 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4687 evt : 251 added weight : 0.999194
+0.999194 : 9
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
-**************************************depth negative! -22.2321
-depth negative! -22.2321
-depth negative! -22.2321
-depth negative! -22.2321
-depth negative! -22.218
-depth negative! -22.218
-depth negative! -22.218
-depth negative! -22.218
-depth negative! -22.2288
-depth negative! -22.2288
-depth negative! -22.2288
-depth negative! -22.2288
-depth negative! -22.2217
-depth negative! -22.2217
-depth negative! -22.2217
-depth negative! -22.2217
-****************************** end loop
+*************************************************Thrown 300
+****************************************************************************************************Thrown 400
+*************************************************************************************************** end loop
 Total Events Thrown: 500
-Total_Global_Pass : 4
-Total_Weight : 2.89898
-Total_Probability : 4542.12
-weight bin values : 1, 0, 0, 0, 0, 0, 0, 0, 0, 3
+Total_Global_Pass : 3
+Total_Weight : 0.999202
+Total_Probability : 4186.5
+weight bin values : 2, 0, 0, 0, 0, 0, 0, 0, 0, 1
 
 Radius: 5000 [m]
 IceVolume : 2.35619e+11
-test Veff(ice) : 1.71671e+10 m3sr, 17.1671 km3sr
-test Veff(water eq.) : 1.57422e+10 m3sr, 15.7422 km3sr
-And Veff(water eq.) error plus : 9.29678 km3sr and error minus : 7.30422 km3sr
+test Veff(ice) : 5.91703e+09 m3sr, 5.91703 km3sr
+test Veff(water eq.) : 5.42592e+09 m3sr, 5.42592 km3sr
+And Veff(water eq.) error plus : 7.31507 km3sr and error minus : 2.42195 km3sr
+Warning in <THistPainter::MakeChopt>: option SCAT is deprecated.
 Info in <TCanvas::Print>: pdf file sigmaCrossSection.pdf has been created
 outputdir= outputs/.
 This AraSim run is complete and will exit with code 0
diff --git a/test/veff_pa/run_veff_pa_test.sh b/test/veff_pa/run_veff_pa_test.sh
index 46d2212c..324b7231 100644
--- a/test/veff_pa/run_veff_pa_test.sh
+++ b/test/veff_pa/run_veff_pa_test.sh
@@ -4,7 +4,7 @@
 
 # then, do a comparison to check for consistency
 
-EXPECTED_GLOBAL_PASS=9
-EXPECTED_TOTAL_WEIGHT=4.96778
+EXPECTED_GLOBAL_PASS=7
+EXPECTED_TOTAL_WEIGHT=2.98723
 EXPECTED_TOTAL_WEIGHT_SIGMA=0.0001
 python3 test/check_sim.py test/veff_pa/veff_pa_test_output.txt $EXPECTED_GLOBAL_PASS $EXPECTED_TOTAL_WEIGHT $EXPECTED_TOTAL_WEIGHT_SIGMA
diff --git a/test/veff_pa/veff_pa_test_output.txt b/test/veff_pa/veff_pa_test_output.txt
index 98da3a04..0626bbf0 100644
--- a/test/veff_pa/veff_pa_test_output.txt
+++ b/test/veff_pa/veff_pa_test_output.txt
@@ -1,4 +1,4 @@
-The Git Commit Hash: 30328d5436ef38ee5ee2ee10ed82ae595675a984
+The Git Commit Hash: b7a845875883896465a00bfbc7c7c0ee83371f9e
 
 	Default values!
 NNU : 100
@@ -28,12 +28,13 @@ DETECTOR : 5
 POSNU_RADIUS : 5000
 EVENT_GENERATION_MODE: 0
 first random : 0.455064
+Loaded ice attenuation percent table from /home/baclark/scratch/ARA/AraSim_old//data/iceattenuation_systematics.txt with 3001 rows.
 
 Surface at (log:0, lat:0) : 6.35973e+06
 SurfaceAboveGeoid at (log:0, lat:0) : 2977
 Simulating realistic ARA05 and Phased Array.
 Using 7 ARA05 vanilla Vpols
-The number of channels: 16
+Index of refraction for antenna gain file /home/baclark/scratch/ARA/AraSim_old//data/antennas/realizedGain/PVA_RealizedGainAndPhase_Copol_Kansas2024.txt could not be determined. Index will adjust to always match local value.
 DATA_BIN_SIZE: 16384 TIMESTEP: 6.67e-10 df_fft: 91507
 number of f bins: 8192
 ThresOffset ch0 : 1
@@ -86,10 +87,10 @@ System temp ch14 : 325
 System temp ch15 : 325
 System temp ch16 : 0
 Reading in situ PA noise for this station and configuration from file: 
-/home/abishop/AraSim-AraSoft//data/noise/sigmavsfreq_PA_config_4.csv
+/home/baclark/scratch/ARA/AraSim_old//data/noise/sigmavsfreq_PA_config_4.csv
 start read elect chain
  Reading standard PA electronics response from file:
-/home/abishop/AraSim-AraSoft//data/gain/PA_Electronics_TotalGainPhase.csv
+/home/baclark/scratch/ARA/AraSim_old//data/gain/PA_Electronics_TotalGainPhase.csv
  Recalculating in situ electronic response amplitude to ensure consistency with antenna model used here.
 done read elect chain
 Sharesurface: 0 : 0 : 0 : 9999.57 : 9999.57 : 6.35945e+06 : 
@@ -182,11 +183,13 @@ station[0].strings[3].antennas[1] no_ch:14
 station[0].strings[4].antennas[0] no_ch:15
 station[0].strings[4].antennas[1] no_ch:16
 *Thrown 0
-************************************************************
-PA trigger ~~~ raySolNum: 0  avgSNR: 2.53092  Event Number : 60  PA efficiency : 0.930428
+********************************************************** efficiency : 1.00161  avg SNR : 4.4043
+ efficiency : 1.00161  avg SNR : 4.33342
 
-Global_Pass : 8609 evt : 60 added weight : 7.84046e-81
-7.84046e-81 : 0
+PA trigger ~~~  Event Number : 58
+
+Global_Pass : 9094 evt : 58 added weight : 1.46284e-06
+1.46284e-06 : 0
 
 Making useful event
 StationID: 6
@@ -198,121 +201,93 @@ OR
 c) Call UsefulAtriStationEvent *realAtriEvPtr = new UsefulAtriStationEvent(rawAtriEvPtr, AraCalType::kLatestCalib); ***BEFORE*** 
 AraGeomTool::getStationInfo(3,2017)
 Opening 2013-2017 SQliteDB using AraGeomTool::getStationInfo()
-Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /home/abishop/analysis/AraRoot_Install//share/araCalib/AntennaInfo.sqlite
+Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /cvmfs/ara.opensciencegrid.org/trunk/RHEL_8_x86_64/ara_build/share/araCalib/AntennaInfo.sqlite
 Opening default 2013-2017 SQLite tables for all stations 
-****************************************Thrown 100
-*********************************
-PA trigger ~~~ raySolNum: 0  avgSNR: 62.1325  Event Number : 133  PA efficiency : 1.00162
-
-Global_Pass : 8555 evt : 133 added weight : 3.35065e-29
-3.35065e-29 : 0
+******************************************Thrown 100
+****************************************** efficiency : 1.00162  avg SNR : 22.4232
+ efficiency : 1.00149  avg SNR : 3.80011
 
-Making useful event
-StationID: 6
-StationID_AraRoot: 3
-*******************************************************************Thrown 200
-*****
-PA trigger ~~~ raySolNum: 0  avgSNR: 8.62804  Event Number : 205  PA efficiency : 1.00162
+PA trigger ~~~  Event Number : 142
 
-Global_Pass : 8474 evt : 205 added weight : 0.995622
-0.995622 : 9
+Global_Pass : 9881 evt : 142 added weight : 3.4195e-80
+3.4195e-80 : 0
 
 Making useful event
 StationID: 6
 StationID_AraRoot: 3
-*******************************
-PA trigger ~~~ raySolNum: 0  avgSNR: 4.65238  Event Number : 236  PA efficiency : 1.00162
-
-Global_Pass : 8446 evt : 236 added weight : 0.998235
-0.998235 : 9
+******************************************************** efficiency : 1.00162  avg SNR : 16.7165
+ efficiency : 1.00162  avg SNR : 4.9953
 
-Making useful event
-StationID: 6
-StationID_AraRoot: 3
-*****************************
-PA trigger ~~~ raySolNum: 0  avgSNR: 9.40231  Event Number : 265  PA efficiency : 1.00162
+PA trigger ~~~  Event Number : 198
 
-Global_Pass : 8545 evt : 265 added weight : 2.20756e-43
-2.20756e-43 : 0
+Global_Pass : 8508 evt : 198 added weight : 0.997365
+0.997365 : 9
 
 Making useful event
 StationID: 6
 StationID_AraRoot: 3
-***********************************Thrown 300
-********************************************
-PA trigger ~~~ raySolNum: 0  avgSNR: 2.38054  Event Number : 344  PA efficiency : 0.863027
-
-Global_Pass : 8474 evt : 344 added weight : 0.995615
-0.995615 : 9
+**Thrown 200
+****************************************************************************************************Thrown 300
+* efficiency : 1.00162  avg SNR : 5.5357
+ efficiency : 1.00161  avg SNR : 4.32862
 
-Making useful event
-StationID: 6
-StationID_AraRoot: 3
-*************************************depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.06175
-depth negative! -7.05293
-depth negative! -7.05293
-depth negative! -7.07004
-depth negative! -7.07155
-depth negative! -7.07155
-depth negative! -7.05425
-depth negative! -7.05425
-*******************Thrown 400
-***************
-PA trigger ~~~ raySolNum: 0  avgSNR: 32.5518  Event Number : 415  PA efficiency : 1.00162
+PA trigger ~~~  Event Number : 301
 
-Global_Pass : 8616 evt : 415 added weight : 0.0243914
-0.0243914 : 4
+Global_Pass : 8377 evt : 301 added weight : 1.17701e-50
+1.17701e-50 : 0
 
 Making useful event
 StationID: 6
 StationID_AraRoot: 3
-******
-PA trigger ~~~ raySolNum: 0  avgSNR: 3.53149  Event Number : 421  PA efficiency : 1.00109
+*************************** efficiency : 0.958449  avg SNR : 2.63821  random number : 0.377613
+******************************** efficiency : 1.00162  avg SNR : 20.8026
+ efficiency : 1.00118  avg SNR : 3.56964
+
+PA trigger ~~~  Event Number : 360
 
-Global_Pass : 8400 evt : 421 added weight : 0.994757
-0.994757 : 9
+Global_Pass : 8528 evt : 360 added weight : 4.09918e-20
+4.09918e-20 : 0
 
 Making useful event
 StationID: 6
 StationID_AraRoot: 3
-***************************************************
-PA trigger ~~~ raySolNum: 0  avgSNR: 4.51783  Event Number : 472  PA efficiency : 1.00162
+************* efficiency : 1.00162  avg SNR : 21.6302
+ efficiency : 1.0016  avg SNR : 4.22
 
-Global_Pass : 8625 evt : 472 added weight : 0.992111
-0.992111 : 9
+PA trigger ~~~  Event Number : 373
+
+Global_Pass : 8361 evt : 373 added weight : 0.996485
+0.996485 : 9
 
 Making useful event
 StationID: 6
 StationID_AraRoot: 3
-*****
-PA trigger ~~~ raySolNum: 0  avgSNR: 7.39083  Event Number : 477  PA efficiency : 1.00162
+***************************Thrown 400
+********** efficiency : 1.00162  avg SNR : 24.0975
+ efficiency : 1.00157  avg SNR : 3.98589
+
+PA trigger ~~~  Event Number : 410
 
-Global_Pass : 8436 evt : 477 added weight : 7.54711e-75
-7.54711e-75 : 0
+Global_Pass : 8486 evt : 410 added weight : 0.993382
+0.993382 : 9
 
 Making useful event
 StationID: 6
 StationID_AraRoot: 3
-********************** end loop
+************************************ efficiency : 0.783637  avg SNR : 2.27016  random number : 0.464669
+***************************************************** end loop
 Total Events Thrown: 500
-Total_Global_Pass : 10
-Total_Weight : 5.00073
-Total_Probability : 4.8271e-315
-weight bin values : 4, 0, 0, 0, 1, 0, 0, 0, 0, 5
+Total_Global_Pass : 7
+Total_Weight : 2.98723
+Total_Probability : 3689.84
+weight bin values : 4, 0, 0, 0, 0, 0, 0, 0, 0, 3
 
 Radius: 5000 [m]
 IceVolume : 2.35619e+11
-test Veff(ice) : 2.96131e+10 m3sr, 29.6131 km3sr
-test Veff(water eq.) : 2.71552e+10 m3sr, 27.1552 km3sr
-And Veff(water eq.) error plus : 11.1788 km3sr and error minus : 8.65008 km3sr
+test Veff(ice) : 1.76897e+10 m3sr, 17.6897 km3sr
+test Veff(water eq.) : 1.62214e+10 m3sr, 16.2214 km3sr
+And Veff(water eq.) error plus : 9.29679 km3sr and error minus : 7.30423 km3sr
+Warning in <THistPainter::MakeChopt>: option SCAT is deprecated.
 Info in <TCanvas::Print>: pdf file sigmaCrossSection.pdf has been created
 outputdir= outputs/.
 This AraSim run is complete and will exit with code 0
diff --git a/test/veff_pulser/veff_pulser_test_output.txt b/test/veff_pulser/veff_pulser_test_output.txt
index 606a1083..53132383 100644
--- a/test/veff_pulser/veff_pulser_test_output.txt
+++ b/test/veff_pulser/veff_pulser_test_output.txt
@@ -1,3 +1,4 @@
+The Git Commit Hash: b7a845875883896465a00bfbc7c7c0ee83371f9e
 
 	Default values!
 NNU : 100
@@ -29,6 +30,7 @@ DETECTOR : 4
 POSNU_RADIUS : 3000
 EVENT_GENERATION_MODE: 0
 first random : 0.809467
+Loaded ice attenuation percent table from /home/baclark/scratch/ARA/AraSim_old//data/iceattenuation_systematics.txt with 3001 rows.
 
 Surface at (log:0, lat:0) : 6.35973e+06
 SurfaceAboveGeoid at (log:0, lat:0) : 2977
@@ -43,7 +45,7 @@ read z_btw
 	Error, station number not match !
 InstalledStations[StationID].nChannels is 20
 Opening 2013-2017 SQliteDB using AraGeomTool::getStationInfo()
-Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /home/abishop/analysis/AraRoot_Install//share/araCalib/AntennaInfo.sqlite
+Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /cvmfs/ara.opensciencegrid.org/trunk/RHEL_8_x86_64/ara_build/share/araCalib/AntennaInfo.sqlite
 Opening default 2013-2017 SQLite tables for all stations 
 Borehole ch: 0 inserted station: 0 station: 2 string: 3 ant: 2 Type: 0
 Borehole ch: 1 inserted station: 0 station: 2 string: 0 ant: 2 Type: 0
@@ -100,17 +102,17 @@ Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 0 : 10010.6 : 10002.
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 1 : 10010.6 : 10002.3 : -186.546 : 14152.5 : 1.58398 : 0.784986 : 6.37813e+06 : 
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 2 : 10010.6 : 10002.3 : -170.247 : 14152.3 : 1.58283 : 0.784986 : 6.37813e+06 : 
 Detector:station:string:antenna:X:Y:Z:R:Theta:Phi:: 0 : 3 : 3 : 10010.6 : 10002.3 : -167.492 : 14152.2 : 1.58263 : 0.784985 : 6.37813e+06 : 
-The number of channels: 16
+Index of refraction for antenna gain file /home/baclark/scratch/ARA/AraSim_old//data/antennas/realizedGain/PVA_RealizedGainAndPhase_Copol_Kansas2024.txt could not be determined. Index will adjust to always match local value.
 DATA_BIN_SIZE: 16384 TIMESTEP: 5e-10 df_fft: 122070
 number of f bins: 8192
 start read elect chain
  Reading standard ARA electronics response from file:
-/home/abishop/AraSim-AraSoft//data/gain/ARA_Electronics_TotalGain_TwoFilters.csv
+/home/baclark/scratch/ARA/AraSim_old//data/gain/ARA_Electronics_TotalGain_TwoFilters.csv
 WARNING - Antenna model used to calculate in situ gain model may not match that used in this simulation!
 	 To ensure consistency load an in-situ noise model with NOISE = 1.
 done read elect chain
  Reading trigger formation values for this station and configuration from file:
-/home/abishop/AraSim-AraSoft//data/trigger/delays_masking_A2_C6.csv
+/home/baclark/scratch/ARA/AraSim_old//data/trigger/delays_masking_A2_C6.csv
 Sharesurface: 0 : 0 : 0 : 10004.4 : 9989.17 : 6.35944e+06 : 
 Sharesurface: 0 : 0 : 1 : 10004.4 : 9989.17 : 6.35945e+06 : 
 Sharesurface: 0 : 0 : 2 : 10004.4 : 9989.2 : 6.35946e+06 : 
@@ -195,8 +197,8 @@ Calculating noise waveforms for RMS
 0% done
 100% done
 From pure noise waveforms, diode responses
-mean, rms diode are -4.65116e-15 9.28804e-14
-rms voltage is 0.0400492
+mean, rms diode are -4.76064e-15 9.45211e-14
+rms voltage is 0.0405195
  DATA_BIN_SIZE : 16384
 powerthreshold : -6.15
 begin looping events!!
@@ -221,10 +223,10 @@ station[0].strings[3].antennas[2] no_ch:15
 station[0].strings[3].antennas[3] no_ch:16
 *Thrown 0
 
-trigger passed at bin 4181  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 13 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4251  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 15 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4284  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 10 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4284 evt : 0 added weight : 1
+trigger passed at bin 4097  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 3 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4167  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 9 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4198  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 7 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4198 evt : 0 added weight : 1
 
 Making useful event
 StationID: 2
@@ -236,86 +238,85 @@ OR
 c) Call UsefulAtriStationEvent *realAtriEvPtr = new UsefulAtriStationEvent(rawAtriEvPtr, AraCalType::kLatestCalib); ***BEFORE*** 
 AraGeomTool::getStationInfo(3,2017)
 Opening 2013-2017 SQliteDB using AraGeomTool::getStationInfo()
-Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /home/abishop/analysis/AraRoot_Install//share/araCalib/AntennaInfo.sqlite
+Pre 2018: AraStationInfo::readChannelMapDbAtri_2(): INFO - /cvmfs/ara.opensciencegrid.org/trunk/RHEL_8_x86_64/ara_build/share/araCalib/AntennaInfo.sqlite
 Opening default 2013-2017 SQLite tables for all stations 
 *
-trigger passed at bin 4181  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 4 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4256  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 9 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4284  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 8 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4284 evt : 1 added weight : 1
+trigger passed at bin 4093  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 3 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4167  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 15 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4198  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 10 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4198 evt : 1 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4181  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 5 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4255  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 6 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4284  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 3 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4284 evt : 2 added weight : 1
+trigger passed at bin 4097  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 1 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4170  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 3 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4201  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 15 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4201 evt : 2 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4181  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 7 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4255  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 15 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4285  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 0 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4285 evt : 3 added weight : 1
+trigger passed at bin 4096  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 3 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4168  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 5 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4059  passed ch : 13 (1type) Direct dist btw posnu : 2576.36 noiseID : 10 ViewAngle : 67.183 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4168 evt : 3 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4182  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 9 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4255  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 15 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4285  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 3 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4285 evt : 4 added weight : 1
+trigger passed at bin 4095  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 3 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4170  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 12 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4157  passed ch : 13 (1type) Direct dist btw posnu : 2576.36 noiseID : 15 ViewAngle : 67.183 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4170 evt : 4 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4182  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 7 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4255  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 5 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4284  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 15 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4284 evt : 5 added weight : 1
+trigger passed at bin 4097  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 6 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4172  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 0 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4062  passed ch : 13 (1type) Direct dist btw posnu : 2576.36 noiseID : 15 ViewAngle : 67.183 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4172 evt : 5 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4181  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 4 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4254  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 13 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4283  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 3 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4283 evt : 6 added weight : 1
+trigger passed at bin 4096  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 10 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4168  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 14 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4198  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 7 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4198 evt : 6 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4181  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 10 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4254  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 12 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4154  passed ch : 6 (0type) Direct dist btw posnu : 2586.84 noiseID : 9 ViewAngle : 66.8717 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4284  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 0 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4284 evt : 7 added weight : 1
+trigger passed at bin 4096  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 15 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4170  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 5 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4200  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 13 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4200 evt : 7 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4179  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 14 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4254  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 12 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4282  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 6 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4282 evt : 8 added weight : 1
+trigger passed at bin 4097  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 0 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4167  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 12 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4050  passed ch : 7 (1type) Direct dist btw posnu : 2588.06 noiseID : 15 ViewAngle : 66.7967 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4167 evt : 8 added weight : 1
 
 Making useful event
 StationID: 2
 StationID_AraRoot: 2
 *
-trigger passed at bin 4182  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 9 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4255  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 13 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
-trigger passed at bin 4283  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 6 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
-Global_Pass : 4283 evt : 9 added weight : 1
+trigger passed at bin 4096  passed ch : 1 (1type) Direct dist btw posnu : 2564.13 noiseID : 9 ViewAngle : 67.0586 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4169  passed ch : 3 (1type) Direct dist btw posnu : 2571.43 noiseID : 8 ViewAngle : 66.6135 LikelyTrigSignal : interaction 0, ray 0
+trigger passed at bin 4198  passed ch : 11 (1type) Direct dist btw posnu : 2573.59 noiseID : 12 ViewAngle : 66.7716 LikelyTrigSignal : interaction 0, ray 0
+Global_Pass : 4198 evt : 9 added weight : 1
 
 Making useful event
 StationID: 2
@@ -324,11 +325,10 @@ StationID_AraRoot: 2
 Total Events Thrown: 10
 Total_Global_Pass : 10
 Total_Weight : 10
-Total_Probability : -88.04
+Total_Probability : -2.85546e+191
 weight bin values : 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
 
-max_dt : 6.0702e-07
-rmsdiode= 9.28804e-14
+Warning in <THistPainter::MakeChopt>: option SCAT is deprecated.
 Info in <TCanvas::Print>: pdf file sigmaCrossSection.pdf has been created
 outputdir= outputs/.
 This AraSim run is complete and will exit with code 0