|
|
[Rivet-svn] r4106 - in branches/2012-06-aidarivet: data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Wed Jan 23 10:38:58 GMT 2013
Author: fsiegert
Date: Wed Jan 23 10:38:58 2013
New Revision: 4106
Log:
Split MC_WJETS into three separate analyses: an inclusive one (MC_WINC), a jet property one (MC_WJETS) and a kt splitting scale one (MC_WKTSPLITTINGS).
Added:
branches/2012-06-aidarivet/data/anainfo/MC_WINC.info
- copied, changed from r4103, branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info
branches/2012-06-aidarivet/data/anainfo/MC_WKTSPLITTINGS.info
- copied, changed from r4103, branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info
branches/2012-06-aidarivet/data/plotinfo/MC_WINC.plot
- copied, changed from r4103, branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot
branches/2012-06-aidarivet/data/plotinfo/MC_WKTSPLITTINGS.plot
- copied, changed from r4103, branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot
branches/2012-06-aidarivet/src/Analyses/MC_WINC.cc
- copied, changed from r4103, branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc
branches/2012-06-aidarivet/src/Analyses/MC_WKTSPLITTINGS.cc
- copied, changed from r4103, branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc
Modified:
branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info
branches/2012-06-aidarivet/data/anainfo/Makefile.am
branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot
branches/2012-06-aidarivet/data/plotinfo/Makefile.am
branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc
branches/2012-06-aidarivet/src/Analyses/Makefile.am
Copied and modified: branches/2012-06-aidarivet/data/anainfo/MC_WINC.info (from r4103, branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info)
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info Tue Jan 22 15:07:30 2013 (r4103, copy source)
+++ branches/2012-06-aidarivet/data/anainfo/MC_WINC.info Wed Jan 23 10:38:58 2013 (r4106)
@@ -1,5 +1,5 @@
-Name: MC_WJETS
-Summary: Monte Carlo validation observables for $W[e \, \nu]$ + jets production
+Name: MC_WINC
+Summary: Monte Carlo validation observables for inclusive $W[e \, \nu]$ production
Status: VALIDATED
Authors:
- Frank Siegert <frank.siegert at cern.ch>
@@ -9,8 +9,5 @@
NumEvents: 1000000
PtCuts: [0]
Description:
- Available observables are W mass, pT of jets 1-4, jet multiplicity,
- $\Delta\eta(W, \text{jet1})$, $\Delta R(\text{jet2}, \text{jet3})$,
- differential jet rates 0->1, 1->2, 2->3, 3->4, integrated 0--4 jet
- rates.
+ Monte Carlo validation observables for inclusive $W[e \, \nu]$ production
NeedCrossSection: yes
Modified: branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info Wed Jan 23 10:37:49 2013 (r4105)
+++ branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info Wed Jan 23 10:38:58 2013 (r4106)
@@ -9,8 +9,5 @@
NumEvents: 1000000
PtCuts: [0]
Description:
- Available observables are W mass, pT of jets 1-4, jet multiplicity,
- $\Delta\eta(W, \text{jet1})$, $\Delta R(\text{jet2}, \text{jet3})$,
- differential jet rates 0->1, 1->2, 2->3, 3->4, integrated 0--4 jet
- rates.
+ Monte Carlo validation observables for $W[e \, \nu]$ + jets production
NeedCrossSection: yes
Copied and modified: branches/2012-06-aidarivet/data/anainfo/MC_WKTSPLITTINGS.info (from r4103, branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info)
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/MC_WJETS.info Tue Jan 22 15:07:30 2013 (r4103, copy source)
+++ branches/2012-06-aidarivet/data/anainfo/MC_WKTSPLITTINGS.info Wed Jan 23 10:38:58 2013 (r4106)
@@ -1,5 +1,5 @@
-Name: MC_WJETS
-Summary: Monte Carlo validation observables for $W[e \, \nu]$ + jets production
+Name: MC_WKTSPLITTINGS
+Summary: Monte Carlo validation observables for $k_\perp$ splitting scales in $W[e \, \nu]$ + jets events
Status: VALIDATED
Authors:
- Frank Siegert <frank.siegert at cern.ch>
@@ -9,8 +9,5 @@
NumEvents: 1000000
PtCuts: [0]
Description:
- Available observables are W mass, pT of jets 1-4, jet multiplicity,
- $\Delta\eta(W, \text{jet1})$, $\Delta R(\text{jet2}, \text{jet3})$,
- differential jet rates 0->1, 1->2, 2->3, 3->4, integrated 0--4 jet
- rates.
+ Monte Carlo validation observables for $k_\perp$ splitting scales in $W[e \, \nu]$ + jets events
NeedCrossSection: yes
Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/Makefile.am Wed Jan 23 10:37:49 2013 (r4105)
+++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Wed Jan 23 10:38:58 2013 (r4106)
@@ -192,7 +192,9 @@
MC_SUSY.info \
MC_TTBAR.info \
MC_VH2BB.info \
+ MC_WINC.info \
MC_WJETS.info \
+ MC_WKTSPLITTINGS.info \
MC_WPOL.info \
MC_WWJETS.info \
MC_XS.info \
Copied and modified: branches/2012-06-aidarivet/data/plotinfo/MC_WINC.plot (from r4103, branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot)
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot Tue Jan 22 15:07:30 2013 (r4103, copy source)
+++ branches/2012-06-aidarivet/data/plotinfo/MC_WINC.plot Wed Jan 23 10:38:58 2013 (r4106)
@@ -1,4 +1,4 @@
-# BEGIN PLOT /MC_WJETS/W_chargeasymm_eta
+# BEGIN PLOT /MC_WINC/W_chargeasymm_eta
Title=Charge asymmetry vs. $\eta(\ell)$
XLabel=$\eta(\ell)$ [GeV]
YLabel=$(\mathrm{d}\sigma(\ell^+)-\mathrm{d}\sigma(\ell^-))\;/\;(\mathrm{d}\sigma(\ell^+)+\mathrm{d}\sigma(\ell^-))$
@@ -6,7 +6,7 @@
Scale=1.0
# END PLOT
-# BEGIN PLOT /MC_WJETS/W_chargeasymm_pT
+# BEGIN PLOT /MC_WINC/W_chargeasymm_pT
Title=Charge asymmetry vs. $p_\perp(W)$
XLabel=$p_\perp(W)$ [GeV]
YLabel=$\mathrm{d}\sigma(W^+)\;/\;\mathrm{d}\sigma(W^-)$
@@ -17,27 +17,13 @@
YMax=3.0
# END PLOT
-# BEGIN PLOT /MC_WJETS/W_jet1_dR
-Title=Separation between W boson and leading jet
-XLabel=$\Delta R(\text{W, 1st jet})$
-YLabel=$\text{d}\sigma/\text{d}\Delta R(\text{W, 1st jet})$ [pb]
-LegendXPos=0.35
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_jet1_deta
-Title=
-XLabel=$\Delta{\eta}(\text{W, 1st jet})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{\eta}(\text{W, 1st jet})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_mass
+# BEGIN PLOT /MC_WINC/W_mass
Title=W mass
XLabel=$m_{\text{W}}$ [GeV]
YLabel=$\text{d}\sigma/\text{d}m_{\text{W}}$ [pb/GeV]
# END PLOT
-# BEGIN PLOT /MC_WJETS/W_pT$
+# BEGIN PLOT /MC_WINC/W_pT$
Title=W boson $p_\perp$
XLabel=$p_\perp^{\text{W}}$ [GeV]
YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
@@ -46,7 +32,7 @@
LegendYPos=0.5
# END PLOT
-# BEGIN PLOT /MC_WJETS/W_pT_peak
+# BEGIN PLOT /MC_WINC/W_pT_peak
Title=W boson $p_\perp$ in peak region
XLabel=$p_\perp^{\text{W}}$ [GeV]
YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
@@ -56,7 +42,7 @@
LegendXPos=0.30
# END PLOT
-# BEGIN PLOT /MC_WJETS/W_y
+# BEGIN PLOT /MC_WINC/W_y
Title=W boson rapidity
XLabel=$y_{\text{W}}$
YLabel=$\text{d}\sigma/\text{d}y_{\text{W}}$ [pb]
@@ -64,13 +50,13 @@
LegendYPos=0.5
# END PLOT
-# BEGIN PLOT /MC_WJETS/W_phi
+# BEGIN PLOT /MC_WINC/W_phi
Title=W azimuthal angle
XLabel=$\phi_{\text{W}}$
YLabel=$\text{d}\sigma/\text{d}\phi_{\text{W}}$ [pb]
# END PLOT
-# BEGIN PLOT /MC_WJETS/Wplus_pT$
+# BEGIN PLOT /MC_WINC/Wplus_pT$
Title=$p_\perp$ of $W^+$ boson
XLabel=$p_\perp^{W^+}$ [GeV]
YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^+}$ [pb/GeV]
@@ -79,7 +65,7 @@
LegendYPos=0.5
# END PLOT
-# BEGIN PLOT /MC_WJETS/Wminus_pT$
+# BEGIN PLOT /MC_WINC/Wminus_pT$
Title=$p_\perp$ of $W^-$ boson
XLabel=$p_\perp^{W^-}$ [GeV]
YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^-}$ [pb/GeV]
@@ -88,353 +74,16 @@
LegendYPos=0.5
# END PLOT
-# BEGIN PLOT /MC_WJETS/lepton_pT
+# BEGIN PLOT /MC_WINC/lepton_pT
Title=Lepton $p_\perp$
XLabel=$p_\perp^{\text{l}}$ [GeV]
YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{l}}$ [pb/GeV]
LogX=1
# END PLOT
-# BEGIN PLOT /MC_WJETS/lepton_eta
+# BEGIN PLOT /MC_WINC/lepton_eta
Title=Lepton $\eta$
XLabel=$\eta_{\text{l}}$
YLabel=$\text{d}\sigma/\text{d}\eta_{\text{l}}$ [pb]
# END PLOT
-# BEGIN PLOT /MC_WJETS/jets_dR_
-Title=$\Delta R$ separation between jets
-LegendXPos=0.10
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dR_12
-XLabel=$\Delta{R}(\text{jet 1, jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dR_13
-XLabel=$\Delta{R}(\text{jet 1, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dR_23
-XLabel=$\Delta{R}(\text{jet 2, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 2, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_
-Title=Pseudorapidity separation between jets
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_12
-XLabel=$\Delta\eta(\text{jet 1, jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_13
-XLabel=$\Delta\eta(\text{jet 1, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_23
-XLabel=$\Delta\eta(\text{jet 2, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 2, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_
-Title=Pseudorapidity separation between jets
-LegendXPos=0.1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_12
-XLabel=$\Delta\phi(\text{jet 1, jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 1, jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_13
-XLabel=$\Delta\phi(\text{jet 1, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 1, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_23
-XLabel=$\Delta\phi(\text{jet 2, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 2, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_1
-Title=Pseudorapidity of leading jet
-XLabel=$\eta(\text{jet 1})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 1})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_2
-Title=Pseudorapidity of second jet
-XLabel=$\eta(\text{jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 2})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_3
-Title=Pseudorapidity of third jet
-XLabel=$\eta(\text{jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 3})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_4
-Title=Pseudorapidity of fourth jet
-XLabel=$\eta(\text{jet 4})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 4})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_1
-Title=Pseudorapidity $+/-$ ratio of first jet
-XLabel=$\eta(\text{jet 1})_+/\eta(\text{jet 1})_-$
-YLabel=$|\eta(\text{jet 1}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_2
-Title=Pseudorapidity $+/-$ ratio of second jet
-XLabel=$\eta(\text{jet 2})_+/\eta(\text{jet 2})_-$
-YLabel=$|\eta(\text{jet 2}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_3
-Title=Pseudorapidity $+/-$ ratio of third jet
-XLabel=$\eta(\text{jet 3})_+/\eta(\text{jet 3})_-$
-YLabel=$|\eta(\text{jet 3}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_4
-Title=Pseudorapidity $+/-$ ratio of fourth jet
-XLabel=$\eta(\text{jet 4})_+/\eta(\text{jet 4})_-$
-YLabel=$|\eta(\text{jet 4}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_1
-Title=Rapidity of first jet
-XLabel=$y(\text{jet 1})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 1})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_2
-Title=Rapidity of second jet
-XLabel=$y(\text{jet 2})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_3
-Title=Rapidity of third jet
-XLabel=$y(\text{jet 3})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_4
-Title=Rapidity of fourth jet
-XLabel=$y(\text{jet 4})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 4})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_1
-Title=Rapidity $+/-$ ratio of first jet
-XLabel=$y(\text{jet 1})_+/y(\text{jet 1})_-$
-YLabel=$|y(\text{jet 1}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_2
-Title=Rapidity $+/-$ ratio of second jet
-XLabel=$y(\text{jet 2})_+/y(\text{jet 2})_-$
-YLabel=$|y(\text{jet 2}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_3
-Title=Rapidity $+/-$ ratio of third jet
-XLabel=$y(\text{jet 3})_+/y(\text{jet 3})_-$
-YLabel=$|y(\text{jet 3}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_4
-Title=Rapidity $+/-$ ratio of fourth jet
-XLabel=$y(\text{jet 4})_+/y(\text{jet 4})_-$
-YLabel=$|y(\text{jet 4}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_1
-Title=Mass of first jet
-XLabel=$m(\text{jet 1})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 1})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_2
-Title=Mass of second jet
-XLabel=$m(\text{jet 2})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 2})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_3
-Title=Mass of third jet
-XLabel=$m(\text{jet 3})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 3})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_4
-Title=Mass of fourth jet
-XLabel=$m(\text{jet 4})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 4})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_multi_exclusive
-Title=Exclusive jet multiplicity
-XLabel=$N_{\text{jet}}$
-YLabel=$\sigma(N_{\text{jet}})$ [pb]
-XMajorTickMarks=10
-XMinorTickMarks=0
-ErrorBands=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_multi_inclusive
-Title=Inclusive jet multiplicity
-XLabel=$N_{\text{jet}}$
-YLabel=$\sigma(\geq N_{\text{jet}})$ [pb]
-XMajorTickMarks=10
-XMinorTickMarks=0
-ErrorBands=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_multi_ratio
-Title=Ratio of jet multiplicity
-XLabel=$N_{\text{jet}}$
-YLabel=$\sigma(\geq N_{\text{jet}})/\sigma(\geq N_{\text{jet}}-1)$
-XMajorTickMarks=10
-XMinorTickMarks=0
-LogY=0
-ErrorBands=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_0
-Title=$\log_{10}$(Integrated $0$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{0}$
-Rebin=2
-LegendYPos=0.8
-LegendXPos=0.75
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_1
-Title=$\log_{10}$(Integrated $1$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{1}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_2
-Title=$\log_{10}$(Integrated $2$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{2}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_3
-Title=$\log_{10}$(Integrated $3$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{3}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_4
-Title=$\log_{10}$(Integrated $4$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{\geq4}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_01
-Title=$\log_{10}$($k_\perp$ jet resolution $0 \to 1$ [GeV])
-XLabel=$\log_{10}(d_{01}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{01}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_12
-Title=$\log_{10}$($k_\perp$ jet resolution $1 \to 2$ [GeV])
-XLabel=$\log_{10}(d_{12}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{12}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_23
-Title=$\log_{10}$($k_\perp$ jet resolution $2 \to 3$ [GeV])
-XLabel=$\log_{10}(d_{23}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{23}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_34
-Title=$\log_{10}$($k_\perp$ jet resolution $3 \to 4$ [GeV])
-XLabel=$\log_{10}(d_{34}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{34}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_1
-Title=Transverse momentum of leading jet
-XLabel=$p_\perp(\text{jet 1})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 1})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_2
-Title=Transverse momentum of second jet
-XLabel=$p_\perp(\text{jet 2})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 2})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_3
-Title=Transverse momentum of third jet
-XLabel=$p_\perp(\text{jet 3})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 3})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_4
-Title=Transverse momentum of fourth jet
-XLabel=$p_\perp(\text{jet 4})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 4})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_HT
-Title=Scalar sum of jet transverse momenta (HT)
-XLabel=$HT$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
-LogX=1
-# END PLOT
Modified: branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot Wed Jan 23 10:37:49 2013 (r4105)
+++ branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot Wed Jan 23 10:38:58 2013 (r4106)
@@ -1,22 +1,3 @@
-# BEGIN PLOT /MC_WJETS/W_chargeasymm_eta
-Title=Charge asymmetry vs. $\eta(\ell)$
-XLabel=$\eta(\ell)$ [GeV]
-YLabel=$(\mathrm{d}\sigma(\ell^+)-\mathrm{d}\sigma(\ell^-))\;/\;(\mathrm{d}\sigma(\ell^+)+\mathrm{d}\sigma(\ell^-))$
-LogY=0
-Scale=1.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_chargeasymm_pT
-Title=Charge asymmetry vs. $p_\perp(W)$
-XLabel=$p_\perp(W)$ [GeV]
-YLabel=$\mathrm{d}\sigma(W^+)\;/\;\mathrm{d}\sigma(W^-)$
-LogX=1
-LogY=0
-Scale=1.0
-YMin=0.0
-YMax=3.0
-# END PLOT
-
# BEGIN PLOT /MC_WJETS/W_jet1_dR
Title=Separation between W boson and leading jet
XLabel=$\Delta R(\text{W, 1st jet})$
@@ -31,76 +12,6 @@
YLabel=$\text{d}\sigma/\text{d}\Delta{\eta}(\text{W, 1st jet})$ [pb]
# END PLOT
-# BEGIN PLOT /MC_WJETS/W_mass
-Title=W mass
-XLabel=$m_{\text{W}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}m_{\text{W}}$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_pT$
-Title=W boson $p_\perp$
-XLabel=$p_\perp^{\text{W}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_pT_peak
-Title=W boson $p_\perp$ in peak region
-XLabel=$p_\perp^{\text{W}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
-LogY=0
-LogX=0
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_y
-Title=W boson rapidity
-XLabel=$y_{\text{W}}$
-YLabel=$\text{d}\sigma/\text{d}y_{\text{W}}$ [pb]
-LegendXPos=0.35
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_phi
-Title=W azimuthal angle
-XLabel=$\phi_{\text{W}}$
-YLabel=$\text{d}\sigma/\text{d}\phi_{\text{W}}$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/Wplus_pT$
-Title=$p_\perp$ of $W^+$ boson
-XLabel=$p_\perp^{W^+}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^+}$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/Wminus_pT$
-Title=$p_\perp$ of $W^-$ boson
-XLabel=$p_\perp^{W^-}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^-}$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/lepton_pT
-Title=Lepton $p_\perp$
-XLabel=$p_\perp^{\text{l}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{l}}$ [pb/GeV]
-LogX=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/lepton_eta
-Title=Lepton $\eta$
-XLabel=$\eta_{\text{l}}$
-YLabel=$\text{d}\sigma/\text{d}\eta_{\text{l}}$ [pb]
-# END PLOT
-
# BEGIN PLOT /MC_WJETS/jets_dR_
Title=$\Delta R$ separation between jets
LegendXPos=0.10
@@ -319,79 +230,6 @@
ErrorBands=1
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_R_0
-Title=$\log_{10}$(Integrated $0$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{0}$
-Rebin=2
-LegendYPos=0.8
-LegendXPos=0.75
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_1
-Title=$\log_{10}$(Integrated $1$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{1}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_2
-Title=$\log_{10}$(Integrated $2$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{2}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_3
-Title=$\log_{10}$(Integrated $3$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{3}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_4
-Title=$\log_{10}$(Integrated $4$ jet rate in $k_\perp$ [GeV])
-XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
-YLabel=$R_{\geq4}$
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_01
-Title=$\log_{10}$($k_\perp$ jet resolution $0 \to 1$ [GeV])
-XLabel=$\log_{10}(d_{01}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{01}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_12
-Title=$\log_{10}$($k_\perp$ jet resolution $1 \to 2$ [GeV])
-XLabel=$\log_{10}(d_{12}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{12}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_23
-Title=$\log_{10}$($k_\perp$ jet resolution $2 \to 3$ [GeV])
-XLabel=$\log_{10}(d_{23}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{23}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_d_34
-Title=$\log_{10}$($k_\perp$ jet resolution $3 \to 4$ [GeV])
-XLabel=$\log_{10}(d_{34}/\text{GeV})$
-YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{34}/\text{GeV})$ [pb]
-LegendXPos=0.15
-LegendYPos=0.5
-Rebin=2
-# END PLOT
-
# BEGIN PLOT /MC_WJETS/jet_pT_1
Title=Transverse momentum of leading jet
XLabel=$p_\perp(\text{jet 1})$ [GeV]
Copied and modified: branches/2012-06-aidarivet/data/plotinfo/MC_WKTSPLITTINGS.plot (from r4103, branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot)
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/MC_WJETS.plot Tue Jan 22 15:07:30 2013 (r4103, copy source)
+++ branches/2012-06-aidarivet/data/plotinfo/MC_WKTSPLITTINGS.plot Wed Jan 23 10:38:58 2013 (r4106)
@@ -1,325 +1,4 @@
-# BEGIN PLOT /MC_WJETS/W_chargeasymm_eta
-Title=Charge asymmetry vs. $\eta(\ell)$
-XLabel=$\eta(\ell)$ [GeV]
-YLabel=$(\mathrm{d}\sigma(\ell^+)-\mathrm{d}\sigma(\ell^-))\;/\;(\mathrm{d}\sigma(\ell^+)+\mathrm{d}\sigma(\ell^-))$
-LogY=0
-Scale=1.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_chargeasymm_pT
-Title=Charge asymmetry vs. $p_\perp(W)$
-XLabel=$p_\perp(W)$ [GeV]
-YLabel=$\mathrm{d}\sigma(W^+)\;/\;\mathrm{d}\sigma(W^-)$
-LogX=1
-LogY=0
-Scale=1.0
-YMin=0.0
-YMax=3.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_jet1_dR
-Title=Separation between W boson and leading jet
-XLabel=$\Delta R(\text{W, 1st jet})$
-YLabel=$\text{d}\sigma/\text{d}\Delta R(\text{W, 1st jet})$ [pb]
-LegendXPos=0.35
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_jet1_deta
-Title=
-XLabel=$\Delta{\eta}(\text{W, 1st jet})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{\eta}(\text{W, 1st jet})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_mass
-Title=W mass
-XLabel=$m_{\text{W}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}m_{\text{W}}$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_pT$
-Title=W boson $p_\perp$
-XLabel=$p_\perp^{\text{W}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_pT_peak
-Title=W boson $p_\perp$ in peak region
-XLabel=$p_\perp^{\text{W}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{W}}$ [pb/GeV]
-LogY=0
-LogX=0
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_y
-Title=W boson rapidity
-XLabel=$y_{\text{W}}$
-YLabel=$\text{d}\sigma/\text{d}y_{\text{W}}$ [pb]
-LegendXPos=0.35
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/W_phi
-Title=W azimuthal angle
-XLabel=$\phi_{\text{W}}$
-YLabel=$\text{d}\sigma/\text{d}\phi_{\text{W}}$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/Wplus_pT$
-Title=$p_\perp$ of $W^+$ boson
-XLabel=$p_\perp^{W^+}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^+}$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/Wminus_pT$
-Title=$p_\perp$ of $W^-$ boson
-XLabel=$p_\perp^{W^-}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{W^-}$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/lepton_pT
-Title=Lepton $p_\perp$
-XLabel=$p_\perp^{\text{l}}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{\text{l}}$ [pb/GeV]
-LogX=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/lepton_eta
-Title=Lepton $\eta$
-XLabel=$\eta_{\text{l}}$
-YLabel=$\text{d}\sigma/\text{d}\eta_{\text{l}}$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dR_
-Title=$\Delta R$ separation between jets
-LegendXPos=0.10
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dR_12
-XLabel=$\Delta{R}(\text{jet 1, jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dR_13
-XLabel=$\Delta{R}(\text{jet 1, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 1, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dR_23
-XLabel=$\Delta{R}(\text{jet 2, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{R}(\text{jet 2, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_
-Title=Pseudorapidity separation between jets
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_12
-XLabel=$\Delta\eta(\text{jet 1, jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_13
-XLabel=$\Delta\eta(\text{jet 1, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 1, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_deta_23
-XLabel=$\Delta\eta(\text{jet 2, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\eta(\text{jet 2, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_
-Title=Pseudorapidity separation between jets
-LegendXPos=0.1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_12
-XLabel=$\Delta\phi(\text{jet 1, jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 1, jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_13
-XLabel=$\Delta\phi(\text{jet 1, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 1, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jets_dphi_23
-XLabel=$\Delta\phi(\text{jet 2, jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\Delta\phi(\text{jet 2, jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_1
-Title=Pseudorapidity of leading jet
-XLabel=$\eta(\text{jet 1})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 1})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_2
-Title=Pseudorapidity of second jet
-XLabel=$\eta(\text{jet 2})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 2})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_3
-Title=Pseudorapidity of third jet
-XLabel=$\eta(\text{jet 3})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 3})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_4
-Title=Pseudorapidity of fourth jet
-XLabel=$\eta(\text{jet 4})$
-YLabel=$\text{d}\sigma/\text{d}\eta(\text{jet 4})$ [pb]
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_1
-Title=Pseudorapidity $+/-$ ratio of first jet
-XLabel=$\eta(\text{jet 1})_+/\eta(\text{jet 1})_-$
-YLabel=$|\eta(\text{jet 1}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_2
-Title=Pseudorapidity $+/-$ ratio of second jet
-XLabel=$\eta(\text{jet 2})_+/\eta(\text{jet 2})_-$
-YLabel=$|\eta(\text{jet 2}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_3
-Title=Pseudorapidity $+/-$ ratio of third jet
-XLabel=$\eta(\text{jet 3})_+/\eta(\text{jet 3})_-$
-YLabel=$|\eta(\text{jet 3}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_eta_pmratio_4
-Title=Pseudorapidity $+/-$ ratio of fourth jet
-XLabel=$\eta(\text{jet 4})_+/\eta(\text{jet 4})_-$
-YLabel=$|\eta(\text{jet 4}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_1
-Title=Rapidity of first jet
-XLabel=$y(\text{jet 1})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 1})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_2
-Title=Rapidity of second jet
-XLabel=$y(\text{jet 2})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 2})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_3
-Title=Rapidity of third jet
-XLabel=$y(\text{jet 3})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 3})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_4
-Title=Rapidity of fourth jet
-XLabel=$y(\text{jet 4})$
-YLabel=$\text{d}\sigma/\text{d}y(\text{jet 4})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_1
-Title=Rapidity $+/-$ ratio of first jet
-XLabel=$y(\text{jet 1})_+/y(\text{jet 1})_-$
-YLabel=$|y(\text{jet 1}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_2
-Title=Rapidity $+/-$ ratio of second jet
-XLabel=$y(\text{jet 2})_+/y(\text{jet 2})_-$
-YLabel=$|y(\text{jet 2}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_3
-Title=Rapidity $+/-$ ratio of third jet
-XLabel=$y(\text{jet 3})_+/y(\text{jet 3})_-$
-YLabel=$|y(\text{jet 3}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_y_pmratio_4
-Title=Rapidity $+/-$ ratio of fourth jet
-XLabel=$y(\text{jet 4})_+/y(\text{jet 4})_-$
-YLabel=$|y(\text{jet 4}|$
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_1
-Title=Mass of first jet
-XLabel=$m(\text{jet 1})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 1})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_2
-Title=Mass of second jet
-XLabel=$m(\text{jet 2})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 2})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_3
-Title=Mass of third jet
-XLabel=$m(\text{jet 3})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 3})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_mass_4
-Title=Mass of fourth jet
-XLabel=$m(\text{jet 4})$
-YLabel=$\text{d}\sigma/\text{d}m(\text{jet 4})$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_multi_exclusive
-Title=Exclusive jet multiplicity
-XLabel=$N_{\text{jet}}$
-YLabel=$\sigma(N_{\text{jet}})$ [pb]
-XMajorTickMarks=10
-XMinorTickMarks=0
-ErrorBands=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_multi_inclusive
-Title=Inclusive jet multiplicity
-XLabel=$N_{\text{jet}}$
-YLabel=$\sigma(\geq N_{\text{jet}})$ [pb]
-XMajorTickMarks=10
-XMinorTickMarks=0
-ErrorBands=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_multi_ratio
-Title=Ratio of jet multiplicity
-XLabel=$N_{\text{jet}}$
-YLabel=$\sigma(\geq N_{\text{jet}})/\sigma(\geq N_{\text{jet}}-1)$
-XMajorTickMarks=10
-XMinorTickMarks=0
-LogY=0
-ErrorBands=1
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/log10_R_0
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_R_0
Title=$\log_{10}$(Integrated $0$ jet rate in $k_\perp$ [GeV])
XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
YLabel=$R_{0}$
@@ -328,35 +7,35 @@
LegendXPos=0.75
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_R_1
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_R_1
Title=$\log_{10}$(Integrated $1$ jet rate in $k_\perp$ [GeV])
XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
YLabel=$R_{1}$
Rebin=2
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_R_2
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_R_2
Title=$\log_{10}$(Integrated $2$ jet rate in $k_\perp$ [GeV])
XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
YLabel=$R_{2}$
Rebin=2
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_R_3
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_R_3
Title=$\log_{10}$(Integrated $3$ jet rate in $k_\perp$ [GeV])
XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
YLabel=$R_{3}$
Rebin=2
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_R_4
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_R_4
Title=$\log_{10}$(Integrated $4$ jet rate in $k_\perp$ [GeV])
XLabel=$\log_{10}(d_{\text{cut}}/\text{GeV})$
YLabel=$R_{\geq4}$
Rebin=2
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_d_01
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_d_01
Title=$\log_{10}$($k_\perp$ jet resolution $0 \to 1$ [GeV])
XLabel=$\log_{10}(d_{01}/\text{GeV})$
YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{01}/\text{GeV})$ [pb]
@@ -365,7 +44,7 @@
Rebin=2
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_d_12
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_d_12
Title=$\log_{10}$($k_\perp$ jet resolution $1 \to 2$ [GeV])
XLabel=$\log_{10}(d_{12}/\text{GeV})$
YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{12}/\text{GeV})$ [pb]
@@ -374,7 +53,7 @@
Rebin=2
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_d_23
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_d_23
Title=$\log_{10}$($k_\perp$ jet resolution $2 \to 3$ [GeV])
XLabel=$\log_{10}(d_{23}/\text{GeV})$
YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{23}/\text{GeV})$ [pb]
@@ -383,7 +62,7 @@
Rebin=2
# END PLOT
-# BEGIN PLOT /MC_WJETS/log10_d_34
+# BEGIN PLOT /MC_WKTSPLITTINGS/log10_d_34
Title=$\log_{10}$($k_\perp$ jet resolution $3 \to 4$ [GeV])
XLabel=$\log_{10}(d_{34}/\text{GeV})$
YLabel=$\text{d}\sigma/\text{d}\log_{10}(d_{34}/\text{GeV})$ [pb]
@@ -391,50 +70,3 @@
LegendYPos=0.5
Rebin=2
# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_1
-Title=Transverse momentum of leading jet
-XLabel=$p_\perp(\text{jet 1})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 1})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_2
-Title=Transverse momentum of second jet
-XLabel=$p_\perp(\text{jet 2})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 2})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_3
-Title=Transverse momentum of third jet
-XLabel=$p_\perp(\text{jet 3})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 3})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_pT_4
-Title=Transverse momentum of fourth jet
-XLabel=$p_\perp(\text{jet 4})$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp(\text{jet 4})$ [pb/GeV]
-LogX=1
-LegendXPos=0.05
-LegendYPos=0.5
-XMin=20.0
-# END PLOT
-
-# BEGIN PLOT /MC_WJETS/jet_HT
-Title=Scalar sum of jet transverse momenta (HT)
-XLabel=$HT$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
-LogX=1
-# END PLOT
Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Wed Jan 23 10:37:49 2013 (r4105)
+++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Wed Jan 23 10:38:58 2013 (r4106)
@@ -188,7 +188,9 @@
MC_VH2BB.plot \
MC_HJETS.plot \
MC_PDFS.plot \
+ MC_WINC.plot \
MC_WJETS.plot \
+ MC_WKTSPLITTINGS.plot \
MC_WPOL.plot \
MC_WWJETS.plot \
MC_XS.plot \
Copied and modified: branches/2012-06-aidarivet/src/Analyses/MC_WINC.cc (from r4103, branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc)
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc Tue Jan 22 15:07:30 2013 (r4103, copy source)
+++ branches/2012-06-aidarivet/src/Analyses/MC_WINC.cc Wed Jan 23 10:38:58 2013 (r4106)
@@ -9,13 +9,13 @@
namespace Rivet {
- /// @brief MC validation analysis for W + jets events
- class MC_WJETS : public MC_JetAnalysis {
+ /// @brief MC validation analysis for inclusive W events
+ class MC_WINC : public Analysis {
public:
/// Default constructor
- MC_WJETS()
- : MC_JetAnalysis("MC_WJETS", 4, "Jets")
+ MC_WINC()
+ : Analysis("MC_WINC")
{ }
@@ -27,35 +27,30 @@
FinalState fs;
WFinder wfinder(fs, -3.5, 3.5, 25.0*GeV, ELECTRON, 60.0*GeV, 100.0*GeV, 25.0*GeV, 0.2);
addProjection(wfinder, "WFinder");
- FastJets jetpro(wfinder.remainingFinalState(), FastJets::KT, 0.7);
- addProjection(jetpro, "Jets");
_h_W_mass = bookHistogram1D("W_mass", 50, 55.0, 105.0);
_h_W_pT = bookHistogram1D("W_pT", logBinEdges(100, 1.0, 0.5*sqrtS()));
_h_W_pT_peak = bookHistogram1D("W_pT_peak", 25, 0.0, 125.0);
_h_W_y = bookHistogram1D("W_y", 40, -4.0, 4.0);
_h_W_phi = bookHistogram1D("W_phi", 25, 0.0, TWOPI);
- _h_W_jet1_deta = bookHistogram1D("W_jet1_deta", 50, -5.0, 5.0);
- _h_W_jet1_dR = bookHistogram1D("W_jet1_dR", 25, 0.5, 7.0);
_h_Wplus_pT = bookHistogram1D("Wplus_pT", logBinEdges(25, 10.0, 0.5*sqrtS()));
_h_Wminus_pT = bookHistogram1D("Wminus_pT", logBinEdges(25, 10.0, 0.5*sqrtS()));
_h_lepton_pT = bookHistogram1D("lepton_pT", logBinEdges(100, 10.0, 0.25*sqrtS()));
_h_lepton_eta = bookHistogram1D("lepton_eta", 40, -4.0, 4.0);
_htmp_dsigminus_deta = bookHistogram1D("lepton_dsigminus_deta", 20, 0.0, 4.0);
_htmp_dsigplus_deta = bookHistogram1D("lepton_dsigplus_deta", 20, 0.0, 4.0);
-
- MC_JetAnalysis::init();
}
/// Do the analysis
void analyze(const Event & e) {
+ const double weight = e.weight();
+
const WFinder& wfinder = applyProjection<WFinder>(e, "WFinder");
if (wfinder.bosons().size() != 1) {
vetoEvent;
}
- const double weight = e.weight();
double charge3_x_eta = 0;
int charge3 = 0;
@@ -88,15 +83,6 @@
} else {
_h_Wplus_pT->fill(wmom.pT(), weight);
}
-
- const FastJets& jetpro = applyProjection<FastJets>(e, "Jets");
- const Jets& jets = jetpro.jetsByPt(20.0*GeV);
- if (jets.size() > 0) {
- _h_W_jet1_deta->fill(wmom.eta()-jets[0].momentum().eta(), weight);
- _h_W_jet1_dR->fill(deltaR(wmom, jets[0].momentum()), weight);
- }
-
- MC_JetAnalysis::analyze(e);
}
@@ -107,8 +93,6 @@
scale(_h_W_pT_peak, crossSection()/sumOfWeights());
scale(_h_W_y, crossSection()/sumOfWeights());
scale(_h_W_phi, crossSection()/sumOfWeights());
- scale(_h_W_jet1_deta, crossSection()/sumOfWeights());
- scale(_h_W_jet1_dR, crossSection()/sumOfWeights());
scale(_h_lepton_pT, crossSection()/sumOfWeights());
scale(_h_lepton_eta, crossSection()/sumOfWeights());
@@ -127,8 +111,6 @@
hf.divide(histoDir() + "/W_chargeasymm_pT", *_h_Wplus_pT, *_h_Wminus_pT);
scale(_h_Wplus_pT, crossSection()/sumOfWeights());
scale(_h_Wminus_pT, crossSection()/sumOfWeights());
-
- MC_JetAnalysis::finalize();
}
//@}
@@ -143,8 +125,6 @@
AIDA::IHistogram1D * _h_W_pT_peak;
AIDA::IHistogram1D * _h_W_y;
AIDA::IHistogram1D * _h_W_phi;
- AIDA::IHistogram1D * _h_W_jet1_deta;
- AIDA::IHistogram1D * _h_W_jet1_dR;
AIDA::IHistogram1D * _h_Wplus_pT;
AIDA::IHistogram1D * _h_Wminus_pT;
AIDA::IHistogram1D * _h_lepton_pT;
@@ -159,6 +139,6 @@
// The hook for the plugin system
- DECLARE_RIVET_PLUGIN(MC_WJETS);
+ DECLARE_RIVET_PLUGIN(MC_WINC);
}
Modified: branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc Wed Jan 23 10:37:49 2013 (r4105)
+++ branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc Wed Jan 23 10:38:58 2013 (r4106)
@@ -27,22 +27,11 @@
FinalState fs;
WFinder wfinder(fs, -3.5, 3.5, 25.0*GeV, ELECTRON, 60.0*GeV, 100.0*GeV, 25.0*GeV, 0.2);
addProjection(wfinder, "WFinder");
- FastJets jetpro(wfinder.remainingFinalState(), FastJets::KT, 0.7);
+ FastJets jetpro(wfinder.remainingFinalState(), FastJets::ANTIKT, 0.4);
addProjection(jetpro, "Jets");
- _h_W_mass = bookHistogram1D("W_mass", 50, 55.0, 105.0);
- _h_W_pT = bookHistogram1D("W_pT", logBinEdges(100, 1.0, 0.5*sqrtS()));
- _h_W_pT_peak = bookHistogram1D("W_pT_peak", 25, 0.0, 125.0);
- _h_W_y = bookHistogram1D("W_y", 40, -4.0, 4.0);
- _h_W_phi = bookHistogram1D("W_phi", 25, 0.0, TWOPI);
_h_W_jet1_deta = bookHistogram1D("W_jet1_deta", 50, -5.0, 5.0);
_h_W_jet1_dR = bookHistogram1D("W_jet1_dR", 25, 0.5, 7.0);
- _h_Wplus_pT = bookHistogram1D("Wplus_pT", logBinEdges(25, 10.0, 0.5*sqrtS()));
- _h_Wminus_pT = bookHistogram1D("Wminus_pT", logBinEdges(25, 10.0, 0.5*sqrtS()));
- _h_lepton_pT = bookHistogram1D("lepton_pT", logBinEdges(100, 10.0, 0.25*sqrtS()));
- _h_lepton_eta = bookHistogram1D("lepton_eta", 40, -4.0, 4.0);
- _htmp_dsigminus_deta = bookHistogram1D("lepton_dsigminus_deta", 20, 0.0, 4.0);
- _htmp_dsigplus_deta = bookHistogram1D("lepton_dsigplus_deta", 20, 0.0, 4.0);
MC_JetAnalysis::init();
}
@@ -51,46 +40,16 @@
/// Do the analysis
void analyze(const Event & e) {
+ const double weight = e.weight();
+
const WFinder& wfinder = applyProjection<WFinder>(e, "WFinder");
if (wfinder.bosons().size() != 1) {
vetoEvent;
}
- const double weight = e.weight();
-
- double charge3_x_eta = 0;
- int charge3 = 0;
- FourMomentum emom;
FourMomentum wmom(wfinder.bosons().front().momentum());
- _h_W_mass->fill(wmom.mass(), weight);
- _h_W_pT->fill(wmom.pT(), weight);
- _h_W_pT_peak->fill(wmom.pT(), weight);
- _h_W_y->fill(wmom.rapidity(), weight);
- _h_W_phi->fill(wmom.azimuthalAngle(), weight);
- Particle l=wfinder.constituentLeptons()[0];
- _h_lepton_pT->fill(l.momentum().pT(), weight);
- _h_lepton_eta->fill(l.momentum().eta(), weight);
- if (PID::threeCharge(l.pdgId()) != 0) {
- emom = l.momentum();
- charge3_x_eta = PID::threeCharge(l.pdgId()) * emom.eta();
- charge3 = PID::threeCharge(l.pdgId());
- }
- assert(charge3_x_eta != 0);
- assert(charge3!=0);
- if (emom.Et() > 30/GeV) {
- if (charge3_x_eta < 0) {
- _htmp_dsigminus_deta->fill(emom.eta(), weight);
- } else {
- _htmp_dsigplus_deta->fill(emom.eta(), weight);
- }
- }
- if (charge3 < 0) {
- _h_Wminus_pT->fill(wmom.pT(), weight);
- } else {
- _h_Wplus_pT->fill(wmom.pT(), weight);
- }
const FastJets& jetpro = applyProjection<FastJets>(e, "Jets");
- const Jets& jets = jetpro.jetsByPt(20.0*GeV);
+ const Jets& jets = jetpro.jetsByPt(m_jetptcut);
if (jets.size() > 0) {
_h_W_jet1_deta->fill(wmom.eta()-jets[0].momentum().eta(), weight);
_h_W_jet1_dR->fill(deltaR(wmom, jets[0].momentum()), weight);
@@ -102,31 +61,8 @@
/// Finalize
void finalize() {
- scale(_h_W_mass, crossSection()/sumOfWeights());
- scale(_h_W_pT, crossSection()/sumOfWeights());
- scale(_h_W_pT_peak, crossSection()/sumOfWeights());
- scale(_h_W_y, crossSection()/sumOfWeights());
- scale(_h_W_phi, crossSection()/sumOfWeights());
scale(_h_W_jet1_deta, crossSection()/sumOfWeights());
scale(_h_W_jet1_dR, crossSection()/sumOfWeights());
- scale(_h_lepton_pT, crossSection()/sumOfWeights());
- scale(_h_lepton_eta, crossSection()/sumOfWeights());
-
- // Construct asymmetry: (dsig+/deta - dsig-/deta) / (dsig+/deta + dsig-/deta) for each Et region
- AIDA::IHistogramFactory& hf = histogramFactory();
- IHistogram1D* numtmp = hf.subtract("/numtmp", *_htmp_dsigplus_deta, *_htmp_dsigminus_deta);
- IHistogram1D* dentmp = hf.add("/dentmp", *_htmp_dsigplus_deta, *_htmp_dsigminus_deta);
- assert(numtmp && dentmp);
- hf.divide(histoDir() + "/W_chargeasymm_eta", *numtmp, *dentmp);
- hf.destroy(numtmp);
- hf.destroy(dentmp);
- hf.destroy(_htmp_dsigminus_deta);
- hf.destroy(_htmp_dsigplus_deta);
-
- // W charge asymmetry vs. pTW: dsig+/dpT / dsig-/dpT
- hf.divide(histoDir() + "/W_chargeasymm_pT", *_h_Wplus_pT, *_h_Wminus_pT);
- scale(_h_Wplus_pT, crossSection()/sumOfWeights());
- scale(_h_Wminus_pT, crossSection()/sumOfWeights());
MC_JetAnalysis::finalize();
}
@@ -138,20 +74,8 @@
/// @name Histograms
//@{
- AIDA::IHistogram1D * _h_W_mass;
- AIDA::IHistogram1D * _h_W_pT;
- AIDA::IHistogram1D * _h_W_pT_peak;
- AIDA::IHistogram1D * _h_W_y;
- AIDA::IHistogram1D * _h_W_phi;
AIDA::IHistogram1D * _h_W_jet1_deta;
AIDA::IHistogram1D * _h_W_jet1_dR;
- AIDA::IHistogram1D * _h_Wplus_pT;
- AIDA::IHistogram1D * _h_Wminus_pT;
- AIDA::IHistogram1D * _h_lepton_pT;
- AIDA::IHistogram1D * _h_lepton_eta;
-
- AIDA::IHistogram1D * _htmp_dsigminus_deta;
- AIDA::IHistogram1D * _htmp_dsigplus_deta;
//@}
};
Copied and modified: branches/2012-06-aidarivet/src/Analyses/MC_WKTSPLITTINGS.cc (from r4103, branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc)
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/MC_WJETS.cc Tue Jan 22 15:07:30 2013 (r4103, copy source)
+++ branches/2012-06-aidarivet/src/Analyses/MC_WKTSPLITTINGS.cc Wed Jan 23 10:38:58 2013 (r4106)
@@ -1,5 +1,5 @@
// -*- C++ -*-
-#include "Rivet/Analyses/MC_JetAnalysis.hh"
+#include "Rivet/Analyses/MC_JetSplittings.hh"
#include "Rivet/Projections/WFinder.hh"
#include "Rivet/Projections/FastJets.hh"
#include "Rivet/Tools/Logging.hh"
@@ -9,13 +9,13 @@
namespace Rivet {
- /// @brief MC validation analysis for W + jets events
- class MC_WJETS : public MC_JetAnalysis {
+ /// @brief MC validation analysis for kt splitting scales in W + jets events
+ class MC_WKTSPLITTINGS : public MC_JetSplittings {
public:
/// Default constructor
- MC_WJETS()
- : MC_JetAnalysis("MC_WJETS", 4, "Jets")
+ MC_WKTSPLITTINGS()
+ : MC_JetSplittings("MC_WKTSPLITTINGS", 4, "Jets")
{ }
@@ -27,138 +27,35 @@
FinalState fs;
WFinder wfinder(fs, -3.5, 3.5, 25.0*GeV, ELECTRON, 60.0*GeV, 100.0*GeV, 25.0*GeV, 0.2);
addProjection(wfinder, "WFinder");
- FastJets jetpro(wfinder.remainingFinalState(), FastJets::KT, 0.7);
+ FastJets jetpro(wfinder.remainingFinalState(), FastJets::KT, 0.6);
addProjection(jetpro, "Jets");
- _h_W_mass = bookHistogram1D("W_mass", 50, 55.0, 105.0);
- _h_W_pT = bookHistogram1D("W_pT", logBinEdges(100, 1.0, 0.5*sqrtS()));
- _h_W_pT_peak = bookHistogram1D("W_pT_peak", 25, 0.0, 125.0);
- _h_W_y = bookHistogram1D("W_y", 40, -4.0, 4.0);
- _h_W_phi = bookHistogram1D("W_phi", 25, 0.0, TWOPI);
- _h_W_jet1_deta = bookHistogram1D("W_jet1_deta", 50, -5.0, 5.0);
- _h_W_jet1_dR = bookHistogram1D("W_jet1_dR", 25, 0.5, 7.0);
- _h_Wplus_pT = bookHistogram1D("Wplus_pT", logBinEdges(25, 10.0, 0.5*sqrtS()));
- _h_Wminus_pT = bookHistogram1D("Wminus_pT", logBinEdges(25, 10.0, 0.5*sqrtS()));
- _h_lepton_pT = bookHistogram1D("lepton_pT", logBinEdges(100, 10.0, 0.25*sqrtS()));
- _h_lepton_eta = bookHistogram1D("lepton_eta", 40, -4.0, 4.0);
- _htmp_dsigminus_deta = bookHistogram1D("lepton_dsigminus_deta", 20, 0.0, 4.0);
- _htmp_dsigplus_deta = bookHistogram1D("lepton_dsigplus_deta", 20, 0.0, 4.0);
-
- MC_JetAnalysis::init();
+ MC_JetSplittings::init();
}
-
/// Do the analysis
void analyze(const Event & e) {
+
const WFinder& wfinder = applyProjection<WFinder>(e, "WFinder");
if (wfinder.bosons().size() != 1) {
vetoEvent;
}
- const double weight = e.weight();
-
- double charge3_x_eta = 0;
- int charge3 = 0;
- FourMomentum emom;
- FourMomentum wmom(wfinder.bosons().front().momentum());
- _h_W_mass->fill(wmom.mass(), weight);
- _h_W_pT->fill(wmom.pT(), weight);
- _h_W_pT_peak->fill(wmom.pT(), weight);
- _h_W_y->fill(wmom.rapidity(), weight);
- _h_W_phi->fill(wmom.azimuthalAngle(), weight);
- Particle l=wfinder.constituentLeptons()[0];
- _h_lepton_pT->fill(l.momentum().pT(), weight);
- _h_lepton_eta->fill(l.momentum().eta(), weight);
- if (PID::threeCharge(l.pdgId()) != 0) {
- emom = l.momentum();
- charge3_x_eta = PID::threeCharge(l.pdgId()) * emom.eta();
- charge3 = PID::threeCharge(l.pdgId());
- }
- assert(charge3_x_eta != 0);
- assert(charge3!=0);
- if (emom.Et() > 30/GeV) {
- if (charge3_x_eta < 0) {
- _htmp_dsigminus_deta->fill(emom.eta(), weight);
- } else {
- _htmp_dsigplus_deta->fill(emom.eta(), weight);
- }
- }
- if (charge3 < 0) {
- _h_Wminus_pT->fill(wmom.pT(), weight);
- } else {
- _h_Wplus_pT->fill(wmom.pT(), weight);
- }
-
- const FastJets& jetpro = applyProjection<FastJets>(e, "Jets");
- const Jets& jets = jetpro.jetsByPt(20.0*GeV);
- if (jets.size() > 0) {
- _h_W_jet1_deta->fill(wmom.eta()-jets[0].momentum().eta(), weight);
- _h_W_jet1_dR->fill(deltaR(wmom, jets[0].momentum()), weight);
- }
- MC_JetAnalysis::analyze(e);
+ MC_JetSplittings::analyze(e);
}
/// Finalize
void finalize() {
- scale(_h_W_mass, crossSection()/sumOfWeights());
- scale(_h_W_pT, crossSection()/sumOfWeights());
- scale(_h_W_pT_peak, crossSection()/sumOfWeights());
- scale(_h_W_y, crossSection()/sumOfWeights());
- scale(_h_W_phi, crossSection()/sumOfWeights());
- scale(_h_W_jet1_deta, crossSection()/sumOfWeights());
- scale(_h_W_jet1_dR, crossSection()/sumOfWeights());
- scale(_h_lepton_pT, crossSection()/sumOfWeights());
- scale(_h_lepton_eta, crossSection()/sumOfWeights());
-
- // Construct asymmetry: (dsig+/deta - dsig-/deta) / (dsig+/deta + dsig-/deta) for each Et region
- AIDA::IHistogramFactory& hf = histogramFactory();
- IHistogram1D* numtmp = hf.subtract("/numtmp", *_htmp_dsigplus_deta, *_htmp_dsigminus_deta);
- IHistogram1D* dentmp = hf.add("/dentmp", *_htmp_dsigplus_deta, *_htmp_dsigminus_deta);
- assert(numtmp && dentmp);
- hf.divide(histoDir() + "/W_chargeasymm_eta", *numtmp, *dentmp);
- hf.destroy(numtmp);
- hf.destroy(dentmp);
- hf.destroy(_htmp_dsigminus_deta);
- hf.destroy(_htmp_dsigplus_deta);
-
- // W charge asymmetry vs. pTW: dsig+/dpT / dsig-/dpT
- hf.divide(histoDir() + "/W_chargeasymm_pT", *_h_Wplus_pT, *_h_Wminus_pT);
- scale(_h_Wplus_pT, crossSection()/sumOfWeights());
- scale(_h_Wminus_pT, crossSection()/sumOfWeights());
-
- MC_JetAnalysis::finalize();
+ MC_JetSplittings::finalize();
}
//@}
-
- private:
-
- /// @name Histograms
- //@{
- AIDA::IHistogram1D * _h_W_mass;
- AIDA::IHistogram1D * _h_W_pT;
- AIDA::IHistogram1D * _h_W_pT_peak;
- AIDA::IHistogram1D * _h_W_y;
- AIDA::IHistogram1D * _h_W_phi;
- AIDA::IHistogram1D * _h_W_jet1_deta;
- AIDA::IHistogram1D * _h_W_jet1_dR;
- AIDA::IHistogram1D * _h_Wplus_pT;
- AIDA::IHistogram1D * _h_Wminus_pT;
- AIDA::IHistogram1D * _h_lepton_pT;
- AIDA::IHistogram1D * _h_lepton_eta;
-
- AIDA::IHistogram1D * _htmp_dsigminus_deta;
- AIDA::IHistogram1D * _htmp_dsigplus_deta;
- //@}
-
};
-
-
// The hook for the plugin system
- DECLARE_RIVET_PLUGIN(MC_WJETS);
+ DECLARE_RIVET_PLUGIN(MC_WKTSPLITTINGS);
}
Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/Makefile.am Wed Jan 23 10:37:49 2013 (r4105)
+++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Wed Jan 23 10:38:58 2013 (r4106)
@@ -360,7 +360,9 @@
MC_SUSY.cc \
MC_TTBAR.cc \
MC_VH2BB.cc \
+ MC_WINC.cc \
MC_WJETS.cc \
+ MC_WKTSPLITTINGS.cc \
MC_WPOL.cc \
MC_WWJETS.cc \
MC_XS.cc \
More information about the Rivet-svn
mailing list
|