"Use Docker" documentation to setup environment and run

To run

tdis
-pnthreads=1
-pjana:nevents=10
-ppodio:output_file=/mnt/data/test_output_v01.root
-ppodio:print_collections=McHits
-pacts:geometry=/mnt/data/
/mnt/data/g4sbsout_EPCEvents_200000.txt

Geometry IDs

The geometry description is gives us that each pad is identified by plane, ring, pad:

• plane The mTPC will consist of 10 separate TPC volumes. Numbered 0 - 9 from upstream to downstream.
• ring The pads are arranged into 21 concentric rings, The rings are numbered from 0 to 20 sequentially, with 0 being the innermost ring.
• pad Each ring has 122 pads, 0 is at or closest to phi=0 and numbering is clockwise

This is convenient to encode in a singe uint32_t ID. We assign 8 bits each to plane and ring, and 16 bits to pad. The encodePadID function packs these three values into a single uint32_t using bitwise shifts and OR operations.

import numpy as np

def encode_pad_id(plane, ring, pad):
    """Encode plane, ring, pad into an uint32_t id.
    Works with both scalar values and numpy arrays."""
    plane = np.asarray(plane, dtype=np.uint32)
    ring = np.asarray(ring, dtype=np.uint32)
    pad = np.asarray(pad, dtype=np.uint32)
    
    pad_id = (plane << 24) | (ring << 16) | pad
    return pad_id

def decode_pad_id(pad_id):
    """Decode pad_id into plane, ring, pad.
    Works with both scalar values and numpy arrays."""
    
    pad_id = np.asarray(pad_id, dtype=np.uint32)
    plane = (pad_id >> 24) & 0xFF  # Upper 8 bits
    ring = (pad_id >> 16) & 0xFF   # Next 8 bits
    pad = pad_id & 0xFFFF          # Lower 16 bits
    return plane, ring, pad