Skip to content

response Module

For simulating the detector response function, mostly to integrate charge

apply_response(response, electrons)

Apply the response function to get a signal amplitude, integral

Scale the response by the number of electrons. Retrieve the amplitude, integral of the response in ADC units. Signals above 11-bit max (4096) are clipped

Parameters:

Name Type Description Default
response ndarray

The response signal

 required
electrons float

The number of electrons to scale by

 required

Returns:

Type Description
tuple[float, float]

The signal amplitude, integral respectively
Source code in src/attpc_engine/detector/response.py 
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
@njit
def apply_response(response: np.ndarray, electrons: float) -> tuple[float, float]:
    """Apply the response function to get a signal amplitude, integral

    Scale the response by the number of electrons. Retrieve the amplitude,
    integral of the response in ADC units. Signals above 11-bit max (4096) are clipped

    Parameters
    ----------
    response: numpy.ndarray
        The response signal
    electrons: float
        The number of electrons to scale by

    Returns
    -------
    tuple[float, float]
        The signal amplitude, integral respectively

    """
    resp_sig = response * electrons
    resp_sig[resp_sig > 4095] = 4095  # Set saturated pads to max ADC
    return (resp_sig.max(), resp_sig.sum())

get_response(config)

Theoretical response function of GET electronics provided by the chip manufacturer. See https://doi.org/10.1016/j.nima.2016.09.018.

Parameters:

Name Type Description Default
config Config

The simulation configuration

 required

Returns:

Type Description
ndarray

Returns 1xNUM_TB array of the response function evaluated at each time bucket
Source code in src/attpc_engine/detector/response.py 
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
def get_response(config: Config) -> np.ndarray:
    """
    Theoretical response function of GET electronics provided by the
    chip manufacturer. See https://doi.org/10.1016/j.nima.2016.09.018.

    Parameters
    ----------
    config: Config
        The simulation configuration

    Returns
    -------
    np.ndarray
        Returns 1xNUM_TB array of the response function evaluated
        at each time bucket
    """
    response = np.zeros(NUM_TB)
    c1 = 4095 * E_CHARGE / config.elec_params.amp_gain / 1e-15
    tbs = np.linspace(0.0, NUM_TB, NUM_TB)
    c2 = tbs / (config.elec_params.shaping_time * config.elec_params.clock_freq * 0.001)
    response = c1 * np.exp(-3.0 * c2) * (c2**3) * np.sin(c2)
    # Cannot have negative ADC values
    response[response < 0] = 0

    return response