1 / 10

The HV1.0 Monte Carlo

The HV1.0 Monte Carlo. Steve Mrenna, Peter Skands, and MJS. Features and Important Facts. The new program can do more things  Final state can be in form of Les Houches accord event record All v-hadrons in final state; user must decay them

bruis
Download Presentation

The HV1.0 Monte Carlo

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. The HV1.0 Monte Carlo Steve Mrenna, Peter Skands, and MJS

  2. Features and Important Facts • The new program can do more things  • Final state can be in form of • Les Houches accord event record • All v-hadrons in final state; user must decay them • Stable v-hadron and unstable-v-hadron standard-model decay products NEW • HEPEVT hadron-level output • Two approaches to v-sector • Scaled-up QCD • User-defined theory with independent fragmentation NEW • More production processes • Z’ • Others in development (not tested yet)

  3. Features and Important Facts • User Input has changed and is still changing  • User defines the process via standard Pythia call NEW • User defines the theory through two upload files and standard Pythia calls NEW • Currently CALL PYUPDA (Pythia < 6.410) • Moving to CALL PYSLHA (Pythia = 6.410) • HV.input.par file in process of adjustment • Some “IMODE”s of operation no longer active, new ones added • New input file needed to control fragmentation, other settings; not present yet

  4. HV.input.par Note the running modes have changed! Use 2,4,5 only '! IMODE 1 lha + hepevt output (UE subtlety, see README) 2 lha output (vhadron level) 3 lha input -> pythia output (not tested) 4 lha output (parton level) 5 hepevt output (hadron level)' 3/15/07 !!!!only use modes 2,4,5 for now '! In modes 2,4 program shuts off fragmentation of proton remnants in initial state, showering/hadronization in final state '! Add 60 to IMODE for A0,H0 bosons to be standins for vpions in lha file' '! IHVTHY0 -- 1 is scaled-up vQCD, 9 is user-defined theory '! number of flavors inferred from input file where vquarks x00000f appear '! scale for scaled-up vQCD inferred from vQCDscale (code x000100) '! '! IFRAG: 0 default, 1 independent vfragmentation, '! 2 scaled-up qcd hadronization, 3 scaled-up qcd showering and fragmentation' '! Default is ifrag=3 for ihvthy0=1, ifrag=1 for ihvthy=9 '! '! The START LINE initiates read-in' '! The next line reads in:' '! IMODE,IDIAGN,NEV,IRANDM,ICLDR,IHVTHY0,IFRAG' 'START' 2,9,3,24181784,0,9,0 ! do not insert a comment line next myinput ! this line must be the input file prefix followed by a space ! myoutput ! this line must be the output file prefix followed by a space ! MSEL=0 ! user control ! MSUB(141)=1 ! Z' production ! MSTP(44)=3 ! no interference ! MSTP(47)=0 ! phase space Z' CKIN(1)=1650.D0 ! minimum partonic collision energy 'END' Temporary – is this needed anymore? Note masses and decays are no longer set here… But the process (Z’ production) is set here

  5. Input v-sector (PYUPDA version) pdg code name conjugate mass width decay info 6000021 gv 0 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 0 6000001 dv dvbar 0 0 1 0.06000 0.00000 0.00000 0.00000E+00 0 0 6000002 uv uvbar 0 0 1 0.03000 0.00000 0.00000 0.00000E+00 0 0 6000091 vcluster 0 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 0 6000092 vstring 0 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 0 6000093 vindep 0 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 0 6000094 CMvshwr 0 0 0 0.00000 0.00000 0.00000 0.00000E+00 0 0 6000100 Lambda_V 0 0 0 55.00000 0.00000 0.00000 0.00000E+00 0 0 6000111 piv0 0 0 0 13.49800 0.00001 0.00000 1.30680E+01 2 1 1 0 1.000000 5 -5 0 0 0 6000113 rhov0 0 0 0 76.85000 0.15100 0.00000 1.30680E-12 2 1 1 0 1.000000 6000211 -6000211 0 0 0 6000115 a_2v0 0 0 0 131.80000 0.10700 0.00000 1.84418E-12 2 1 1 0 0.414132 6000213 -6000211 0 0 0 1 0 0.414132 -6000213 6000211 0 0 0 1 0 0.171735 6000221 6000111 0 0 0 6000211 piv+ piv- 0 0 1 13.95700 0.00000 0.00000 0.00000E+00 0 0 6000213 rhov+ rhov- 0 0 1 76.69000 0.14900 0.00000 1.32434E-12 2 1 1 0 1.000000 6000211 6000111 0 0 0 6000215 a_2v+ a_2v- 0 0 1 131.80000 0.10700 0.00000 1.84418E-12 2 1 1 0 1.000000 6000221 6000211 0 0 0 6000221 etav 0 0 0 54.74500 0.10000 0.00000 1.97327E-12 2 1 1 0 1.000000 6000111 6000111 6000111 0 0 Decays

  6. Input Z’ Br’s (PYUPDA version) 32 Z'0 0 0 0 3200.00000 14.54029 145.40294 0.00000E+00 1 1 0 32 0.145835 1 -1 0 0 0 0 32 0.113276 2 -2 0 0 0 0 32 0.145835 3 -3 0 0 0 0 32 0.113271 4 -4 0 0 0 0 32 0.145781 5 -5 0 0 0 0 32 0.029002 6 -6 0 0 0 0 32 0.020000 7 -7 0 0 0 0 32 0.000000 8 -8 0 0 0 0 0 0.032025 11 -11 0 0 0 0 0 0.063642 12 -12 0 0 0 0 0 0.032025 13 -13 0 0 0 0 0 0.063642 14 -14 0 0 0 1 0 0.032022 6000001 -6000001 0 0 0 1 0 0.063642 6000002 -6000002 0 0 0 0 0 0.000000 17 -17 0 0 0 0 0 0.000000 18 -18 0 0 0 0 0 0.000000 24 -24 0 0 0 0 0 0.000000 37 -37 0 0 0 0 0 0.000000 23 22 0 0 0 0 0 0.000000 23 25 0 0 0 0 0 0.000000 25 36 0 0 0 0 0 0.000000 35 36 0 0 0 Note all decays are disabled except those to v-quarks

  7. Output MODE 2: V-hadrons <event> 17 661 1. 25.5996552 0.00781999994 0.115999997 -1 -1 0 0 0 501 -0.17 -0.96 867.46 867.46 0.00 0. 9. 1 -1 0 0 501 0 -2.41 -2.29 -2903.10 2903.10 0.00 0. 9. 32 2 1 2 0 0 -2.58 -3.25 -2035.64 3770.56 3173.85 0.0000 9. 6000002 2 3 0 0 0 -1263.33 939.33 -1256.00 2013.92 0.03 0.0000 9. -6000002 2 3 0 0 0 1260.75 -942.58 -779.64 1756.64 0.03 0.0000 9. 6000092 2 4 0 0 0 -2.58 -3.25 -2035.64 3770.56 3173.85 0.0000 9. 6000111 1 6 0 0 0 81.47 -4.92 41.97 95.28 25.60 0.0000 9. -6000211 1 6 0 0 0 19.69 -40.24 -29.20 59.66 26.47 0.0000 9. -6000211 1 6 0 0 0 -494.84 353.82 -432.07 746.62 26.47 0.0000 9. 6000211 1 6 0 0 0 -250.89 244.01 -309.06 467.66 26.47 0.0000 9. 6000111 1 6 0 0 0 -81.67 48.28 -78.96 126.06 25.60 0.0000 9. 6000211 1 6 0 0 0 -281.16 187.12 -365.39 498.27 26.47 0.0000 9. 6000111 1 6 0 0 0 -100.93 55.15 -75.22 139.79 25.60 0.0000 9. -6000211 1 6 0 0 0 157.95 -55.25 -37.65 173.55 26.47 0.0000 9. 6000211 1 6 0 0 0 82.56 -69.14 -129.98 170.85 26.47 0.0000 9. -6000211 1 6 0 0 0 167.64 -199.98 -157.34 305.86 26.47 0.0000 9. 6000211 1 6 0 0 0 697.61 -522.10 -462.75 986.96 26.47 0.0000 9. </event> 6000111 is pizero, 6000211 is piplus; no “neutrinos” anymore

  8. Output MODE 4: <event> 23 661 1. 25.5996552 0.00781999994 0.115999997 -1 -1 0 0 0 501 -0.17 -0.96 867.46 867.46 0.00 0. 9. 1 -1 0 0 501 0 -2.41 -2.29 -2903.10 2903.10 0.00 0. 9. 32 2 1 2 0 0 -2.58 -3.25 -2035.64 3770.56 3173.85 0.0000 9. 6000002 2 3 0 0 0 -1263.33 939.33 -1256.00 2013.92 0.03 0.0000 9. -6000002 2 3 0 0 0 1260.75 -942.58 -779.64 1756.64 0.03 0.0000 9. 6000092 2 4 0 0 0 -2.58 -3.25 -2035.64 3770.56 3173.85 0.0000 9. 6000111 2 6 0 0 0 81.47 -4.92 41.97 95.28 25.60 4.7274 9. -6000211 1 6 0 0 0 19.69 -40.24 -29.20 59.66 26.47 0.0000 9. -6000211 1 6 0 0 0 -494.84 353.82 -432.07 746.62 26.47 0.0000 9. 6000211 1 6 0 0 0 -250.89 244.01 -309.06 467.66 26.47 0.0000 9. 6000111 2 6 0 0 0 -81.67 48.28 -78.96 126.06 25.60 4.3217 9. 6000211 1 6 0 0 0 -281.16 187.12 -365.39 498.27 26.47 0.0000 9. 6000111 2 6 0 0 0 -100.93 55.15 -75.22 139.79 25.60 12.1460 9. -6000211 1 6 0 0 0 157.95 -55.25 -37.65 173.55 26.47 0.0000 9. 6000211 1 6 0 0 0 82.56 -69.14 -129.98 170.85 26.47 0.0000 9. -6000211 1 6 0 0 0 167.64 -199.98 -157.34 305.86 26.47 0.0000 9. 6000211 1 6 0 0 0 697.61 -522.10 -462.75 986.96 26.47 0.0000 9. 5 1 7 0 511 0 2.18 -2.26 0.86 5.80 4.80 0.0000 9. -5 1 7 0 0 511 79.28 -2.66 41.10 89.47 4.80 0.0000 9. 5 1 11 0 515 0 -54.30 41.01 -64.59 93.95 4.80 0.0000 9. -5 1 11 0 0 515 -27.37 7.27 -14.37 32.11 4.80 0.0000 9. 5 1 13 0 517 0 -7.16 3.01 -0.88 9.17 4.80 0.0000 9. -5 1 13 0 0 517 -93.76 52.14 -74.34 130.62 4.80 0.0000 9. </event> Can you use this? How much information about the new v-particles do you need?

  9. Output Mode 5: Hadrons • Same as before, except now the new 6000xxx particles appear in the HEPEVT record • This will cause problems for the event display program • However we can either • Keep a mode in the current HV program so that 6000xxx particles are replaced with H0, A0 (this is currently the case, add 60 to IMODE • Write a script to convert all 6000xxx particles to something the programs recognize • Convert the programs to recognize 6000xxx particles

  10. Summary • Many advances but not stable yet • Currently must use a non-standard Pythia version because of problems and bugs • Usable I think! But will need to adjust input details as time goes on • Need to decide on output modes; is current arrangement for output mode 4 ideal, or must it be altered? • How should we move away from using H0/A0 in our code? • How are Pythia/GEANT/Athena to be taught about the new particles and their properties, which are needed even in mode 4? • New theories can also be studied now, but that’s for the future • PYUPDA will move to PYSLHA when Atlas moves to Pythia 6.411

More Related