All direct photon production in bins of PhaseSpace:mHatMin, in the range pT=50-25,000 GeV. Cross sections depend on pT and are not given. 2000 events per file, 50 files per pT bin. Cross sections are included in the log files.

How to decode name

Name: tev100_pythia8_ptPTMIN_PTMAX_NUMBER.promc
where PTMIN is min pT for PhaseSpace:ptHatMin,
PTMAX - PhaseSpace:pTHatMax*2,
NUMBER corresponds to the cluster production node. When you run, use "PTMIN-PTMAX" ranges. The actual cut on PhaseSpace:pTHatMax=2*PTMAX to reduce fluctuations in the tails.

How to use

To get a sample with a given pT, use "glob" regular expressions (in Jython, Python etc). Example for the code: pT>100 corresponds to "glob.glob(*pt100_*)"

Slimmed as

Particle records are slimmed (all stable with pT>0.3 GeV) and (PID=5 || PID=6) or PID>22 && PID<38) or PID>10 && PID<17).

Input card

. input =Random:seed = 0
.. input =Beams:idA = 2212
.. input =Beams:idB = 2212
.. input =Beams:eCM = 100000.
.. input =HardQCD:all = off
PromptPhoton:all = on
.. input =Tune:pp = 5
.. input =PDF:pSet = LHAPDF6:MSTW2008lo68cl.LHgrid
.. input =MultipartonInteractions:bProfile = 4
.. input =MultipartonInteractions:a1 = 0.03
.. input =MultipartonInteractions:pT0Ref = 1.90
.. input =MultipartonInteractions:ecmPow = 0.30
.. input =SpaceShower:rapidityOrder=0
.. input =ParticleDecays:limitTau0 = on
.. input =ParticleDecays:tau0Max = 10
.. input =EventsNumber=10
.. input =PhaseSpace:pTHatMin = PTMIN
.. input =PhaseSpace:mHatMax = PTMIN