|
[Rivet-svn] r1994 - in trunk: data/plotinfo src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgWed Nov 4 10:58:20 GMT 2009
Author: fsiegert Date: Wed Nov 4 10:58:19 2009 New Revision: 1994 Log: Revert r1953 such that histograms of the affected analyses are normalised in the finalize method again instead of in the plotting tools. Modified: trunk/data/plotinfo/D0_1996_S3214044.plot trunk/data/plotinfo/D0_1996_S3324664.plot trunk/data/plotinfo/D0_2004_S5992206.plot trunk/data/plotinfo/D0_2007_S7075677.plot trunk/data/plotinfo/D0_2008_S6879055.plot trunk/data/plotinfo/D0_2008_S7554427.plot trunk/data/plotinfo/D0_2009_S8320160.plot trunk/src/Analyses/D0_1996_S3214044.cc trunk/src/Analyses/D0_1996_S3324664.cc trunk/src/Analyses/D0_2004_S5992206.cc trunk/src/Analyses/D0_2007_S7075677.cc trunk/src/Analyses/D0_2008_S6879055.cc trunk/src/Analyses/D0_2008_S7554427.cc trunk/src/Analyses/D0_2009_S8320160.cc Modified: trunk/data/plotinfo/D0_1996_S3214044.plot ============================================================================== --- trunk/data/plotinfo/D0_1996_S3214044.plot Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/data/plotinfo/D0_1996_S3214044.plot Wed Nov 4 10:58:19 2009 (r1994) @@ -3,7 +3,6 @@ XLabel=$x_3 = \frac{2E_3}{\sqrt{s}}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}x_3$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d02-x01-y01 @@ -11,7 +10,6 @@ XLabel=$x_5 = \frac{2E_5}{\sqrt{s}}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}x_5$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d03-x01-y01 @@ -19,7 +17,6 @@ XLabel=$\cos(\theta_3^*)$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\theta_3^*)$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d04-x01-y01 @@ -27,7 +24,6 @@ XLabel=$\psi^*$ (degree) YLabel=$1/N \mathrm{d}N/\mathrm{d}\psi^*$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d05-x01-y01 @@ -35,7 +31,6 @@ XLabel=$\mu_{34}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{34}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d06-x01-y01 @@ -43,7 +38,6 @@ XLabel=$\mu_{35}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{35}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d07-x01-y01 @@ -51,7 +45,6 @@ XLabel=$\mu_{45}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{45}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d08-x01-y01 @@ -59,7 +52,6 @@ XLabel=$x_3 = \frac{2E_3}{\sqrt{s}}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}x_3$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d09-x01-y01 @@ -67,7 +59,6 @@ XLabel=$x_4 = \frac{2E_4}{\sqrt{s}}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}x_4$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d10-x01-y01 @@ -75,7 +66,6 @@ XLabel=$x_5 = \frac{2E_5}{\sqrt{s}}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}x_5$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d11-x01-y01 @@ -83,7 +73,6 @@ XLabel=$x_6 = \frac{2E_6}{\sqrt{s}}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}x_6$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d12-x01-y01 @@ -91,7 +80,6 @@ XLabel=$\cos(\theta_3^*)$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\theta_3^*)$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d13-x01-y01 @@ -99,7 +87,6 @@ XLabel=$\cos(\theta_4^*)$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\theta_4^*)$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d14-x01-y01 @@ -107,7 +94,6 @@ XLabel=$\cos(\theta_5^*)$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\theta_5^*)$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d15-x01-y01 @@ -115,7 +101,6 @@ XLabel=$\cos(\theta_6^*)$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\theta_6^*)$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d16-x01-y01 @@ -123,7 +108,6 @@ XLabel=$\cos(\omega_{34})$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\omega_{34}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d17-x01-y01 @@ -131,7 +115,6 @@ XLabel=$\cos(\omega_{35})$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\omega_{35})$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d18-x01-y01 @@ -139,7 +122,6 @@ XLabel=$\cos(\omega_{36})$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\omega_{36})$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d19-x01-y01 @@ -147,7 +129,6 @@ XLabel=$\cos(\omega_{45})$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\omega_{45})$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d20-x01-y01 @@ -155,7 +136,6 @@ XLabel=$\cos(\omega_{46})$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\omega_{46})$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d21-x01-y01 @@ -163,7 +143,6 @@ XLabel=$\cos(\omega_{56})$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\omega_{56})$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d22-x01-y01 @@ -171,7 +150,6 @@ XLabel=$\mu_{34}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{34}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d23-x01-y01 @@ -179,7 +157,6 @@ XLabel=$\mu_{35}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{35}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d24-x01-y01 @@ -187,7 +164,6 @@ XLabel=$\mu_{36}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{36}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d25-x01-y01 @@ -195,7 +171,6 @@ XLabel=$\mu_{45}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{45}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d26-x01-y01 @@ -203,7 +178,6 @@ XLabel=$\mu_{46}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{46}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d27-x01-y01 @@ -211,7 +185,6 @@ XLabel=$\mu_{56}$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\mu_{56}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d28-x01-y01 @@ -219,7 +192,6 @@ XLabel=$\chi_\mathrm{BZ}$ (degree) YLabel=$1/N \mathrm{d}N/\mathrm{d}\chi_\mathrm{BZ}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_1996_S3214044/d29-x01-y01 @@ -227,5 +199,4 @@ XLabel=$\cos(\theta_\mathrm{NR})$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\cos(\theta_\mathrm{NR})$ LogY=0 -NormalizeToIntegral=1 # END PLOT Modified: trunk/data/plotinfo/D0_1996_S3324664.plot ============================================================================== --- trunk/data/plotinfo/D0_1996_S3324664.plot Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/data/plotinfo/D0_1996_S3324664.plot Wed Nov 4 10:58:19 2009 (r1994) @@ -2,39 +2,24 @@ Title=Pseudorapidity difference of the two opposite jets XLabel=$\Delta \eta = |\eta_1 - \eta_2|$ YLabel=$\mathrm{d}N/\mathrm{d}\Delta\eta$ -NormalizeToIntegral=1 -# Normalised to #events -Scale=8830.0 # END PLOT # BEGIN PLOT /D0_1996_S3324664/d02-x01-y01 Title=Azimuthal angle difference for $0 < \Delta \eta < 2$ XLabel=$1 - \Delta \phi/\pi$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\Delta\phi$ -# I have no idea what this is normalised to... in the paper it says unity! -# @todo Understand this! -NormalizeToIntegral=1 -Scale=0.0798 # END PLOT # BEGIN PLOT /D0_1996_S3324664/d02-x01-y02 Title=Azimuthal angle difference for $2 < \Delta \eta < 4$ XLabel=$1 - \Delta \phi/\pi$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\Delta\phi$ -# I have no idea what this is normalised to... in the paper it says unity! -# @todo Understand this! -NormalizeToIntegral=1 -Scale=0.0798 # END PLOT # BEGIN PLOT /D0_1996_S3324664/d02-x01-y03 Title=Azimuthal angle difference for $4 < \Delta \eta < 6$ XLabel=$1 - \Delta \phi/\pi$ YLabel=$1/N \mathrm{d}N/\mathrm{d}\Delta\phi$ -# I have no idea what this is normalised to... in the paper it says unity! -# @todo Understand this! -NormalizeToIntegral=1 -Scale=0.0798 # END PLOT # BEGIN PLOT /D0_1996_S3324664/d03-x01-y01 Modified: trunk/data/plotinfo/D0_2004_S5992206.plot ============================================================================== --- trunk/data/plotinfo/D0_2004_S5992206.plot Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/data/plotinfo/D0_2004_S5992206.plot Wed Nov 4 10:58:19 2009 (r1994) @@ -2,27 +2,23 @@ Title=Jet--jet azimuthal angle, $p_\perp^\text{max} \in [75,100]\text{ GeV}$ XLabel=$\Delta \phi_\text{dijet}$ / rad YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2004_S5992206/d02-x02-y01 Title=Jet--jet azimuthal angle, $p_\perp^\text{max} \in [100..130]\text{ GeV}$ XLabel=$\Delta \phi_\text{dijet}$ / rad YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2004_S5992206/d03-x02-y01 Title=Jet--jet azimuthal angle, $p_\perp^\text{max} \in [130..180]\text{ GeV}$ XLabel=$\Delta \phi_\text{dijet}$ / rad YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2004_S5992206/d04-x02-y01 Title=Jet--jet azimuthal angle, $p_\perp^\text{max} > 180\text{ GeV}$ XLabel=$\Delta \phi_\text{dijet}$ / rad YLabel=$1/\sigma_\text{dijet} \, \text{d}\sigma_\text{dijet}/\text{d}\Delta\phi_\text{dijet}$ -NormalizeToIntegral=1 # END PLOT Modified: trunk/data/plotinfo/D0_2007_S7075677.plot ============================================================================== --- trunk/data/plotinfo/D0_2007_S7075677.plot Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/data/plotinfo/D0_2007_S7075677.plot Wed Nov 4 10:58:19 2009 (r1994) @@ -4,9 +4,5 @@ YLabel=$1/\sigma \; \text{d}\sigma/\text{d}|y|(Z)$ LegendXPos=0.5 LegendYPos=0.6 -NormalizeToIntegral=1 -# Data seems to have been normalized for the avg of the two sides -# (+ve & -ve rapidity) rather than the sum, hence the 0.5: -Scale=0.5 # END PLOT Modified: trunk/data/plotinfo/D0_2008_S6879055.plot ============================================================================== --- trunk/data/plotinfo/D0_2008_S6879055.plot Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/data/plotinfo/D0_2008_S6879055.plot Wed Nov 4 10:58:19 2009 (r1994) @@ -10,32 +10,17 @@ Title=$p_\perp$ of 1st jet (not detector-corrected!) XLabel=$p_\perp$(1st jet) [GeV] YLabel=$1/\sigma \; \text{d}\sigma/\text{d}p_\perp$ -# Normalise jet pT's to integral of data -# there is no other way to do this, because these quantities are not -# detector corrected -NormalizeToIntegral=1 -Scale=10439.0 # END PLOT # BEGIN PLOT /D0_2008_S6879055/d03-x01-y01 Title=$p_\perp$ of 2nd jet (not detector-corrected!) XLabel=$p_\perp$(2nd jet) [GeV] YLabel=$1/\sigma \; \text{d}\sigma/\text{d}p_\perp$ -# Normalise jet pT's to integral of data -# there is no other way to do this, because these quantities are not -# detector corrected -NormalizeToIntegral=1 -Scale=1461.5 # END PLOT # BEGIN PLOT /D0_2008_S6879055/d04-x01-y01 Title=$p_\perp$ of 3rd jet (not detector-corrected!) XLabel=$p_\perp$(3rd jet) [GeV] YLabel=$1/\sigma \; \text{d}\sigma/\text{d}p_\perp$ -# Normalise jet pT's to integral of data -# there is no other way to do this, because these quantities are not -# detector corrected -NormalizeToIntegral=1 -Scale=217.0 # END PLOT Modified: trunk/data/plotinfo/D0_2008_S7554427.plot ============================================================================== --- trunk/data/plotinfo/D0_2008_S7554427.plot Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/data/plotinfo/D0_2008_S7554427.plot Wed Nov 4 10:58:19 2009 (r1994) @@ -2,13 +2,11 @@ Title=Z boson pT XLabel=$p_{\perp}$(Z) [GeV] YLabel=$1/\sigma \; \text{d}\sigma/\text{d}p_\perp(Z)$ -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2008_S7554427/d03-x01-y01 Title=Z boson pT (forward region only) XLabel=$p_{\perp}$(Z) [GeV] YLabel=$1/\sigma \; \text{d}\sigma/\text{d}p_\perp(Z)$ -NormalizeToIntegral=1 # END PLOT Modified: trunk/data/plotinfo/D0_2009_S8320160.plot ============================================================================== --- trunk/data/plotinfo/D0_2009_S8320160.plot Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/data/plotinfo/D0_2009_S8320160.plot Wed Nov 4 10:58:19 2009 (r1994) @@ -3,7 +3,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d02-x01-y01 @@ -11,7 +10,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d03-x01-y01 @@ -19,7 +17,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d04-x01-y01 @@ -27,7 +24,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d05-x01-y01 @@ -35,7 +31,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d06-x01-y01 @@ -43,7 +38,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d07-x01-y01 @@ -51,7 +45,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d08-x01-y01 @@ -59,7 +52,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d09-x01-y01 @@ -67,7 +59,6 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT # BEGIN PLOT /D0_2009_S8320160/d10-x01-y01 @@ -75,5 +66,4 @@ XLabel=$\chi_{\mathrm{dijet}} = \exp(|y_1-y_2|)$ YLabel=$1/\sigma_{\mathrm{dijet}} \, \mathrm{d}\sigma/\mathrm{d}\chi_{\mathrm{dijet}}$ LogY=0 -NormalizeToIntegral=1 # END PLOT Modified: trunk/src/Analyses/D0_1996_S3214044.cc ============================================================================== --- trunk/src/Analyses/D0_1996_S3214044.cc Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/src/Analyses/D0_1996_S3214044.cc Wed Nov 4 10:58:19 2009 (r1994) @@ -21,7 +21,7 @@ D0_1996_S3214044() : Analysis("D0_1996_S3214044") { setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); + setNeedsCrossSection(false); } @@ -103,36 +103,35 @@ void finalize() { - const double factor = crossSection()/sumOfWeights(); - scale(_h_3j_x3, factor); - scale(_h_3j_x5, factor); - scale(_h_3j_costheta3, factor); - scale(_h_3j_psi, factor); - scale(_h_3j_mu34, factor); - scale(_h_3j_mu35, factor); - scale(_h_3j_mu45, factor); - scale(_h_4j_x3, factor); - scale(_h_4j_x4, factor); - scale(_h_4j_x5, factor); - scale(_h_4j_x6, factor); - scale(_h_4j_costheta3, factor); - scale(_h_4j_costheta4, factor); - scale(_h_4j_costheta5, factor); - scale(_h_4j_costheta6, factor); - scale(_h_4j_cosomega34, factor); - scale(_h_4j_cosomega35, factor); - scale(_h_4j_cosomega36, factor); - scale(_h_4j_cosomega45, factor); - scale(_h_4j_cosomega46, factor); - scale(_h_4j_cosomega56, factor); - scale(_h_4j_mu34, factor); - scale(_h_4j_mu35, factor); - scale(_h_4j_mu36, factor); - scale(_h_4j_mu45, factor); - scale(_h_4j_mu46, factor); - scale(_h_4j_mu56, factor); - scale(_h_4j_theta_BZ, factor); - scale(_h_4j_costheta_NR, factor); + normalize(_h_3j_x3, 1.0); + normalize(_h_3j_x5, 1.0); + normalize(_h_3j_costheta3, 1.0); + normalize(_h_3j_psi, 1.0); + normalize(_h_3j_mu34, 1.0); + normalize(_h_3j_mu35, 1.0); + normalize(_h_3j_mu45, 1.0); + normalize(_h_4j_x3, 1.0); + normalize(_h_4j_x4, 1.0); + normalize(_h_4j_x5, 1.0); + normalize(_h_4j_x6, 1.0); + normalize(_h_4j_costheta3, 1.0); + normalize(_h_4j_costheta4, 1.0); + normalize(_h_4j_costheta5, 1.0); + normalize(_h_4j_costheta6, 1.0); + normalize(_h_4j_cosomega34, 1.0); + normalize(_h_4j_cosomega35, 1.0); + normalize(_h_4j_cosomega36, 1.0); + normalize(_h_4j_cosomega45, 1.0); + normalize(_h_4j_cosomega46, 1.0); + normalize(_h_4j_cosomega56, 1.0); + normalize(_h_4j_mu34, 1.0); + normalize(_h_4j_mu35, 1.0); + normalize(_h_4j_mu36, 1.0); + normalize(_h_4j_mu45, 1.0); + normalize(_h_4j_mu46, 1.0); + normalize(_h_4j_mu56, 1.0); + normalize(_h_4j_theta_BZ, 1.0); + normalize(_h_4j_costheta_NR, 1.0); } //@} Modified: trunk/src/Analyses/D0_1996_S3324664.cc ============================================================================== --- trunk/src/Analyses/D0_1996_S3324664.cc Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/src/Analyses/D0_1996_S3324664.cc Wed Nov 4 10:58:19 2009 (r1994) @@ -18,7 +18,7 @@ /// Constructor D0_1996_S3324664() : Analysis("D0_1996_S3324664") { setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); + setNeedsCrossSection(false); } @@ -85,10 +85,13 @@ void finalize() { - scale(_h_deta, crossSection()/sumOfWeights()); + // Normalised to #events + normalize(_h_deta, 8830.0); + // I have no idea what this is normalised to... in the paper it says unity! + /// @todo Understand this! foreach (IHistogram1D* histo, _h_dphi.getHistograms()) { - scale(histo, crossSection()/sumOfWeights()); + normalize(histo, 0.0798); } } Modified: trunk/src/Analyses/D0_2004_S5992206.cc ============================================================================== --- trunk/src/Analyses/D0_2004_S5992206.cc Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/src/Analyses/D0_2004_S5992206.cc Wed Nov 4 10:58:19 2009 (r1994) @@ -33,7 +33,6 @@ D0_2004_S5992206() : Analysis("D0_2004_S5992206") { setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); } //@} @@ -109,10 +108,11 @@ // Finalize void finalize() { - scale(_histJetAzimuth_pTmax75_100, crossSection()/sumOfWeights()); - scale(_histJetAzimuth_pTmax100_130, crossSection()/sumOfWeights()); - scale(_histJetAzimuth_pTmax130_180, crossSection()/sumOfWeights()); - scale(_histJetAzimuth_pTmax180_, crossSection()/sumOfWeights()); + // Normalize histograms to unit area + normalize(_histJetAzimuth_pTmax75_100); + normalize(_histJetAzimuth_pTmax100_130); + normalize(_histJetAzimuth_pTmax130_180); + normalize(_histJetAzimuth_pTmax180_); } //@} Modified: trunk/src/Analyses/D0_2007_S7075677.cc ============================================================================== --- trunk/src/Analyses/D0_2007_S7075677.cc Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/src/Analyses/D0_2007_S7075677.cc Wed Nov 4 10:58:19 2009 (r1994) @@ -20,7 +20,6 @@ { // Run II Z rapidity setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); } @@ -63,7 +62,9 @@ // Finalize void finalize() { - scale(_h_yZ, crossSection()/sumOfWeights()); + // Data seems to have been normalized for the avg of the two sides + // (+ve & -ve rapidity) rather than the sum, hence the 0.5: + normalize(_h_yZ, 0.5); } //@} Modified: trunk/src/Analyses/D0_2008_S6879055.cc ============================================================================== --- trunk/src/Analyses/D0_2008_S6879055.cc Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/src/Analyses/D0_2008_S6879055.cc Wed Nov 4 10:58:19 2009 (r1994) @@ -20,7 +20,6 @@ D0_2008_S6879055() : Analysis("D0_2008_S6879055") { setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); } @@ -145,9 +144,12 @@ // Now divide by the inclusive result _crossSectionRatio->scale(1.0/_crossSectionRatio->binHeight(0)); - scale(_pTjet1, crossSection()/sumOfWeights()); - scale(_pTjet2, crossSection()/sumOfWeights()); - scale(_pTjet3, crossSection()/sumOfWeights()); + // Normalise jet pT's to integral of data + // there is no other way to do this, because these quantities are not + // detector corrected + normalize(_pTjet1, 10439.0); + normalize(_pTjet2, 1461.5); + normalize(_pTjet3, 217.0); } //@} Modified: trunk/src/Analyses/D0_2008_S7554427.cc ============================================================================== --- trunk/src/Analyses/D0_2008_S7554427.cc Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/src/Analyses/D0_2008_S7554427.cc Wed Nov 4 10:58:19 2009 (r1994) @@ -21,7 +21,6 @@ { // Run II Z pT setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); } @@ -63,8 +62,8 @@ // Finalize void finalize() { - scale(_h_ZpT, crossSection()/sumOfWeights()); - scale(_h_forward_ZpT, crossSection()/sumOfWeights()); + normalize(_h_ZpT); + normalize(_h_forward_ZpT); } //@} Modified: trunk/src/Analyses/D0_2009_S8320160.cc ============================================================================== --- trunk/src/Analyses/D0_2009_S8320160.cc Wed Nov 4 01:41:01 2009 (r1993) +++ trunk/src/Analyses/D0_2009_S8320160.cc Wed Nov 4 10:58:19 2009 (r1994) @@ -21,7 +21,6 @@ D0_2009_S8320160() : Analysis("D0_2009_S8320160") { setBeams(PROTON, ANTIPROTON); - setNeedsCrossSection(true); } //@} @@ -74,7 +73,7 @@ /// Finalize void finalize() { foreach (AIDA::IHistogram1D* hist, _h_chi_dijet.getHistograms()) { - scale(hist, crossSection()/sumOfWeights()); + normalize(hist); } }
More information about the Rivet-svn mailing list |