Ticket #10: ttb_dec_wspin.in
File ttb_dec_wspin.in, 4.9 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:DebugLevel 0 |
17 | set EventGenerator:PrintEvent 10 |
18 | set EventGenerator:MaxErrors 10000 |
19 | |
20 | ############################################################# |
21 | # Create a LesHouches File Handler |
22 | ############################################################# |
23 | cd /Herwig/EventHandlers |
24 | create ThePEG::LesHouchesEventHandler _lhefHandler |
25 | set _lhefHandler:PartonExtractor /Herwig/Partons/PPExtractor |
26 | set _lhefHandler:CascadeHandler /Herwig/Shower/ShowerHandler |
27 | set _lhefHandler:DecayHandler /Herwig/Decays/DecayHandler |
28 | set _lhefHandler:HadronizationHandler /Herwig/Hadronization/ClusterHadHandler |
29 | |
30 | set _lhefHandler:WeightOption VarNegWeight #for MC@NLO |
31 | #set _lhefHandler:WeightWarnings 0 |
32 | |
33 | set /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 | ########################################################################### |
39 | create 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 | ################################################# |
51 | create ThePEG::LHAPDF /Herwig/Partons/LHAPDF ThePEGLHAPDF.so |
52 | set /Herwig/Partons/LHAPDF:PDFName MSTW2008nlo68cl.LHgrid |
53 | set /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 | ############################################################# |
64 | create ThePEG::LesHouchesFileReader _lhefReader |
65 | set _lhefReader:FileName pwhgevent_1.lhe |
66 | set _lhefReader:AllowedToReOpen No |
67 | set _lhefReader:CacheFileName cacheevents.tmp |
68 | |
69 | #### |
70 | set _lhefReader:Cuts /Herwig/Cuts/NoCuts |
71 | set _lhefReader:InitPDFs 0 |
72 | set _lhefReader:PDFA /Herwig/Partons/LHAPDF |
73 | set _lhefReader:PDFB /Herwig/Partons/LHAPDF |
74 | #### |
75 | |
76 | insert _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 |
85 | set /Herwig/Shower/ShowerHandler:MaxPtIsMuF Yes |
86 | set /Herwig/Shower/ShowerHandler:RestrictPhasespace Yes |
87 | # Shower parameters |
88 | # treatment of wide angle radiation |
89 | set /Herwig/Shower/PartnerFinder:PartnerMethod Random |
90 | set /Herwig/Shower/PartnerFinder:ScaleChoice Partner |
91 | # fix issue before 7.0.5 (not needed after this) |
92 | set /Herwig/Shower/GtoQQbarSplitFn:AngularOrdered Yes |
93 | set /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 | ################################# |
106 | set /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 |
115 | create ThePEG::HepMCFile HepMC HepMCAnalysis.so |
116 | insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMC |
117 | |
118 | |
119 | ################ |
120 | # Save the run # |
121 | ################ |
122 | cd /Herwig/Generators |
123 | saverun ttb_dec_wspin EventGenerator |
124 | |
125 |