[Brahms-dev-l] Dst problem

From: Radoslaw Karabowicz <r.karabowicz@if.uj.edu.pl>
Date: Thu Jan 12 2006 - 06:34:12 EST
As some of you have already heard it, I have problems
with the dst3 version of auau run04 200.
The problem has been localized and it is connected
with the lack of G.fNoFsTrack, which makes it impossible
to nicely loop over FS tracks.
It shows up like this:
A. Loop over Ffs tracks:
EVENT 9080
track 0 out of 1 tracks, momentum 34.1937 richRadius = 8.47939
EVENT 9081
track 0 out of 3 tracks, momentum 24.3179 richRadius = 8.89382
track 1 out of 3 tracks, momentum 2.89575e+32 richRadius = 0
track 2 out of 3 tracks, momentum 1.33417e+31 richRadius = 0
EVENT 9082
track 0 out of 2 tracks, momentum 22.971 richRadius = 9.34088
track 1 out of 2 tracks, momentum 2.89575e+32 richRadius = 0
EVENT 9083
EVENT 9084
EVENT 9085
track 0 out of 3 tracks, momentum 22.971 richRadius = 9.34088
track 1 out of 3 tracks, momentum 2.89575e+32 richRadius = 0
track 2 out of 3 tracks, momentum 1.33417e+31 richRadius = 0
EVENT 9086
EVENT 9087
track 0 out of 1 tracks, momentum 22.971 richRadius = 9.34088
EVENT 9088
track 0 out of 1 tracks, momentum 22.971 richRadius = 9.34088
EVENT 9089
track 0 out of 1 tracks, momentum 22.971 richRadius = 9.34088

B. Loop over Bfs tracks
EVENT 9080
track 0 out of 1 tracks, momentum 34.1937 richRadius = 8.47939
EVENT 9081
track 0 out of 1 tracks, momentum 24.3179 richRadius = 8.89382
EVENT 9082
track 0 out of 1 tracks, momentum 22.971 richRadius = 9.34088
EVENT 9083
track 0 out of 1 tracks, momentum 22.971 richRadius = 9.34088
EVENT 9084
track 0 out of 1 tracks, momentum 22.971 richRadius = 9.34088
EVENT 9085
EVENT 9086
EVENT 9087
EVENT 9088
EVENT 9089

So, it can be only solved by putting:
       if ( (TMath::Abs(bdst->FS_fP[t]-oldMomentum)<0.001) &&
            (TMath::Abs(bdst->FS_fRich_fRadius[t]-oldRichRadius)<0.001) )
        continue;
      oldMomentum = bdst->FS_fP[t];
      oldRichRadius = bdst->FS_fRich_fRadius[t];
in the track loop -- but it's not elegant;)

Is it a bug in DST? Or is there some nice way to get the number of fs
tracks?

And another problem that seems strange to us. Why are in the loop over
Bfs tracks so many Bfs tracks that doesn't have FFS? Pawel thinks that:
- number of such tracks should be something like 10% at the beginning (and
it seems that is
as high as 50%),
- they shouldn't be put into Dst's anyway (they are, or at least information
about them).

Greetings, radek


_______________________________________________
Brahms-dev-l mailing list
Brahms-dev-l@lists.bnl.gov
http://lists.bnl.gov/mailman/listinfo/brahms-dev-l
Received on Thu Jan 12 06:35:33 2006

This archive was generated by hypermail 2.1.8 : Thu Jan 12 2006 - 06:35:40 EST