/tibet/work01/kasahara/CosmosData/VFluxmu/icrmcf26/muons etc.
(6 stations).  By NewLund (>9GeV).

 CosZenith = (1.0,1.00),
 CutOffFile = '/tibet/home/kasahara/Cosmos/Contrib/Data/CutoffNew/Kamioka',
 Ddelta = 5.00,
 DeadLine = ' ',
 DepthList = 10, 20, 40, 50, 100, 200, 500, 1000, 2000, 4000, 5000,
        6000, 7000, 8000, 10300,
 DestEventNo = -100000000, -1000000,

 
 primary > 6 GeV. BESS(p+He) + CNO

# -----------------# of Component defined= 3----------------------
# Seq.# Code  Comp.  E_ Unit  E_type Diff/Integ  Emin    Emax   # of seg. inte.flux 
#   1     6   p        gev     ke/n       d      6.00    0.440E+09  53  273.    
#   2    10   he       gev     ke/n       d      6.00    0.800E+09  53  15.8    
#   3    12   cno      gev     ke/n       d      6.00    0.300E+10  16 0.955    
  getST3 with fai=0~2pi may give you a different result for cos=1
  both should give the same result but due to the
  table usage, the results are different. 

           cos value=   1.00000000000000     
      X-axis From South=   89.9999655705780      deg

  No rigidity cut is assumed
  primary  Int(dI/dE)  sum(cumlative);
  p                272.5682         272.5682    
  he               15.78990         288.3581    
  cno             0.9552743         289.3134    
 If N primaries are generated in simulation,
  STdOmega= N/sum  

N=6x10^6
 1/STdOmega = 289.31/N =289.31/6x10^6 = 4.822E-05

The data:  where code subcode Ek (no charge)

-----------------------------

500000 x6 events for Old Lund (file name ByOldLund)
(up to 1000 GeV).

 1/STdOmega = 289.31/N =289.31/3x10^6 =  9.6447E-05

data:  where code subcode charge Ek


ByNewLund2:
>9 GeV ; NewLund.  500000 x 6 = 3x10^6  events


