 |
 |
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) > >