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