herwig is hosted by Hepforge, IPPP Durham
close Warning: Can't synchronize with repository "(default)" (Repository path '/hepforge/hg/herwig/public/herwig' does not exist.). Look in the Trac log for more information.

Ticket #10: ttb_dec_wspin.in

File ttb_dec_wspin.in, 4.9 KB (added by new, 6 years ago)
Line 
1##################################################
2# Powheg Input to Herwig
3##################################################
4read snippets/PPCollider.in
5
6#######################################################
7# Load libraries needed to read LesHouces Event Files #
8######################################################
9library LesHouches.so
10
11#####################################
12# Technical parameters for this run #
13#####################################
14cd /Herwig/Generators
15set EventGenerator:NumberOfEvents 50000
16set EventGenerator:DebugLevel 0
17set EventGenerator:PrintEvent 10
18set EventGenerator:MaxErrors 10000
19
20#############################################################
21# Create a LesHouches File Handler
22#############################################################
23cd /Herwig/EventHandlers
24create ThePEG::LesHouchesEventHandler _lhefHandler
25set _lhefHandler:PartonExtractor /Herwig/Partons/PPExtractor
26set _lhefHandler:CascadeHandler /Herwig/Shower/ShowerHandler
27set _lhefHandler:DecayHandler /Herwig/Decays/DecayHandler
28set _lhefHandler:HadronizationHandler /Herwig/Hadronization/ClusterHadHandler
29
30set _lhefHandler:WeightOption VarNegWeight #for MC@NLO
31#set _lhefHandler:WeightWarnings 0
32
33set /Herwig/Generators/EventGenerator:EventHandler /Herwig/EventHandlers/_lhefHandler
34
35###########################################################################
36# Set cuts to match those used when generating the events with Powheg-Box #
37# Set up an EMPTY CUTS object in particular for POWHEG and MC@NLO
38###########################################################################
39create ThePEG::Cuts /Herwig/Cuts/NoCuts
40#set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
41#set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
42#set /Herwig/Cuts/PhotonKtCut:MinKT 0.0*GeV
43#set /Herwig/Cuts/LeptonKtCut:MinKT 0.0*GeV
44#set /Herwig/Cuts/WBosonKtCut:MinKT 0.0*GeV
45#set /Herwig/Cuts/MassCut:MinM 0.0*GeV
46#set lhefReader:Cuts /Herwig/Cuts/QCDCuts
47
48#################################################
49# Set the same PDFs as those used in Powheg-Box #
50#################################################
51create ThePEG::LHAPDF /Herwig/Partons/LHAPDF ThePEGLHAPDF.so
52set /Herwig/Partons/LHAPDF:PDFName MSTW2008nlo68cl.LHgrid
53set /Herwig/Partons/LHAPDF:RemnantHandler /Herwig/Partons/HadronRemnants
54#set /Herwig/Particles/p+:PDF /Herwig/Partons/LHAPDF
55#set /Herwig/Particles/pbar-:PDF /Herwig/Partons/LHAPDF
56#set /Herwig/Partons/PPExtractor:FirstPDF  /Herwig/Partons/LHAPDF
57#set /Herwig/Partons/PPExtractor:SecondPDF /Herwig/Partons/LHAPDF
58#set /Herwig/Shower/ShowerHandler:PDFA /Herwig/Partons/LHAPDF
59#set /Herwig/Shower/ShowerHandler:PDFB /Herwig/Partons/LHAPDF
60
61#############################################################
62# Create a LesHouches reader object and open the input file #
63#############################################################
64create ThePEG::LesHouchesFileReader _lhefReader
65set _lhefReader:FileName pwhgevent_1.lhe
66set _lhefReader:AllowedToReOpen No
67set _lhefReader:CacheFileName cacheevents.tmp
68
69####
70set _lhefReader:Cuts /Herwig/Cuts/NoCuts
71set _lhefReader:InitPDFs 0
72set _lhefReader:PDFA /Herwig/Partons/LHAPDF
73set _lhefReader:PDFB /Herwig/Partons/LHAPDF
74####
75
76insert _lhefHandler:LesHouchesReaders 0 _lhefReader
77
78############################
79# Set up the POWHEG shower #
80############################
81#cd /Herwig/Shower
82# normally, especially for POWHEG, you want
83# the scale supplied in the event files (SCALUP)
84# to be used as a pT veto scale in the parton shower
85set /Herwig/Shower/ShowerHandler:MaxPtIsMuF Yes
86set /Herwig/Shower/ShowerHandler:RestrictPhasespace Yes
87# Shower parameters
88# treatment of wide angle radiation
89set /Herwig/Shower/PartnerFinder:PartnerMethod Random
90set /Herwig/Shower/PartnerFinder:ScaleChoice Partner
91# fix issue before 7.0.5 (not needed after this)
92set /Herwig/Shower/GtoQQbarSplitFn:AngularOrdered Yes
93set /Herwig/Shower/GammatoQQbarSplitFn:AngularOrdered Yes
94# with MC@NLO these parameters are required for consistency of the subtraction terms
95# suggested parameters (give worse physics results with POWHEG)
96#set /Herwig/Shower/KinematicsReconstructor:InitialInitialBoostOption LongTransBoost
97#set /Herwig/Shower/KinematicsReconstructor:ReconstructionOption General
98#set /Herwig/Shower/KinematicsReconstructor:FinalStateReconOption Default
99#set /Herwig/Shower/KinematicsReconstructor:InitialStateReconOption Rapidity
100#set /Herwig/Shower/ShowerHandler:SpinCorrelations No
101
102
103#################################
104# Set up LHC physics parameters #
105#################################
106set /Herwig/Particles/t:NominalMass 173.2
107#cd /Herwig/Generators
108
109
110#############################
111# Set up the RIVET analysis #
112#############################
113#create ThePEG::RivetAnalysis /Herwig/Analysis/RivetAnalysis RivetAnalysis.so
114#insert /Herwig/Analysis/RivetAnalysis:Analyses 0 W_ANALYSIS
115create ThePEG::HepMCFile HepMC HepMCAnalysis.so
116insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMC
117
118
119################
120# Save the run #
121################
122cd /Herwig/Generators
123saverun ttb_dec_wspin EventGenerator
124
125