Hi,
I believe I have uncovered the problem...and there IS a problem,
regardless of
Hiro's message.
When the spot array was redefined to have a length of 9 in gcuser.inc, the
corresponding change was not made in gbr2c.f This results in the NEVPAR
parameter being overwritten at the end of an event, being set to 0.
Unfortunately,
gbr2c looks for NEVPAR being returned with a negative value to write the
-802 end-of-event flag to the cdat file. This never happens. The
consequence in the
BrGeantInput only sees one event.
I will update gbr2c.f later this afternoon after I am sure that I have
not introduced any
new errors.
Regards, steve
Hironori Ito wrote:
> 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 : Wed Aug 01 2001 - 15:28:28 EDT