Back to home page

EIC code displayed by LXR

 
 

    


Warning, /detector_benchmarks/benchmarks/nhcal_dimuon_photoproduction/scripts/ep_DiMuon.cmnd is written in an unsupported language. File is not indexed.

0001 #==============================================================================
0002 # STAR Heavy Flavor Tune 1.1
0003 #
0004 # PYTHIA Version 8.1.08
0005 # Date:  November 7, 2008       
0006 # Last updated by: Thomas Ullrich
0007 #
0008 # This file contains commands to be read in for a Pythia8 run. 
0009 # Lines not beginning with a letter or digit are comments.
0010 # Names are case-insensitive  -  but spellings-sensitive!
0011 #
0012 # Important notes:
0013 # ================
0014 # For hard 2->2 processes, such as most heavy flavor producing one,
0015 # there's  divergence in PYTHIA that needs to be suppressed. Use
0016 # these two lines in the code before you initialize PYTHIA to avoid
0017 # any trouble: 
0018 # UserHooks *oniumUserHook = new SuppressSmallPT();
0019 # pythia.setUserHooksPtr(oniumUserHook);
0020 #
0021 #==============================================================================
0022 
0023 #------------------------------------------------------------------------------
0024 # Parameters that need to be set by whoever runs it.
0025 # Note that they have no meaning unless restored and used
0026 # in the user provided code (main program).
0027 # This is not part of the star_hf tune (just convenient)
0028 # Documentation: <pyhiadir>/htmldoc/MainProgramSettings.html
0029 #------------------------------------------------------------------------------
0030 Main:numberOfEvents = 1000       ! number of events to generate
0031 Next:numberShowEvent = 20             ! show event record
0032 Next:numberCount = 0              ! number of events to print
0033 Next:numberShowInfo = 100             ! show how far along run is this many times
0034 Main:timesAllowErrors = 100        ! abort run after this many flawed events
0035 Init:showChangedSettings = on      ! print changed flags/modes/parameters
0036 Init:showAllSettings = on          ! print all flags/modes/parameters
0037 #Init:showAllStatistics = on        ! print statistics at the end
0038 
0039 #------------------------------------------------------------------------------
0040 # Colliding beams and center-of-mass energy
0041 # Documentation: <pyhiadir>/htmldoc/BeamParameters.html
0042 #------------------------------------------------------------------------------
0043 Beams:frameType = 2
0044 Beams:idA = 2212                   ! proton
0045 Beams:idB = 11                 ! e-
0046 Beams:eA = 275.                  ! RHIC nominal (GeV)
0047 Beams:eB = 18.                  ! RHIC nominal (GeV)
0048 
0049 #------------------------------------------------------------------------------
0050 # Process Selection
0051 #------------------------------------------------------------------------------
0052 
0053 
0054 # Enable equivalent photon approximation (EPA) for both beams
0055 Photon:Q2Max = 1.0                ! Upper Q^2 limit for EPA photons (in GeV^2)
0056 Photon:ProcessType = 4    ! 4 = direct-direct photons
0057 #Photon:EPA = on
0058 
0059 PDF:beamA2gamma = on    ! EPA photon flux from beam A
0060 PDF:beamB2gamma = on    ! EPA photon flux from beam B
0061 
0062 # Enable gamma-gamma -> mu+ mu−
0063 PhotonCollision:gmgm2mumu = on
0064 
0065 # Optional: Turn off other QED or QCD backgrounds if you want exclusivity
0066 PartonLevel:ISR = off
0067 PartonLevel:FSR = off
0068 HadronLevel:all = off
0069 
0070 
0071 #------------------------------------------------------------------------------
0072 # K factor
0073 # Multiply almost all cross sections by this common fix factor.
0074 # This is usually no very useful. The data can be shifted up and down
0075 # later anyhow as we please. 
0076 # Documentation: <pyhiadir>/htmldoc/CouplingsAndScales.html
0077 #------------------------------------------------------------------------------
0078 # SigmaProcess:Kfactor = 3
0079 
0080 #------------------------------------------------------------------------------
0081 # Scales (Ramona's suggestions)
0082 # This sets the scale to settings typically for hard probes:
0083 # mu_F = mu_R = 2*mT
0084 # Documentation: <pyhiadir>/htmldoc/CouplingsAndScales.html
0085 #------------------------------------------------------------------------------
0086 # SigmaProcess:renormScale2 = 3
0087 # SigmaProcess:factorScale2 = 3
0088 # SigmaProcess:renormMultFac = 2   ! 2mT
0089 # SigmaProcess:factorMultFac = 2   ! 2mT
0090 
0091 #------------------------------------------------------------------------------
0092 # To limit particle production to a certain pthat range uncomment
0093 # these lines. Use only when you 100% know what you are doing.
0094 # It is extremely useful to split runs up in ptHat bins to generate
0095 # statistics evenly in pt. Book keeping is important then (cross-sections,
0096 # number of events) to compile the final complete spectra.
0097 # Documentation: <pyhiadir>/htmldoc/PhaseSpaceCuts.html
0098 #------------------------------------------------------------------------------
0099 # PhaseSpace:pTHatMin = 4               ! Minimal pT cut
0100 # PhaseSpace:pTHatMax = 2
0101 
0102 #------------------------------------------------------------------------------
0103 # Random Number
0104 # Initialize random generator according to time. Otherwise multiple jobs
0105 # will produce the same sequence (unless you pass a different seed every
0106 # time which is not practical).
0107 # Documentation: <pyhiadir>/htmldoc/RandomNumberSeed.html
0108 #------------------------------------------------------------------------------
0109 Random:setSeed = on
0110 Random:seed = 0
0111 
0112 #------------------------------------------------------------------------------
0113 # PDF Selection:
0114 # Note: you need LHAPDF to be installed. Pythia 8 only provides a 
0115 # minimal set to get started.
0116 # The choice of PDF here is greatly motivated by:
0117 # A.~Sherstnev and R.~S.~Thorne, arXiv:0807.2132 and arXiv:0711.2473v3
0118 # and W. Vogelsang (private communication)
0119 # These are PDFs especially made for LO Monte-Carlo generators such
0120 # as PYTHIA.
0121 # The state-of-the-art NLO PDF is cteq66.LHgrid which can be used
0122 # as an alternative (also from LHAPDF) but with the known problems
0123 # that arise when using a NLO PDF in an LO simulator.
0124 # Documentation: <pyhiadir>/htmldoc/PDFSelection.html
0125 #------------------------------------------------------------------------------
0126 #PDF:useLHAPDF = off
0127 #PDF:LHAPDFset = MRSTMCal.LHgrid
0128 #PDF:extrapolateLHAPDF = on
0129 #PDF:pSet = 8
0130 #PDF:beamA2gamma = on
0131 
0132 #------------------------------------------------------------------------------
0133 # Settings for the event generation process in the Pythia8 library
0134 # Effect/Relevance of MI, ISR and FSR need to be checked. For sure
0135 # the use more CPU and hence less events/s.
0136 # If HadronLevel:Hadronize = off we end up with the pure c, b spectra
0137 # (which might be useful at times)
0138 # Documentation: <pyhiadir>/htmldoc/MasterSwitches.html
0139 # Documentation: <pyhiadir>/htmldoc/MultipleInteractions.html
0140 #------------------------------------------------------------------------------
0141 #PartonLevel:MPI = on              ! multiple interactions
0142 #PartonLevel:ISR = on             ! initial-state radiation 
0143 #BeamRemnants:primordialKT = on   ! primordial kt
0144 #PartonLevel:FSR = on             ! final-state radiation
0145 #HadronLevel:Hadronize = off     ! no hadronization use
0146 
0147 #------------------------------------------------------------------------------
0148 # Relative production ratio vector/pseudoscalar for charm and bottom mesons
0149 # This was originally PARJ(13) where PARJ(13) = V/(PS+V) that is the 
0150 # vector meson  fraction of primary charm+bottom mesons. 
0151 # Andre David (CERN/NA60) made an exhaustive study and found that the
0152 # world data supports 0.6 while PYTHIA default was PARJ(13) = 3/4 = 0.75
0153 # from simple spin counting.
0154 # In PYTHIA8 we now use V/PS not V/(PS+V)
0155 # Documentation: <pyhiadir>/htmldoc/FlavourSelection.html
0156 #------------------------------------------------------------------------------
0157 #StringFlav:mesonCvector = 1.5    ! same as PARJ(13)=0.6
0158 #StringFlav:mesonBvector = 3      ! leave at 0.75 
0159 
0160 #------------------------------------------------------------------------------
0161 # Heavy quark masses.
0162 # Note that this should match with the ones used in the PDF.
0163 # The masses are listed in the header of the refering PDF file.
0164 # Documentation: <pyhiadir>/htmldoc/ParticleDataScheme.html
0165 # Documentation: <pyhiadir>/htmldoc/ParticleData.html
0166 #------------------------------------------------------------------------------
0167 #4:m0 = 1.43
0168 #5:m0 = 4.30
0169 
0170 #------------------------------------------------------------------------------
0171 # Particle Decay limits
0172 # When on, only particles with a decay within a volume limited by 
0173 # rho = sqrt(x^2 + y^2) < xyMax and |z| < zMax are decayed. 
0174 # The above xyMax, expressed in mm/c.
0175 #------------------------------------------------------------------------------
0176 # ParticleDecays:limitCylinder = on
0177 # ParticleDecays:xyMax = 600
0178 # ParticleDecays:zMax = 1000
0179 
0180 # EOF