NuPath.h

Go to the documentation of this file.

 
 
 
#ifndef NUPATH_H
#define NUPATH_H
 
#include <vector>
 
namespace OscProb {
 
  struct NuPath {
      NuPath(double l = 0, double d = 0, double z = 0.5, int ly = 0)
      {
        SetPath(l, d, z, ly);
      }
 
      void SetPath(double l = 0, double d = 0, double z = 0.5, int ly = 0)
      {
        length  = l;
        density = d;
        zoa     = z;
        layer   = ly;
      }
 
      double length;  
      double density; 
      double zoa;     
      int    layer;   
  };
 
  struct PremLayer {
      PremLayer(double r = 0, double d = 0, double z = 0.5, int ly = 0)
      {
        SetLayer(r, d, z, ly);
      }
 
      void SetLayer(double r = 0, double d = 0, double z = 0.5, int ly = 0)
      {
        radius  = r;
        density = d;
        zoa     = z;
        layer   = ly;
      }
 
      double radius;  
      double density; 
      double zoa;     
      int    layer;   
 
      bool operator==(PremLayer& rhs)
      {
        return rhs.density == density && rhs.zoa == zoa && rhs.layer == layer;
      }
  };
 
  NuPath AvgPath(NuPath& p1, NuPath& p2); 
  NuPath AvgPath(
      std::vector<NuPath>& pv); 
  std::vector<NuPath> MergePaths(std::vector<NuPath>& inputPath, int j,
                                 int k); 
 
} // namespace OscProb
 
#endif