FWD Software

[Akio Ogawa <akio AT bnl.gov>]

Hello

Please see modified /star/u/akio/forWeibin/test.C

It seems

   int nl = hits->GetLast();

gives "correct" number of hits minus 1. It looks like GetSize() and 
GetEntries() give current size of TRefArray which incldes allocated 
but not currently used elements.

Note that loop need to be

  int nl = hits->GetLast();
  for (int i=0; i<=nl; i++){}  //include i==nl!!!

to loop over all hits.

There is TRefArray::Compress(). But it makes GetLast() point to same
place -1 as GetSize/GetEntries, not otherway around.

Removing each by hand
  if(nl+1 < nh){ for(i=nl+1; i<nh; i++) {hits->RemoveAt(i);} }
makes GetEntries() consistent with GetLast() but GetSize() remains same.

I'm not sure how TRefArray are added/reset/reused/cleared in Mudst and
if this is what Daniel expects to see... Or if this is optimize for file 
size?

OGAWA, akio                 Email : akio AT bnl.gov
BrookHaven National Lab. Physics Dep. 510A 2-235
Upton, NY, 11973-5000          Tel. 631-344-5293

On Fri, 17 Nov 2023, Weibin Zhang wrote:

> Hi, Akio and Daniel:
> I found that the MuDst FCS cluster is not properly filled with corresponding 
> hits, specifically, one cluster has
> hits that have different detector Ids.  Here is a simple macro to reproduce 
> this problem:
>
> macro:  /gpfs/mnt/gpfs01/star/pwg/weibin/DY/test.C
> and the MuDst file: /star/data01/pwg/weibin/DY/condor/pi-/pi_0.MuDst.root
>
> Please check it.
>
> Best,
> Weibin