adding HWP support to OpSimConviqt.

src/toast/pipeline_tools/sky_signal.py

@@ -296,6 +296,55 @@ def apply_conviqt(args, comm, data, cache_prefix="signal", verbose=True):
     return cache_prefix
 
 
+@function_timer
+def apply_weighted_conviqt(args, comm, data, cache_prefix="signal", verbose=True):
+    if (
+        args.conviqt_sky_file is None
+        or args.conviqt_beam_file is None
+        or not args.simulate_sky
+    ):
+        return None
+
+    log = Logger.get()
+    timer = Timer()
+    timer.start()
+
+    if comm.world_rank == 0 and verbose:
+        log.info("Running Conviqt")
+
+    verbosity = 0
+    if verbose:
+        verbosity = 1
+    if args.debug:
+        verbosity = 10
+
+    conviqt = OpSimWeightedConviqt(
+        getattr(comm, "comm_" + args.conviqt_mpi_comm),
+        args.conviqt_sky_file,
+        args.conviqt_beam_file,
+        lmax=args.conviqt_lmax,
+        beammmax=args.conviqt_beam_mmax,
+        pol=True,
+        fwhm=args.conviqt_fwhm,
+        order=args.conviqt_order,
+        calibrate=True,
+        dxx=args.conviqt_dxx,
+        out=cache_prefix,
+        remove_monopole=args.conviqt_remove_monopole,
+        remove_dipole=args.conviqt_remove_dipole,
+        normalize_beam=args.conviqt_normalize_beam,
+        verbosity=verbosity,
+    )
+    conviqt.exec(data)
+
+    if comm.comm_world is not None:
+        comm.comm_world.barrier()
+    if comm.world_rank == 0 and verbose:
+        timer.report_clear("Read and sample map")
+
+    return cache_prefix
+
+
 @function_timer
 def scan_sky_signal(
     args, comm, data, cache_prefix="signal", verbose=True, pixels="pixels", mc=0