GMACRO-generated file in FZEX zebra exchange format are available with
interactions generated by the NEUGEN (MINOS) neutrino generator.
There are 16 files, each containing 10,000 events, for a total of
100,000 numu/numubar events and 60,000 nue/nuebar events. They can be
found on cithe504.cithep.caltech.edu (use the macro account with the
usual password; cithe504 runs ssh so you can use scp to avoid
transmitting an unencrypted password). The files are
numu/numubar:
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.0.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.1.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.2.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.3.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.4.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.5.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.6.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.7.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.8.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.mu.10k.9.zeb
nue/nuebar:
/macro7/nolty/gmacro-output/gmacro.nugn.e.10k.0.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.e.10k.1.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.e.10k.2.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.e.10k.3.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.e.10k.4.zeb
/macro7/nolty/gmacro-output/gmacro.nugn.e.10k.5.zeb
My generation volume is slightly bigger than MACRO, so about 2% of
these events interact in the rock around MACRO. In that enlarged
volume, the calculated event rates are
numu/numubar: 1259.611328 ev/yr
neu/nuebar: 650.1838989 ev/yr
so if you use all 100,000 mu events and 51,618 e events, you are
simulating 79.39 years.
If you want to decode and use the Monte Carlo truth information in the
SIMU banks, contact me -- I use a different format than others.
Here is some of the log file from the generation for numu/numubar
events:
###### MACRO GEANT3 SIMULATION (VERSION 7.21) June 1995 ######
###### FOR COMMENTS PLEASE SEND MAIL TO lamanna@roma1.infn.it ######
1 USER'S DIRECTIVES TO RUN THIS JOB
----------------------------------
***** DATA CARD CONTENT LIST
***** DATA CARD CONTENT TRIG 10000 selects number of events to be analysed
***** DATA CARD CONTENT RNDM 2125420028 31992246637
***** DATA CARD CONTENT NUSE 19780501
***** DATA CARD CONTENT *FOUT '/macro7/nolty/gmacro-output/gmacro.nugn.10k.1.zeb'
***** DATA CARD CONTENT KEYS
***** DATA CARD CONTENT RUNG 8000
***** DATA CARD CONTENT DATE 1650893216 9985603
***** DATA CARD CONTENT TRAK 1 to track the particle
***** DATA CARD CONTENT HITS 1 to evaluate the hits
***** DATA CARD CONTENT DIGI 1 to evaluate the digitisings
***** DATA CARD CONTENT RAWD 1 to evaluate raw data
***** DATA CARD CONTENT NSMOD 6
***** DATA CARD CONTENT NHMOD 2
***** DATA CARD CONTENT BACK 'NONE'
***** DATA CARD CONTENT MULS 1
***** DATA CARD CONTENT PFIS 0
***** DATA CARD CONTENT LOSS 1
***** DATA CARD CONTENT ANNI 1
***** DATA CARD CONTENT BREM 1
***** DATA CARD CONTENT COMP 1
***** DATA CARD CONTENT DRAY 1
***** DATA CARD CONTENT PAIR 1
***** DATA CARD CONTENT PHOT 0
***** DATA CARD CONTENT HADR 1
***** DATA CARD CONTENT MUNU 1
***** DATA CARD CONTENT RAYL 0
***** DATA CARD CONTENT DCAY 1
***** DATA CARD CONTENT C cutgam cutele cutneu cuthad cutmuo bcute bcutm dcute dcutm ppcutm
***** DATA CARD CONTENT CUTS .001 .001 .001 .001 .001 .001 .001 .001 .001 .001
***** DATA CARD CONTENT C
***** DATA CARD CONTENT C SET DETECTORS
***** DATA CARD CONTENT SETS 1='CPLS' 'LPLS' 'FPLS' 'CPLT' 'LPLT' 'FPLT'
***** DATA CARD CONTENT C The following sets must be used for the attico
***** DATA CARD CONTENT SETS 7='HTPT' 'TPT1' 'TPT2'
***** DATA CARD CONTENT SETS 10='HTPS' 'TPS1' 'TPS2'
***** DATA CARD CONTENT SETS 13='SQTP'
***** DATA CARD CONTENT C WNOI 1
***** DATA CARD CONTENT C SNOI 1
***** DATA CARD CONTENT TROT 1000
***** DATA CARD CONTENT TRST 5.
***** DATA CARD CONTENT
***** DATA CARD CONTENT C internal neutrino event generator
***** DATA CARD CONTENT KINE 1=7
***** DATA CARD CONTENT UPDN 0
***** DATA CARD CONTENT NUTY 2 2
***** DATA CARD CONTENT ECUT 0.100 999999.
***** DATA CARD CONTENT FLMO 6 4
***** DATA CARD CONTENT C-SOLP 0.5
***** DATA CARD CONTENT SOLP 3
***** DATA CARD CONTENT NBOX 1=8000.0 2=1400.0 3=1050.0
***** DATA CARD CONTENT C-NBOX 1=9000.0 2=1500.0 3=1300.0
***** DATA CARD CONTENT C-PTRK 1=5 2=6 3=13 4=14
***** DATA CARD CONTENT
***** DATA CARD CONTENT SAVE 'RAWZ'
***** DATA CARD CONTENT TIME 0. 0. 0
***** DATA CARD CONTENT
***** DATA CARD CONTENT ERPP 1=1 2=10 3=5 4=10 5=5
***** DATA CARD CONTENT END
FFREAD KEYWORDS C LIST READ END STOP NOLI WRIT HOLL KEYS *
USER KEYWORDS ABAN ANNI AUTO BACK BREM CKOV COMP CUTS DART DATE DCAY DEBU DIGI DRAY ECUT ELAT ERAN
ERPP ESTR EWIR FDIS FINP FLMO FOUT FRWV GEOM GET HADR HITS HSTA KINE LABS LOSS MTRG
MULS MUNU NBOX NHMO NSMO NUSE NUTY OPTI PAIR PCUT PFIS PHOT PLOT PNUM PRIN PTRK RAWD
RAYL RGET RNDM RSAV RUNG SAVE SCAL SCAN SCAP SETS SNOI SOLP SORD STAT STRA STRF SWIT
FZFILE. LUN= 2 initialize for OPT= TLO
FZFILE. Use LREC= 5400, options= TLO
* *
*----------------------------------------------------------*
* *
* Standard TPAR for this run are *
* ------------------------------ *
* *
* CUTGAM= 1.00 MeV CUTELE= 1.00 MeV CUTNEU= 1.00 MeV *
* CUTHAD= 1.00 MeV CUTMUO= 1.00 MeV *
* BCUTE = 1.00 MeV BCUTM = 1.00 MeV *
* DCUTE = 1.00 MeV DCUTM = 1.00 MeV PPCUTM= 2.04 MeV *
* IPAIR = 1. ICOMP = 1. IPHOT = 0. *
* IPFIS = 0. IDRAY = 1. IANNI = 1. *
* IBREM = 1. IHADR = 1. IMUNU = 1. *
* IDCAY = 1. ILOSS = 1. IMULS = 1. *
* IRAYL = 0. ILABS = 0. ISYNC = 0. *
* ISTRA = 0. *
* *
* *
************************************************************
###################################################################################################################
# XNORTH,XSOUTH,YEAST,YWEST,ZBOTTOM,ZTOP OF MACRO IN MARS : 119.77 7810.18 49.00 1266.00 69.00 1023.70 #
###################################################################################################################
####################################################################################################################
# XMUTIN,XMUTIS,YMUTIE,YMUTIW,ZMUTIB,ZMUTIT OF MACRO IN MARS 130.99 7808.96 48.75 1266.75 74.49 1019.91 #
####################################################################################################################
Initializing NEUGEN generator
Warning: using only last 2 digits of NUSEed 19780501
Possible neutrino fluxes :
1 = Volkova
2 = Mitsui
3 = Butkevitch
4 = Bartol (including Bartol-low)
5 = Bartol (old) (including Bartol-low)
6 = Bartol "NEW" (low energy)
(Note that only 4, 5 and 6 are valid choices below 1 GeV)
Give index JFLUX : ENTER 1 FOR SOLAR MINIMUM
2 FOR SOLAR MAXIMUM
3 FOR SOLAR MID (= (MIN+MAX)/2)
TAUOLA enabled for Tau decays
mc_fill_rsqdat: read arrays from rsq_xsects.dat
# Rate of nu 3 (kton year)**-1 .1709165E+03
# Rate of nu 4 (kton year)**-1 .6302871E+02
Calculating interaction rates
------------------------------
| nu | p/n | rate |
| -- | --- | ---- |
| numu | p | 1.519E-31 |
| numu | n | 4.157E-31 |
| numubar | p | 1.422E-31 |
| numubar | n | .671E-31 |
------------------------------
Total rate: p: 0.2940941340E-30 n: 0.4827468478E-30
Generation volume:
-35.02490234 < x < 7964.975098
-42.50000000 < y < 1357.500000
21.34996033 < z < 1071.349976
Please wait while I throw 1000000 darts at the generation volume
100000 ...
200000 ...
300000 ...
400000 ...
500000 ...
600000 ...
700000 ...
800000 ...
900000 ...
1000000 ...
3145 events changed to absorber
---------------------------------------------------------------
| # | medium | %ofvol | density| rate/nucl | ratinvol |
|----|--------------|--------|--------|------------|----------|
| 1 | AIR | 74.39% | .0012 | .3885E-30 | 2.467 | 8748.06 +/- 10.14
| 2 | GAS | 3.89% | .0010 | .3884E-30 | .107 | 456.99 +/- 2.32
| 3 | IRON | 1.16% | 7.8700 | .3949E-30 | 254.773 | 136.10 +/- 1.27
| 4 | SCINTILLATOR | 5.47% | .8700 | .3872E-30 | 130.536 | 643.32 +/- 2.75
| 5 | PVC | .00% | 2.2650 | .3885E-30 | .000 | .00 +/- .00
| 6 | ABSORBER | 12.59% | 2.1700 | .3886E-30 | 752.208 | 1480.84 +/- 4.17
| 7 | PLASTIC | 1.79% | 1.6300 | .3894E-30 | 80.647 | 210.97 +/- 1.58
| 8 | | .00% | .0000 | .0000E+00 | .000 | .00 +/- .00
| 9 | | .00% | .0000 | .0000E+00 | .000 | .00 +/- .00
| 10 | GS_ROCK | .43% | 2.7100 | .3886E-30 | 32.183 | 50.73 +/- .77
| 11 | MINERAL_OIL | .28% | .8700 | .3872E-30 | 6.691 | 32.98 +/- .62
---------------------------------------------------------------
Total rate 1259.611328 ev/yr
10000 events = 7.938956738 years
------------------------------------------------------------------------
For the nue events, here are the relevant cross section calculations:
# Rate of nu 1 (kton year)**-1 .9118724E+02
# Rate of nu 2 (kton year)**-1 .2897952E+02
Calculating interaction rates
------------------------------
| nu | p/n | rate |
| -- | --- | ---- |
| nue | p | .777E-31 |
| nue | n | 2.251E-31 |
| nuebar | p | .657E-31 |
| nuebar | n | .306E-31 |
------------------------------
Total rate: p: 0.1433841208E-30 n: 0.2556417389E-30
Generation volume:
-35.02490234 < x < 7964.975098
-42.50000000 < y < 1357.500000
21.34996033 < z < 1071.349976
Please wait while I throw 1000000 darts at the generation volume
100000 ...
200000 ...
300000 ...
400000 ...
500000 ...
600000 ...
700000 ...
800000 ...
900000 ...
1000000 ...
3153 events changed to absorber
---------------------------------------------------------------
| # | medium | %ofvol | density| rate/nucl | ratinvol |
|----|--------------|--------|--------|------------|----------|
| 1 | AIR | 74.35% | .0012 | .1996E-30 | 1.266 | 8743.74 +/- 10.14
| 2 | GAS | 3.90% | .0010 | .1995E-30 | .055 | 458.62 +/- 2.32
| 3 | IRON | 1.17% | 7.8700 | .2034E-30 | 132.816 | 137.77 +/- 1.27
| 4 | SCINTILLATOR | 5.43% | .8700 | .1988E-30 | 66.503 | 638.38 +/- 2.74
| 5 | PVC | .00% | 2.2650 | .1996E-30 | .000 | .00 +/- .00
| 6 | ABSORBER | 12.68% | 2.1700 | .1996E-30 | 389.143 | 1491.32 +/- 4.19
| 7 | PLASTIC | 1.77% | 1.6300 | .2001E-30 | 40.999 | 208.73 +/- 1.57
| 8 | | .00% | .0000 | .0000E+00 | .000 | .00 +/- .00
| 9 | | .00% | .0000 | .0000E+00 | .000 | .00 +/- .00
| 10 | GS_ROCK | .42% | 2.7100 | .1996E-30 | 16.046 | 49.24 +/- .76
| 11 | MINERAL_OIL | .27% | .8700 | .1988E-30 | 3.356 | 32.21 +/- .62
---------------------------------------------------------------
Total rate 650.1838989 ev/yr
10000 events = 15.38026428 years
Bob