[Rivet-svn] r2051 - in trunk: . data/anainfo src/Analyses

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Tue Nov 10 12:46:37 2009
New Revision: 2051

Log:
STAR_2006_S6500200: Fixing normalisation after adding the
trigger requirement. I consider this analysis now validated.

Modified:
   trunk/ChangeLog
   trunk/data/anainfo/STAR_2006_S6500200.info
   trunk/src/Analyses/STAR_2006_S6500200.cc

Modified: trunk/ChangeLog
==============================================================================
--- trunk/ChangeLog	Tue Nov 10 01:13:37 2009	(r2050)
+++ trunk/ChangeLog	Tue Nov 10 12:46:37 2009	(r2051)
@@ -1,3 +1,8 @@
+2009-11-10  Hendrik Hoeth <hendrik.hoeth at cern.ch>
+
+	* Fixing last issues in STAR_2006_S6500200 and setting it to
+	VALIDATED.
+
 2009-11-09  Andy Buckley  <andy at insectnation.org>
 
 	* Overhaul of jet caching and ParticleBase interface.

Modified: trunk/data/anainfo/STAR_2006_S6500200.info
==============================================================================
--- trunk/data/anainfo/STAR_2006_S6500200.info	Tue Nov 10 01:13:37 2009	(r2050)
+++ trunk/data/anainfo/STAR_2006_S6500200.info	Tue Nov 10 12:46:37 2009	(r2051)
@@ -4,7 +4,7 @@
 Experiment: STAR
 Collider: RHIC pp 200 GeV
 SpiresID: 6500200
-Status: UNVALIDATED
+Status: VALIDATED
 Authors:
  - Bedanga Mohanty <bedanga at rcf.bnl.gov>
  - Hendrik Hoeth <hendrik.hoeth at cern.ch>

Modified: trunk/src/Analyses/STAR_2006_S6500200.cc
==============================================================================
--- trunk/src/Analyses/STAR_2006_S6500200.cc	Tue Nov 10 01:13:37 2009	(r2050)
+++ trunk/src/Analyses/STAR_2006_S6500200.cc	Tue Nov 10 12:46:37 2009	(r2051)
@@ -13,7 +13,8 @@
 
     /// Constructor
     STAR_2006_S6500200()
-      : Analysis("STAR_2006_S6500200")
+      : Analysis("STAR_2006_S6500200"),
+        _sumWeightSelected(0.0)
     {
       setBeams(PROTON, PROTON);
     }
@@ -73,6 +74,7 @@
           }
         }
       }
+      _sumWeightSelected += event.weight();
     }
 
 
@@ -86,14 +88,18 @@
       hf.divide(dir + "/d02-x03-y01", *_h_pT_proton, *_h_pT_piplus);
       hf.divide(dir + "/d02-x04-y01", *_h_pT_antiproton, *_h_pT_piminus);
 
-      scale(_h_pT_piplus,     1./(2*M_PI*sumOfWeights()));
-      scale(_h_pT_piminus,    1./(2*M_PI*sumOfWeights()));
-      scale(_h_pT_proton,     1./(2*M_PI*sumOfWeights()));
-      scale(_h_pT_antiproton, 1./(2*M_PI*sumOfWeights()));
+      scale(_h_pT_piplus,     1./(2*M_PI*_sumWeightSelected));
+      scale(_h_pT_piminus,    1./(2*M_PI*_sumWeightSelected));
+      scale(_h_pT_proton,     1./(2*M_PI*_sumWeightSelected));
+      scale(_h_pT_antiproton, 1./(2*M_PI*_sumWeightSelected));
+      getLog() << Log::DEBUG << "sumOfWeights()     = " << sumOfWeights() << std::endl;
+      getLog() << Log::DEBUG << "_sumWeightSelected = " << _sumWeightSelected << std::endl;
     }
 
   private:
 
+    double _sumWeightSelected;
+
     AIDA::IHistogram1D * _h_pT_piplus;
     AIDA::IHistogram1D * _h_pT_piminus;
     AIDA::IHistogram1D * _h_pT_proton;