|
|
[Rivet-svn] r3829 - in schools/2009-Lund/setups: . herwig pythia sherpa sherpa/ttbar sherpa/wbbjj sherpa/wjets sherpa/zjets ttbar_rivetplugin
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Tue Jul 17 14:05:03 BST 2012
Author: hoeth
Date: Tue Jul 17 14:05:03 2012
New Revision: 3829
Log:
Lund setups
Added:
schools/2009-Lund/setups/
schools/2009-Lund/setups/herwig/
schools/2009-Lund/setups/herwig/Graph.in
schools/2009-Lund/setups/herwig/TopBackground.in
schools/2009-Lund/setups/herwig/TopSignal.in
schools/2009-Lund/setups/herwig/WJetsRivet.in
schools/2009-Lund/setups/pythia/
schools/2009-Lund/setups/pythia/topbackground.cmnd
schools/2009-Lund/setups/pythia/topsignal.cmnd
schools/2009-Lund/setups/sherpa/
schools/2009-Lund/setups/sherpa/ttbar/
schools/2009-Lund/setups/sherpa/ttbar/Analysis.dat
schools/2009-Lund/setups/sherpa/ttbar/Run.dat
schools/2009-Lund/setups/sherpa/wbbjj/
schools/2009-Lund/setups/sherpa/wbbjj/Analysis.dat
schools/2009-Lund/setups/sherpa/wbbjj/Run.dat
schools/2009-Lund/setups/sherpa/wjets/
schools/2009-Lund/setups/sherpa/wjets/Analysis.dat
schools/2009-Lund/setups/sherpa/wjets/Run.dat
schools/2009-Lund/setups/sherpa/zjets/
schools/2009-Lund/setups/sherpa/zjets/Analysis.dat
schools/2009-Lund/setups/sherpa/zjets/Run.dat
schools/2009-Lund/setups/ttbar_rivetplugin/
schools/2009-Lund/setups/ttbar_rivetplugin/Makefile
schools/2009-Lund/setups/ttbar_rivetplugin/TopAnalysis.cc
schools/2009-Lund/setups/ttbar_rivetplugin/TopAnalysis.hh
Added: schools/2009-Lund/setups/herwig/Graph.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/herwig/Graph.in Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,98 @@
+##################################################
+# Example generator based on LHC parameters
+# usage: Herwig++ read LHC.in
+##################################################
+
+##################################################
+# Technical parameters for this run
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 1
+set LHCGenerator:RandomNumberGenerator:Seed 31101
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+
+##################################################
+# LHC physics parameters (override defaults here)
+##################################################
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 14000.0
+
+# Intrinsic pT tune extrapolated to LHC energy
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 2.2*GeV
+
+##################################################
+# Matrix Elements for hadron-hadron collisions
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements/
+
+# Drell-Yan Z/gamma
+insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff
+# Drell-Yan W
+# insert SimpleQCD:MatrixElements[0] MEqq2W2ff
+# gamma-gamma
+# insert SimpleQCD:MatrixElements[0] MEGammaGamma
+# gamma+jet
+# insert SimpleQCD:MatrixElements[0] MEGammaJet
+# gg/qqbar -> Higgs (recommend including q qbar->Hg as not in ME correction)
+# insert SimpleQCD:MatrixElements[0] MEHiggs
+# insert SimpleQCD:MatrixElements[0] MEHiggsJet
+# set MEHiggsJet:Process qqbar
+# set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+# higgs+jet
+# insert SimpleQCD:MatrixElements[0] MEHiggsJet
+# QCD 2-2 scattering
+# insert SimpleQCD:MatrixElements[0] MEQCD2to2
+# top-antitop production
+# insert SimpleQCD:MatrixElements[0] MEHeavyQuark
+# W+jet
+#insert SimpleQCD:MatrixElements[0] MEWJet
+# Z+jet
+#insert SimpleQCD:MatrixElements[0] MEZJet
+# higgs + W (N.B. if considering all W decay modes useful to set )
+# (jet pT cut to zero so no cut on W decay products )
+# insert SimpleQCD:MatrixElements[0] MEPP2WH
+# set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+# higgs + Z (N.B. if considering all Z decay modes useful to set )
+# (jet pT cut to zero so no cut on Z decay products )
+# insert SimpleQCD:MatrixElements[0] MEPP2ZH
+# set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+#
+cd /Herwig/Generators
+
+##################################################
+# Useful analysis handlers for hadron-hadron physics
+##################################################
+# analysis of W/Z events
+# insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/DrellYan
+# analysis of top-antitop events
+# insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/TTbar
+# analysis of gamma+jet events
+# insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/GammaJet
+# analysis of gamma-gamma events
+# insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/GammaGamma
+# analysis of higgs-jet events
+# insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HiggsJet
+
+##################################################
+# Useful analysis handlers for HepMC related output
+##################################################
+# Schematic overview of an event (requires --with-hepmc to be set at configure time
+# and the graphviz program 'dot' to produce a plot)
+insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/Plot
+# A HepMC dump file (requires --with-hepmc to be set at configure time)
+# insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+# set /Herwig/Analysis/HepMCFile:PrintEvent 100
+# set /Herwig/Analysis/HepMCFile:Format GenEvent
+# set /Herwig/Analysis/HepMCFile:Units GeV_mm
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+
+set /Herwig/EventHandlers/LHCHandler:HadronizationHandler NULL
+set /Herwig/EventHandlers/LHCHandler:DecayHandler NULL
+set /Herwig/Analysis/Basics:CheckQuark No
+
+run Graph LHCGenerator
Added: schools/2009-Lund/setups/herwig/TopBackground.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/herwig/TopBackground.in Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,89 @@
+##################################################
+# Example generator based on LHC parameters
+# usage: Herwig++ read LHC.in
+##################################################
+cd /Herwig/EventHandlers
+# create the object ot read to madgraph sample
+create ThePEG::MadGraphReader TopReader LesHouches.so MadGraphReader.so
+set TopReader:PDFA /Herwig/Partons/MRST
+set TopReader:PDFB /Herwig/Partons/MRST
+set TopReader:EBeamA 7000.*GeV
+set TopReader:EBeamB 7000.*GeV
+set TopReader:BeamA /Herwig/Particles/p+
+set TopReader:BeamB /Herwig/Particles/p+
+set TopReader:Cuts /Herwig/Cuts/QCDCuts
+set TopReader:PartonExtractor /Herwig/Partons/QCDExtractor
+
+# Different prepared samples
+#
+# Alternative 0: W b bbar jet jet with both W+ and W-.
+# Here jet = d, dbar, u, ubar, s, sbar, c, cbar, b, bbar, g.
+# /usr/local/mcnet/lhef/wbbjj10k.lhe
+
+# In subsequent files jet is not b or bbar.
+
+# Alternative 1: W b bbar jet jet background with W+ -> e+.
+#/usr/local/mcnet/lhef/wposbbjj50k.lhe
+
+# Alternative 2: W b bbar jet jet background with W- -> e-.
+#/usr/local/mcnet/lhef/wnegbbjj50k.lhe
+
+# Alternative 3: W b bbar jet background with W+ -> e+.
+#/usr/local/mcnet/lhef/wposbbj50k.lhe
+
+# Alternative 4: W b bbar jet background with W- -> e-.
+#/usr/local/mcnet/lhef/wnegbbj50k.lhe
+
+# Alternative 5: W b bbar background with W+ -> e+.
+#/usr/local/mcnet/lhef/wposbb36k.lhe
+
+# Alternative 6: W b bbar background with W- -> e-.
+#/usr/local/mcnet/lhef/wnegbb28k.lhe
+
+# Use Alt. 1
+set TopReader:FileName /usr/local/mcnet/lhef/wposbbjj50k.lhe
+
+
+##################################################
+# Connect the LesHouches reader to the generator
+##################################################
+# create the event handler
+create ThePEG::LesHouchesEventHandler TopEventHandler LesHouches.so
+insert TopEventHandler:LesHouchesReaders[0] TopReader
+set TopEventHandler:CascadeHandler /Herwig/Shower/ShowerHandler
+set TopEventHandler:HadronizationHandler /Herwig/Hadronization/ClusterHadHandler
+set TopEventHandler:DecayHandler /Herwig/Decays/DecayHandler
+cd /Herwig/Generators
+set LHCGenerator:EventHandler /Herwig/EventHandlers/TopEventHandler
+
+##################################################
+# Technical parameters for this run
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+
+# Intrinsic pT tune extrapolated to LHC energy
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 2.2*GeV
+
+##################################################
+# Enable HepMC format output
+##################################################
+cd /Herwig/Generators
+insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+set /Herwig/Analysis/HepMCFile:PrintEvent 5000
+set /Herwig/Analysis/HepMCFile:Format GenEvent
+set /Herwig/Analysis/HepMCFile:Units GeV_mm
+set /Herwig/Analysis/HepMCFile:Filename hw-top-bkgrd.hepmc
+
+# For the long run, use the pipe! Uncomment the lines below.
+# set /Herwig/Analysis/HepMCFile:Filename hepmc.fifo
+# set /Herwig/Analysis/HepMCFile:PrintEvent 1000000
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+saverun TopBackground LHCGenerator
Added: schools/2009-Lund/setups/herwig/TopSignal.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/herwig/TopSignal.in Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,78 @@
+##################################################
+# Example generator based on LHC parameters
+# usage: Herwig++ read LHC.in
+##################################################
+
+##################################################
+# Technical parameters for this run
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+
+##################################################
+# Changes to top decays to get semi-leptonic top
+# and hadronic tbar
+##################################################
+# change so t and tbar properties not automatically
+# charge conjugates of each other
+set /Herwig/Particles/t:Synchronized Not_synchronized
+set /Herwig/Particles/tbar:Synchronized Not_synchronized
+# top decay modes, only semi-leptonic on
+set /Herwig/Particles/t/t->b,bbar,c;:OnOff Off
+set /Herwig/Particles/t/t->b,c,dbar;:OnOff Off
+set /Herwig/Particles/t/t->b,c,sbar;:OnOff Off
+set /Herwig/Particles/t/t->b,sbar,u;:OnOff Off
+set /Herwig/Particles/t/t->b,u,dbar;:OnOff Off
+set /Herwig/Particles/t/t->nu_e,e+,b;:OnOff On
+set /Herwig/Particles/t/t->nu_mu,mu+,b;:OnOff On
+set /Herwig/Particles/t/t->nu_tau,tau+,b;:OnOff Off
+# tbar decay modes, only hadronic on
+set /Herwig/Particles/tbar/tbar->b,bbar,cbar;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,cbar,d;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,cbar,s;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,s,ubar;:OnOff On
+set /Herwig/Particles/tbar/tbar->bbar,ubar,d;:OnOff On
+set /Herwig/Particles/tbar/tbar->nu_ebar,e-,bbar;:OnOff Off
+set /Herwig/Particles/tbar/tbar->nu_mubar,mu-,bbar;:OnOff Off
+set /Herwig/Particles/tbar/tbar->nu_taubar,tau-,bbar;:OnOff Off
+
+##################################################
+# LHC physics parameters (override defaults here)
+##################################################
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 14000.0
+
+# Intrinsic pT tune extrapolated to LHC energy
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 2.2*GeV
+
+##################################################
+# Matrix Elements for hadron-hadron collisions
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements/
+
+# top-antitop production
+insert SimpleQCD:MatrixElements[0] MEHeavyQuark
+
+##################################################
+# Enable HepMC format output
+##################################################
+cd /Herwig/Generators
+insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+set /Herwig/Analysis/HepMCFile:PrintEvent 5000
+set /Herwig/Analysis/HepMCFile:Format GenEvent
+set /Herwig/Analysis/HepMCFile:Units GeV_mm
+set /Herwig/Analysis/HepMCFile:Filename hw-top-signal.hepmc
+
+# For the long run, use the pipe! Uncomment the lines below.
+# set /Herwig/Analysis/HepMCFile:Filename hepmc.fifo
+# set /Herwig/Analysis/HepMCFile:PrintEvent 1000000
+
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+saverun TopSignal LHCGenerator
Added: schools/2009-Lund/setups/herwig/WJetsRivet.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/herwig/WJetsRivet.in Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,91 @@
+##################################################
+# Example generator based on Tevatron parameters
+# using NLO matrix elements and matching in
+# the Powheg formalism
+# usage: Herwig++ read TVT.in
+#
+# Since most parameters are identical to LHC,
+# we use the default LHCGenerator and adapt only
+# for the differences
+##################################################
+
+##################################################
+# Technical parameters for this run
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+set LHCGenerator:EventHandler:Sampler:Ntry 10000
+
+##################################################
+# Need to use an NLO PDF
+##################################################
+set /Herwig/Particles/p+:PDF /Herwig/Partons/MRST-NLO
+set /Herwig/Particles/pbar-:PDF /Herwig/Partons/MRST-NLO
+create Herwig::O2AlphaS /Herwig/O2AlphaS
+set /Herwig/Model:QCD/AlphaS /Herwig/O2AlphaS
+
+##################################################
+# Setup the POWHEG shower
+##################################################
+cd /Herwig/Shower
+# use the general recon for now
+set KinematicsReconstructor:ReconstructionOption General
+# create the Powheg evolver and use it instead of the default one
+create Herwig::PowhegEvolver PowhegEvolver HwPowhegShower.so
+set ShowerHandler:Evolver PowhegEvolver
+set PowhegEvolver:ShowerModel ShowerModel
+set PowhegEvolver:SplittingGenerator SplittingGenerator
+set PowhegEvolver:MECorrMode 0
+# create and use the Drell-yan hard emission generator
+create Herwig::DrellYanHardGenerator DrellYanHardGenerator
+set DrellYanHardGenerator:ShowerAlpha AlphaQCD
+insert PowhegEvolver:HardGenerator 0 DrellYanHardGenerator
+# create and use the gg->H hard emission generator
+create Herwig::GGtoHHardGenerator GGtoHHardGenerator
+set GGtoHHardGenerator:ShowerAlpha AlphaQCD
+insert PowhegEvolver:HardGenerator 0 GGtoHHardGenerator
+
+##################################################
+# Tevatron physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 1960.0
+set LHCGenerator:EventHandler:BeamB /Herwig/Particles/pbar-
+
+##################################################
+# Matrix Elements for hadron-hadron collisions
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+
+# Drell-Yan W
+insert SimpleQCD:MatrixElements[0] PowhegMEqq2W2ff
+# electrons only
+set PowhegMEqq2W2ff:Process 3
+
+##################################################
+# Enable HepMC format output
+##################################################
+cd /Herwig/Generators
+insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+set /Herwig/Analysis/HepMCFile:PrintEvent 1000000
+set /Herwig/Analysis/HepMCFile:Format GenEvent
+set /Herwig/Analysis/HepMCFile:Units GeV_mm
+set /Herwig/Analysis/HepMCFile:Filename hepmc.fifo
+
+##################################################
+# optionally, we can also call Rivet from inside Herwig
+##################################################
+#create Herwig::RivetAnalysis /Herwig/Analysis/RivetAnalysis RivetAnalysis.so
+#insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CDF_2008_S7541902
+#insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/RivetAnalysis
+#set LHCGenerator:EventHandler:StatLevel Full
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+saverun WJetsRivet LHCGenerator
Added: schools/2009-Lund/setups/pythia/topbackground.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/pythia/topbackground.cmnd Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,23 @@
+! File: topbackground.cmnd
+! The intent is to study the top mass in top pair production,
+! with one W decaying leptonically and the other hadronically.
+! This file generates the background of W + 4 jets from a MadGraph LHEF,
+! alternatively fewer jets, with rest to come from showers etc.
+! At parton-level the lepton and all jets have pT > 20 and Delta-R > 0.7.
+
+! 1) Settings that will be used in a main program.
+Main:numberOfEvents = 4000 ! number of events to generate
+Main:timesToShow = 0 ! show how far along run is this many times
+Main:timesAllowErrors = 10 ! abort run after this many flawed events
+Main:showChangedSettings = on ! print changed flags/modes/parameters
+Main:showAllSettings = off ! print all flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+Main:showAllParticleData = off ! print all particle and decay data
+
+! 2) Beam and process selection from a Les Houches Event File.
+Beams:frameType = 4 ! read info from a LHEF
+
+! W b bbar jet jet with W+
+! Here jet = d, dbar, u, ubar, s, sbar, c, cbar, g.
+Beams:LHEF = /home/rivet/school/day3/shared/Wposbbjj.lhe
+
Added: schools/2009-Lund/setups/pythia/topsignal.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/pythia/topsignal.cmnd Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,35 @@
+! File: topsignal.cmnd
+! The intent is to study the top mass in top pair production,
+! with one W decaying leptonically and the other hadronically.
+
+! 1) Settings that will be used in a main program.
+Main:numberOfEvents = 4000 ! number of events to generate
+Main:timesToShow = 0 ! show how far along run is this many times
+Main:timesAllowErrors = 10 ! abort run after this many flawed events
+Main:showChangedSettings = on ! print changed flags/modes/parameters
+Main:showAllSettings = off ! print all flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+Main:showAllParticleData = off ! print all particle and decay data
+
+! 2) Beam parameter settings. Values below agree with default ones.
+Beams:idA = 2212 ! first beam, p = 2212, pbar = -2212
+Beams:idB = 2212 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 7000. ! CM energy of collision
+
+! 3) Pick processes.
+Top:gg2ttbar = on ! g g -> t tbar
+Top:qqbar2ttbar = on ! q qbar -> t tbar
+
+! 4) Pick sign of lepton by forcing W decay modes.
+24:onMode = off ! all W decays are switched off
+
+! 4.1) Alternative 1: W+ -> e+, W- -> jets.
+24:onPosIfAny = 11 ! W+ can decay to e+
+24:onNegIfAny = 1 2 3 4 5 ! W- can decay to quark (d u s c b)
+
+! 4.2) Alternative 2: W+ -> jets, W- -> e-.
+#24:onPosIfAny = 1 2 3 4 5 ! W+ can decay to quark (d u s c b)
+#24:onNegIfAny = 11 ! W- can decay to e-
+
+! 5) Other settings.
+#6:m0 = 174.3 ! change the top mass
Added: schools/2009-Lund/setups/sherpa/ttbar/Analysis.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/ttbar/Analysis.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,101 @@
+BEGIN_ANALYSIS {
+
+ LEVEL ME Shower
+
+ // jets for control sample
+ Trigger {
+ InList FinalState;
+ OutList KTJetsWT;
+ Finder 93 20 -4.5 4.5 1.0 -1;
+ Qual !KF(24);
+ }
+
+ // jets for analysis
+ Trigger {
+ InList FinalState;
+ OutList KTJetsW;
+ Finder 93 20 -4.5 4.5 1.0 0;
+ Qual !KF(24);
+ }
+
+ // global observables
+ PT4 -24 24 97 97 0 1000 200 Lin KTJetsW
+ PT4 -24 24 97 97 0.1 1000 200 Log KTJetsW
+
+ // boson observables
+ DPhi -24 24 0 3.2 64 Lin FinalState
+ DEta -24 24 0 5 100 Lin FinalState
+ DR -24 24 0 6 100 Lin FinalState
+
+ PT2 -24 24 0 1000 250 Lin FinalState
+ Mass -24 24 0 1000 200 Lin FinalState
+
+ Eta -24 -5 5 200 Lin FinalState
+ PT -24 0 1000 200 Lin FinalState
+
+ // jet rates
+ JetDRate 1.0e-1 1.0e3 400 1 1 5 Log KTJetsW
+
+ // jet observables
+ OneEta 93 0 -5 5 200 Lin KTJetsW
+ OneEta 93 1 -5 5 200 Lin KTJetsW
+ OnePT 93 0 0 250 250 Lin KTJetsW
+ OnePT 93 1 0 250 250 Lin KTJetsW
+ OnePhi 93 0 0 3.2 128 Lin KTJetsW
+ OnePhi 93 1 0 3.2 128 Lin KTJetsW
+ TwoDPhi 93 0 93 1 0 3.2 200 Lin KTJetsW
+ TwoDR 93 0 93 1 0 6 150 Lin KTJetsW
+
+ // jet observables control sample
+ OneEta 93 0 -5 5 200 Lin KTJetsWT
+ OneEta 93 1 -5 5 200 Lin KTJetsWT
+ OnePT 93 0 0 250 250 Lin KTJetsWT
+ OnePT 93 1 0 250 250 Lin KTJetsWT
+ OnePhi 93 0 0 3.2 128 Lin KTJetsWT
+ OnePhi 93 1 0 3.2 128 Lin KTJetsWT
+ TwoDPhi 93 0 93 1 0 3.2 200 Lin KTJetsWT
+ TwoDR 93 0 93 1 0 6 150 Lin KTJetsWT
+
+ // bjet observables
+ OneEta 97 0 -5 5 200 Lin KTJetsW
+ OnePT 97 0 0 500 250 Lin KTJetsW
+ TwoDPhi 97 0 97 1 0 3.2 64 Lin KTJetsW
+ TwoDEta 97 0 97 1 0 5 100 Lin KTJetsW
+ TwoDR 97 0 97 1 0 6 100 Lin KTJetsW
+ TwoPT 97 0 97 1 0 1000 250 Lin KTJetsW
+ TwoMass 97 0 97 1 0 1000 200 Lin KTJetsW
+
+ // bjet observables control sample
+ OneEta 97 0 -5 5 200 Lin KTJetsWT
+ OnePT 97 0 0 500 250 Lin KTJetsWT
+ TwoDPhi 97 0 -97 0 0 3.2 64 Lin KTJetsWT
+ TwoDEta 97 0 -97 0 0 5 100 Lin KTJetsWT
+ TwoDR 97 0 -97 0 0 6 100 Lin KTJetsWT
+ TwoPT 97 0 -97 0 0 1000 250 Lin KTJetsWT
+ TwoMass 97 0 -97 0 0 1000 200 Lin KTJetsWT
+
+ // boson-jet observables
+ TwoDPhi 24 0 97 0 0 3.2 64 Lin KTJetsW
+ TwoDEta 24 0 97 0 0 5 100 Lin KTJetsW
+ TwoDR 24 0 97 0 0 6 100 Lin KTJetsW
+ TwoPT 24 0 97 0 0 1000 250 Lin KTJetsW
+ TwoMass 24 0 97 0 0 600 150 Lin KTJetsW
+
+ // boson-jet observables control sample
+ TwoDPhi 24 0 97 0 0 3.2 64 Lin KTJetsWT
+ TwoDEta 24 0 97 0 0 5 100 Lin KTJetsWT
+ TwoDR 24 0 97 0 0 6 100 Lin KTJetsWT
+ TwoPT 24 0 97 0 0 1000 250 Lin KTJetsWT
+ TwoMass 24 0 97 0 0 600 150 Lin KTJetsWT
+
+ // angles
+ FourPlaneAngle 24 0 -24 0 97 0 97 1 0 3.2 128 Lin KTJetsW
+ FourPlaneAngle 24 0 97 0 -24 0 97 1 0 3.2 128 Lin KTJetsW
+ FourPlaneAngle 24 0 97 1 -24 0 97 0 0 3.2 128 Lin KTJetsW
+
+ // angles control sample
+ FourPlaneAngle 24 0 -24 0 97 0 -97 0 0 3.2 128 Lin KTJetsWT
+ FourPlaneAngle 24 0 97 0 -24 0 -97 0 0 3.2 128 Lin KTJetsWT
+ FourPlaneAngle 24 0 -97 0 -24 0 97 0 0 3.2 128 Lin KTJetsWT
+
+} END_ANALYSIS
Added: schools/2009-Lund/setups/sherpa/ttbar/Run.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/ttbar/Run.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,117 @@
+(run){
+ EVENTS = 1000 ! Number of Events
+ OUTPUT = 2 ! Output level:
+ ! 0 = errors, 1 = events, 2 = run info,
+ ! 3 = 1+2 = events + run info
+ ANALYSIS = 0 ! Built-in analysis (see Analysis.dat): 0=off, 1=on
+ FILE_SIZE = 10000000 ! Events per file (if written out)
+ NUM_ACCURACY = 1.e-5
+}(run)
+
+(beam){
+ BEAM_1 = 2212 ! beam particle: 2212 = P+, -2212 = P-, 11 = e-
+ BEAM_ENERGY_1 = 7000. ! in GeV
+ K_PERP_MEAN_1 = 0.33 ! Settings for intrinsic kperp model
+ K_PERP_SIGMA_1 = 0.8 !
+
+ BEAM_2 = 2212
+ BEAM_ENERGY_2 = 7000.
+ K_PERP_MEAN_2 = 0.33
+ K_PERP_SIGMA_2 = 0.8
+}(beam)
+
+(processes){
+ %-------------------------------------------------------
+ %-- Processes to calculate -----------------------------
+ %-------------------------------------------------------
+
+ Process : 93 93 -> 6[a] -6[b] 93{1}
+ Order electroweak : 0
+ Decay : 6[a] -> -11 12 5
+ Decay : -6[b] -> 93 93 -5
+ Integration_Error 0.0017 {6}
+ Integration_Error 0.013 {7}
+ End process
+}(processes)
+
+(selector){
+ !-------------------------------------------------------
+ !-- Parton level cuts ----------------------------------
+ !-------------------------------------------------------
+ JetFinder sqr(30/E_CMS)[sqr(30/E_CMS)][sqr(30/E_CMS)] 1.
+ "m" -11,12,5 150.0,200.0
+ "m" 93,93,-5 150.0,200.0
+}(selector)
+
+(model){
+ !-------------------------------------------------------
+ !-- Model parameters -----------------------------------
+ !-------------------------------------------------------
+ MODEL = SM ! Model
+ STABLE[15] = 0 ! Set tau lepton (PDG-Id 15) unstable
+ MASSIVE[5] = 1
+}(model)
+
+(isr){
+ PDF_SET = cteq6l.LHpdf ! Choose PDF set to use with LHAPDF
+}(isr)
+
+(me){
+ !-------------------------------------------------------
+ !-- ME generators --------------------------------------
+ !-------------------------------------------------------
+ EVENT_GENERATION_MODE = Unweighted
+ SCALE_SCHEME = CKKW ! see 'http://www.sherpa-mc.de/scales.html'
+}(me)
+
+(integration){
+ FINISH_OPTIMIZATION = Off ! don't wait for complete phase space generator optimization
+}(integration)
+
+(shower){
+ !-------------------------------------------------------
+ !-- Parton showers -------------------------------------
+ !-------------------------------------------------------
+ SHOWER_GENERATOR = Apacic ! Apacic
+ !
+ FSR_SHOWER = 1 ! 1=On,0=Off
+ ISR_SHOWER = 1 ! 1=On,0=Off
+
+ SHOWER_MODE = 3 ! include shower from decays
+}(shower)
+
+(fragmentation){
+ !-------------------------------------------------------
+ !-- Fragmentation parameters----------------------------
+ !-------------------------------------------------------
+ FRAGMENTATION = Ahadic
+ DECAYMODEL = Hadrons
+ YFS_MODE = 2 ! Soft QED on hadron decays
+}(fragmentation)
+
+(mi){
+ !-------------------------------------------------------
+ ! Underlying Event Setup
+ !-------------------------------------------------------
+ !
+ ! general parameters
+ !
+ MI_HANDLER = None ! Amisic / None
+ !
+ ! hard underlying event parameters
+ !
+ CREATE_GRID 93 93 -> 93 93 ! processes to generate
+ PS_ERROR = 1.0e-2 ! error for integration
+ MI_SCALE_SCHEME = G_MEAN_PT2 ! Scale scheme
+ MI_K_FACTOR_SCHEME = 1 ! K-factor scheme
+ REGULATE_XS = 0 ! regulate cross section
+ XS_REGULATOR = QCD_Trivial ! regulation scheme
+ XS_REGULATION = 2.225 ! regulation parameter
+ SCALE_MIN = 2.225 ! minimum scale
+ JET_VETO = 1 ! apply jet veto
+ RESCALE_EXPONENT = 0.25 ! rescaling exponent
+ REFERENCE_SCALE = 1800.0 ! reference energy scale
+ PROFILE_FUNCTION = Double_Gaussian ! hadron profile function
+ PROFILE_PARAMETERS = 1.0 0.5 0.5 ! size (must be 1), coresize,
+ ! matter fraction
+}(mi)
Added: schools/2009-Lund/setups/sherpa/wbbjj/Analysis.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/wbbjj/Analysis.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,42 @@
+BEGIN_ANALYSIS {
+
+ LEVEL Hadron;
+
+ // Total cross section
+ Statistics FinalState
+
+ // Kt-Jets
+ Trigger {
+ InList FinalState
+ OutList kt_lepton_jets
+ JetMode 0
+ Finder 93 10. -5. 5. 0.4 1
+ }
+ Trigger {
+ InList kt_lepton_jets
+ OutList kt_jets
+ Keep 93
+ }
+
+ // W- variables
+ Mass 11 -12 50. 110. 20 LinErr FinalState
+ PT2 11 -12 0. 200. 50 LinErr FinalState
+ PT2 11 -12 0. 20. 20 LinErr FinalState
+ Eta2 11 -12 -6. 6. 48 LinErr FinalState
+
+ // W+ variables
+ Mass -11 12 50. 110. 20 LinErr FinalState
+ PT2 -11 12 0. 200. 50 LinErr FinalState
+ PT2 -11 12 0. 20. 20 LinErr FinalState
+ Eta2 -11 12 -6. 6. 48 LinErr FinalState
+
+ // jet variables
+ Multi -0.5 19.5 20 LinErr kt_jets
+ JetPT 0. 200. 50 1 1 4 LinErr kt_jets
+ JetEta -6. 6. 48 1 1 4 LinErr kt_jets
+ JetPhi -3.2 3.2 16 1 1 4 LinErr kt_jets
+ JetDR 0. 5. 25 1 2 4 LinErr kt_jets
+ JetDEta 0. 5. 25 1 2 4 LinErr kt_jets
+ JetDPhi 0. 3.2 16 1 2 4 LinErr kt_jets
+
+} END_ANALYSIS;
Added: schools/2009-Lund/setups/sherpa/wbbjj/Run.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/wbbjj/Run.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,124 @@
+(run){
+ EVENTS = 1000 ! Number of Events
+ OUTPUT = 2 ! Output level:
+ ! 0 = errors, 1 = events, 2 = run info,
+ ! 3 = 1+2 = events + run info
+ ANALYSIS = 0 ! Built-in analysis (see Analysis.dat): 0=off, 1=on
+ FILE_SIZE = 10000000 ! Events per file (if written out)
+ SHERPA_CPP_PATH = /usr/local/mcnet/sherpa113/setups/wbbjj
+ RESULT_DIRECTORY = /usr/local/mcnet/sherpa113/setups/wbbjj/Results
+}(run)
+
+(beam){
+ BEAM_1 = 2212 ! beam particle: 2212 = P+, -2212 = P-, 11 = e-
+ BEAM_ENERGY_1 = 7000. ! in GeV
+ K_PERP_MEAN_1 = 0.33 ! Settings for intrinsic kperp model
+ K_PERP_SIGMA_1 = 0.8 !
+
+ BEAM_2 = 2212
+ BEAM_ENERGY_2 = 7000.
+ K_PERP_MEAN_2 = 0.33
+ K_PERP_SIGMA_2 = 0.8
+}(beam)
+
+(processes){
+ %-------------------------------------------------------
+ %-- Processes to calculate -----------------------------
+ %-------------------------------------------------------
+
+ %
+ % jet jet -> electron neutrino + b bar + 2 jets
+ %
+ Process : 93 93 -> 11 -12 5 -5 93 93
+ Order electroweak : 2
+ Max_Epsilon : 1e-4
+ End process
+
+ %
+ % jet jet -> positron neutrino + b bbar + 2 jets
+ %
+ Process : 93 93 -> -11 12 5 -5 93 93
+ Order electroweak : 2
+ Max_Epsilon : 1e-4
+ End process
+}(processes)
+
+(selector){
+ !-------------------------------------------------------
+ !-- Parton level cuts ----------------------------------
+ !-------------------------------------------------------
+ JetFinder sqr(30/E_CMS) 1
+}(selector)
+
+(model){
+ !-------------------------------------------------------
+ !-- Model parameters -----------------------------------
+ !-------------------------------------------------------
+ MODEL = SM ! Model
+ MASSIVE[5] = 1
+ STABLE[15] = 0 ! Set tau lepton (PDG-Id 15) unstable
+}(model)
+
+(isr){
+ PDF_SET = cteq6l.LHpdf ! Choose PDF set to use with LHAPDF
+}(isr)
+
+(me){
+ !-------------------------------------------------------
+ !-- ME generators --------------------------------------
+ !-------------------------------------------------------
+ EVENT_GENERATION_MODE = Unweighted
+ SUDAKOV_WEIGHT = 0
+ SCALE_SCHEME = FIX_SCALE ! see 'http://www.sherpa-mc.de/scales.html'
+ FIXED_SCALE = 6467.215561
+}(me)
+
+(integration){
+ FINISH_OPTIMIZATION = On ! don't wait for complete phase space generator optimization
+}(integration)
+
+(shower){
+ !-------------------------------------------------------
+ !-- Parton showers -------------------------------------
+ !-------------------------------------------------------
+ SHOWER_GENERATOR = Apacic ! Apacic
+ !
+ FSR_SHOWER = 1 ! 1=On,0=Off
+ ISR_SHOWER = 1 ! 1=On,0=Off
+}(shower)
+
+(fragmentation){
+ !-------------------------------------------------------
+ !-- Fragmentation parameters----------------------------
+ !-------------------------------------------------------
+ FRAGMENTATION = Ahadic
+ DECAYMODEL = Hadrons
+ YFS_MODE = 2 ! Soft QED on hadron decays
+}(fragmentation)
+
+(mi){
+ !-------------------------------------------------------
+ ! Underlying Event Setup
+ !-------------------------------------------------------
+ !
+ ! general parameters
+ !
+ MI_HANDLER = None ! Amisic / None
+ !
+ ! hard underlying event parameters
+ !
+ CREATE_GRID 93 93 -> 93 93 ! processes to generate
+ PS_ERROR = 1.0e-2 ! error for integration
+ MI_SCALE_SCHEME = G_MEAN_PT2 ! Scale scheme
+ MI_K_FACTOR_SCHEME = 1 ! K-factor scheme
+ REGULATE_XS = 0 ! regulate cross section
+ XS_REGULATOR = QCD_Trivial ! regulation scheme
+ XS_REGULATION = 2.225 ! regulation parameter
+ SCALE_MIN = 2.225 ! minimum scale
+ JET_VETO = 1 ! apply jet veto
+ RESCALE_EXPONENT = 0.25 ! rescaling exponent
+ REFERENCE_SCALE = 1800.0 ! reference energy scale
+ PROFILE_FUNCTION = Double_Gaussian ! hadron profile function
+ PROFILE_PARAMETERS = 1.0 0.5 0.5 ! size (must be 1), coresize,
+ ! matter fraction
+}(mi)
Added: schools/2009-Lund/setups/sherpa/wjets/Analysis.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/wjets/Analysis.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,8 @@
+BEGIN_ANALYSIS {
+
+ LEVEL Hadron;
+
+ // Total cross section
+ Statistics FinalState
+
+} END_ANALYSIS;
Added: schools/2009-Lund/setups/sherpa/wjets/Run.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/wjets/Run.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,123 @@
+(run){
+ EVENTS = 1000 ! Number of Events
+ OUTPUT = 2 ! Output level:
+ ! 0 = errors, 1 = events, 2 = run info,
+ ! 3 = 1+2 = events + run info
+ ANALYSIS = 0 ! Built-in analysis (see Analysis.dat): 0=off, 1=on
+ FILE_SIZE = 10000000 ! Events per file (if written out)
+}(run)
+
+(beam){
+ BEAM_1 = 2212 ! beam particle: 2212 = P+, -2212 = P-, 11 = e-
+ BEAM_ENERGY_1 = 980. ! in GeV
+ K_PERP_MEAN_1 = 0.33 ! Settings for intrinsic kperp model
+ K_PERP_SIGMA_1 = 0.8 !
+
+ BEAM_2 = -2212
+ BEAM_ENERGY_2 = 980.
+ K_PERP_MEAN_2 = 0.33
+ K_PERP_SIGMA_2 = 0.8
+}(beam)
+
+(processes){
+ %-------------------------------------------------------
+ %-- Processes to calculate -----------------------------
+ %-------------------------------------------------------
+
+ %
+ % jet jet -> electron neutrino + 0,1,2,3 jets
+ %
+ Process : 93 93 -> 11 -12 93{3}
+ Integration_Error 0.0025 {3} ! Accuracy to which the different final
+ Integration_Error 0.005 {4} ! state multiplicities are integrated
+ Integration_Error 0.01 {5}
+ Order electroweak : 2
+ End process
+
+ %
+ % jet jet -> positron neutrino + 0,1,2,3 jets
+ %
+ Process : 93 93 -> -11 12 93{3}
+ Integration_Error 0.0025 {3} ! Accuracy to which the different final
+ Integration_Error 0.005 {4} ! state multiplicities are integrated
+ Integration_Error 0.01 {5}
+ Order electroweak : 2
+ End process
+}(processes)
+
+(selector){
+ !-------------------------------------------------------
+ !-- Parton level cuts ----------------------------------
+ !-------------------------------------------------------
+ JetFinder sqr(15/E_CMS) 1
+}(selector)
+
+(model){
+ !-------------------------------------------------------
+ !-- Model parameters -----------------------------------
+ !-------------------------------------------------------
+ MODEL = SM ! Model
+ STABLE[15] = 0 ! Set tau lepton (PDG-Id 15) unstable
+}(model)
+
+(isr){
+ PDF_SET = cteq6l.LHpdf ! Choose PDF set to use with LHAPDF
+}(isr)
+
+(me){
+ !-------------------------------------------------------
+ !-- ME generators --------------------------------------
+ !-------------------------------------------------------
+ EVENT_GENERATION_MODE = Unweighted
+ SCALE_SCHEME = CKKW ! see 'http://www.sherpa-mc.de/scales.html'
+}(me)
+
+(integration){
+ FINISH_OPTIMIZATION = Off ! don't wait for complete phase space generator optimization
+}(integration)
+
+(shower){
+ !-------------------------------------------------------
+ !-- Parton showers -------------------------------------
+ !-------------------------------------------------------
+ SHOWER_GENERATOR = Apacic ! Apacic
+ !
+ FSR_SHOWER = 1 ! 1=On,0=Off
+ ISR_SHOWER = 1 ! 1=On,0=Off
+}(shower)
+
+(fragmentation){
+ !-------------------------------------------------------
+ !-- Fragmentation parameters----------------------------
+ !-------------------------------------------------------
+ FRAGMENTATION = Ahadic
+ DECAYMODEL = Hadrons
+ YFS_MODE = 2 ! Soft QED on hadron decays
+}(fragmentation)
+
+(mi){
+ !-------------------------------------------------------
+ ! Underlying Event Setup
+ !-------------------------------------------------------
+ !
+ ! general parameters
+ !
+ MI_HANDLER = None ! Amisic / None
+ !
+ ! hard underlying event parameters
+ !
+ CREATE_GRID 93 93 -> 93 93 ! processes to generate
+ PS_ERROR = 1.0e-2 ! error for integration
+ MI_SCALE_SCHEME = G_MEAN_PT2 ! Scale scheme
+ MI_K_FACTOR_SCHEME = 1 ! K-factor scheme
+ REGULATE_XS = 0 ! regulate cross section
+ XS_REGULATOR = QCD_Trivial ! regulation scheme
+ XS_REGULATION = 2.225 ! regulation parameter
+ SCALE_MIN = 2.225 ! minimum scale
+ JET_VETO = 1 ! apply jet veto
+ RESCALE_EXPONENT = 0.25 ! rescaling exponent
+ REFERENCE_SCALE = 1800.0 ! reference energy scale
+ PROFILE_FUNCTION = Double_Gaussian ! hadron profile function
+ PROFILE_PARAMETERS = 1.0 0.5 0.5 ! size (must be 1), coresize,
+ ! matter fraction
+}(mi)
Added: schools/2009-Lund/setups/sherpa/zjets/Analysis.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/zjets/Analysis.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,46 @@
+BEGIN_ANALYSIS {
+
+ LEVEL Hadron;
+
+ // Total cross section
+ Statistics FinalState
+
+ // Kt-Jets
+ Trigger {
+ InList FinalState
+ OutList kt_lepton_jets
+ JetMode 0
+ Finder 93 10. -5. 5. 0.4 1
+ }
+ Trigger {
+ InList kt_lepton_jets
+ OutList kt_jets
+ Keep 93
+ }
+
+ // Z variables
+ Mass 11 -11 50. 130. 20 LinErr FinalState
+ PT2 11 -11 0. 200. 50 LinErr FinalState
+ PT2 11 -11 0. 20. 20 LinErr FinalState
+ Eta2 11 -11 -6. 6. 48 LinErr FinalState
+
+ // electron variables
+ PT 11 0. 100. 50 LinErr FinalState
+ PT -11 0. 100. 50 LinErr FinalState
+ Eta 11 -6. 6. 48 LinErr FinalState
+ Eta -11 -6. 6. 48 LinErr FinalState
+
+ // electron-jet separation
+ DR 11 93 0. 5. 25 LinErr kt_lepton_jets
+ DR -11 93 0. 5. 25 LinErr kt_lepton_jets
+
+ // jet variables
+ Multi -0.5 19.5 20 LinErr kt_jets
+ JetPT 0. 200. 50 1 1 2 LinErr kt_jets
+ JetEta -6. 6. 48 1 1 2 LinErr kt_jets
+ JetPhi -3.2 3.2 16 1 1 2 LinErr kt_jets
+ JetDR 0. 5. 25 1 2 2 LinErr kt_jets
+ JetDEta 0. 5. 25 1 2 2 LinErr kt_jets
+ JetDPhi 0. 3.2 16 1 2 2 LinErr kt_jets
+
+} END_ANALYSIS;
Added: schools/2009-Lund/setups/sherpa/zjets/Run.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/sherpa/zjets/Run.dat Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,111 @@
+(run){
+ EVENTS = 1000 ! Number of Events
+ OUTPUT = 2 ! Output level:
+ ! 0 = errors, 1 = events, 2 = run info,
+ ! 3 = 1+2 = events + run info
+ ANALYSIS = 0 ! Built-in analysis (see Analysis.dat): 0=off, 1=on
+}(run)
+
+(beam){
+ BEAM_1 = 2212 ! beam particle: 2212 = P+, -2212 = P-, 11 = e-
+ BEAM_ENERGY_1 = 7000. ! in GeV
+ K_PERP_MEAN_1 = 0.33 ! Settings for intrinsic kperp model
+ K_PERP_SIGMA_1 = 0.8 !
+
+ BEAM_2 = 2212
+ BEAM_ENERGY_2 = 7000.
+ K_PERP_MEAN_2 = 0.33
+ K_PERP_SIGMA_2 = 0.8
+}(beam)
+
+(processes){
+ %-------------------------------------------------------
+ %-- Hard scattering process(es) ------------------------
+ %-------------------------------------------------------
+ %
+ % jet jet -> electron positron + 0,1 jets
+ %
+ Process : 93 93 -> 11 -11 93{1}
+ Order electroweak : 2
+ Integration_Error 0.0025 {3} ! Accuracy to which the different final
+ Integration_Error 0.005 {4} ! state multiplicities are integrated
+ End process
+}(processes)
+
+(selector){
+ !-------------------------------------------------------
+ !-- Parton level cuts ----------------------------------
+ !-------------------------------------------------------
+ JetFinder sqr(20/E_CMS) 1
+ Mass 11 -11 66 116
+}(selector)
+
+(model){
+ !-------------------------------------------------------
+ !-- Model parameters -----------------------------------
+ !-------------------------------------------------------
+ MODEL = SM ! Model
+ STABLE[15] = 0 ! Set tau lepton (PDG-Id 15) unstable
+}(model)
+
+(isr){
+ PDF_SET = cteq6l.LHpdf ! Choose PDF set to use with LHAPDF
+}(isr)
+
+(me){
+ !-------------------------------------------------------
+ !-- ME generators --------------------------------------
+ !-------------------------------------------------------
+ EVENT_GENERATION_MODE = Unweighted
+ SCALE_SCHEME = CKKW ! see 'http://www.sherpa-mc.de/scales.html'
+}(me)
+
+(integration){
+ FINISH_OPTIMIZATION = Off ! don't wait for complete phase space generator optimization
+}(integration)
+
+(shower){
+ !-------------------------------------------------------
+ !-- Parton showers -------------------------------------
+ !-------------------------------------------------------
+ SHOWER_GENERATOR = Apacic ! Apacic
+ !
+ FSR_SHOWER = 1 ! 1=On,0=Off
+ ISR_SHOWER = 1 ! 1=On,0=Off
+}(shower)
+
+(fragmentation){
+ !-------------------------------------------------------
+ !-- Fragmentation parameters----------------------------
+ !-------------------------------------------------------
+ FRAGMENTATION = Ahadic
+ DECAYMODEL = Hadrons
+ YFS_MODE = 2 ! Soft QED on hadron decays
+}(fragmentation)
+
+(mi){
+ !-------------------------------------------------------
+ ! Underlying Event Setup
+ !-------------------------------------------------------
+ !
+ ! general parameters
+ !
+ MI_HANDLER = None ! Amisic / None
+ !
+ ! hard underlying event parameters
+ !
+ CREATE_GRID 93 93 -> 93 93 ! processes to generate
+ PS_ERROR = 1.0e-2 ! error for integration
+ MI_SCALE_SCHEME = G_MEAN_PT2 ! Scale scheme
+ MI_K_FACTOR_SCHEME = 1 ! K-factor scheme
+ REGULATE_XS = 0 ! regulate cross section
+ XS_REGULATOR = QCD_Trivial ! regulation scheme
+ XS_REGULATION = 2.225 ! regulation parameter
+ SCALE_MIN = 2.225 ! minimum scale
+ JET_VETO = 1 ! apply jet veto
+ RESCALE_EXPONENT = 0.25 ! rescaling exponent
+ REFERENCE_SCALE = 1800.0 ! reference energy scale
+ PROFILE_FUNCTION = Double_Gaussian ! hadron profile function
+ PROFILE_PARAMETERS = 1.0 0.5 0.5 ! size (must be 1), coresize,
+ ! matter fraction
+}(mi)
Added: schools/2009-Lund/setups/ttbar_rivetplugin/Makefile
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/ttbar_rivetplugin/Makefile Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,2 @@
+libRivetTopAnalysis.so: TopAnalysis.cc TopAnalysis.hh
+ g++ -o libRivetTopAnalysis.so -shared -I /usr/local/mcnet/local/include TopAnalysis.cc
Added: schools/2009-Lund/setups/ttbar_rivetplugin/TopAnalysis.cc
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/ttbar_rivetplugin/TopAnalysis.cc Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,183 @@
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/VetoedFinalState.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Projections/ChargedLeptons.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/AnalysisLoader.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "TopAnalysis.hh"
+
+
+namespace Rivet {
+
+ /// Constructor:
+ //
+ // Projections are a way of selecting final states or calculating observables,
+ // and they provide automatic result caching. We use two projections (well,
+ // acutally four):
+ //
+ // - The ChargedFinalState is used to select charged particles within
+ // |eta|<3.5 and pT>25 GeV out of which the ChargedLeptons projection picks
+ // only the leptons. The ChargedLeptons projection then is registered as "LFS".
+ //
+ // - The FinalState projection selects all stable particles (charged and
+ // neutral) within |eta|<2.5, any pT. We then use this in a "vetoed" final
+ // state to exclude the neutrinos, since they are not seen by the
+ // experiment. The remaining visible particles are clustered into jets
+ // using Fastjet's kT cluster algorithm. The jet projection is then
+ // registered as "JETS".
+ //
+ TopAnalysis::TopAnalysis() {
+ addProjection(ChargedLeptons(ChargedFinalState(-3.5, 3.5, 25*GeV)), "LFS");
+ VetoedFinalState visiblefs = VetoedFinalState(FinalState(-2.5, 2.5, 0*GeV));
+ visiblefs.addVetoPairId(NU_E)
+ .addVetoPairId(NU_MU)
+ .addVetoPairId(NU_TAU);
+ addProjection(FastJets(visiblefs, FastJets::KT, 0.5), "JETS");
+ }
+
+ /// Initialisation. This method is called exactly once at the
+ /// beginning of your run. That's where you should book histograms.
+ /// Don't forget to declare the histograms in the header file.
+ void TopAnalysis::init() {
+ // bookHistogram1D expects a unique histogram name, a histogram title, x-axis label,
+ // y-axis label, number of bins, lower edge and upper edge.
+ _h_jet_0_pT = bookHistogram1D("jet_0_pT", "Jet 0 pT","jet $p_T$/GeV","", 50, 0, 250);
+ _h_jet_1_pT = bookHistogram1D("jet_1_pT", "Jet 1 pT","jet $p_T$/GeV","", 50, 0, 250);
+ _h_jet_2_pT = bookHistogram1D("jet_2_pT", "Jet 2 pT","jet $p_T$/GeV","", 50, 0, 250);
+ _h_jet_3_pT = bookHistogram1D("jet_3_pT", "Jet 3 pT","jet $p_T$/GeV","", 50, 0, 250);
+
+ _h_bjet_0_pT = bookHistogram1D("bjet_0_pT", "b-jet 0 pT","jet $p_T$/GeV","", 50, 0, 250);
+ _h_bjet_1_pT = bookHistogram1D("bjet_1_pT", "b-jet 1 pT","jet $p_T$/GeV","", 50, 0, 250);
+
+ _h_ljet_0_pT = bookHistogram1D("ljet_0_pT", "light jet 0 pT","jet $p_T$/GeV","", 50, 0, 250);
+ _h_ljet_1_pT = bookHistogram1D("ljet_1_pT", "light jet 1 pT","jet $p_T$/GeV","", 50, 0, 250);
+
+ _h_W_mass = bookHistogram1D("W_mass", "$W$ mass","$m_W$/GeV","", 75, 30, 180);
+ _h_t_mass = bookHistogram1D("t_mass", "$t$ mass","$m_t$/GeV","", 150, 120, 420);
+ _h_t_mass_W_cut = bookHistogram1D("t_mass_W_cut", "$t$ mass after $m_W$ cut","$m_t$/GeV","", 150, 120, 420);
+ _h_W_comb_mass = bookHistogram1D("W_comb_mass", "combinatoric di-jet mass","$m_W$/GeV","", 75, 30, 180);
+ _h_t_comb_mass = bookHistogram1D("t_comb_mass", "combinatoric tri-jet mass","$m_t$/GeV","", 150, 120, 420);
+ }
+
+ /// Analysis method. This method is called for each event.
+ void TopAnalysis::analyze(const Event& event) {
+ // Rivet supports weighted events, so let's get the event weight (even though
+ // in this tutorial all events are unweighted):
+ double weight = event.weight();
+
+ // First we want to apply the lepton cut, so we need to select electrons by
+ // applying the ChargedLeptons projection and then only take the electrons:
+ const ChargedLeptons& lfs = applyProjection<ChargedLeptons>(event, "LFS");
+ ParticleVector electrons;
+ foreach (Particle lepton, lfs.chargedLeptons()) { // loop over leptons
+ if ( abs(lepton.pdgId()) == ELECTRON ) { // select electrons
+ electrons.push_back(lepton); // store electrons for later use
+ }
+ }
+ if (electrons.size() != 1) { // we want exactly one electron
+ vetoEvent(event);
+ }
+
+ // Now it's time to look at the jets. First we run the jet finding:
+ const FastJets& jetpro = applyProjection<FastJets>(event, "JETS");
+
+ // We want to have jets sorted by pT:
+ const Jets jets = jetpro.jetsByPt();
+
+ // If we have less than four jets, we reject the event right away:
+ if (jets.size() < 4) {
+ vetoEvent(event);
+ }
+
+ // Fill jet pT histograms for the first four leading jets
+ _h_jet_0_pT->fill(jets[0].momentum().pT(), weight);
+ _h_jet_1_pT->fill(jets[1].momentum().pT(), weight);
+ _h_jet_2_pT->fill(jets[2].momentum().pT(), weight);
+ _h_jet_3_pT->fill(jets[3].momentum().pT(), weight);
+
+ // Require the first jet to be >60GeV, the second >50GeV, the 3rd and 4th >30GeV:
+ if (jets[0].momentum().pT() < 60 || jets[1].momentum().pT() < 50 || jets[3].momentum().pT() < 30) {
+ vetoEvent(event);
+ }
+
+ // Check if jets with pT>30GeV contain a b-quark
+ // Here we are cheating a bit -- containsBottom() looks at the history
+ // and checks if there is a b-meson inside the jet. That's because real
+ // b-tagging on MC is not trivial.
+ Jets bjets, ljets;
+ foreach (Jet jet, jets) {
+ if (jet.containsBottom() && jet.momentum().pT() > 30)
+ bjets.push_back(jet); // that's a b-jet
+ else
+ ljets.push_back(jet); // that's a light jet
+ }
+
+
+ // We want exactly two b-jets:
+ if (bjets.size() !=2) {
+ vetoEvent(event);
+ }
+
+ // Fill the jet pT of the b-jets and of the two leading light jets into histograms
+ _h_bjet_0_pT->fill(bjets[0].momentum().pT(), weight);
+ _h_bjet_1_pT->fill(bjets[1].momentum().pT(), weight);
+ _h_ljet_0_pT->fill(ljets[0].momentum().pT(), weight);
+ _h_ljet_1_pT->fill(ljets[1].momentum().pT(), weight);
+
+ // Calculate the invariant masses of the W (m(jj)) and the two possible combinations
+ // of the top (m(bjj))
+ FourVector W = ljets[0].momentum() + ljets[1].momentum();
+ FourVector t1 = W + bjets[0].momentum();
+ FourVector t2 = W + bjets[1].momentum();
+
+ // Plot those masses
+ _h_W_mass->fill(mass(W), weight);
+ _h_t_mass->fill(mass(t1), weight);
+ _h_t_mass->fill(mass(t2), weight);
+
+ // Apply a W mass window cut and plot the top mass after that cut
+ if (mass(W) > 65 && mass(W) < 95) {
+ _h_t_mass_W_cut->fill(mass(t1), weight);
+ _h_t_mass_W_cut->fill(mass(t2), weight);
+ }
+
+ // Look at the combinatoric di-jet masses of the four leading jets
+ _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[2].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[3].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[1].momentum() + jets[2].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[1].momentum() + jets[3].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[2].momentum() + jets[3].momentum()), weight);
+
+ // Look at the combinatoric three jet masses of the four leading jets
+ _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum() + jets[2].momentum()), weight);
+ _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum() + jets[3].momentum()), weight);
+ _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[2].momentum() + jets[3].momentum()), weight);
+ _h_t_comb_mass->fill(mass(jets[1].momentum() + jets[2].momentum() + jets[3].momentum()), weight);
+ }
+
+ /// The finalize() method is called exactly once at the end
+ /// of the run. Good place to normalise histograms if you need to.
+ void TopAnalysis::finalize() {
+ scale(_h_jet_0_pT, 1.0/sumOfWeights());
+ scale(_h_jet_1_pT, 1.0/sumOfWeights());
+ scale(_h_jet_2_pT, 1.0/sumOfWeights());
+ scale(_h_jet_3_pT, 1.0/sumOfWeights());
+
+ scale(_h_bjet_0_pT, 1.0/sumOfWeights());
+ scale(_h_bjet_1_pT, 1.0/sumOfWeights());
+
+ scale(_h_ljet_0_pT, 1.0/sumOfWeights());
+ scale(_h_ljet_1_pT, 1.0/sumOfWeights());
+
+ scale(_h_W_mass, 1.0/sumOfWeights());
+ scale(_h_t_mass, 1.0/sumOfWeights());
+ scale(_h_W_comb_mass, 1.0/sumOfWeights());
+ scale(_h_t_comb_mass, 1.0/sumOfWeights());
+ scale(_h_t_mass_W_cut, 1.0/sumOfWeights());
+ }
+
+
+} // end of Rivet namespace
Added: schools/2009-Lund/setups/ttbar_rivetplugin/TopAnalysis.hh
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2009-Lund/setups/ttbar_rivetplugin/TopAnalysis.hh Tue Jul 17 14:05:03 2012 (r3829)
@@ -0,0 +1,112 @@
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Projections/ChargedLeptons.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/AnalysisLoader.hh"
+#include "Rivet/RivetAIDA.hh"
+
+
+namespace Rivet {
+
+ class TopAnalysis : public Analysis {
+
+ // Here you can declare your histograms:
+ private:
+ AIDA::IHistogram1D * _h_jet_0_pT;
+ AIDA::IHistogram1D * _h_jet_1_pT;
+ AIDA::IHistogram1D * _h_jet_2_pT;
+ AIDA::IHistogram1D * _h_jet_3_pT;
+
+ AIDA::IHistogram1D * _h_bjet_0_pT;
+ AIDA::IHistogram1D * _h_bjet_1_pT;
+
+ AIDA::IHistogram1D * _h_ljet_0_pT;
+ AIDA::IHistogram1D * _h_ljet_1_pT;
+
+ AIDA::IHistogram1D * _h_W_mass;
+ AIDA::IHistogram1D * _h_t_mass;
+ AIDA::IHistogram1D * _h_W_comb_mass;
+ AIDA::IHistogram1D * _h_t_comb_mass;
+ AIDA::IHistogram1D * _h_t_mass_W_cut;
+
+ /*************************************************/
+ // Anything below this line is just meta-data or needed
+ // for the Rivet plugin mechanism. You don't need to
+ // change anything here.
+ public:
+
+ /// Default constructor
+ TopAnalysis();
+
+
+ /// Factory method
+ static Analysis* create() {
+ return new TopAnalysis();
+ }
+
+
+ /// Return the name of this analysis
+ string name() const {
+ return "TOPANALYSIS";
+ }
+
+ /// Get the SPIRES ID code
+ string spiresId() const {
+ return "NONE";
+ }
+
+ /// Get a description of the analysis.
+ string description() const {
+ return "MCnet school tutorial analysis";
+ }
+
+ /// Experiment which performed and published this analysis.
+ string experiment() const {
+ return "NONE";
+ }
+
+ /// When published (preprint year according to SPIRES).
+ string year() const {
+ return "NONE";
+ }
+
+ /// Summary of analysis
+ string summary() const{
+ return "Sample analysis for the Lund MCnet school tutorials.";
+ }
+
+ /// Beam conditions for this analysis
+ string runInfo() const{
+ return "inclusive ttbar or W+4jets";
+ }
+
+ string collider() const{
+ return "LHC 14 TeV";
+ }
+
+ vector<string> authors() const{
+ return vector<string>();
+ }
+
+ vector<string> references() const{
+ return vector<string>();
+ }
+
+ // Analysis methods
+ void init();
+ void analyze(const Event& event);
+ void finalize();
+
+ };
+
+ extern "C" {
+ AnalysisBuilders getAnalysisBuilders() {
+ AnalysisBuilders fns;
+ fns["TOPANALYSIS"] = Rivet::TopAnalysis::create;
+ return fns;
+ }
+ }
+
+
+}
More information about the Rivet-svn
mailing list
|