#Fe XMCD
off_plus = 0.0
off_minus = 2.0

Estart = 630
Eend   = 670

sample_name = 'Mn'
nr_loops=1

set_range('k1', 10)
set_range('k2', 10)
set_au_mesh(125.0)
set_har(1)

set_slit(3.0)
set_energy(Estart)

open_valve()

for i in range(nr_loops):
    wait_temp()
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
 
close_valve()

#Fe hyst
off_plus = 0.0
off_minus = 2.0

sample_name = 'Mn'
rrate = 0.7        # ramping speed of mag field
e1 = 639        # energy 1 (usually on edge)
e2 = 637        # energy 2 (pre edge)

B1 = 6.8            # field span
B2 = -B1
nr_loops = 1        # number of loops to acquire

set_range('k1', 10)
set_range('k2', 10)
set_au_mesh(125.0)
set_har(1)

set_slit(3.0)
set_energy(e1)

open_valve()

for i in range(nr_loops):

    set_pol_cplus(offset = off_plus)
    set_file(sample_name + '_hyst_plus')
    hyst_cont(field_x, B1, B2, rrate, [e1, e2])
    
    set_pol_cminus(offset = off_minus)
    set_file(sample_name + '_hyst_minus')
    hyst_cont(field_x, B2, B1, rrate, [e1, e2])

    set_pol_cminus(offset = off_minus)
    set_file(sample_name + '_hyst_minus')
    hyst_cont(field_x, B1, B2, rrate, [e1,e2])
    
    set_pol_cplus(offset = off_plus)
    set_file(sample_name + '_hyst_plus')
    hyst_cont(field_x, B2, B1, rrate, [e1,e2])
    

    print "Loop # ",i+1," finished."

close_valve()

#Fe XMCD
off_plus = 0.0
off_minus = 2.0

Estart = 630
Eend   = 670

sample_name = 'Mn'
nr_loops=1

set_range('k1', 10)
set_range('k2', 10)
set_au_mesh(125.0)
set_har(1)

set_slit(3.0)
set_energy(Estart)

open_valve()

for i in range(nr_loops):
    wait_temp()
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    otf(start=Estart, end=Eend, time=2, delay=5, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
 
close_valve()