Лого Сделано у нас
43

СуН на «Эльбрусе» :)

читать полностью

Читайте также...

Вступайте в наши группы и добавляйте нас в друзья :)

Подпишитесь на наш канал в Яндекс.Дзен и сделайте вашу ленту объективнее!
  • 2
    vx8400 . vx8400 .
    15.12.1513:54:59

    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. как здесь сохранить пробелы в начале строки?

    Отредактировано: vx8400 .~15:08 15.12.15
    • 1
      Нет аватара nesk
      15.12.1514:08:38

      спасибо     жду когда соберется pythia-6.4.28.f

      CPU TTY PID USERNAME PRI NI SIZE RES STATE TIME %WCPU %CPU COMMAND

      1 pts/0 18315 root 241 20 255M 184M run 67:41 100.68 100.51 f90com32

      пока еще компилируется, уже более 67 минут чистого времени на CPU

    • 1
      Нет аватара nesk
      15.12.1514:55:38

      там по семь по семь пробелов перед инструкциями и комментарии (С…) с первой позиции? Да?

      я расставлю.

      Отредактировано: nesk~15:56 15.12.15
    • 1
      Нет аватара nesk
      15.12.1514:56:59

      а можешь на почту мне кинуть?

Написать комментарий
Отмена
Для комментирования вам необходимо зарегистрироваться и войти на сайт,