scan,
sky_in,
CESnumber=CESnumber,
projection=params.projection,
nside_out=params.nside_out,
pixel_size=params.pixel_size,
width=params.width,
array_noise_level=params.array_noise_level,
array_noise_seed=seeds_for_noise[CESnumber],
mapping_perpair=False)
## Initialise map containers for each processor
if pos_CES == 0:
sky_out_tot = OutputSkyMap(projection=tod.projection,
nside=tod.nside_out,
obspix=tod.obspix,
npixsky=tod.npixsky,
pixel_size=tod.pixel_size)
## Scan input map to get TODs
d = np.array(
[tod.map2tod(det) for det in range(inst.focal_plane.nbolometer)])
## Inject crosstalk
inject_crosstalk_inside_SQUID(d,
squid_ids,
bolo_ids,
radius=args.radius,
mu=args.mu,
sigma=args.sigma,
beta=args.beta,
scan,
sky_in,
CESnumber=CESnumber,
projection=params.projection,
nside_out=params.nside_out,
pixel_size=params.pixel_size,
width=params.width,
array_noise_level=params.array_noise_level,
array_noise_seed=seeds_for_noise[CESnumber],
mapping_perpair=params.mapping_perpair)
## Initialise map containers for each processor
if pos_CES == 0:
sky_out_tot = OutputSkyMap(projection=tod.projection,
nside=tod.nside_out,
obspix=tod.obspix,
npixsky=tod.npixsky,
pixel_size=tod.pixel_size)
## Scan input map to get TODs
for pair in tod.pair_list:
d = np.array([tod.map2tod(det) for det in pair])
## Project TOD to maps
tod.tod2map(d, sky_out_tot)
MPI.COMM_WORLD.barrier()
## Coaddition over all processors.
## Note that all processors will then have the coadded data.
## If you want informations at the level of each CES (or group of),
tod = TimeOrderedDataPairDiff(
inst, scan, sky_in,
CESnumber=CESnumber,
projection=params.projection,
nside_out=params.nside_out,
pixel_size=params.pixel_size,
width=params.width,
array_noise_level=params.array_noise_level,
array_noise_seed=seeds_for_noise[CESnumber],
mapping_perpair=params.mapping_perpair)
## Initialise map containers for each processor
if pos_CES == 0:
sky_out_tot = OutputSkyMap(projection=tod.projection,
nside=tod.nside_out,
obspix=tod.obspix,
npixsky=tod.npixsky,
pixel_size=tod.pixel_size)
## Scan input map to get TODs
d = []
for det in tqdm(range(inst.focal_plane.nbolometer)):
d.append(tod.map2tod(det))
## Project TOD to maps
tod.tod2map(np.array(d), sky_out_tot)
MPI.COMM_WORLD.barrier()
## Coaddition over all processors.
## Note that all processors will then have the coadded data.