hermes. ihep. su: 8001/pool/atlas/p6/main79.f

Все там открывается. Ловите код сюда:

PROGRAM MAIN79

C…Test program for new gcc4 compiler.

C…Top pair events are generated at the LHC, 14 TeV.

C…Final charged multiplicity is histogrammed.

C…The final total cross section should be 4.9E-07 mb,

C…the average multiplicity 243, and the rms 70.2,

C…all within statistical fluctuations.

C******************************************************************

C…All real arithmetic in double precision.

IMPLICIT DOUBLE PRECISION(A-H, O-Z)

C…Three Pythia functions return integers, so need declaring.

INTEGER PYK,PYCHGE,PYCOMP

C…EXTERNAL statement links PYDATA on most machines.

EXTERNAL PYDATA

C…Commonblocks.

C…The event record.

COMMON/PYJETS/N,NPAD,K(4000,5),P(4000,5),V(4000,5)

C…Selection of hard scattering subprocesses.

COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)

C…Number of events.

NEV=1000

C…Processes qqbar, gg -> ttbar.

MSEL=6

C…Initialize for the LHC.

CALL PYINIT('CMS','p','p',14000D0)

C…Histogram.

CALL PYBOOK(1,'Final charged multiplicity',100,-1D0,599D0)

C…Event generation loop.

DO 200 IEV=1,NEV

if (MOD(IEV,100) .eq. 0) write(*,*) 'evt: ',IEV

CALL PYEVNT

C…Remove all but charged particles and histogram multiplicity.

CALL PYEDIT(3)

CALL PYFILL(1,DBLE(N),1D0)

C…End of event generation loop.

200 CONTINUE

C…Cross section. Histogram.

CALL PYSTAT(1)

CALL PYHIST

END

P. S. как здесь сохранить пробелы в начале строки?