[Rivet-svn] r2597 - in trunk: data/plotinfo src/Analyses

[blackhole at projects.hepforge.org]

Author: richardn
Date: Thu Jul 15 08:41:41 2010
New Revision: 2597

Log:
fixed to D0_2008_S7837160 for #431

Modified:
   trunk/data/plotinfo/D0_2008_S7837160.plot
   trunk/src/Analyses/D0_2008_S7837160.cc

Modified: trunk/data/plotinfo/D0_2008_S7837160.plot
==============================================================================
--- trunk/data/plotinfo/D0_2008_S7837160.plot	Wed Jul 14 23:01:22 2010	(r2596)
+++ trunk/data/plotinfo/D0_2008_S7837160.plot	Thu Jul 15 08:41:41 2010	(r2597)
@@ -3,6 +3,8 @@
 XLabel=$|\eta|$ of charged lepton
 YLabel=$\mathcal{A}(\mathrm{d}{\sigma^+}\!/\!\mathrm{d}{|\eta|},\, \mathrm{d}{\sigma^-}\!/\!\mathrm{d}{|\eta|})$
 LogY=0
+LegendXPos=0.5
+LegendYPos=0.4
 # END PLOT
 
 # BEGIN PLOT /D0_2008_S7837160/d01-x01-y02
@@ -10,6 +12,8 @@
 XLabel=$|\eta|$ of charged lepton
 YLabel=$\mathcal{A}(\mathrm{d}{\sigma^+}\!/\!\mathrm{d}{|\eta|},\, \mathrm{d}{\sigma^-}\!/\!\mathrm{d}{|\eta|})$
 LogY=0
+LegendXPos=0.5
+LegendYPos=0.4
 # END PLOT
 
 # BEGIN PLOT /D0_2008_S7837160/d01-x01-y03
@@ -17,4 +21,6 @@
 XLabel=$|\eta|$ of charged lepton
 YLabel=$\mathcal{A}(\mathrm{d}{\sigma^+}\!/\!\mathrm{d}{|\eta|},\, \mathrm{d}{\sigma^-}\!/\!\mathrm{d}{|\eta|})$
 LogY=0
+LegendXPos=0.5
+LegendYPos=0.4
 # END PLOT

Modified: trunk/src/Analyses/D0_2008_S7837160.cc
==============================================================================
--- trunk/src/Analyses/D0_2008_S7837160.cc	Wed Jul 14 23:01:22 2010	(r2596)
+++ trunk/src/Analyses/D0_2008_S7837160.cc	Thu Jul 15 08:41:41 2010	(r2597)
@@ -78,7 +78,9 @@
       assert(chg_e != 0);
 
       const double weight = event.weight();
-      const double eta_e = fabs(p_e.pseudorapidity());
+      double eta_e = p_e.pseudorapidity();
+      if(eta_e<0) chg_e *= -1;
+      eta_e = fabs(eta_e);
       const double et_e = p_e.Et();
       if (et_e < 35*GeV) {
         // 25 <= ET < 35
@@ -110,6 +112,7 @@
       AIDA::IHistogramFactory& hf = histogramFactory();
 
       IHistogram1D* num25_35 = hf.subtract("/num25_35", *_h_dsigplus_deta_25_35, *_h_dsigminus_deta_25_35);
+      num25_35->scale(100.);
       IHistogram1D* denom25_35 = hf.add("/denom25_35", *_h_dsigplus_deta_25_35, *_h_dsigminus_deta_25_35);
       assert(num25_35 && denom25_35);
       hf.divide(histoDir() + "/d01-x01-y01", *num25_35, *denom25_35);
@@ -117,6 +120,7 @@
       hf.destroy(denom25_35);
       //
       IHistogram1D* num35 = hf.subtract("/num35", *_h_dsigplus_deta_35, *_h_dsigminus_deta_35);
+      num35->scale(100.);
       IHistogram1D* denom35 = hf.add("/denom35", *_h_dsigplus_deta_35, *_h_dsigminus_deta_35);
       assert(num35 && denom35);
       hf.divide(histoDir() + "/d01-x01-y02", *num35, *denom35);
@@ -124,6 +128,7 @@
       hf.destroy(denom35);
       //
       IHistogram1D* num25 = hf.subtract("/num25", *_h_dsigplus_deta_25, *_h_dsigminus_deta_25);
+      num25->scale(100.);
       IHistogram1D* denom25 = hf.add("/denom25", *_h_dsigplus_deta_25, *_h_dsigminus_deta_25);
       assert(num25 && denom25);
       hf.divide(histoDir() + "/d01-x01-y03", *num25, *denom25);