Are brag and brat2 compatible?

From: Stephen J. Sanders (ssanders@ku.edu)
Date: Mon Jul 30 2001 - 22:39:05 EDT

  • Next message: Christian Holm Christensen: "Re: Are brag and brat2 compatible?"

    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 : Mon Jul 30 2001 - 22:40:39 EDT