стань автором. присоединяйся к сообществу!
Добавить новость можно всем, без премодерации, только регистрация
Лого Сделано у нас
shigorin 10 декабря 2015, 20:28 37 Михаил Шигорин

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

читать полностью
Источник: sdelanounas.ru
МЦСТ, Эльбрус, Эльбрус-4С, Эльбрус-401, фотофакт
все комментарии
  • 2
    Нет аватара guest
    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
    #723556 ↑
    • 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

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

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

      я расставлю.

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

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

      #723600 ↑
Для комментирования вам необходимо зарегистрироваться и войти на сайт,