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,
seed=args.seed)
## Project TOD to maps
tod.tod2map(d, sky_out_tot)
MPI.COMM_WORLD.barrier()
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),
## use instead:
## final_map = OutputSkyMap(nside=nside_out, obspix=tod.obspix)
## final_map.coadd_MPI(sky_out_tot, MPI=MPI)
sky_out_tot.coadd_MPI(sky_out_tot, MPI=MPI)
if rank == 0:
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.
## If you want informations at the level of each CES (or group of),
## use instead:
## final_map = OutputSkyMap(nside=nside_out, obspix=tod.obspix)
## final_map.coadd_MPI(sky_out_tot, MPI=MPI)
sky_out_tot.coadd_MPI(sky_out_tot, MPI=MPI)
## Check that output = input