On Thu, Sep 18, 2014 at 6:35 AM, Rafayel Paremuzyan <[log in to unmask]> wrote:

Thank you Omar,

After that fixes I don't have a segmentation.

The recon LCIO was deleted, and I repeated readout, reconstruction and dst_maker again,
so these new files are different.

I put them in this directory
/w/hallb-1/hps/rafopar/tests/Omar_check

Rafayel

On 09/18/2014 12:53 AM, Omar Moreno wrote:

Hi Rafayel,

Can you send me or point me to the recon LCIO file that you used to create the DST?

There are a few reasons your code isn't working:

1) You don't instantiate the HpsEvent pointer so it's pointing to junk. So when you set the branch address, it's being set to some random spot in memory. Instead do the following:

HpsEvent *ev1 = new HpsEvent();

before you set the branch address.

2) The track index starts at 0. So getting a track from an event as follows:

for(int j = 0; j < n_svt_trk; j++){

svt_tr = ev1->getTrack(j+1)

}

Will eventually cause the code to crash because the track j+1 doesn't exist.

3) Most of the tracks in the event point to NULL. If I check for NULL tracks and skip them, everything runs fine (after the modification above were made). I'm not sure if this is a problem with the DST maker or the recon file you used to create the DST. I'll need to look at the recon file and generate the DST myself so I can debug it. I haven't seen this before so I need to dig a bit to find the issue.

Another thing, since you are linking against the HpsEvent library, you don't need to load it inside your code. Specifically, the following line of code is unnecessary:

gSystem->Load("../hps-dst/build/lib/libHpsEvent.so");

You do need to make sure that the library path is added to LD_LIBRARY_PATH.

--Omar Moreno

On Wed, Sep 17, 2014 at 8:01 PM, Rafayel Paremuzyan <[log in to unmask]> wrote:

Dear Omar,

I am getting a segmentation fault error when I am reading dst files, and I didn't manage to figure out
the cause.

I shrank the code to the minimum to reproduce the problem (it is attached).
Could you please take a look into it?

Some additional info that might help to find the cause of error:

The way this file produced:
A slic output file from /mss/hallb/hps/production/slic/tritrig/2pt2/
was read out with HPS_Java using HPS2014TruthReadoutToLcio.lcsim and then
reconstructed using HPS2014OfflineTruthRecon.lcsim, later with the latest version of dst (downloaded yesterday)
is converted into root file (this file is also attached).

The compile command:
g++ test3.cc -o test3.exe -Wl,--no-as-needed -L/usr/local/hps_dst/hps-dst/build/lib -lHpsEvent -I/usr/local/hps_dst/hps-dst/include/hps_event/ -I/usr/local/hps_dst/hps-dst/include/dst -I/usr/local/hps_dst/hps-dst/include/hps_event/utils `root-config --cflags` `root-config --libs`

Rafo

########################################################################
Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1