Hi Richard,

Many thanks for the offer of help. I'm transferring the generator
output files to my unix SLAC account (/a/juno/u18/ey/homer) so that they can be
tested. Their names are
-rw-r-----   1 neal     ppg       298172 Feb  2 09:18 isa_susy.out
-rw-r-----   1 neal     ppg       345788 Feb  2 09:26 isa_susy2.out
-rw-r-----   1 neal     ppg        37436 Feb  2 12:52 isa_sugra.out

I generated them using ij737. As an example, I used the following
to create the first file:
isajet <</EOF
../isadecay.dat
susy.input
10
isa_susy.out
/EOF

Where susy.input contains:
cat susy.input
SAMPLE SUSY JOB
500.,100,1,10/
SUPERSYM
PT
50,100,50,100/
JETTYPE1
'GLSS','SQUARKS'/
JETTYPE2
'GLSS','SQUARKS'/
GAUGINO
60,1,40,40/
SQUARK
80.3,80.3,80.5,81.6,85,110/
FORCE
29,30,1,-1/
FORCE
21,29,1/
FORCE
22,29,2/
FORCE
23,29,3/
FORCE
24,29,4/
FORCE
25,29,5/
FORCE
26,29,6/
END
STOP

I'm in the process of producing some output from PYTHIA
to see if that works. By the way, I tried to build PYBMS
using the suggested pybms_57.tar file code. I quickly discovered
the AIX xlf is necessary. xlf is not on my OSF machine.
Is PYBMS available in standard fortran.

By the way, I did try providing all the arguments for nld_bat
and I got the same results.

The tar that I extracted the nld_ binaries from is:
lcd-working19990112-mask-binosf1.tar.gz and I'm
using the latest Small.ini and Large.ini detector definitions.

I can put particles into Beam.ini and nld_bat/nld_win seem to
successfully simulate them. EXCEPT ... unless I change the material
of the calorimeter to air, anything other than muons seems to take
an enormous time to simulate.

I suspect that there may be an incompatibility with the stdhep I'm
using. It is version V3_02.

The traceback for nld_bat, which I'm not sure is reliable for the reason that
David G. mentions on the web, is as follows:
 1 event read
signal Floating point exception at   [annih__3EGS:6 +0x23c,0x12020c1dc]  Source
not available
(dbx) where
>  0 annih__3EGS()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//egs/annih.cxx":6, 0x12020c1dc]
   1 interact__3EGSRC8PEGSData()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//egs/EGS.cxx":6, 0x1201fe338]
   2 interact__13EGSInteractorP9GParticlePC6Medium()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//egs/EGSInteractor.cxx":6,
0x1201ff6f4]
   3 endStep__10MCParticle()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/MCParticle.cxx":6,
0x120174094]
   4 propagate__10MCParticlePC6Medium()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/MCParticle.cxx":6,
0x1201716e8]
   5 propagate__10MCParticlePC6Medium()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/MCParticle.cxx":6,
0x1201719f4]
   6 propagate__10MCParticlePC6Medium()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/MCParticle.cxx":6,
0x1201719f4]
   7 propagate__10MCParticlePC6Medium()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/MCParticle.cxx":6,
0x1201719f4]
   8 propagate__10MCParticlePC6Medium()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/MCParticle.cxx":6,
0x1201719f4]
   9 propagate__6HepEvtPC6Medium()
["/cdf/data/gerdes/lcd/lcd-working/source//nld/HepEvt.cxx":6, 0x12007da9c]
  10 processEvent__5WorldR9Generator()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/World.cxx":6,
0x1201649d0]
  11 generate__5Gismo()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//gismo/Gismo.cxx":6,
0x120150ec0]
  12 execute__14GismoInterface()
["/cdf/data/gerdes/lcd/lcd-working/source//nld/main.cxx":6, 0x1200c551c]
  13 execute__10CheckPoint()
["/cdf/data/gerdes/lcd/lcd-working/source//nld/CheckPoint.cxx":6, 0x12011ef80]
  14 execute__11DataControl()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//arve/DataControl.cxx":6,
0x1201236b0]
  15 execute__4Arve()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//arve/Arve.cxx":6, 0x12012a398]
  16 execute__t13SimpleCommand1Z4Arve()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//arve/Arve.cxx":6, 0x12012d10c]
  17 for_each__H2ZPP7CommandZ8Executor_X01T0X11_X11()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//control/CompoundCommand.cxx":6,
0x12013dfe0]
  18 execute__15CompoundCommand()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//control/CompoundCommand.cxx":6,
0x12013d298]
  19 resume__11LoopCommand()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//control/LoopCommand.cxx":6,
0x12013bf6c]
  20 execute__11LoopCommand()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//control/LoopCommand.cxx":6,
0x12013bee8]
  21 execute__9EventLoop()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//control/EventLoop.cxx":6,
0x12013aaec]
  22 strike__4Menuc()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//control/Menu.cxx":6,
0x1201354bc]
  23 run__5Batch()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//batch/Batch.cxx":6,
0x120270dfc]
  24 run__4Menu()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//control/Menu.cxx":6,
0x120134254]
  25 run__4Arve()
["/cdf/data/gerdes/lcd/lcd-working/source/arve//arve/Arve.cxx":6, 0x12012a7a4]
  26 main() ["/cdf/data/gerdes/lcd/lcd-working/source//nld/main.cxx":6,
0x1200bac38]

Once again, thanks! By the way I guess that I'll see you at SLAC next week
after Tuesday. Let me know if there is anything you need from the CERN area.
I'll be departing for Yale on Saturday.

Cheers,
   Homer


On Wed, 3 Feb 1999, Richard Dubois wrote:

>  Hi Homer,
> 
>  Where does your input generator file come from? Is there a chance it is sick?
> Why don't you send me a copy and I can try it here. Have you run single pions,
> electrons and muons through NBeam?
> 
>  You should subscribe to the lcd-sim mailing list so others can be aware of
> pitfalls and solutions...
> 
> Richard
> 
> David Gerdes wrote:
> 
> > Hi Homer,
> >
> > I always get the mcfio_Block errors at the beginning, but have never seen
> > the crash you report. I've run through several thousand ttbar events
> > without an error. The exe is statically linked, so I don't think the
> > problem is a difference of runtime libraries. Frankly, I'm a bit stumped.
> > Are you using the most recent exe (1/24) from my web page? And the latest
> > and greatest detector.ini files from SLAC?
> >
> > A dumb question... are you specifying all the command-line parameters when
> > you run nld_bat? The full systax is:
> >
> > nld_bat Detector.ini HepEvt Infile NumEvents Outfile RunNum Randomseed
> >
> > Eg. when I run I use something like
> >
> > nld_bat Small.ini HepEvt pybms-ttbar500-01.dat 100
> >            sim/pybms-ttbar500-01.asc 1 52342
> >
> > Good luck,
> >
> >    Dave
> >
> > --------------------------------------
> > David Gerdes, University of Michigan
> > (734) 647-3807 / (734) 936-1817 FAX
> > [log in to unmask]
> > http://umaxp1.physics.lsa.umich.edu/~gerdes/
> >
> > On Wed, 3 Feb 1999, Homer Neal wrote:
> >
> > > Dear David,
> > >
> > > I'm using your executable on by osf1 system to simulate some LCD
> > > events. It works fine for BEAM inputs but when I try to input
> > > a hepevt file it crashes on the first event. I've tried many
> > > different events with low and high multiplicities and with
> > > eliminating the calorimeter material etc... and it always
> > > crashes. Do you have any suggestions?
> > >
> > > ylws41> nld_bat Small.ini hepevt isa_sugra.out
> > >  StdHepXdrReadOpen: successfully opened input stream 1
> > >           title: ISAJET stdHEP output file
> > >           date: Tue Feb  2 12:52:33 1999
> > >
> > >                     12 events
> > >                     4 blocks per event
> > >  mcfio_Block: Unable to find block i.d. 101 in Stream 1
> > >  mcfio_Block: Unable to find block i.d. 105 in Stream 1
> > >  Reading StdHep begin run record.
> > >  at event 1 with 10 particles
> > >  1 event read
> > > Floating exception
> > >
> > > Many thanks for any help you may provide and also for
> > > making your executables available. I've been having many
> > > difficulties building the executables on my system.
> > >
> > > Best regards,
> > >     Homer
> > >
> > >
> 
> --
> Richard Dubois
> SLD, Stanford Linear Accelerator Center
> [log in to unmask]
> http://www.slac.stanford.edu/~richard/
> 650-926-3824
> 650-926-2923 (FAX)
> 
>