[Rivet-svn] r2414 - trunk/src/Analyses

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Fri Apr 23 21:11:02 2010
New Revision: 2414

Log:
Forgot to scale the deta and dR histograms in MC_JetAnalysis.

Modified:
   trunk/src/Analyses/MC_JetAnalysis.cc

Modified: trunk/src/Analyses/MC_JetAnalysis.cc
==============================================================================
--- trunk/src/Analyses/MC_JetAnalysis.cc	Thu Apr 22 16:06:46 2010	(r2413)
+++ trunk/src/Analyses/MC_JetAnalysis.cc	Fri Apr 23 21:11:02 2010	(r2414)
@@ -141,14 +141,21 @@
       scale(_h_pT_jet[i], crossSection()/sumOfWeights());
       scale(_h_eta_jet[i], crossSection()/sumOfWeights());
 
-      for (size_t j=i+1; j<m_njet; ++j) {
-      }
     }
     for (int ibin=0; ibin<_h_log10_R[m_njet]->size(); ++ibin) {
       IDataPoint* dp=_h_log10_R[m_njet]->point(ibin);
       dp->coordinate(1)->setValue(dp->coordinate(1)->value()*crossSection()/sumOfWeights());
     }
 
+    // scale the d{eta,R} histograms
+    map<pair<size_t, size_t>, AIDA::IHistogram1D*>::iterator it;
+    for (it=_h_deta_jets.begin(); it!=_h_deta_jets.begin(); ++it) {
+      scale(it->second, crossSection()/sumOfWeights());
+    }
+    for (it=_h_dR_jets.begin(); it!=_h_dR_jets.begin(); ++it) {
+      scale(it->second, crossSection()/sumOfWeights());
+    }
+    
     // fill inclusive jet multi ratio
     int Nbins=_h_jet_multi_inclusive->axis().bins();
     std::vector<double> ratio(Nbins-1, 0.0);