|
[Rivet-svn] r4114 - branches/2012-06-aidarivet/data/anainfo branches/2012-06-aidarivet/data/plotinfo branches/2012-06-aidarivet/src/Analyses trunk/data/anainfo trunk/data/plotinfo trunk/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Jan 31 14:42:29 GMT 2013
Author: fsiegert Date: Thu Jan 31 14:42:29 2013 New Revision: 4114 Log: Add generic QCD validation analysis at parton level: MC_QCD_PARTONS Added: branches/2012-06-aidarivet/data/anainfo/MC_QCD_PARTONS.info branches/2012-06-aidarivet/data/plotinfo/MC_QCD_PARTONS.plot branches/2012-06-aidarivet/src/Analyses/MC_QCD_PARTONS.cc trunk/data/anainfo/MC_QCD_PARTONS.info trunk/data/plotinfo/MC_QCD_PARTONS.plot trunk/src/Analyses/MC_QCD_PARTONS.cc Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am branches/2012-06-aidarivet/data/plotinfo/Makefile.am branches/2012-06-aidarivet/src/Analyses/Makefile.am trunk/data/anainfo/Makefile.am trunk/data/plotinfo/Makefile.am trunk/src/Analyses/Makefile.am Added: branches/2012-06-aidarivet/data/anainfo/MC_QCD_PARTONS.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/anainfo/MC_QCD_PARTONS.info Thu Jan 31 14:42:29 2013 (r4114) @@ -0,0 +1,15 @@ +Name: MC_QCD_PARTONS +Summary: Generic parton-level Monte Carlo validation analysis for <X> + jets. +Status: VALIDATED +Authors: + - Frank Siegert <frank.siegert at cern.ch> +References: +RunInfo: + Any <X> + jets. +NumEvents: 1000000 +PtCuts: [0] +Description: + Only partons (excluding top quarks) are taken into account to construct a + kt cluster sequence. Thus this analysis can be used as a generic validation + tool for QCD activity. +NeedCrossSection: yes Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/anainfo/Makefile.am Wed Jan 30 14:15:47 2013 (r4113) +++ branches/2012-06-aidarivet/data/anainfo/Makefile.am Thu Jan 31 14:42:29 2013 (r4114) @@ -194,6 +194,7 @@ MC_PHOTONJETUE.info \ MC_PHOTONS.info \ MC_PRINTEVENT.info \ + MC_QCD_PARTONS.info \ MC_SUSY.info \ MC_TTBAR.info \ MC_VH2BB.info \ Added: branches/2012-06-aidarivet/data/plotinfo/MC_QCD_PARTONS.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/data/plotinfo/MC_QCD_PARTONS.plot Thu Jan 31 14:42:29 2013 (r4114) @@ -0,0 +1,57 @@ +# BEGIN PLOT /MC_QCD_PARTONS/.* +LegendAlign=r +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/data/plotinfo/Makefile.am Wed Jan 30 14:15:47 2013 (r4113) +++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am Thu Jan 31 14:42:29 2013 (r4114) @@ -188,6 +188,7 @@ MC_PHOTONJETS.plot \ MC_PHOTONKTSPLITTINGS.plot \ MC_PHOTONS.plot \ + MC_QCD_PARTONS.plot \ MC_VH2BB.plot \ MC_HINC.plot \ MC_HJETS.plot \ Added: branches/2012-06-aidarivet/src/Analyses/MC_QCD_PARTONS.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ branches/2012-06-aidarivet/src/Analyses/MC_QCD_PARTONS.cc Thu Jan 31 14:42:29 2013 (r4114) @@ -0,0 +1,67 @@ +// -*- C++ -*- +#include "Rivet/Analyses/MC_JetSplittings.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" +#include "Rivet/Projections/FastJets.hh" +#include "Rivet/RivetAIDA.hh" + +namespace Rivet { + + + /// Generic analysis looking at kt splitting scales of partons + class MC_QCD_PARTONS : public MC_JetSplittings { + public: + + /// Constructor + MC_QCD_PARTONS() + : MC_JetSplittings("MC_QCD_PARTONS", 4, "Jets") + { } + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Projections + IdentifiedFinalState partonfs; + for (int i=1; i<6; ++i) partonfs.acceptIdPair(i); + partonfs.acceptId(GLUON); + addProjection(FastJets(partonfs, FastJets::KT, 0.6), "Jets"); + + MC_JetSplittings::init(); + } + + + + /// Perform the per-event analysis + void analyze(const Event& event) { + MC_JetSplittings::analyze(event); + } + + + + /// Finalize + void finalize() { + MC_JetSplittings::finalize(); + } + + //@} + + + private: + + /// @name Histograms + //@{ + //@} + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(MC_QCD_PARTONS); + +} Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/Makefile.am Wed Jan 30 14:15:47 2013 (r4113) +++ branches/2012-06-aidarivet/src/Analyses/Makefile.am Thu Jan 31 14:42:29 2013 (r4114) @@ -347,6 +347,7 @@ lib_LTLIBRARIES += RivetMCAnalyses.la RivetMCAnalyses_la_SOURCES = \ EXAMPLE.cc \ + MC_QCD_PARTONS.cc \ MC_DIPHOTON.cc \ MC_GENERIC.cc \ MC_HINC.cc \ Added: trunk/data/anainfo/MC_QCD_PARTONS.info ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/anainfo/MC_QCD_PARTONS.info Thu Jan 31 14:42:29 2013 (r4114) @@ -0,0 +1,15 @@ +Name: MC_QCD_PARTONS +Summary: Generic parton-level Monte Carlo validation analysis for <X> + jets. +Status: VALIDATED +Authors: + - Frank Siegert <frank.siegert at cern.ch> +References: +RunInfo: + Any <X> + jets. +NumEvents: 1000000 +PtCuts: [0] +Description: + Only partons (excluding top quarks) are taken into account to construct a + kt cluster sequence. Thus this analysis can be used as a generic validation + tool for QCD activity. +NeedCrossSection: yes Modified: trunk/data/anainfo/Makefile.am ============================================================================== --- trunk/data/anainfo/Makefile.am Wed Jan 30 14:15:47 2013 (r4113) +++ trunk/data/anainfo/Makefile.am Thu Jan 31 14:42:29 2013 (r4114) @@ -183,6 +183,7 @@ MC_PHOTONJETUE.info \ MC_PHOTONS.info \ MC_PRINTEVENT.info \ + MC_QCD_PARTONS.info \ MC_SUSY.info \ MC_TTBAR.info \ MC_VH2BB.info \ Added: trunk/data/plotinfo/MC_QCD_PARTONS.plot ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/data/plotinfo/MC_QCD_PARTONS.plot Thu Jan 31 14:42:29 2013 (r4114) @@ -0,0 +1,57 @@ +# BEGIN PLOT /MC_QCD_PARTONS/.* +LegendAlign=r +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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}$ +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT + +# BEGIN PLOT /MC_QCD_PARTONS/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] +# END PLOT Modified: trunk/data/plotinfo/Makefile.am ============================================================================== --- trunk/data/plotinfo/Makefile.am Wed Jan 30 14:15:47 2013 (r4113) +++ trunk/data/plotinfo/Makefile.am Thu Jan 31 14:42:29 2013 (r4114) @@ -177,6 +177,7 @@ MC_PHOTONJETS.plot \ MC_PHOTONKTSPLITTINGS.plot \ MC_PHOTONS.plot \ + MC_QCD_PARTONS.plot \ MC_VH2BB.plot \ MC_HINC.plot \ MC_HJETS.plot \ Added: trunk/src/Analyses/MC_QCD_PARTONS.cc ============================================================================== --- /dev/null 00:00:00 1970 (empty, because file is newly added) +++ trunk/src/Analyses/MC_QCD_PARTONS.cc Thu Jan 31 14:42:29 2013 (r4114) @@ -0,0 +1,67 @@ +// -*- C++ -*- +#include "Rivet/Analyses/MC_JetSplittings.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/IdentifiedFinalState.hh" +#include "Rivet/Projections/FastJets.hh" +#include "Rivet/RivetYODA.hh" + +namespace Rivet { + + + /// Generic analysis looking at kt splitting scales of partons + class MC_QCD_PARTONS : public MC_JetSplittings { + public: + + /// Constructor + MC_QCD_PARTONS() + : MC_JetSplittings("MC_QCD_PARTONS", 4, "Jets") + { } + + + public: + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Projections + IdentifiedFinalState partonfs; + for (int i=1; i<6; ++i) partonfs.acceptIdPair(i); + partonfs.acceptId(GLUON); + addProjection(FastJets(partonfs, FastJets::KT, 0.6), "Jets"); + + MC_JetSplittings::init(); + } + + + + /// Perform the per-event analysis + void analyze(const Event& event) { + MC_JetSplittings::analyze(event); + } + + + + /// Finalize + void finalize() { + MC_JetSplittings::finalize(); + } + + //@} + + + private: + + /// @name Histograms + //@{ + //@} + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(MC_QCD_PARTONS); + +} Modified: trunk/src/Analyses/Makefile.am ============================================================================== --- trunk/src/Analyses/Makefile.am Wed Jan 30 14:15:47 2013 (r4113) +++ trunk/src/Analyses/Makefile.am Thu Jan 31 14:42:29 2013 (r4114) @@ -335,6 +335,7 @@ lib_LTLIBRARIES += RivetMCAnalyses.la RivetMCAnalyses_la_SOURCES = \ EXAMPLE.cc \ + MC_QCD_PARTONS.cc \ MC_DIPHOTON.cc \ MC_GENERIC.cc \ MC_HINC.cc \
More information about the Rivet-svn mailing list |