Home
The MC@ NLO 3.1 Event Generator
Contents
1. and it is only used by the sample analysis routines mcatnlo_hwanzzz f As such the user will not need it when linking to a self contained analysis code We also remind the reader that the HERWIG version must be 6 5 or higher since the Les Houches interface is used A 2 From MC NLO version 2 0 to version 2 1 In this appendix we list the changes that occurred in the package from version 2 0 to version 2 1 e Higgs production has been added which implies new process specific files mcatnlo_hgmain f mcatnlo_hgxsec f hgscblks h mcatnlo_hwanhgg f anda modification to mcatnlo_hwlhin f e Post 1999 PDF sets have been added to the MCQ NLO PDF library e Script variables have been added to MCatNLO inputs Most of them are only relevant to Higgs production and don t affect processes implemented in version 2 0 One of them LAMBDAHERW may affect all processes in version 2 1 the variables LAMBDAFIVE and LAMBDAHERW are used to set the value of Agcp in NLO and MC runs respectively whereas in version 2 0 LAMBDAFIVE controlled both The new setup is necessary since modern PDF sets have Agcp values which are too large to be supported by HERWIG Recall that the effect of using LAMBDAHERW different from LAMBDAFIVE is beyond NLO e The new script variable PDFPATH should be set equal to the name of the direc tory where the PDF grid files which can be downloaded from the MC NLO web page are stored At run time when executing runNLO or runMC or
2. which does include leading order spin correlations Particular decay modes of vector bosons may be forced in MCQ NLO in the same way as in standard HERWIG using the MODBOS variables see sect 3 4 of ref 20 However top decays cannot be forced in this way because the decay is treated as a three body process the W boson entry in HEPEVT is for information only In stead the top branching ratios can be altered using the HWMODK subroutine see sect 7 of ref 20 This is done separately for the t and t For example CALL HWMODK 6 1 D0 100 12 11 5 0 0 forces the decay t v etb while leaving t de cays unaffected Note that the order of the decay products is important for the decay matrix element NME 100 to be applied correctly The relevant statements should be inserted in the HERWIG main program corresponding to mcatnlo_hwdriver f in this package after the statement CALL HWUINC and before the loop over events A separate run with CALL HWMODK 6 1 D0 100 12 11 5 0 0 should be per formed if one wishes to symmetrize the forcing of t and t decays since calls to HWMODK from within the event loop do not produce the desired result 3 3 Results As in the case of standard HERWIG the form of the results will be determined by the user s analysis routines However in addition to any files written by the user s analysis routines the MCQNLO writes the following files 3Non factorizable spin correlations of virtual origi
3. corrections Phys Lett B 264 1991 440 G Altarelli R K Ellis and G Martinelli Large Perturbative Corrections To The Drell Yan Process In QCD Nucl Phys B 157 1979 461 P Aurenche and J Lindfors QCD Corrections To Direct Lepton Production In Hadronic Collisions Nucl Phys B 185 1981 274 C Oleari V Del Duca S Frixione and B R Webber MC NLO for Higgs Boson Production in preparation S Frixione and B R Webber The MC NLO event generator arXiv hep ph 0207182 S Frixione and B R Webber The MC NLO 2 0 event generator arXiv hep ph 0307146 S Frixione and B R Webber The MC NLO 2 2 event generator arXiv hep ph 0309186 S Frixione and B R Webber The MC NLO 2 3 event generator arXiv hep ph 0402116 E Boos et al Generic user process interface for event generators arXiv hep ph 0109068 G Marchesini B R Webber G Abbiendi I G Knowles M H Seymour and L Stanco HERWIG A Monte Carlo event generator for simulating hadron emission reactions with interfering gluons Version 5 1 April 1991 Comput Phys Commun 67 1992 465 G Corcella I G Knowles G Marchesini S Moretti K Odagiri P Richardson M H Seymour and B R Webber HERWIG 6 An event generator for hadron emis sion reactions with interfering gluons including supersymmetric processes JHEP 0101 2001 010 hep ph 0011363 20
4. on shell scheme was previously used for single vector boson production e The script variables FRENMC and FFACTMC have been eliminated e The structure of pseudo random number generation in heavy flavour produc tion has been changed to avoid a correlation that affected the azimuthal angle dis tribution for the products of the hard partonic subprocesses e A few minor bugs have been corrected which affected the rapidity of the vector bosons in single vector boson production a 2 3 effect and the assignment of Agop for the LO and NLO PDF sets of Alekhin A 5 From MC NLO version 2 3 to version 3 1 In this appendix we list the changes that occurred in the package from version 2 3 to version 3 1 e Associated Higgs production has been added which implies new process specific files mcatnlo_vhmain f mcatnlo_vhxsec f vhgcblks h mcatnlo_hwanvhg f and modifications to mcatnlo_hwlhin f and mcatnlo_hwdriver f e Spin correlations in WtW production and leptonic decay have been added the relevant codes mcatnlo_vpmain f mcatnlo_vhxsec f have been modified 15 the sample analysis routines mcatnlo_hwanvbp f have also been changed Tree level matrix elements have been computed with MadGraph MadEvent 23 24 which uses HELAS 25 the relevant routines and common blocks are included in mcatnlo_helas2 f and MEcoupl inc e The format of the event file has changed in several respects the most relevant of which is that the four moment
5. the user wants to create another event file to increase the statistics The user simply sets BASES OFF since the integration step is not necessary any longer however the data files must not be removed the information stored there is still used by the NLO code changes the string EVPREFIX failure to do so overwrites the existing event file while keeping FPREFIX at the same value as before and changes the value of RNDEVSEED the random number seed used in the event generation step failure to do so results in an event file identical to the previous one the number NEVENTS generated may or may not be equal to the one chosen in generating the former event file s We point out that data and event files may be very large If the user wants to store them in a scratch area this can be done by setting the script variable SCRTCH equal to the physical address of the scratch area see sect 3 3 3 2 Decays MC NLO is intended primarily for the study of NLO corrections to production cross sections and distributions NLO corrections to the decays of produced particles are not included As for spin correlations in decays the situation in version 3 1 is summarized in table 1 they are included for all processes except tt ZZ and WZ production For the latter processes quantities sensitive to the polarisation of produced particles are not given correctly even to leading order For such quantities it may be preferable to use the standard HERWIG MC
6. to fully determine the process type and new variables IV IL and IL have been introduced see table 1 The variables IL and IL can take the same values as IL relevant to lepton pair production notice however that in the latter case IL is not an independent variable and its value is included via IPROC in addition IL 7 implies that lepton spin correlations for the decay products of the corresponding vector boson are not taken into account as indicated in table 1 Apart from the above differences MCQ NLO and HERWIG behave in exactly the same way Thus the available user s analysis routines can be used in the case of MC NLO One should recall however that MC NLO always generates some events with negative weights see refs 1 2 therefore the correct distributions are obtained by summing weights with their signs i e the absolute values of the weights must NOT be used when filling the histograms With such a structure it is natural to create two separate executables which we improperly denote as NLO and MC The former has the sole scope of creating the event file the latter is just HERWIG augmented by the capability of reading the event file 1 2 Package files The package consists of the following files e Shell utilities MCatNLO Script MCatNLO inputs Makefile e Utility codes MEcoupl inc alpha f dummies f linux f mcatnlo_date f mcatnlo_hbook f mcatnlo_helas2 f mcatnlo_hwdummy f mcatnlo_int f mcatnlo_libo
7. 21 22 23 24 25 G Corcella et al HERWIG 6 5 release note arXiv hep ph 0210213 S Kawabata A New version of the multidimensional integration and event gener ation package BASES SPRING Comput Phys Commun 88 1995 309 T Stelzer and W F Long Automatic generation of tree level helicity amplitudes Comput Phys Commun 81 1994 357 arXiv hep ph 9401258 F Maltoni and T Stelzer MadEvent Automatic event generation with Mad Graph JHEP 0302 2003 027 arXiv hep ph 0208156 H Murayama I Watanabe and K Hagiwara HELAS HELicity amplitude subrou tines for Feynman diagram evaluations KEK 91 11 21
8. 4 15 16 17 respectively The new processes implemented since version 3 1 are Higgs production associated with a vector boson and WtW production with spin correlations for the decay into leptons For precise details of version changes see app A 1 A 5 1 1 Mode of operation In the case of standard MC a hard kinematic configuration is generated on a event by event basis and it is subsequently showered and hadronized In the case of MC NLO all of the hard kinematic configurations are generated in advance and stored in a file which we call event file see sect 3 1 the event file is then read by HERWIG which showers and hadronizes each hard configuration Since version 2 0 the events are handled by the Les Houches generic user process interface 18 see ref 2 for more details Therefore in MC NLO the reading of a hard configuration from the event file is equivalent to the generation of such a configuration in the standard MC The signal to HERWIG that configurations should be read from an event file using the Les Houches interface is a negative value of the process code IPROC this accounts for the negative values in table 1 In the case of heavy quark pair Higgs Higgs in association with a W or Z and lepton pair through Z 7y exchange production the codes are simply the negative of those for the corresponding standard HERWIG MC processes Where possible this convention will be adopted for additional MCQNLO processe
9. ER have been described in sect 4 Finally CERNLIB must be set in order to link the local version of CERN PDFLIB In order to create the object files eventually linked static compilation is always recommended for example g77 Wall fno automatic on Linux B 2 The input files In this appendix we describe the inputs to be given to the NLO and MC executables in the case of vector boson pair production The case of other production processes is completely analogous When the shell scripts are used to run the MCQNLO two files are created FPREFIXNLOinput and FPREFIXMCinput which are read by the NLO and MC executable respectively We start by considering the inputs for the NLO executable presented in table 2 The variables whose name is in uppercase characters have been described in sect 4 The other variables are assigned by the shell script Their default values are given in table 3 Users who run the package without the script should use the values given in table 3 The variable zi controls to a certain extent the number of negative weight events generated by the MC NLO see ref 1 Therefore the user may want to tune this parameter in order to reduce as much as possible the number of negative weight events We stress that the MC code will not change this number thus the tuning can and must be done only by running the NLO code The variables nitn control the integration step see sect 3 1 which can be skipped by setting nitn 0 I
10. WGTTYPE Valid entries are 0 and 1 When set to 0 the weights in the event file are 1 When set to 1 they are w with w a constant such that the sum of the weights gives the total NLO cross section N B These weights are redefined by HERWIG at MC run time according to its own convention see HERWIG manual RNDEVSEED The seed for the random number generation in the event generation step must be changed in order to obtain statistically equivalent but different event files BASES Controls the integration step valid entries are ON and OFF At least one run with BASES ON must be performed see sect 3 1 PDFLIBRARY Valid entries are PDFLIB and THISLIB In the former case PDFLIB is used to compute the parton densities whereas in the latter case the densities are obtained from our self contained faster package HERPDF If set to DEFAULT HERWIG uses its internal PDF set controlled by NSTRU regardless of the densities adopted at the NLO level If set to EXTPDF HERWIG uses the same PDFs as the NLO code see sect 2 1 HWPATH The physical address of the directory where the user s preferred version of HERWIG is stored SCRTCH The physical address of the directory where the user wants to store the data and event files If left blank these files are stored in the running directory HWUTI This variables must be set equal to a list of object files needed by the analysis routines of the user for example HWUTI obj1 0 obj2 0 obj3 o is a val
11. a are now given as Px Py Pz m up to version 2 3 we had pz py pz E Event files generated with version 2 3 or lower must not be used with version 3 1 or higher the code will prevent the user from doing so e The script variables GAMMAX MASSINF and MASSSUP have been replaced with xGAMMAX xMASSINF and xMASSSUP with x H V1 V2 e New script variables IVCODE IL1CODE and IL2CODE have been introduced e Minor changes have been made to the routines that put the partons on the HERWIG mass shell for lepton pair heavy quark and vector boson pair production effects are beyond the fourth digit e The default electroweak parameters have been changed for vector boson pair production in order to make them consistent with those used in other processes The cross sections are generally smaller in version 3 1 wrt previous versions the dominant effect being the value of sin w we have now sin w 0 2311 in lower versions sin Ow 1 m m The cross sections are inversely proportional to sinf Ow B Running the package without the shell scripts In this appendix we describe the actions that the user needs to take in order to run the package without using the shell scripts and the Makefile Examples are given for vector boson pair production but only trivial modifications are necessary in order to treat other production processes B 1 Creating the executables An MC NLO run requires the creation of two executables for the NLO an
12. arXiv hep ph 0506182v1 17 Jun 2005 Preprint typeset in JHEP style PAPER VERSION Cavendish HEP 05 09 GEF TH 6 2005 The MC NLO 3 1 Event Generator Stefano Frixione INFN Sezione di Genova Via Dodecaneso 33 16146 Genova Italy E mail Stefano Frixione cern ch Bryan R Webber Cavendish Laboratory Madingley Road Cambridge CB3 OHE U K E mail webber hep phy cam ac uk ABSTRACT This is the user s manual of MCQ NLO 3 1 This package is a practical implementation based upon the HERWIG event generator of the MC NLO formal ism which allows one to incorporate NLO QCD matrix elements consistently into a parton shower framework Processes available in this version include the hadropro duction of single vector and Higgs bosons vector boson pairs heavy quark pairs lepton pairs and Higgs bosons in association with a W or Z Spin correlations in decays are included for all processes except tt ZZ and WZ production This document is self contained but we emphasise the main differences with respect to previous versions KEYWORDS QCD Monte Carlo NLO Computations Resummation Hadronic Colliders Work supported in part by the UK Particle Physics and Astronomy Research Council and by the EU Fourth Framework Programme Training and Mobility of Researchers Network Quantum Chromodynamics and the Deep Structure of Elementary Particles contract FMRX CT98 0194 DG 12 MIHT Contents 1 Generalities 1 1 1 M
13. by calls to compileNLO compileMC runNLO runMC which have obvious meanings We point out that the command runMC may be used with IPROC 1350 IL 1450 IL 1600 ID 1699 1705 1706 2600 ID 2699 2700 ID 2799 to generate Z 7 W Higgs bb tt H W or H Z events with stan dard HERWIG see the HERWIG manual for more details We stress that the input parameters are not solely related to physics masses CM energy and so on there are a few of them which control other things such as the number of events generated These must also be set by the user according to his her needs see sect 4 Two such variables are HERWIGVER and HWUTI which were moved in version 2 0 from the Makefile to MCatNLO inputs The former variable must be set equal to the object file name of the version of HERWIG currently adopted matching the one whose common blocks are included in the files mentioned in sect 2 The variable HWUTI must be set equal to the list of object files that the user needs in the analysis routines If the shell scripts are not used to run the codes the inputs are given to the NLO or MC codes during an interactive talk to phase the complete sets of inputs for our codes are reported in app B 2 for vector boson pair production 3 1 Event file The NLO code creates the event file In order to do so it goes through two steps first it integrates the cross sections integration step and then using the information gathered in the int
14. c tion set it to NO to use Qem 1 137 0359895 Process number that identifies the hard subprocess see table 1 for valid entries Identifies the nature of the vector boson in associated Higgs production It corresponds to variable IV of table 1 Identify the nature of the leptons emerging from vector boson decays x 1 2 They correspond to variables IL and IL of table 1 The type of the incoming particle n with n 1 2 HERWIG naming conven tions are used P PBAR N NBAR The name of the group fitting the parton densities used the labeling conven tions of PDFLIB are adopted The number of the parton density set according to PFDLIB the pair PDFGROUP PDFSET identifies the densities for a given particle type The value of Agon for five flavours and in the MS scheme used in the compu tation of NLO cross sections The value of Agcp used in MC runs this parameter has the same meaning as Agcy in HERWIG The subtraction scheme in which the parton densities are defined Our integration routine creates files with name beginning by the string FPREFIX These files are not directly accessed by the user for more details see sect 3 1 The name of the event file begins with this string for more details see sect 3 1 The names of the NLO and MC executables begin with this string this is useful in the case of simultaneous runs The number of events stored in the event file eventually processed by HERWIG 12
15. cally at run time As stressed before consistent inputs must be given to the NLO and MC codes However in ref 1 we found that the dependence upon the PDFs used by the MC is rather weak So one may want to run the NLO and MC adopting a regular NLL evolved set in the former case and the default HERWIG set in the latter the advantage is that this option reduces the amount of running time of the MC In order to do so the user must set the variable HERPDF equal to DEFAULT in the file MCatNLO inputs setting HERPDF EXTPDF will force the MC to use the same PDF set as the NLO code Regardless of the PDFs used in the MC run users must delete the dummy PDFLIB routines PDFSET and STRUCTM from HERWIG as explained earlier In MC NLO 3 1 the PDF library LHAPDF is not supported 3 Running It is straightforward to run the MCQNLO First edit MCatNLO inputs and write there all the input parameters for the complete list of the input parame ters see sect 4 As the last line of the file MCatNLO inputs write runMCatNLO Finally execute MCatNLO inputs from the bash shell This procedure will create the NLO and MC executables and run them using the inputs given in MCatNLO inputs which guarantees that the parameters used in the NLO and MC runs are consistent Should the user only need to create the executables without running them or to run the NLO or the MC only he she should replace the call to runMCatNLO in the last line of MCatNLO inputs
16. d MC codes respectively The files to link depend on whether one uses PDFLIB or the PDF library provided with this package we list them below e NLO without PDFLIB mcatnlo_vbmain o mcatnlo_vbxsec o mcatnlo_date o mcatnlo_int o mcatnlo_uxdate o mcatnlo_uti o mcatnlo_str o mcatnlo_pdftomlm o mcatnlo_libofpdf o dummies o SYSFILE e NLO with PDFLIB mcatnlo_vbmain o mcatnlo_vbxsec o mcatnlo_date o mcatnlo_int o mcatnlo_uxdate o mcatnlo_uti o mcatnlo_str o mcatnlo_mlmtopdf o dummies o SYSFILE CERNLIB 16 e MC without PDFLIB mcatnlo_hwdriver o mcatnlo_hwlhin o mcatnlo_hwanvbp o mcatnlo_hbook o mcatnlo_str o mcatnlo_pdftomlm o mcatnlo_libofpdf o dummies o HWUTI HERWIGVER e MC with PDFLIB mcatnlo_hwdriver o mcatnlo_hwlhin o mcatnlo_hwanvbp o mcatnlo_hbook o mcatnlo_str o mcatnlo_mlmtopdf o dummies o HWUTI HERWIGVER CERNLIB The process specific codes mcatnlo_vbmain o and mcatnlo_vbxsec o for the NLO executable and mcatnlo_hwanvbp o the HERWIG analysis routines in the MC exe cutable need to be replaced by their analogues for other production processes which can be easily read from the list given in sect 1 2 The variable SYSFILE must be set either equal to alpha o or to linux o or to sun o according to the architecture of the machine on which the run is performed For any other architecture the user should provide a file corresponding to alpha f etc which he she will easily obtain by modifying alpha f The variables HWUTI and HERWIGV
17. egration step produces a set of hard events event generation step Integration and event generation are performed with a modified version of the SPRING BASES package 22 We stress that the events stored in the event file just contain the partons in volved in the hard suprocesses Owing to the modified subtraction introduced in the MC NLO formalism see ref 1 they do not correspond to pure NLO configura tions and should not be used to plot physical observables Parton level observables must be reconstructed using the fully showered events The event generation step necessarily follows the integration step however for each integration step one can have an arbitrary number of event generation steps i e an arbitrary number of event files This is useful in the case in which the statistics accumulated with a given event file is not sufficient Suppose the user wants to create an event file editing MCatNLO inputs the user sets BASES ON to enable the integration step sets the parameter NEVENTS equal to the number of events wanted on tape and runs the code the information on the integration step unreadable to the user but needed by the code in the event generation step is written on files whose name begin with FPREFIX a string the user sets in MCatNLO inputs these files which we denotes as data files have extensions data The name of the event file is EVPREFIX events where EVPREFIX is again a string set by the user Now suppose
18. f one needs to perform the integration step we suggest setting these variables as indicated in table 3 We now turn to the inputs for the MC executable presented in table 4 The variables whose names are in uppercase characters have been described in sect 4 The other variables are assigned by the shell script Their default values are given 17 FPREFIX EVPREFIX ECM FFACT FREN FFACTMC FRENMC IPROC WMASS ZMASS UMASS DMASS SMASS CMASS BMASS GMASS PART1 PART2 PDFGROUP PDFSET LAMBDAFIVE SCHEMEOFPDF NEVENTS WGTTYPE RNDEVSEED zi nitn nitns prefix for BASES files prefix for event files energy scalefactors 2850 60 70 80 WW ZZ ZW ZW M_W M_Z quark and gluon masses hadron types PDF group and id number Lambda_5 lt 0 for default scheme number of events 0 gt wgt 1 1 1 gt wgt w w seed for rnd numbers zi itmx1 itmx2 Table 2 Sample input file for the NLO code for vector boson pair production FPREFIX and EVPREFIX must be understood with SCRTCH in front see sect 4 Variable Default value Zi 0 2 nitn 10 0 BASES 0N OFF Table 3 Default values for script generated variables in FPREFIXNLOinput in table 5 The user can freely change the values of esctype and pdftype on the other hand the value of beammom must always be equal to half of the hadronic CM energy In the case of y Z W Higgs or heavy quark production the MC executable can be run with t
19. fpdf f mcatnlo_mlmtopdf f mcatnlo_pdftomlm f mcatnlo_str f mcatnlo_uti f mcatnlo_uxdate c sun f trapfpe c e General HERWIG routines mcatnlo_hwdriver f mcatnlo_hwlhin f e Process specific codes mcatnlo_hwanbtm mcatnlo_hwanhgg mcatnlo_hwanllp mcatnlo_hwantop mcatnlo_hwansvb mcatnlo_hwanvbp Hh Hh Hh hh Hh hh Fh mcatnlo_hwanvhg mcatnlo_hgmain mcatnlo_hgxsec mcatnlo_llmain mcatnlo_llxsec mcatnlo_qqmain Hh Hh Hh Fh OF Fb mcatnlo_qqxsec mcatnlo_sbmain mcatnlo_sbxsec mcatnlo_vbmain mcatnlo_vbxsec mcatnlo_vhmain Hh Hh Hh Hh Fh Fh mcatnlo_vhxsec hgscblks h hvqcblks h llpcblks h svbcblks h vhgcblks h These files can be downloaded from the web page http www hep phy cam ac uk theory webber MCatNLO The files mcatnlo_hwangzz f which appear in the list of the process specific codes are sample HERWIG analysis routines They are provided here to give the user a ready to run package but they should be replaced with appropriate codes according to the user s needs In addition to the files listed above the user will need a version of the HERWIG code 19 20 21 As stressed in ref 1 for the MC NLO we do not modify the existing LL shower algorithm However since MCQ NLO versions 2 0 and higher make use of the Les Houches interface first implemented in HERWIG 6 5 the version must be 6 500 or higher On most systems users will need to delete the dummy subroutines UPEVNT UPINIT PDFSET a
20. he corresponding positive input process codes IPROC 1350 1399 1499 1600 ID 1705 1706 2600 ID or 2700 ID to generate a standard HERWIG run for comparison purposes Then the input event file will not be read instead parton configurations will be generated by HERWIG according to the LO matrix elements 4For vector boson pair production for historical reasons the different process codes 2800 2825 must be used 18 EVPREFIX events NEVENTS pdftype PART1 PART2 beammom beammom IPROC PDFGROUP PDFSET PDFGROUP PDFSET LAMBDAHERW WMASS WMASS ZMASS event file number of events 0 gt Herwig PDFs 1 otherwise hadron types beam momenta 2850 60 70 80 WW ZZ ZW ZW PDF group 1 PDF id number 1 PDF group 2 PDF id number 2 Lambda_5 lt 0 for default M_W M_W MZ UMASS DMASS SMASS CMASS BMASS GMASS quark and gluon masses Table 4 Sample input file for the MC code for vector boson pair production resulting from setting HERPDF EXTPDF which implies pdftype 1 Setting HERPDF DEFAULT results in an analogous file with pdftype 0 and without the lines concerning PDFGROUP and PDFSET EVPREFIX must be understood with SCRTCH in front see sect 4 The negative sign of IPROC tells HERWIG to use Les Houches interface routines Variable esctype pdftype beammom Default value 0 0 1 HERPDF DEFAULT EXTPDF EMC 2 Table 5 Default values for script generated variables in MCinpu
21. id assignment HERWIGVER This variable must to be set equal to the name of the object file corresponding to the version of HERWIG linked to the package for example HERWIGVER herwig65 o is a valid assignment PDFPATH The physical address of the directory where the PDF grids are stored Acknowledgement Many thanks to Paolo Nason for contributions to the heavy quark code and valuable discussions on all aspects of the MCQ NLO project We also thank V Drollinger and B Quayle for testing a preliminary version of the WW code with spin correlations 13 Appendices A Version changes A 1 From MC NLO version 1 0 to version 2 0 In this appendix we list the changes that occurred in the package from version 1 0 to version 2 0 e The Les Houches generic user process interface has been adopted e As a result the convention for process codes has been changed MC NLO process codes IPROC are negative e The code mcatnlo_hwhvvj f which was specific to vector boson pair produc tion in version 1 0 has been replaced by mcatnlo_hwlhin f which reads the event file according to the Les Houches prescription and works for all the production processes implemented e The Makefile need not be edited since the variables HERWIGVER and HWUTI have been moved to MCatNLO inputs where they must be set by the user e A code mcatnlo_hbook f has been added to the list of utility codes It contains a simplified version written by M Mangano of HBOOK
22. n are not included in WtW production 10 FPREFIXNLOinput the input file for the NLO executable created according to the set of input parameters defined in MCatNLO inputs where the user also sets the string FPREFIX See table 2 FPREFIXNLO 1log the log file relevant to the NLO run FPREFIXxxx data xxx can assume several different values These are the data files created by the NLO code They can be removed only if no further event generation step is foreseen with the current choice of parameters FPREFIXMCinput analogous to FPREFIXNLOinput but for the MC executable See table 4 FPREFIXMC 1log analogous to FPREFIXNLO 1log but for the MC run EVPREFIX events the event file where EVPREFIX is the string set by the user in MCatNLO inputs EVPREFIXxxx events xxx can assume several different values These files are temporary event files which are used by the NLO code and eventually removed by the shell scripts They MUST NOT be removed by the user during the run the program will crash or give meaningless results By default all the files produced by the MCQNLO are written in the running directory However if the variable SCRTCH to be set in MCatNLO inputs is not blank the data and event files will be written in the directory whose address is stored in SCRTCH such a directory is not created by the scripts and must already exist at run time 4 Script variables In the following we list all the variables a
23. nd STRUCTM from the standard HERWIG package to permit linkage of the corresponding routines from the MCQ NLO package As a general rule the user is strongly advised to use the most recent version of HERWIG currently 6 507 with versions lower than 6 504 problems can be found in attempting to specify the decay modes of single vector bosons through the variable MODBOS Also crashes in the shower phase have been reported when using HERWIG 6 505 and we therefore recommend not to use that version 1 3 Working environment We have written a number of shell scripts and a Makefile all listed under Shell utilities above which will simplify the use of the package considerably In order to use them the computing system must support bash shell and gmake Should they be unavailable on the user s computing system the compilation and running of MCQ NLO requires more detailed instructions in this case we refer the reader to app B This appendix will serve also as a reference for a more advanced use of the package For Macs running under OSX v10 or higher make can be used instead of gmake 1 4 Source and running directories We assume that all the files of the package sit in the same directory which we call the source directory When creating the executable our shell scripts determine the type of operating system and create a subdirectory of the source directory which we call the running directory whose name is Alpha Sun Linux or Dar
24. ngle vector boson production processes IPROC 1396 1397 1497 1498 can be used in which case the vector boson decay products are distributed according to phase space There are a number of other differences between the lepton pair and single vector boson processes The latter do not feature the y Z interference terms Also their cross sections are fully inclusive in the final state fermions resulting from y Z or W The user can still select a definite decay mode using the variable MODBOS see sect 3 2 but the relevant branching ratio will not be included automatically by MC NLO In the case of 7 production the branching ratios are C q 20 3 qi being the electric charge in units of the positron charge of the fermion i selected through MODBOS and C 1 for leptons and 3 for quarks Notice that 20 3 X Cig the sum including all leptons and quarks except the top Thus the total rate predicted by MC NLO in the case of lepton pair production can also be recovered by multiplying the corresponding single vector boson total rate by the relevant branching ratio In the case of vector boson pair production the process codes are the negative of those adopted in MC NLO 1 0 for which the Les Houches interface was not yet available rather than those of standard HERWIG Furthermore in the case of Higgs production in association with a W or Z as well as vector boson pair production the value of IPROC alone is not sufficient
25. ode of operation 3 1 2 Package files 4 1 3 Working environment 6 1 4 Source and running directories 7 2 Prior to running T 2 1 Parton densities 7 3 Running 8 3 1 Event file 9 3 2 Decays 10 3 3 Results 10 4 Script variables 11 A Version changes 14 A 1 From MC NLO version 1 0 to version 2 0 14 A 2 From MC NLO version 2 0 to version 2 1 14 A 3 From MC NLO version 2 1 to version 2 2 15 A 4 From MC NLO version 2 2 to version 2 3 15 A 5 From MC NLO version 2 3 to version 3 1 15 B Running the package without the shell scripts 16 B 1 Creating the executables 16 B 2 The input files 17 1 Generalities In this documentation file we briefly describe how to run the MCQ NLO package im plemented according to the formalism introduced in ref 1 When using MC NLO please cite refs 1 2 The production processes now available are listed in table 1 The process codes IPROC will be explained below Hj represent hadrons in practice por p The treatment of undecayed vector boson pair production within MC NLO has been described in ref 1 that of heavy quark pair production in ref 2 The NLO matrix elements for these processes have been taken from refs 3 4 5 6 The information given in refs 1 2 allows the implementation in MC NLO of any production process provided that the formalism of refs 7 8 is used for the compu tation of cross sections to NLO accuracy The matrix elements for Standard Model Higgs single vector bo
26. ppearing in MCatNLO inputs these can be changed by the user to suit his her needs This must be done by editing MCatNLO inputs For fuller details see the comments in MCatNLO inputs ECM The CM energy of the colliding particles FREN The ratio between the renormalization scale and a reference mass scale FFACT As FREN for the factorization scale HVQMASS The mass in GeV of the top quark except when IPROC 1 1705 when it is the mass of the bottom quark In this case HVQMASS must coincide with BMASS xMASS The mass in GeV of the particle x with x HGG W Z U D S C B G xWIDTH The physical Breit Wigner width in GeV of the particle x with x HGG W Z IBORNHGG Valid entries are 1 and 2 If set to 1 the exact top mass dependence is retained at the Born level in Higgs production If set to 2 the m oo limit is used 11 xGAMMAX xMASSINF xMASSSUP AEMRUN IPROC IVCODE ILxCODE PARTn PDFGROUP PDFSET LAMBDAFIVE LAMBDAHERW SCHEMEOFPDF FPREFIX EVPREF IX EXEPREF IX NEVENTS If xGAMMAX gt 0 controls the width of the mass range for Higgs x H and vector bosons x V1 V2 the range is MASS GAMMAX x WIDTH Lower limit of the Higgs x H or vector boson x V1 V2 mass range used only when xGAMMAX lt 0 Upper limit of the Higgs x H or vector boson x V1 V2 mass range used only when xGAMMAX lt 0 Set it to YES to use running Qem in lepton pair and single vector boson produ
27. runMCatNLO 14 logical links to these files will be created in the running directory in version 2 0 this operation had to be performed by the user manually e Minor bugs corrected in mcatnlo_hbook f and sample analysis routines A 3 From MC NLO version 2 1 to version 2 2 In this appendix we list the changes that occurred in the package from version 2 1 to version 2 2 e Single vector boson production has been added which implies new process specific files mcatnlo_sbmain f mcatnlo_sbxsec f svbcblks h mcatnlo_hwansvb f and a modification to mcatnlo_hwlhin f e The script variables WWIDTH and ZWIDTH have been added to MCatNLO inputs These denote the physical widths of the W and Z bosons used to generate the mass distributions of the vector bosons according to the Breit Wigner function in the case of single vector boson production vector boson pair production is still implemented only in the zero width approximation A 4 From MC NLO version 2 2 to version 2 3 In this appendix we list the changes that occurred in the package from version 2 2 to version 2 3 e Lepton pair production has been added which implies new process specific files mcatnlo_llmain f mcatnlo_llxsec f 1lpcblks h mcatnlo_hwanllp f and mod ifications to mcatnlo_hwlhin f and mcatnlo_hwdriver f e The script variable AEMRUN has been added since the computation of single vector boson and lepton pair cross sections is performed in the MS scheme the
28. s Consistently with what happens in standard HERWIG by subtracting 10000 from IPROC one generates the same processes as in table 1 but eliminates the underlying event Higgs decays are controlled in the same way as in HERWIG that is by adding ID to the process code The conventions for ID are the same as in HERWIG namely ID 1 6 for u tt 7 9 for ete T r 10 11 for WWT ZZ and 12 for yy Furthermore ID 0 gives quarks of all flavours and ID 99 gives all decays Process codes IPROC 1360 IL and 1370 IL do not have an analogue in HERWIG they are the same as 1350 IL except for the fact that only a Z or a y respectively is exchanged The value of IL determines the lepton identities and the same convention as in HERWIG is adopted IL 1 6 for li e Ve H Vu T Vr respectively At variance with HERWIG IL cannot be set equal to zero Process codes IPROC 1460 IL and 1470 IL are the analogue of HERWIG 1450 IL in HERWIG either Wt or W can be produced whereas MC NLO treats the two vector bosons separately For these processes as in HERWIG IL 1 2 3 for Jy e u T but again the choice IL 0 is not allowed The same effect can be achieved by setting the HERWIG parameter PRSOF 0 The lepton pair processes IPROC 1350 IL 1470 IL include spin corre lations when generating the angular distributions of the produced leptons However if spin correlations are not an issue the si
29. son lepton pair and associated Higgs production have been taken from refs 9 10 ref 11 ref 12 and ref 13 respectively oat v mmo in eX mson v amu iso mm gt 0 nt xX _ cm v amo WF im rao mm WwW nat a e mrX a x amx aw a SX Pe dx amw CS om l amoa Pm amea is pms x amon 200 m li 7 x hho wWwetX O cn i SW E 2600 rf e aS eWe x eom a i i SW oX sro o 7 x mh WZ EX omo fa v mm gt WO it x 2850 7 7 x Mih gt Wtw x P2880 i i i gt Wt AW yp X P2860 r 7 x amo OO so x a SW P2880 iW E Table 1 Processes implemented in MC NLO 3 1 H denotes the Standard Model Higgs boson and the value of ID controls its decay as described in the HERWIG manual and below The values of IV IL ILj and IL control the identities of vector bosons and leptons as described below IPROC 10000 generates the same processes as IPROC but eliminates the underlying event A void entry indicates that the corresponding variable is unused The Spin column indicates whether spin correlations in vector boson or top decays are included v neglected x or absent void entry Spin correlations in Higgs decays are included by HERWIG e g in H WWT ItvI D This documentation refers to MC NLO version 3 1 previous versions 1 0 2 0 2 2 and 2 3 are described in refs 1
30. t References 1 S Frixione and B R Webber Matching NLO QCD computations and parton shower simulations JHEP 0206 2002 029 hep ph 0204244 2 S Frixione P Nason and B R Webber Matching NLO QCD and parton showers in heavy flavour production JHEP 0308 2003 007 arXiv hep ph 0305252 3 B Mele P Nason and G Ridolfi QCD Radiative Corrections To Z Boson Pair Production In Hadronic Collisions Nucl Phys B 357 1991 409 4 S Frixione P Nason and G Ridolfi Strong corrections to W Z production at hadron colliders Nucl Phys B 383 1992 3 19 10 11 12 13 14 15 16 17 18 19 S Frixione A Next to leading order calculation of the cross section for the produc tion of W W pairs in hadronic collisions Nucl Phys B 410 1993 280 M L Mangano P Nason and G Ridolfi Heavy quark correlations in hadron col lisions at next to leading order Nucl Phys B 373 1992 295 S Frixione Z Kunszt and A Signer Three jet cross sections to next to leading order Nucl Phys B 467 1996 399 arXiv hep ph 9512328 S Frixione A general approach to jet cross sections in QCD Nucl Phys B 507 1997 295 arXiv hep ph 9706545 S Dawson Radiative Corrections To Higgs Boson Production Nucl Phys B 359 1991 283 A Djouadi M Spira and P M Zerwas Production of Higgs bosons in proton colliders QCD
31. which contain the user s analysis routines We point out that since version 2 0 the Makefile need not be edited any longer since the corresponding operations are now performed by setting script variables see sect 4 2 1 Parton densities Since the knowledge of the parton densities PDF is necessary in order to get the physical cross section a PDF library must be linked The possibility exists to link the CERNLIB PDF library PDFLIB however we also provide a self contained PDF library with this package which is faster than PDFLIB and contains PDF sets released after the last and final PDFLIB version 8 04 A complete list of the PDFs available in our PDF library can be downloaded from the MC NLO web page The user may link either PDF library all that is necessary is to set the variable PDFLIBRARY in the file MCatNLO inputs equal to THISLIB if one wants to link to our PDF library and equal to PDFLIB if one wants to link to PDFLIB Our PDF library collects the original codes written by the authors of the PDF fits as such for most of the densities it needs to read the files which contain the grids that initialize the PDFs These files which can be also downloaded from the MCQNLO web page must either be copied into the running directory or defined in the running directory as logical links to the physical files by using 1n sn We stress that if the user runs MC NLO with the shell scripts the logical links will be created automati
32. win depending on the operating system If the operating system is not known by our scripts the name of the working directory is Run The running directory contains all the object files and executable files and in general all the files produced by the MCQNLO while running It must also contain the relevant grid files see sect 2 1 or links to them if the library of parton densities provided with the MCQ NLO package is used 2 Prior to running Before running the code the user needs to edit the following files mcatnlo_hwangzz f mcatnlo_hwdriver f mcatnlo_hwlhin f We do not assume that the user will adopt the latest release of HERWIG although as explained above it must be version 6 500 or higher For this reason the files mcatnlo_hwdriver f and mcatnlo_hwlhin f must be edited in order to modify the INCLUDE HERWIGXX INC command to correspond to the version of HERWIG the user is going to adopt mcatnlo_hwdriver f contains a set of read state ments which are necessary for the MC to get the input parameters see sect 3 for the input procedure these read statements must not be modified or eliminated Also mcatnlo_hwdriver f calls the HERWIG routines which perform showering hadronization decays see sect 3 2 for more details on this issue and so forth the user can freely modify this part as customary in MC runs Finally the sample codes mcatnlo_hwangezz f contain analysis related routines these files must be replaced by files
Download Pdf Manuals
Related Search
Related Contents
Innerstaatliche Bauartzulassung - baer Avaya Business Communications Manager Release 5.0 Configuration Guide FT1A Controller Service Manual - Amalgamated Pumping Supplies FC-138/418 Peerless MOD-APC mounting kit The Anti-Tip Bracket - Pdfstream.manualsonline.com 生体の力学実験 はりの曲げたわみと応力の計測 119AS45_FROAG 230_ver 6.indd Copyright © All rights reserved.