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 #11: ttb_dec_wspin.in

File ttb_dec_wspin.in, 5.1 KB (added by new, 6 years ago)

Herwig input file

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