Calibration events in NSB only simulations skipped

[kpfrang]

When reading NSB simulation file it correctly recognizes that calibration events and prints RDLR: CALIBEVENT!. However, I cannot access these events which makes me think that they are skipped.

---

Example for gamma_20deg_180deg_run4___cta-prod3-demo_desert-2150m-Paranal-baseline-5x.simtel.gz

As you can see I don't have any event read with event_source for the NSB file, even though 100 events should be contained in the file (which equals the number for times RDLR: CALIBEVENT! is printed out).

Should they not be readable in exactly the same way as other simtel files as well?

[kosack]

We've never tried reading calib events, so it may require some modification. That means we should also support getting the event trigger type flag so we can cut on that in ctapipe.

This probably just comes from the fact that pyhessio is based on read_hess.c, which was written to analyze normal events only, so it may skip others.

[GernotMaier]

I am using in my hessio reader used for eventdisplay something like:

                   if( find_io_block( iobuf, &item_header ) != 0 )
                    {
                            break;
                    }       

                    // check header types
                    switch( ( int ) item_header.type )
                    {            
                            // ... obviously lot's of other 'case' statements here

                            case IO_TYPE_HESS_CALIBEVENT:
                            {
                                    int type = -1;
                                    rc = read_hess_calib_event( iobuf, &hsdata->event, -1, &type );
                                    if( verbose || rc != 0 )
                                    {
                                            printf( "read_hess_calib_event(), rc = %d, type=%d\n", rc, type );
                                    }
                                   
                                   //.. do here whatever you need to do to convert / expose the data
                            }
                 }

......