Create workflow to replicate purpose of temp_file_renaming about termfactory HOT 4 OPEN

    # DEBUG - temporary override
    maximum_h_rank, maximum_cc_rank, maximum_T_rank, eT_taylor_max_order, omega_max_order = truncations
    path = "./eT_zhz_eqs"
    path += f"_H({maximum_h_rank})"
    path += f"_P({omega_max_order})"
    path += f"_T({maximum_T_rank})"
    path += f"_exp({eT_taylor_max_order})"
    path += f"_Z({maximum_cc_rank})"
    path += ".py"

for p, e, z in it.product(
    [2, 4],
    [1, 2, 3, 4],
    [1, 2, 3],
):

    # for the 'eT_z_t ansatz' only
    maximum_h_rank = 2
    maximum_cc_rank = z  # max order of Z
    maximum_T_rank = 1  # max order of T
    eT_taylor_max_order = e  # Taylor expansion of eT
    omega_max_order = p  # max order of projection operator

    # need to have truncation of e^T
    eT_z_t_truncations = maximum_h_rank, maximum_cc_rank, maximum_T_rank, eT_taylor_max_order, omega_max_order

    switch = 4

    if switch == 1:
        generate_latex_files(
            eT_z_t_truncations,
            only_ground_state=True,
            remove_f_terms=False,
            thermal=False,
            file='eT_z_t ansatz'
        )
    if switch == 2:
        generate_latex_files(
            truncations,
            only_ground_state=False,
            remove_f_terms=False,
            thermal=False,
            file='z_t ansatz'
        )
    if switch == 3:
        generate_latex_files(
            truncations,
            only_ground_state=True,
            remove_f_terms=False,
            thermal=False,
            file='full cc'
        )

    if switch == 4:
        generate_eT_zhz_python(eT_z_t_truncations, only_ground_state=True)
        # generate_full_cc_python(truncations, only_ground_state=True)

    print("We reached the end of main")