ctapipe Bi-Weekly Digest — Activity since 2026-03-19

» cta-observatory/ctapipe | cta-observatory/ctapipe

• cta-observatory/ctapipe repo issues [back to top]

  • Expose CORSIKA photons in the SimulatedEventContainer #2698 [vmharvey] enhancement

    Please describe the use case that requires this feature. I'm doing some studies that use the CORSIKA photon bunches...

    • vmharvey: I'm returning to this issue now because, with the pyeventio v2 rewrite, my workaround of...

    • maxnoe: Do you strongly prefer one way or the other between event.meta['simtel_event'] Yes, I'd just...

  • Reference metadata being checked before its is written in CrossValidator #2973 [kosack] bug

    Describe the bug We get a warning like: ctapipe/core/provenance.py:280: MissingReferenceMetadata: Could not read...

    • maxnoe: Yes, should be moved to the finish method.

  • Add a chunking strategy that adjust the user given bin size to create equal-size chunks #2972 [maxnoe] enhancement

    Please describe the use case that requires this feature. From a discussion with @mgaug, the chunking should work in...

    • kosack: Which chunking is this referencing? TableLoader?

    • maxnoe: The grouping used for aggregation in the statistics tool

  • Apply DVR methods to R1 #2975 [mescobargodoy] enhancement

    Please describe the use case that requires this feature. For the medium sized SCT (11328 pixels), we will require...

    • maxnoe: Naively I would think it is possible to add an extra flag/option to apply it to R1 instead of...

    • maxnoe: Another question: I assume we are talking about simulations here?

    • mescobargodoy: Okay thanks for the clarification. Yes, adding the option to write out DL0 works. Eventually we...

    • maxnoe: Yes, if the data is read-in as R1. Maybe the option to re-do DVR on DL0 would make sense again.

    • mescobargodoy: Yes, the plan is to read-in as R1.

  • Add `allowed_tels` option to `ApplyModels`-Tool #2978 [Hckjs] enhancement discussion

    If you currently want to do e.g. some benchmark studies on different array configurations from preprocessed DL1 MC's...

  • Implement simulated histogram shower-distribution event weighting #2977 [kosack] enhancement

    Please describe the use case that requires this feature. For upcoming CTAO simulation campaigns, they will start to...

    • maxnoe: @orelgueta what importance sampling are you going to use? Only the one resulting in non unity...

    • kosack: For context, here is what is in the simulated histogram for a merged DL2 file: Which if you...

    • orelgueta: At the moment the tests are only for the non-unity area weights. @mariakherlakian has only used...

    • maxnoe: Ok thanks. In this case we only need the area weights.

• cta-observatory/ctapipe repo pull_requests [back to top]

  • allow gain selection when skipping R1 calibration in SimTelEventSourve #2943 [TjarkMiener]

    • kosack

  • Add data quality-related subgroup definitions to HDF5 data format #2965 [mexanick]

    • kosack

  • Prepare release 0.30 #2968 [maxnoe] no-changelog-needed

    • mexanick

    • kosack

  • Implement CTAO CR56: Add field pixel_time_shift to R1 and DL0 #2964 [moralejo]

    New field which contains event-and-pixel-wise time corrections in ns (to be subtracted from reconstructed pixel times)

    • maxnoe

    • kosack

    • maxnoe

  • Introduce new `needs_atmosphere_profile` property to reconstructors #2958 [Hckjs] fix maintenance

    Adding also small test

    • kosack

  • Update psutil intersphinx url #2970 [maxnoe] no-changelog-needed

    The old one started to return a 404, there only seems to be "latest" now.

    • kosack

    • maxnoe: Issue fixed upstream

  • Fix TypeError in HDF5EventSource #2963 [yaochengchen]

    Replace dict_values indexing with next(iter(...)) in _is_hillas_in_camera_frame to avoid a Python 3 TypeError.

    • Hckjs: PyTables group._v_children actually returns a ProxyDict which implements its own .values()...

    • maxnoe: Given that we also want to completely remove the support for this (#2396), I will close this PR....

  • Fix incorrect warning in CrossValidator #2974 [kosack]

    Add output file to provenance only after it is closed. Fixes #2973

  • Store original simtel event in ArrayEventContainer #2976 [vmharvey]

    Closes #2698. Note that the telescope numbering scheme inside the simtel event is different to the final ArrayEvent.

  • add polygon_chord #2969 [burmist-git]

    Polygonal description of the chord. 096_arbitrary_shaped_mirror_and_polygon_chord.pdf Uploading...

    • kosack: Maybe we should consider adding shapely as a dependency? It might simplify some of these...

    • maxnoe: Maybe we should consider adding shapely as a dependency? I already looked into using shapely...

    • burmist-git: This is possible, but the function itself is very small, so using an external library would be...

    • mexanick

      • mexanick: I think we should deduce these (if needed) from the vertices. You can write a helper for this.

        • maxnoe: Yes, I'd also prefer the vertices as input. But maybe the function as it is here is fine and...

      • mexanick: These should be quantities, as we discussed offline, the best would be to use telescope...

        • mexanick: Also, perhaps tuples of (x,y) or arrays of them would be more straightforward to define points...

      • mexanick: Should be also a unit of length (or angular, equivalent in the telescope frame)

      • mexanick: why not using try: ... except ZeroDivisionError: # determinant is zero, the cord goes on the...

    • maxnoe: This function should be njit, as it will be called. If possible, it would also be good if it...

  • Polar PSFModel #2945 [TjarkMiener]

    • kosack: I think there are still a few bugs: seems to be clipping or quadrant error at negative lats...

    • maxnoe: If you print the individual psf.sum() like this: psf = psf_model.pdf( tel_id=1, ...

    • maxnoe: Looking at this, I also realize that the output unit of the PSF should be 1 / deg² or 1 / sr

    • maxnoe

    • maxnoe: Thanks @orelgueta

    • orelgueta: Just FYI, here is how it looks from raytracing the telescope: photon_distributions_all.pdf

    • kosack: For the normalization, I would also check your units, particularly the use of deg and rad in...

    • maxnoe: The issue with the normalization is definitely the missing the jacobian. You just multiply two...

    • maxnoe: Dividing by r works, but needs special treatment at r=0, which I didn't figure out yet.

    • kosack: Dividing by r works, but needs special treatment at r=0, which I didn't figure out yet. It...

    • maxnoe: Yes. That computation is what I didn't invest time into

    • kosack: I guess Tjark will do it, but it's pretty easy with SymPy. Here I did it with the definition...

    • mexanick: r==0 is not a problem, we have a lower limit for r, the whole model is not valid at low r, and...

    • maxnoe: the whole model is not valid at low r, and there's a condition on that. Also the other model...

    • kosack: the whole model is not valid at low r, and there's a condition on that. Also the other model...

  • Waveform noise adder #2847 [moralejo]

    First version of the component to add noise in MC waveforms, for adapting it to different observation conditions. It...

    • kosack: Let's merge this now, and open a second PR to add it to the ctapipe-process tool. (as agreed in...

    • maxnoe: Suggested change @njit ...

    • maxnoe: I find this name a bit confusing. Since "Shuffling the camera" for me sounds like we do mix...

      • moralejo: A pixel never uses a noise waveform from another pixel (we do not want that if we ever use e.g....

      • maxnoe: sample_pixels_independently?

    • maxnoe: Changelog is missing and two minor comments inline, otherwise good to go!

    • moralejo: Back to draft, I messed up when rebasing. I think instructions in...

    • maxnoe: Git seems to suggest doing a pull after the rebase, and this brings into the PR unrelated...

    • moralejo: Git seems to suggest doing a pull after the rebase, and this brings into the PR unrelated...

    • moralejo: Back to draft, I messed up when rebasing. I think instructions in...

    • maxnoe

    • maxnoe: Coverage is only missing for the numba njitted function, which are in fact executed

    • kosack

created at 2026-04-02 07:20 by dinghy 1.4.1.