Ticket #11: ttb_dec_wspin.in
File ttb_dec_wspin.in, 5.1 KB (added by , 6 years ago) |
---|
Line | |
---|---|
1 | ################################################## |
2 | # Powheg Input to Herwig |
3 | ################################################## |
4 | read snippets/PPCollider.in |
5 | |
6 | ####################################################### |
7 | # Load libraries needed to read LesHouces Event Files # |
8 | ###################################################### |
9 | library LesHouches.so |
10 | |
11 | ##################################### |
12 | # Technical parameters for this run # |
13 | ##################################### |
14 | cd /Herwig/Generators |
15 | set EventGenerator:NumberOfEvents 50000 |
16 | set EventGenerator:RandomNumberGenerator:Seed 508680 |
17 | set EventGenerator:DebugLevel 0 |
18 | set EventGenerator:PrintEvent 50000 |
19 | set EventGenerator:MaxErrors 1000 |
20 | |
21 | |
22 | ############################################################# |
23 | # Create a LesHouches File Handler |
24 | ############################################################# |
25 | cd /Herwig/EventHandlers |
26 | create ThePEG::LesHouchesEventHandler _lhefHandler |
27 | set _lhefHandler:PartonExtractor /Herwig/Partons/PPExtractor |
28 | set _lhefHandler:CascadeHandler /Herwig/Shower/ShowerHandler |
29 | set _lhefHandler:DecayHandler /Herwig/Decays/DecayHandler |
30 | set _lhefHandler:HadronizationHandler /Herwig/Hadronization/ClusterHadHandler |
31 | |
32 | set _lhefHandler:WeightOption VarNegWeight #for MC@NLO |
33 | #set _lhefHandler:WeightWarnings 0 |
34 | |
35 | set /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 | ########################################################################### |
41 | create 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 | ################################################# |
53 | create ThePEG::LHAPDF /Herwig/Partons/LHAPDF ThePEGLHAPDF.so |
54 | set /Herwig/Partons/LHAPDF:PDFName MSTW2008nlo68cl.LHgrid |
55 | set /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 | ############################################################# |
66 | create ThePEG::LesHouchesFileReader _lhefReader |
67 | set _lhefReader:FileName pwhgevent_1.lhe |
68 | set _lhefReader:AllowedToReOpen No |
69 | set _lhefReader:CacheFileName cacheevents.tmp |
70 | |
71 | #### |
72 | set _lhefReader:Cuts /Herwig/Cuts/NoCuts |
73 | set _lhefReader:InitPDFs 0 |
74 | set _lhefReader:PDFA /Herwig/Partons/LHAPDF |
75 | set _lhefReader:PDFB /Herwig/Partons/LHAPDF |
76 | #### |
77 | |
78 | insert _lhefHandler:LesHouchesReaders 0 _lhefReader |
79 | |
80 | ############################ |
81 | # Set up the POWHEG shower # |
82 | ############################ |
83 | cd /Herwig/Shower |
84 | set 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 |
89 | set /Herwig/Shower/ShowerHandler:MaxPtIsMuF Yes |
90 | set /Herwig/Shower/ShowerHandler:RestrictPhasespace Yes |
91 | # Shower parameters |
92 | # treatment of wide angle radiation |
93 | set /Herwig/Shower/PartnerFinder:PartnerMethod Random |
94 | set /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 | ################################# |
110 | set /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 |
119 | create ThePEG::HepMCFile /Herwig/Analysis/HepMC HepMCAnalysis.so |
120 | set /Herwig/Analysis/HepMC:PrintEvent 50000 |
121 | set /Herwig/Analysis/HepMC:Format GenEvent |
122 | set /Herwig/Analysis/HepMC:Units GeV_mm |
123 | set /Herwig/Analysis/HepMCFile:Filename ttb_dec_wspin_1.hepmc |
124 | insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMC |
125 | |
126 | |
127 | ################ |
128 | # Save the run # |
129 | ################ |
130 | cd /Herwig/Generators |
131 | saverun ttb_dec_wspin EventGenerator |
132 |