Hello. I was able to run brag (Note: we should change this
name) and look at it by brat2 without any problems. Therefore, it is
something else (Platform, compiler, your machine is on strike etc...)
Hiro
On Mon, 30 Jul 2001, Stephen J. Sanders wrote:
> Hi, I'm having trouble analyzing a recently generated .cdat file using
> brat2. I'm hoping someone else might have done this so that I can rule out
> some of the possible sources for the problem.
>
> I generated the .cdat file using the "standard" gbr2c.f routine. I only had
> the mult detectors set on.
>
> Because the programs that I had used previously were not working with
> brat2, I started
> issuing some of the basic commands interactively to see how far I could
> get (I've edited out lines where I made typing mistakes...):
>
> brat [0] BrGeantInput * gInput = new BrGeantInput("GeantInput")
> brat [1] gInput->Open("simu/sim_104.cdat")
> (unsigned char)1
> brat [4] BrEvent* gEvent =new BrEvent("Event",0,0)
> brat [5] gInput->Event(gEvent)
> brat [6] gEvent->ListObjects()
> 0 - 0x1086a588 - GeantHits TILE - GeantHits TILE
> 1 - 0x1086a698 - GeantHits STRA - GeantHits STRA
> 2 - 0x1086a7a8 - GeantHits STRB - GeantHits STRB
> 3 - 0x1086a8b8 - GeantHits STRC - GeantHits STRC
> 4 - 0x1086a9c8 - GeantHits STRD - GeantHits STRD
> 5 - 0x1086aad8 - GeantHits STRE - GeantHits STRE
> 6 - 0x1086abe8 - GeantHits STRF - GeantHits STRF
> 7 - 0x1086acf8 - GeantHits STRG - GeantHits STRG
> 8 - 0x10865ee8 - GeantTracks - GeantTracks
> 9 - 0x1086ae08 - GeantHeader - GeantInput
> brat [7] BrGeantHeader
> *gHeader=(BrGeantHeader*)gEvent->GetObject("GeantHeader$
> brat [8] b=gHeader->GetImpactParameter()
> Warning: Automatic variable b allocated in global scope
> FILE:/tmp/fileFr1vHR_cint LINE:1
> (const Double_t)1.11246051788330078e+01
> brat [9] vertex=gHeader->GetVertex()
> Warning: Automatic variable vertex allocated in global scope
> FILE:/tmp/file4lGYm0_cint LINE:1
> (const Double_t)(-1.81825294494628906e+01)
> brat [13] BrDataTable
> *GeantTileHit=(BrDataTable*)gEvent->GetObject("GeantHits$
> brat [14] GeantTileHit->Entries()
> (Int_t)1196378
> brat [15] BrDataTable * GeantTrack =
> (BrDataTable*)gEvent->GetObject("GeantTr$
> brat [16] GeantTrack->Entries()
> (Int_t)1390198
>
> The problem is that this is the ONLY event found (the brag run was for
> 10000 events,
> Au+Au at 200AGeV with 10fm<=b<=20fm), and the number of hits and tracks
> look
> much more consistent with 10000 peripheral events than 1 event at
> 11.12fm. It also
> took several minutes for the computer to analyze this one event, during
> which time
> my system "froze"--a situation I usually encounter when the physical memory
> is exhausted and I start using the swap area.
>
> This looks like a possible maping problem between brag and brat. I will
> try to
> track it down tomorrow,
> but since I am doing pretty basic stuff I wanted to check first
> if someone else has actually seen this work with brat2.
>
> Regards, Steve
>
>
>
This archive was generated by hypermail 2b30 : Tue Jul 31 2001 - 14:06:25 EDT