def get_jeol_mqmas_phasecycling(transition):
    """ Returns string for MQMAS phase cycling corresponding to the given transition.
    """
    if transition == 3:
        phase_cycling = """    obs_phs_MQexc = {0,30,60,90,120,150,180,210,240,270,300,330};
	obs_phs_CTsel = {12(0),12(45) ,12(90),12 (135),12(180),12(225), 12(270),12(315)};
	obs_phs_phase_acq = {3(0,270,180,90),3(90,0,270,180),3(180,90,0,270),3(270,180,90,0)};"""
        required_scan = int(96)
    elif transition == 5:
        phase_cycling = """    obs_phs_MQexc = {0,18,36,54,72,90,108,126,144,162,180,198,216,234,252,270,288,306,324,342};
	obs_phs_CTsel = {20(0),20(45) ,20(90),20 (135),20(180),20(225), 20(270),20(315)};
	obs_phs_phase_acq = {5(0,270,180,90),5(90,0,270,180),splitt1 = numpy.full(3,'')5(180,90,0,270),5(270,180,90,0)};"""
        required_scan = int(160)
    elif transition == 7:
        phase_cycling = """    obs_phs_MQexc =  {0,90/7,2*90/7,3*90/7,4*90/7,5*90/7,6*90/7,90,8*90/7,9*90/7,10*90/7,11*90/7,12*90/7,13*90/7,180,15*90/7,16*90/7,17*90/7,18*90/7,19*90/7,20*90/7,270,22*90/7,23*90/7,24*90/7,25*90/7,26*90/7,27*90/7};
	obs_phs_CTsel = {28(0),28(45) ,28(90),28(135),28(180),28(225),28(270),28(315)};
	obs_phs_phase_acq = {7(0,270,180,90),7(90,0,270,180),7(180,90,0,270),7(270,180,90,0)};
"""
        required_scan = int(224)
    elif transition == 9:
        phase_cycling = """    obs_phs_MQexc = {0,10,20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240,250,260,270,280,290,300,310,320,330,340,350};
	obs_phs_CTsel = {36(0),36(45) ,36(90),36(135),36(180),36(225),36(270),36(315)};
	obs_phs_phase_acq = {9(0,270,180,90),9(90,0,270,180),9(180,90,0,270),9(270,180,90,0)};
"""
        required_scan = int(288)
    return phase_cycling, required_scan
