[Rivet-svn] r3318 - in branches/2011-07-aida2yoda: bin data/plotinfo include/Rivet/Math src/Analyses

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Tue Aug 23 11:57:46 2011
New Revision: 3318

Log:
merge r3314-3317 from trunk

Modified:
   branches/2011-07-aida2yoda/bin/make-plots
   branches/2011-07-aida2yoda/bin/rivet-mkanalysis
   branches/2011-07-aida2yoda/data/plotinfo/D0_2010_S8671338.plot
   branches/2011-07-aida2yoda/include/Rivet/Math/MathUtils.hh
   branches/2011-07-aida2yoda/src/Analyses/LHCB_2010_S8758301.cc

Modified: branches/2011-07-aida2yoda/bin/make-plots
==============================================================================
--- branches/2011-07-aida2yoda/bin/make-plots	Tue Aug 23 11:52:06 2011	(r3317)
+++ branches/2011-07-aida2yoda/bin/make-plots	Tue Aug 23 11:57:46 2011	(r3318)
@@ -183,7 +183,7 @@
                                         ## Use strip here to deal with DOS newlines containing \r
                                         obj.description.update({prop.strip() : value.strip()})
                 cf.close()
-        
+
 
 
 

Modified: branches/2011-07-aida2yoda/bin/rivet-mkanalysis
==============================================================================
--- branches/2011-07-aida2yoda/bin/rivet-mkanalysis	Tue Aug 23 11:52:06 2011	(r3317)
+++ branches/2011-07-aida2yoda/bin/rivet-mkanalysis	Tue Aug 23 11:57:46 2011	(r3318)
@@ -232,7 +232,6 @@
 
     /// @name Histograms
     //@{
-
     AIDA::IProfile1D *_h_XXXX;
     AIDA::IHistogram1D *_h_YYYY;
     //@}

Modified: branches/2011-07-aida2yoda/data/plotinfo/D0_2010_S8671338.plot
==============================================================================
--- branches/2011-07-aida2yoda/data/plotinfo/D0_2010_S8671338.plot	Tue Aug 23 11:52:06 2011	(r3317)
+++ branches/2011-07-aida2yoda/data/plotinfo/D0_2010_S8671338.plot	Tue Aug 23 11:57:46 2011	(r3318)
@@ -1,11 +1,17 @@
+# BEGIN PLOT /D0_2010_S8671338/d0
+LegendXPos=0.4
+LegendYPos=0.9
+FullRange=1
+# END PLOT
+
 # BEGIN PLOT /D0_2010_S8671338/d01-x01-y01
-Title=Z boson pT
+Title=$Z$ boson $p_\perp$ (normalised)
 XLabel=$p_{\perp}$(Z) [GeV]
 YLabel=$1/\sigma \; \text{d}\sigma/\text{d}p_\perp(Z)$ [1/GeV]
 # END PLOT
 
 # BEGIN PLOT /D0_2010_S8671338/d02-x01-y01
-Title=Z boson pT
+Title=$Z$ boson $p_\perp$ (unnormalised)
 XLabel=$p_{\perp}$(Z) [GeV]
 YLabel=$\text{d}\sigma/\text{d}p_\perp(Z)$ [pb/GeV]
 # END PLOT

Modified: branches/2011-07-aida2yoda/include/Rivet/Math/MathUtils.hh
==============================================================================
--- branches/2011-07-aida2yoda/include/Rivet/Math/MathUtils.hh	Tue Aug 23 11:52:06 2011	(r3317)
+++ branches/2011-07-aida2yoda/include/Rivet/Math/MathUtils.hh	Tue Aug 23 11:57:46 2011	(r3318)
@@ -18,9 +18,10 @@
     return (fabs(val) < tolerance);
   }
 
-  /// Compare an integral-type number to zero. Since there is no
-  /// risk of floating point error, this function just exists in
-  /// case @c isZero is accidentally used on an integer type, to avoid
+  /// Compare an integral-type number to zero.
+  ///
+  /// Since there is no risk of floating point error, this function just exists
+  /// in case @c isZero is accidentally used on an integer type, to avoid
   /// implicit type conversion. The @a tolerance parameter is ignored.
   inline bool isZero(long val, double UNUSED(tolerance)=1E-8) {
     return val == 0;
@@ -28,15 +29,18 @@
 
 
   /// @brief Compare two floating point numbers for equality with a degree of fuzziness
-  /// The @a tolerance parameter is fractional.
+  ///
+  /// The @a tolerance parameter is fractional, based on absolute values of the args.
   inline bool fuzzyEquals(double a, double b, double tolerance=1E-5) {
-    const double absavg = fabs(a + b)/2.0;
+    const double absavg = (fabs(a) + fabs(b))/2.0;
     const double absdiff = fabs(a - b);
-    const bool rtn = (absavg == 0.0 && absdiff == 0.0) || absdiff < tolerance*absavg;
+    const bool rtn = (isZero(a) && isZero(b)) || absdiff < tolerance*absavg;
+    // cout << a << " == " << b << "? " << rtn << endl;
     return rtn;
   }
 
   /// @brief Compare two integral-type numbers for equality with a degree of fuzziness.
+  ///
   /// Since there is no risk of floating point error with integral types,
   /// this function just exists in case @c fuzzyEquals is accidentally
   /// used on an integer type, to avoid implicit type conversion. The @a
@@ -48,12 +52,14 @@
 
 
   /// @brief Compare two floating point numbers for >= with a degree of fuzziness
+  ///
   /// The @a tolerance parameter on the equality test is as for @c fuzzyEquals.
   inline bool fuzzyGtrEquals(double a, double b, double tolerance=1E-5) {
     return a > b || fuzzyEquals(a, b, tolerance);
   }
 
   /// @brief Compare two integral-type numbers for >= with a degree of fuzziness.
+  ///
   /// Since there is no risk of floating point error with integral types,
   /// this function just exists in case @c fuzzyGtrEquals is accidentally
   /// used on an integer type, to avoid implicit type conversion. The @a
@@ -65,12 +71,14 @@
 
 
   /// @brief Compare two floating point numbers for <= with a degree of fuzziness
+  ///
   /// The @a tolerance parameter on the equality test is as for @c fuzzyEquals.
   inline bool fuzzyLessEquals(double a, double b, double tolerance=1E-5) {
     return a < b || fuzzyEquals(a, b, tolerance);
   }
 
   /// @brief Compare two integral-type numbers for <= with a degree of fuzziness.
+  ///
   /// Since there is no risk of floating point error with integral types,
   /// this function just exists in case @c fuzzyLessEquals is accidentally
   /// used on an integer type, to avoid implicit type conversion. The @a
@@ -86,13 +94,15 @@
   /// @name Ranges and intervals
   //@{
 
-  /// Represents whether an interval is open (non-inclusive) or closed
-  /// (inclusive). For example, the interval \f$ [0, \pi) \f$ is closed (an inclusive
+  /// Represents whether an interval is open (non-inclusive) or closed (inclusive).
+  ///
+  /// For example, the interval \f$ [0, \pi) \f$ is closed (an inclusive
   /// boundary) at 0, and open (a non-inclusive boundary) at \f$ \pi \f$.
   enum RangeBoundary { OPEN=0, SOFT=0, CLOSED=1, HARD=1 };
 
 
   /// @brief Determine if @a value is in the range @a low to @a high, for floating point numbers
+  ///
   /// Interval boundary types are defined by @a lowbound and @a highbound.
   /// @todo Optimise to one-line at compile time?
   template<typename NUM>
@@ -101,11 +111,11 @@
     if (lowbound == OPEN && highbound == OPEN) {
       return (value > low && value < high);
     } else if (lowbound == OPEN && highbound == CLOSED) {
-      return (value > low && (value < high || fuzzyEquals(value, high)));
+      return (value > low && fuzzyLessEquals(value, high));
     } else if (lowbound == CLOSED && highbound == OPEN) {
-      return ((value > low || fuzzyEquals(value, low)) && value < high);
+      return (fuzzyGtrEquals(value, low) && value < high);
     } else { // if (lowbound == CLOSED && highbound == CLOSED) {
-      return ((value > low || fuzzyEquals(value, low)) && (value < high || fuzzyEquals(value, high)));
+      return (fuzzyGtrEquals(value, low) && fuzzyLessEquals(value, high));
     }
   }
 
@@ -118,6 +128,7 @@
 
 
   /// @brief Determine if @a value is in the range @a low to @a high, for integer types
+  ///
   /// Interval boundary types are defined by @a lowbound and @a highbound.
   /// @todo Optimise to one-line at compile time?
   inline bool inRange(int value, int low, int high,
@@ -221,6 +232,7 @@
     const double logstart = std::log(start);
     const double logend = std::log(end);
     const vector<double> logvals = linspace(logstart, logend, nbins);
+    assert(logvals.size() == nbins+1);
     vector<double> rtn;
     foreach (double logval, logvals) {
       rtn.push_back(std::exp(logval));
@@ -231,6 +243,7 @@
 
 
   /// @brief Return the bin index of the given value, @a val, given a vector of bin edges
+  ///
   /// NB. The @a binedges vector must be sorted
   template <typename NUM>
   inline int index_between(const NUM& val, const vector<NUM>& binedges) {
@@ -329,7 +342,7 @@
 
 
     // Calculate the error on the correlation strength
-    const double corr_strength_err = correlation_err*sqrt(var2/var1) + 
+    const double corr_strength_err = correlation_err*sqrt(var2/var1) +
       correlation/(2*sqrt(var2/var1)) * (var2_e/var1 - var2*var1_e/pow(2, var2));
 
     return corr_strength_err;
@@ -340,8 +353,10 @@
   /// @name Angle range mappings
   //@{
 
-  /// Reduce any number to the range [-2PI, 2PI] by repeated addition or
-  /// subtraction of 2PI as required. Used to normalise angular measures.
+  /// @brief Reduce any number to the range [-2PI, 2PI]
+  ///
+  /// Achieved by repeated addition or subtraction of 2PI as required. Used to
+  /// normalise angular measures.
   inline double _mapAngleM2PITo2Pi(double angle) {
     double rtn = fmod(angle, TWOPI);
     if (isZero(rtn)) return 0;
@@ -383,8 +398,9 @@
   /// @name Phase space measure helpers
   //@{
 
-  /// Calculate the difference between two angles in radians,
-  /// returning in the range [0, PI].
+  /// @brief Calculate the difference between two angles in radians
+  ///
+  /// Returns in the range [0, PI].
   inline double deltaPhi(double phi1, double phi2) {
     return mapAngle0ToPi(phi1 - phi2);
   }

Modified: branches/2011-07-aida2yoda/src/Analyses/LHCB_2010_S8758301.cc
==============================================================================
--- branches/2011-07-aida2yoda/src/Analyses/LHCB_2010_S8758301.cc	Tue Aug 23 11:52:06 2011	(r3317)
+++ branches/2011-07-aida2yoda/src/Analyses/LHCB_2010_S8758301.cc	Tue Aug 23 11:57:46 2011	(r3318)
@@ -20,7 +20,6 @@
       sumKs0_30(0.0), sumKs0_35(0.0), sumKs0_40(0.0)
 
     {
-      /// @todo Set whether your finalize method needs the generator cross section
       setNeedsCrossSection(false);
     }
 
@@ -72,22 +71,17 @@
           }
         }
       }
-      
+
 
     }
 
 
     /// Normalise histograms etc., after the run
     void finalize() {
-
-
-      // here we normalise to reference data... 
+      /// @todo We have to normalise to reference data... for now
       scale(_h_K0s_pt_y_30, 4.880000e+02/sumKs0_30);
       scale(_h_K0s_pt_y_35, 4.442000e+02/sumKs0_35);
       scale(_h_K0s_pt_y_40, 3.868000e+02/sumKs0_40);
-
-
-
     }
 
     //@}
@@ -95,15 +89,9 @@
 
   private:
 
-    // Data members like post-cuts event weight counters go here
-
-
-  private:
-
     /// @name Histograms
     //@{
 
-
     Histo1DPtr _h_K0s_pt_y_30;
     Histo1DPtr _h_K0s_pt_y_35;
     Histo1DPtr _h_K0s_pt_y_40;