|
|
[Rivet-svn] r2432 - in trunk: doc src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Sat May 1 12:42:11 BST 2010
Author: sonne
Date: Sat May 1 13:21:50 2010
New Revision: 2432
Log:
fill correction factors CDF/Ideal into hist + corresponding weight te get the right errors
Modified:
trunk/doc/analyses.tex
trunk/src/Analyses/CDF_1994_S2952106.cc
Modified: trunk/doc/analyses.tex
==============================================================================
--- trunk/doc/analyses.tex Fri Apr 30 22:16:16 2010 (r2431)
+++ trunk/doc/analyses.tex Sat May 1 13:21:50 2010 (r2432)
@@ -822,7 +822,7 @@
\typeout{Handling analysis CDF_2004_S5839831}
\subsection[CDF\_2004\_S5839831]{CDF\_2004\_S5839831\,\cite{Acosta:2004wqa}}
\textbf{Transverse cone and `Swiss cheese' underlying event studies}\newline
-\textbf{Experiment:} CDF (Tevatron Run 2) \newline
+\textbf{Experiment:} CDF (Tevatron Run 1) \newline
\textbf{Spires ID:} \href{http://www.slac.stanford.edu/spires/find/hep/www?rawcmd=key+5839831}{5839831}\newline
\textbf{Status:} VALIDATED\newline
\textbf{Authors:}
@@ -2179,7 +2179,7 @@
\penalty 100
\begin{itemize}
- \item $e^+ p$ deep inelastic scattering with $p$ at 820 GeV, $e^+$ at 27.5 GeV \ensuremath{\to} \ensuremath{\sqrt{s}} = 300 GeV\end{itemize}
+ \item $e^+ p$ deep inelastic scattering with $p$ at 820 GeV, $e^+$ at 27.5 GeV \ensuremath{\to} \ensuremath{\sqrt{s}} = 300 GeV\end{itemize}
\noindent Measurements of transverse energy flow for neutral current deep- inelastic scattering events produced in positron-proton collisions at HERA. The kinematic range covers squared momentum transfers $Q^2$ from 3.2 to 2200 GeV$^2$; the Bjorken scaling variable $x$ from $8 \times 10^{-5}$ to 0.11 and the hadronic mass $W$ from 66 to 233 GeV. The transverse energy flow is measured in the hadronic centre of mass frame and is studied as a function of $Q^2$, $x$, $W$ and pseudorapidity. The behaviour of the mean transverse energy in the central pseudorapidity region and an interval corresponding to the photon fragmentation region are analysed as a function of $Q^2$ and $W$. This analysis is useful for exploring the effect of photon PDFs and for tuning models of parton evolution and treatment of fragmentation and the proton remnant in DIS.
@@ -2212,7 +2212,7 @@
\item 820 GeV protons colliding with 27.5 GeV positrons; Direct and resolved photoproduction of dijets; Leading jet $\pT > 14$ GeV, second jet $\pT > 11$ GeV; Jet pseudorapidity $-1 < |\eta| < 2.4$\end{itemize}
-\noindent ZEUS photoproduction of jets from proton--positron collisions at beam energies of 820~GeV on 27.5~GeV. Photoproduction can either be direct, in which case the photon interacts directly with the parton, or resolved, in which case the photon acts as a source of quarks and gluons. A photon-proton centre of mass energy of between 134~GeV and 227~GeV is probed, with values of $x_P$, the fractional momentum of the partons inside the proton, predominantly in the region between 0.01 and 0.1. The fractional momentum of the partons from the photon, $x\gamma$, is in the region 0.1 to 1. Jets are reconstructed in the range $-1 < |\eta| < 2.4$ using the $k_\perp$ algorithm with an $R$ parameter of 1.0. The minimum \pT of the leading jet should be greater than 14~GeV, and at least one other jet must have $\pT > 11$~GeV.
+\noindent ZEUS photoproduction of jets from proton--positron collisions at beam energies of 820~GeV on 27.5~GeV. Photoproduction can either be direct, in which case the photon interacts directly with the parton, or resolved, in which case the photon acts as a source of quarks and gluons. A photon-proton centre of mass energy of between 134~GeV and 227~GeV is probed, with values of $x_P$, the fractional momentum of the partons inside the proton, predominantly in the region between 0.01 and 0.1. The fractional momentum of the partons from the photon, $x\gamma$, is in the region 0.1 to 1. Jets are reconstructed in the range $-1 < |\eta| < 2.4$ using the $k_\perp$ algorithm with an $R$ parameter of 1.0. The minimum \pT of the leading jet should be greater than 14~GeV, and at least one other jet must have $\pT > 11$~GeV.
\clearpage
@@ -2241,7 +2241,7 @@
\item pp at 200 GeV\end{itemize}
-\noindent \pT distributions of charged pions and (anti)protons in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC in non-single-diffractive minbias events.
+\noindent \pT distributions of charged pions and (anti)protons in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC in non-single-diffractive minbias events.
\clearpage
@@ -2271,7 +2271,7 @@
\item pp at 200 GeV\end{itemize}
-\noindent \pT distributions of identified strange particles in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC in non-single-diffractive minbias events. WARNING The $\langle \pT \rangle$ vs. particle mass plot is not validated yet and might be wrong.
+\noindent \pT distributions of identified strange particles in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC in non-single-diffractive minbias events. WARNING The $\langle \pT \rangle$ vs. particle mass plot is not validated yet and might be wrong.
\clearpage
@@ -2301,7 +2301,7 @@
\item pp at 200 GeV\end{itemize}
-\noindent Inclusive jet cross section as a function of \pT in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC.
+\noindent Inclusive jet cross section as a function of \pT in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC.
\clearpage
@@ -2363,7 +2363,7 @@
\item pp at 200 GeV\end{itemize}
-\noindent UE analysis similar to Rick Field's leading jet analysis. SIScone with radius/resolution parameter R=0.7 is used. Particles with $\pT > 0.2~\text{GeV}$ and $|\eta| < 1$ are included in the analysis. All particles are assumed to have zero mass. Only jets with neutral energy $< 0.7$ are included. For the transMIN and transMAX $\Delta(\phi)$ is between $\pi/3$ and $2\pi/3$, and $\Delta(\eta) < 2.0$. For the jet region the area of the jet is used for the normalization, i.e. the scaling factor is $\pi R^2$ and not $\mathrm{d}\phi\mathrm{d}\eta$ (this is different from what Rick Field does!). The tracking efficiency is $\sim 0.8$, but that is an approximation, as below $\pT \sim 0.6~\text{GeV}$ it is falling quite steeply.
+\noindent UE analysis similar to Rick Field's leading jet analysis. SIScone with radius/resolution parameter R=0.7 is used. Particles with $\pT > 0.2~\text{GeV}$ and $|\eta| < 1$ are included in the analysis. All particles are assumed to have zero mass. Only jets with neutral energy $< 0.7$ are included. For the transMIN and transMAX $\Delta(\phi)$ is between $\pi/3$ and $2\pi/3$, and $\Delta(\eta) < 2.0$. For the jet region the area of the jet is used for the normalization, i.e. the scaling factor is $\pi R^2$ and not $\mathrm{d}\phi\mathrm{d}\eta$ (this is different from what Rick Field does!). The tracking efficiency is $\sim 0.8$, but that is an approximation, as below $\pT \sim 0.6~\text{GeV}$ it is falling quite steeply.
\clearpage
@@ -2725,7 +2725,7 @@
\item QCD events, double-diffractive events should be turned on as well.\end{itemize}
-\noindent Charged multiplicities are measured at \ensuremath{\sqrt{s}} = 30.4, 44.5, 52.2 and 62.2 GeV using a minimum-bias trigger. The data is sub-divided into inleastic as well as non-single-diffractive events. However, the implementation of the diffractive events will require some work.
+\noindent Charged multiplicities are measured at \ensuremath{\sqrt{s}} = 30.4, 44.5, 52.2 and 62.2 GeV using a minimum-bias trigger. The data is sub-divided into inleastic as well as non-single-diffractive events. However, the implementation of the diffractive events will require some work.
\clearpage
Modified: trunk/src/Analyses/CDF_1994_S2952106.cc
==============================================================================
--- trunk/src/Analyses/CDF_1994_S2952106.cc Fri Apr 30 22:16:16 2010 (r2431)
+++ trunk/src/Analyses/CDF_1994_S2952106.cc Sat May 1 13:21:50 2010 (r2432)
@@ -56,8 +56,6 @@
_histJet3eta = bookHistogram1D(3,1,1);
//_histR23 = bookHistogram1D("R23", 50, 0., 5.);
_histR23 = bookHistogram1D(4,1,1);
-
- /// @todo Need better title
//_histAlpha = bookHistogram1D("alpha", 42, -PI/2., PI/2.);
_histAlpha = bookHistogram1D(5,1,1);
@@ -90,6 +88,22 @@
/// @todo Need better title
//_histJet3etaIdeal = bookHistogram1D("Jet3etaIdeal", 42, -4., 4.);
_histJet3etaIdeal = bookHistogram1D(8,1,2);
+
+
+ //_histEta3Corr = bookHistogram1D(3,1,1);
+ const string hnameEta3 = "Eta3Corr";
+ //const string htitleEta3 = "Eta3Corr";
+ _histEta3Corr = bookHistogram1D(hnameEta3, 40, -4., 4.);
+
+ //_histR23Corr = bookHistogram1D(4,1,1);
+ const string hnameR23 = "R23Corr";
+ //const string htitleR23 = "R23Corr";
+ _histR23Corr = bookHistogram1D(hnameR23, 35, 0., 4.375);
+
+ //_histAlphaCorr = bookHistogram1D(5,1,1);
+ const string hnameAlpha = "AlphaCorr";
+ //const string htitleAlpha = "AlphaCorr";
+ _histAlphaCorr = bookHistogram1D(hnameAlpha, 40, -90., 90.);
}
@@ -140,7 +154,75 @@
/// Finalize
void finalize() {
-
+
+ const double eta3_CDF_sim[] = {0.0013, 0.0037, 0.0047, 0.0071, 0.0093,
+ 0.0117, 0.0151, 0.0149, 0.0197, 0.0257,
+ 0.0344, 0.0409, 0.0481, 0.0454, 0.0394,
+ 0.0409, 0.0387, 0.0387, 0.0322, 0.0313,
+ 0.029, 0.0309, 0.0412, 0.0417, 0.0412,
+ 0.0397, 0.0417, 0.0414, 0.0376, 0.0316,
+ 0.027, 0.0186, 0.0186, 0.0132, 0.0127,
+ 0.0106, 0.0071, 0.004, 0.002, 0.0013};
+
+ const double eta3_CDF_sim_err[] = {0.0009, 0.0009, 0.0007, 0.0007, 0.0007,
+ 0.001, 0.0012, 0.0012, 0.0013, 0.0016,
+ 0.0017, 0.002, 0.002, 0.0022, 0.002,
+ 0.002, 0.0018, 0.0018, 0.0016, 0.0017,
+ 0.0017, 0.0019, 0.002, 0.0021, 0.002,
+ 0.002, 0.0019, 0.002, 0.0018, 0.0017,
+ 0.0017, 0.0014, 0.0014, 0.0009, 0.001,
+ 0.0009, 0.0009, 0.0008, 0.0008, 0.0009};
+
+ const double eta3_Ideal_sim[] = {0.0017, 0.003, 0.0033, 0.0062, 0.0062,
+ 0.0112, 0.0177, 0.0164, 0.0196, 0.0274,
+ 0.0351, 0.0413, 0.052, 0.0497, 0.0448,
+ 0.0446, 0.0375, 0.0329, 0.0291, 0.0272,
+ 0.0233, 0.0288, 0.0384, 0.0396, 0.0468,
+ 0.0419, 0.0459, 0.0399, 0.0355, 0.0329,
+ 0.0274, 0.023, 0.0201, 0.012, 0.01,
+ 0.008, 0.0051, 0.0051, 0.001, 0.001};
+
+ for (int ibin = 0; ibin < 40; ++ibin) //x-value + weight=err^2/y-value
+ _histEta3Corr->fill(-4. + (ibin+0.5)*8./40., //x-value
+ eta3_CDF_sim_err[ibin]/eta3_Ideal_sim[ibin] //error
+ / eta3_CDF_sim[ibin] * eta3_Ideal_sim[ibin] //norm. weight to value
+ ); //(Ideal error identical zero)
+
+
+
+ const double R23_CDF_sim[] = {0.0005, 0.0161, 0.057, 0.0762, 0.0723,
+ 0.0705, 0.0598, 0.0563, 0.0557, 0.0579,
+ 0.0538, 0.0522, 0.0486, 0.0449, 0.0418,
+ 0.0361, 0.0326, 0.0304, 0.0252, 0.0212,
+ 0.0173, 0.0176, 0.0145, 0.0127, 0.0103,
+ 0.0065, 0.0049, 0.0045, 0.0035, 0.0029,
+ 0.0024, 0.0014, 0.0011, 0.001, 0.0009};
+
+ const double R23_CDF_sim_err[] = {0.0013, 0.0009, 0.0022, 0.0029, 0.0026,
+ 0.0024, 0.0022, 0.0025, 0.0023, 0.0024,
+ 0.0021, 0.0021, 0.0021, 0.0021, 0.0021,
+ 0.0019, 0.0019, 0.0016, 0.0017, 0.0014,
+ 0.001, 0.0014, 0.0012, 0.0013, 0.001,
+ 0.0011, 0.001, 0.001, 0.001, 0.0011,
+ 0.0011, 0.0009, 0.0008, 0.0008, 0.0009};
+
+ const double R23_Ideal_sim[] = {0.0005, 0.0176, 0.0585, 0.0862, 0.0843,
+ 0.0756, 0.0673, 0.0635, 0.0586, 0.0619,
+ 0.0565, 0.0515, 0.0466, 0.0472, 0.0349,
+ 0.0349, 0.0266, 0.0254, 0.0204, 0.0179,
+ 0.0142, 0.0134, 0.0101, 0.009, 0.008,
+ 0.0034, 0.003, 0.0033, 0.0027, 0.0021,
+ 0.0012, 0.0006, 0.0004, 0.0005, 0.0003};
+
+ for (int ibin = 0; ibin < 35; ++ibin) //x-value + weight=err^2/y-value
+ _histR23Corr->fill((ibin+0.5)*4.375/35., //x-value
+ R23_CDF_sim_err[ibin]/R23_Ideal_sim[ibin] //error
+ / R23_CDF_sim[ibin] * R23_Ideal_sim[ibin] //norm. weight to y-value
+ ); //(Ideal error identical zero)
+
+
+
+
const double alpha_CDF_sim[] = {0.0517, 0.0461, 0.049, 0.0452, 0.0451,
0.0435, 0.0317, 0.0287, 0.0294, 0.0261,
0.0231, 0.022, 0.0233, 0.0192, 0.0213,
@@ -151,110 +233,48 @@
0.0269, 0.0257, 0.0276, 0.0246, 0.0286};
const double alpha_CDF_sim_err[] = {0.0024, 0.0025, 0.0024, 0.0024, 0.0024,
- 0.0022, 0.0019, 0.0018, 0.0019, 0.0016,
- 0.0017, 0.0017, 0.0019, 0.0013, 0.0017,
- 0.0014, 0.0016, 0.0013, 0.0012, 0.0009,
- 0.0014, 0.0014, 0.0014, 0.0014, 0.0014,
- 0.0015, 0.0014, 0.0016, 0.0016, 0.0015,
- 0.0016, 0.0016, 0.0019, 0.0017, 0.0019,
- 0.0018, 0.0018, 0.0018, 0.0018, 0.0019};
-
+ 0.0022, 0.0019, 0.0018, 0.0019, 0.0016,
+ 0.0017, 0.0017, 0.0019, 0.0013, 0.0017,
+ 0.0014, 0.0016, 0.0013, 0.0012, 0.0009,
+ 0.0014, 0.0014, 0.0014, 0.0014, 0.0014,
+ 0.0015, 0.0014, 0.0016, 0.0016, 0.0015,
+ 0.0016, 0.0016, 0.0019, 0.0017, 0.0019,
+ 0.0018, 0.0018, 0.0018, 0.0018, 0.0019};
+
const double alpha_Ideal_sim[] = {0.0552, 0.0558, 0.0583, 0.055, 0.0495,
- 0.0433, 0.0393, 0.0346, 0.0331, 0.0296,
- 0.0258, 0.0196, 0.0171, 0.0179, 0.0174,
- 0.0141, 0.0114, 0.0096, 0.0076, 0.0087,
- 0.0099, 0.0079, 0.0102, 0.0114, 0.0124,
- 0.013, 0.0165, 0.016, 0.0177, 0.019,
- 0.0232, 0.0243, 0.0238, 0.0248, 0.0235,
- 0.0298, 0.0292, 0.0291, 0.0268, 0.0316};
-
- for (int ibin = 0; ibin < 40; ++ibin) {
- _histAlphaCDF->fill(alpha_CDF_sim[ibin],
- alpha_CDF_sim_err[ibin]*alpha_CDF_sim_err[ibin]/alpha_CDF_sim[ibin]);
-
- }
-
- /*
- for (int ibin = 0; ibin < 40; ++ibin) {
- double raw = _histAlpha->GetBinContent(ibin);
- double corr = alpha_CDF_sim[ibin]/alpha_Ideal_sim[ibin];
- _histAlpha->SetBinContent(raw*corr);
- double rawerr = _histAlpha->GetBinError(ibin);
- double CDF_sim_err = alpha_CDF_sim_err[ibin];
- _histAlpha->SetBinError(ibin, sqrt(corr*corr*rawerr*rawerr +
- raw*raw/(alpha_Ideal_sim[ibin]*alpha_Ideal_sim[ibin])*
- alpha_CDF_sim_err[ibin]*alpha_CDF_sim_err[ibin]) );
+ 0.0433, 0.0393, 0.0346, 0.0331, 0.0296,
+ 0.0258, 0.0196, 0.0171, 0.0179, 0.0174,
+ 0.0141, 0.0114, 0.0096, 0.0076, 0.0087,
+ 0.0099, 0.0079, 0.0102, 0.0114, 0.0124,
+ 0.013, 0.0165, 0.016, 0.0177, 0.019,
+ 0.0232, 0.0243, 0.0238, 0.0248, 0.0235,
+ 0.0298, 0.0292, 0.0291, 0.0268, 0.0316};
+
+ for (int ibin = 0; ibin < 40; ++ibin) {//x-value + weight=err^2/y-value
+ _histAlphaCorr->fill(-90. + (ibin+0.5)*180./40., //y-value
+ alpha_CDF_sim_err[ibin]/alpha_Ideal_sim[ibin] //error
+ / alpha_CDF_sim[ibin] * alpha_Ideal_sim[ibin] //norm. weight to y-value
+ ); //(Ideal error identical zero)
+ cout << "bin i=" << ibin << " val=" << alpha_CDF_sim[ibin]/alpha_Ideal_sim[ibin]
+ << " err^2/value=" << alpha_CDF_sim_err[ibin]/alpha_Ideal_sim[ibin]
+ / alpha_CDF_sim[ibin] * alpha_Ideal_sim[ibin] << endl;
}
- */
-
- const double R23_CDF_sim[] = {0.0005, 0.0161, 0.057, 0.0762, 0.0723,
- 0.0705, 0.0598, 0.0563, 0.0557, 0.0579,
- 0.0538, 0.0522, 0.0486, 0.0449, 0.0418,
- 0.0361, 0.0326, 0.0304, 0.0252, 0.0212,
- 0.0173, 0.0176, 0.0145, 0.0127, 0.0103,
- 0.0065, 0.0049, 0.0045, 0.0035, 0.0029,
- 0.0024, 0.0014, 0.0011, 0.001, 0.0009};
-
- const double R23_CDF_sim_err[] = {0.0013, 0.0009, 0.0022, 0.0029, 0.0026,
- 0.0024, 0.0022, 0.0025, 0.0023, 0.0024,
- 0.0021, 0.0021, 0.0021, 0.0021, 0.0021,
- 0.0019, 0.0019, 0.0016, 0.0017, 0.0014,
- 0.001, 0.0014, 0.0012, 0.0013, 0.001,
- 0.0011, 0.001, 0.001, 0.001, 0.0011,
- 0.0011, 0.0009, 0.0008, 0.0008, 0.0009};
-
- const double R23_Ideal_sim[] = {0.0005, 0.0176, 0.0585, 0.0862, 0.0843,
- 0.0756, 0.0673, 0.0635, 0.0586, 0.0619,
- 0.0565, 0.0515, 0.0466, 0.0472, 0.0349,
- 0.0349, 0.0266, 0.0254, 0.0204, 0.0179,
- 0.0142, 0.0134, 0.0101, 0.009, 0.008,
- 0.0034, 0.003, 0.0033, 0.0027, 0.0021,
- 0.0012, 0.0006, 0.0004, 0.0005, 0.0003};
- const double eta3_CDF_sim[] = {0.0013, 0.0037, 0.0047, 0.0071, 0.0093,
- 0.0117, 0.0151, 0.0149, 0.0197, 0.0257,
- 0.0344, 0.0409, 0.0481, 0.0454, 0.0394,
- 0.0409, 0.0387, 0.0387, 0.0322, 0.0313,
- 0.029, 0.0309, 0.0412, 0.0417, 0.0412,
- 0.0397, 0.0417, 0.0414, 0.0376, 0.0316,
- 0.027, 0.0186, 0.0186, 0.0132, 0.0127,
- 0.0106, 0.0071, 0.004, 0.002, 0.0013};
-
- const double eta3_CDF_sim_err[] = {0.0009, 0.0009, 0.0007, 0.0007, 0.0007,
- 0.001, 0.0012, 0.0012, 0.0013, 0.0016,
- 0.0017, 0.002, 0.002, 0.0022, 0.002,
- 0.002, 0.0018, 0.0018, 0.0016, 0.0017,
- 0.0017, 0.0019, 0.002, 0.0021, 0.002,
- 0.002, 0.0019, 0.002, 0.0018, 0.0017,
- 0.0017, 0.0014, 0.0014, 0.0009, 0.001,
- 0.0009, 0.0009, 0.0008, 0.0008, 0.0009};
-
- const double eta3_Ideal_sim[] = {0.0017, 0.003, 0.0033, 0.0062, 0.0062,
- 0.0112, 0.0177, 0.0164, 0.0196, 0.0274,
- 0.0351, 0.0413, 0.052, 0.0497, 0.0448,
- 0.0446, 0.0375, 0.0329, 0.0291, 0.0272,
- 0.0233, 0.0288, 0.0384, 0.0396, 0.0468,
- 0.0419, 0.0459, 0.0399, 0.0355, 0.0329,
- 0.0274, 0.023, 0.0201, 0.012, 0.01,
- 0.008, 0.0051, 0.0051, 0.001, 0.001};
-
-
-
-
- /// @todo Apply correction
- // double a, b, c, erra, errb, errc;
- // for (int ibin = 0; ibin < _histAlpha->getNbins(); ++ibin) {
- // a = _histAlpha->GetBinContent(ibin);
- // erra = _histAlpha->GetBinError(ibin);
- // b = _histAlphaIdeal->GetBinContent(ibin);
- // errb = _histAlphaIdeal->GetBinError(ibin);
- // c = _histAlphaCDF->GetBinContent(ibin);
- // errc = _histAlphaCDF->GetBinError(ibin);
- // _histAlpha->SetBinContent(ibin, b/c);
- // _histAlpha->SetBinError(ibin, sqrt(sqr(b)/sqr(c)*sqr(erra) + sqr(a)/sqr(c)*sqr(errb) +
- // sqr(a*b/(sqr(c)))*sqr(errc) ) );
- // }
- /// @todo Same correction to be applied for _hisR23 and _histJet3eta histograms
+ AIDA::IHistogramFactory& hf = histogramFactory();
+
+ //hf.multiply(histoDir() + "/d03-x01-y01", *_histJet3eta, *_histEta3Corr);
+ hf.multiply("/_histJet3eta", *_histJet3eta, *_histEta3Corr);
+ //_histJet3eta = hf.multiply("/_histJet3eta", *_histJet3eta, *_histEta3Corr);
+ hf.destroy(_histEta3Corr);
+
+ //hf.multiply(histoDir() + "/d04-x01-y01", *_histR23, *_histR23Corr);
+ hf.multiply("/_histR23", *_histR23, *_histR23Corr);
+ hf.destroy(_histR23Corr);
+
+ //hf.multiply(histoDir() + "/d05-x01-y01", *_histAlpha, *_histAlphaCorr);
+ hf.multiply("/_histAlpha", *_histAlpha, *_histAlphaCorr);
+ hf.destroy(_histAlphaCorr);
+
//getLog() << Log::INFO << "Cross-section = " << crossSection()/picobarn << " pb" << endl;
normalize(_histJet1Et, 12.3025); //norm to integral of Ref data
@@ -278,13 +298,17 @@
AIDA::IHistogram1D* _histR23;
AIDA::IHistogram1D* _histJet3eta;
AIDA::IHistogram1D* _histAlpha;
- // AIDA::IHistogram1D* _histAlphaMCvsDat;
+
AIDA::IHistogram1D* _histAlphaIdeal;
AIDA::IHistogram1D* _histAlphaCDF;
AIDA::IHistogram1D* _histR23Ideal;
AIDA::IHistogram1D* _histR23CDF;
AIDA::IHistogram1D* _histJet3etaIdeal;
AIDA::IHistogram1D* _histJet3etaCDF;
+
+ AIDA::IHistogram1D* _histEta3Corr;
+ AIDA::IHistogram1D* _histR23Corr;
+ AIDA::IHistogram1D* _histAlphaCorr;
//@}
};
More information about the Rivet-svn
mailing list
|