[Rivet-svn] r3720 - branches/2011-07-aida2yoda/src/Analyses

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Wed May  2 13:59:05 2012
New Revision: 3720

Log:
First analysis with YODA histogram division

Modified:
   branches/2011-07-aida2yoda/src/Analyses/MC_GENERIC.cc

Modified: branches/2011-07-aida2yoda/src/Analyses/MC_GENERIC.cc
==============================================================================
--- branches/2011-07-aida2yoda/src/Analyses/MC_GENERIC.cc	Wed May  2 12:25:20 2012	(r3719)
+++ branches/2011-07-aida2yoda/src/Analyses/MC_GENERIC.cc	Wed May  2 13:59:05 2012	(r3720)
@@ -62,6 +62,11 @@
 
       _histPhi    = bookHisto1D("Phi", 50, 0, TWOPI);
       _histPhiCh  = bookHisto1D("PhiCh", 50, 0, TWOPI);
+
+      _histEtaPMRatio = bookScatter2D("EtaPMRatio");
+      _histEtaChPMRatio = bookScatter2D("EtaChPMRatio");
+      _histRapidityPMRatio = bookScatter2D("RapidityPMRatio");
+      _histRapidityChPMRatio = bookScatter2D("RapidityChPMRatio");
     }
 
 
@@ -144,11 +149,10 @@
       scale(_histPhi, 1/sumOfWeights());
       scale(_histPhiCh, 1/sumOfWeights());
 
-      // \todo YODA divide
-      // histogramFactory().divide(histoPath("EtaPMRatio"), *_tmphistEtaPlus, *_tmphistEtaMinus);
-      // histogramFactory().divide(histoPath("EtaChPMRatio"), *_tmphistEtaChPlus, *_tmphistEtaChMinus);
-      // histogramFactory().divide(histoPath("RapidityPMRatio"), *_tmphistRapPlus, *_tmphistRapMinus);
-      // histogramFactory().divide(histoPath("RapidityChPMRatio"), *_tmphistRapChPlus, *_tmphistRapChMinus);
+      *_histEtaPMRatio = *_tmphistEtaPlus/ *_tmphistEtaMinus;
+      *_histEtaChPMRatio = *_tmphistEtaChPlus/ *_tmphistEtaChMinus;
+      *_histRapidityPMRatio = *_tmphistRapPlus/ *_tmphistRapMinus;
+      *_histRapidityChPMRatio = *_tmphistRapChPlus/ *_tmphistRapChMinus;
     }
 
     //@}
@@ -171,6 +175,10 @@
     Histo1DPtr _histPt, _histPtCh;
     Histo1DPtr _histE, _histECh;
     Histo1DPtr _histPhi, _histPhiCh;
+    Scatter2DPtr _histEtaPMRatio;
+    Scatter2DPtr _histEtaChPMRatio;
+    Scatter2DPtr _histRapidityPMRatio;
+    Scatter2DPtr _histRapidityChPMRatio;
     //@}
 
   };