path = os.path.dirname("/Users/hamilton/idl/pro/Qubic/People/pierre/qubic/script/script_ga.py") nside = 256 input_map = np.zeros(12 * nside ** 2) input_map[0] = 1 #### QUBIC Instrument kmax = 2 qubic = QubicInstrument("monochromatic,nopol", nside=nside) pointings = create_random_pointings(1, 20) pointings[0, :] = [0.1, 0, 0] #### configure observation obs = QubicConfiguration(qubic, pointings) C = obs.get_convolution_peak_operator() P = obs.get_projection_peak_operator(kmax=kmax) H = P * C # Produce the Time-Ordered data tod = H(input_map) plot(tod) corner = qubic.pack(qubic.detector.corner).view(float).reshape((-1, 4, 2)) clf() x = corner[:, 0, 0] y = corner[:, 0, 1] scatter(x, y, c=tod.ravel(), marker="s", s=60) colorbar()