|
|
[Rivet-svn] r4033 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo src/Analyses
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Sat Nov 17 22:09:29 GMT 2012
Author: buckley
Date: Sat Nov 17 22:09:29 2012
New Revision: 4033
Log:
Adding MC_PHOTONS by Steve Lloyd and AB, for testing general unisolated photon properties, especially those associated with charged leptons (e and mu).
Added:
branches/2012-06-aidarivet/data/anainfo/MC_PHOTONS.info
- copied, changed from r4032, branches/2012-06-aidarivet/data/anainfo/MC_PHOTONJETS.info
branches/2012-06-aidarivet/data/plotinfo/MC_PHOTONS.plot
- copied, changed from r4032, branches/2012-06-aidarivet/data/plotinfo/MC_PHOTONJETS.plot
branches/2012-06-aidarivet/src/Analyses/MC_PHOTONS.cc
Modified:
branches/2012-06-aidarivet/ChangeLog
branches/2012-06-aidarivet/data/anainfo/Makefile.am
branches/2012-06-aidarivet/data/plotinfo/Makefile.am
branches/2012-06-aidarivet/src/Analyses/Makefile.am
Modified: branches/2012-06-aidarivet/ChangeLog
==============================================================================
--- branches/2012-06-aidarivet/ChangeLog Fri Nov 16 21:51:44 2012 (r4032)
+++ branches/2012-06-aidarivet/ChangeLog Sat Nov 17 22:09:29 2012 (r4033)
@@ -1,3 +1,9 @@
+2012-11-17 Andy Buckley <andy.buckley at cern.ch>
+
+ * Adding MC_PHOTONS by Steve Lloyd and AB, for testing general
+ unisolated photon properties, especially those associated with
+ charged leptons (e and mu).
+
2012-11-16 Andy Buckley <andy.buckley at cern.ch>
* Adding MC_PRINTEVENT, a convenient (but verbose!) analysis for
Copied and modified: branches/2012-06-aidarivet/data/anainfo/MC_PHOTONS.info (from r4032, branches/2012-06-aidarivet/data/anainfo/MC_PHOTONJETS.info)
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/MC_PHOTONJETS.info Fri Nov 16 21:51:44 2012 (r4032, copy source)
+++ branches/2012-06-aidarivet/data/anainfo/MC_PHOTONS.info Sat Nov 17 22:09:29 2012 (r4033)
@@ -1,13 +1,18 @@
-Name: MC_PHOTONJETS
-Summary: Monte Carlo validation observables for photon + jets production
-Status: VALIDATED
+Name: MC_PHOTONS
+Summary: Monte Carlo validation observables for general photons
+Status: UNVALIDATED
Authors:
- - Frank Siegert <frank.siegert at cern.ch>
-References:
+ - Steve Lloyd
+ - Andy Buckley <andy.buckley at cern.ch>
RunInfo:
- Tevatron Run II ppbar -> gamma + jets.
-NumEvents: 10000000
+ Any event type, but there are many observables for photons associated to
+ (semi-)hard leptons.
+NumEvents: 10000
PtCuts: [0]
Description:
- Different observables related to the photon and extra jets.
-NeedCrossSection: yes
+ Observables for testing general unisolated photon properties, especially those
+ associated with charged leptons (e and mu).
+NeedCrossSection: no
+ToDo:
+ - Add a cone size cut on the photon--lepton association.
+ - Add FastJet-based clustering of photons on to leptons.
Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/Makefile.am Fri Nov 16 21:51:44 2012 (r4032)
+++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Sat Nov 17 22:09:29 2012 (r4033)
@@ -175,6 +175,7 @@
MC_PDFS.info \
MC_PHOTONJETS.info \
MC_PHOTONJETUE.info \
+ MC_PHOTONS.info \
MC_PRINTEVENT.info \
MC_SUSY.info \
MC_TTBAR.info \
Copied and modified: branches/2012-06-aidarivet/data/plotinfo/MC_PHOTONS.plot (from r4032, branches/2012-06-aidarivet/data/plotinfo/MC_PHOTONJETS.plot)
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/MC_PHOTONJETS.plot Fri Nov 16 21:51:44 2012 (r4032, copy source)
+++ branches/2012-06-aidarivet/data/plotinfo/MC_PHOTONS.plot Sat Nov 17 22:09:29 2012 (r4033)
@@ -1,374 +1,66 @@
-# BEGIN PLOT /MC_PHOTONJETS/photon_jet1_dR
-Title=
-XLabel=$\Delta R(\gamma\text{, 1st jet})$
-YLabel=$\text{d}\sigma/\text{d}\Delta R(\gamma\text{, 1st jet})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/photon_jet1_deta
-Title=
-XLabel=$\Delta{\eta}(\gamma\text{, 1st jet})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{\eta}(\gamma\text{, 1st jet})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/photon_jet1_dphi
-Title=
-XLabel=$\Delta{\phi}(\gamma\text{, 1st jet})$
-YLabel=$\text{d}\sigma/\text{d}\Delta{\phi}(\gamma\text{, 1st jet})$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/photon_pT$
+# BEGIN PLOT /MC_PHOTONS/Ptgamma
Title=Photon $p_\perp$
-XLabel=$p_\perp^{\gamma}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_\perp^{\gamma}$ [pb/GeV]
-LogX=1
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/photon_pT_lin$
-Title=Photon $p_\perp$
-XLabel=$p_\perp^{\gamma}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}y_{\gamma}$ [pb/GeV]
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/photon_y
-Title=Photon rapidity
-XLabel=$y_{\gamma}$
-YLabel=$\text{d}\sigma/\text{d}y_{\gamma}$ [pb]
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/jets_dR_
-Title=$\Delta R$ separation between jets
-LegendXPos=0.10
-LegendYPos=0.5
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/jets_deta_
-Title=Pseudorapidity separation between jets
-LegendYPos=0.5
-LegendXPos=0.30
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/jets_dphi_
-Title=Pseudorapidity separation between jets
-LegendXPos=0.1
-# END PLOT
-
-# BEGIN PLOT /MC_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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_PHOTONJETS/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
+XLabel=$p_\perp^\gamma$ [GeV]
+#YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}p_\perp^\gamma$ [GeV$^{-1}$]
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/Egamma
+Title=Photon energy
+XLabel=$E_\gamma$ [GeV]
+#YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}E_\gamma$ [GeV$^{-1}$]
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/sumPtgamma
+Title=Scalar sum of photon $p_\perp$
+XLabel=$\sum{p_\perp^\gamma}$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}\sum{p_\perp^\gamma}$ [GeV$^{-1}$]
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/sumEgamma
+Title=Sum of photon energies
+XLabel=$\sum{E_\gamma}$ [GeV]
+YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}\sum{E_\gamma}$ [GeV$^{-1}$]
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/deltaR
+Title=$\Delta{R}$ between photons and associated leptons
+XLabel=$\Delta{R}$
+#YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}\Delta{R}$
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/deltaR_weighted
+Title=$p_\perp$-weighted $\Delta{R}$ between photons and associated leptons
+XLabel=$\Delta{R}$
+#YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}\Delta{R}$
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/deltaR_R
+Title=Phase-space weighted $\Delta{R}$ between photons and associated leptons
+XLabel=$\Delta{R}$
+#YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}\Delta{R}$
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/deltaR_R_weighted
+Title=Phase-space \& $p_\perp$-weighted $\Delta{R}$ between photons and associated leptons
+XLabel=$\Delta{R}$
+#YLabel=$1/\sigma \, \mathrm{d}\sigma/\mathrm{d}\Delta{R}$
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/.*_vs_pTlep
+XLabel=Lepton $p_\perp$ [GeV]
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/DeltaR_vs_pTlep
+Title=Mean $\Delta{R}$ between photons and associated leptons vs. lepton $p_\perp$
+YLabel=$\angle\Delta{R}\rangle$
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/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
+# BEGIN PLOT /MC_PHOTONS/DeltaR_weighted_vs_pTlep
+Title=Mean $p_\perp$-weighted $\Delta{R}$ between photons and associated leptons vs. lepton $p_\perp$
+YLabel=$\angle\text{Weighted}\quad\Delta{R}\rangle$
# END PLOT
-# BEGIN PLOT /MC_PHOTONJETS/jet_HT
-Title=Scalar sum of jet transverse momenta (HT)
-XLabel=$HT$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}HT$ [pb/GeV]
-LogX=1
+# BEGIN PLOT /MC_PHOTONS/sumPtGamma_vs_pTlep
+Title=Mean scalar $\sum{p_\perp}$ of photons on associated leptons vs. lepton $p_\perp$
+YLabel=$\angle\sum{p_\perp}\rangle$
# END PLOT
Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Fri Nov 16 21:51:44 2012 (r4032)
+++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Sat Nov 17 22:09:29 2012 (r4033)
@@ -172,6 +172,7 @@
MC_TTBAR.plot \
MC_JETS.plot \
MC_PHOTONJETS.plot \
+ MC_PHOTONS.plot \
MC_VH2BB.plot \
MC_HJETS.plot \
MC_PDFS.plot \
Added: branches/2012-06-aidarivet/src/Analyses/MC_PHOTONS.cc
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2012-06-aidarivet/src/Analyses/MC_PHOTONS.cc Sat Nov 17 22:09:29 2012 (r4033)
@@ -0,0 +1,144 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/IdentifiedFinalState.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+
+namespace Rivet {
+
+
+ class MC_PHOTONS : public Analysis {
+ public:
+
+ /// @name Constructors etc.
+ //@{
+
+ /// Constructor
+ MC_PHOTONS()
+ : Analysis("MC_PHOTONS")
+ { }
+
+ //@}
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Book histograms and initialise projections before the run
+ void init() {
+ IdentifiedFinalState leptons(-5.0, 5.0, 10*GeV);
+ IdentifiedFinalState photons(-5.0, 5.0);
+ leptons.acceptChLeptons();
+ photons.acceptId(PHOTON);
+ addProjection(leptons, "lFS");
+ addProjection(photons, "gammaFS");
+
+ _h_Ptgamma = bookHistogram1D("Ptgamma", logspace(50, 0.01, 30));
+ _h_Egamma = bookHistogram1D("Egamma", logspace(50, 0.01, 200));
+ _h_sumPtgamma = bookHistogram1D("sumPtgamma", logspace(50, 0.1, 100));
+ _h_sumEgamma = bookHistogram1D("sumEgamma", logspace(50, 0.1, 500));
+ _h_DelR = bookHistogram1D("DeltaR", 50, 0, 2);
+ _h_DelR_weighted = bookHistogram1D("DeltaR_weighted", 50, 0, 2);
+ _h_DelR_R = bookHistogram1D("DeltaR_R", 50, 0, 2);
+ _h_DelR_R_weighted = bookHistogram1D("DeltaR_R_weighted", 50, 0, 2);
+ _p_DelR_vs_pTl = bookProfile1D("DeltaR_vs_pTlep", logspace(50, 10, 50*log10(sqrtS())));
+ _p_DelR_weighted_vs_pTl = bookProfile1D("DeltaR_weighted_vs_pTlep", logspace(50, 10, 50*log10(sqrtS())));
+ _p_sumPtgamma_vs_pTl = bookProfile1D("sumPtGamma_vs_pTlep", logspace(50, 10, 50*log10(sqrtS())));
+ }
+
+
+ /// Perform the per-event analysis
+ void analyze(const Event& event) {
+ const double weight = event.weight();
+
+ /// Get photons and leptons
+ const ParticleVector& photons = applyProjection<FinalState>(event, "gammaFS").particles();
+ MSG_DEBUG("Photon multiplicity = " << photons.size());
+ const ParticleVector& leptons = applyProjection<FinalState>(event, "lFS").particles();
+ MSG_DEBUG("Photon multiplicity = " << leptons.size());
+
+ // Initialise a map of sumPtgamma for each lepton
+ map<size_t, double> sumpT_per_lep;
+ for (size_t il = 0; il < leptons.size(); ++il) sumpT_per_lep[il] = 0;
+
+ // Calculate photon energies and transverse momenta
+ double sumPtgamma(0), sumEgamma(0);
+ foreach (const Particle& p, photons) {
+ // Individual and summed pTs and energies
+ double pTgamma = p.momentum().pT()/GeV;
+ double Egamma = p.momentum().E()/GeV;
+ _h_Ptgamma->fill(pTgamma, weight);
+ _h_Egamma->fill(Egamma, weight);
+ sumPtgamma += pTgamma;
+ sumEgamma += Egamma;
+
+ // Calculate delta R with respect to the nearest lepton
+ int ilep = -1;
+ double delR = 10000;
+ for (size_t il = 0; il < leptons.size(); ++il) {
+ const double tmpdelR = deltaR(leptons[il].momentum(), p.momentum());
+ if (tmpdelR < delR) {
+ ilep = il;
+ delR = tmpdelR;
+ }
+ }
+ if (ilep != -1) {
+ _h_DelR->fill(delR, weight);
+ _h_DelR_weighted->fill(delR*pTgamma/GeV, weight);
+ _h_DelR_R->fill(delR, weight/(delR+1e-5));
+ _h_DelR_R_weighted->fill(delR*pTgamma/GeV, weight/(delR+1e-5));
+ _p_DelR_vs_pTl->fill(leptons[ilep].momentum().pT()/GeV, delR, weight);
+ _p_DelR_weighted_vs_pTl->fill(leptons[ilep].momentum().pT()/GeV, delR*pTgamma, weight);
+ sumpT_per_lep[ilep] += pTgamma;
+ }
+ }
+
+ // Histogram whole-event photon HT/energy
+ _h_sumPtgamma->fill(sumPtgamma/GeV, weight);
+ _h_sumEgamma->fill(sumEgamma/GeV, weight);
+
+ // Histogram per-lepton sum(pT)
+ for (size_t il = 0; il < leptons.size(); ++il) {
+ _p_sumPtgamma_vs_pTl->fill(leptons[il].momentum().pT()/GeV, sumpT_per_lep[il]/GeV, weight);
+ }
+
+ }
+
+
+ /// Normalise histograms etc., after the run
+ void finalize() {
+ normalize(_h_Ptgamma);
+ normalize(_h_Egamma);
+ normalize(_h_sumPtgamma);
+ normalize(_h_sumEgamma);
+ normalize(_h_DelR);
+ normalize(_h_DelR_weighted);
+ normalize(_h_DelR_R);
+ normalize(_h_DelR_R_weighted);
+ }
+
+ //@}
+
+
+ private:
+
+ /// @name Histograms
+ //@{
+ AIDA::IHistogram1D *_h_Ptgamma, *_h_Egamma;
+ AIDA::IHistogram1D *_h_sumPtgamma, *_h_sumEgamma;
+ AIDA::IHistogram1D *_h_DelR, *_h_DelR_weighted;
+ AIDA::IHistogram1D *_h_DelR_R, *_h_DelR_R_weighted;
+ AIDA::IProfile1D *_p_DelR_vs_pTl, *_p_DelR_weighted_vs_pTl;
+ AIDA::IProfile1D *_p_sumPtgamma_vs_pTl;
+ //@}
+
+ };
+
+
+ // The hook for the plugin system
+ DECLARE_RIVET_PLUGIN(MC_PHOTONS);
+
+
+}
Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/Makefile.am Fri Nov 16 21:51:44 2012 (r4032)
+++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Sat Nov 17 22:09:29 2012 (r4033)
@@ -341,6 +341,7 @@
MC_LEADJETUE.cc \
MC_PDFS.cc \
MC_PHOTONJETS.cc \
+ MC_PHOTONS.cc \
MC_PRINTEVENT.cc \
MC_SUSY.cc \
MC_TTBAR.cc \
More information about the Rivet-svn
mailing list
|