#Tb_hyst

off_cplus = -6.5
off_cminus = -2.5
rrate = 1.0        # ramping speed of mag field


e1 = 1235.8        # energy 1 (usually on edge)
e2 = 1229       # energy 2 (pre edge)

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

set_slit(50.0)
set_har(3.0)
set_range('k1', 10)
set_range('k2', 10)


open_valve()

loop_count = 0
while loop_count < nr_loops:

    set_pol_cplus(offset = off_cplus)
    set_file('Tb_hyst_plus')
    hyst_cont(field_x, B1, B2, rrate, [e1,e2])
    
    set_pol_cminus(offset = off_cminus)
    set_file('Tb_hyst_minus')
    hyst_cont(field_x, B2, B1, rrate, [e1,e2])
    
    set_pol_cminus(offset = off_cminus)
    set_file('Tb_hyst_minus')
    hyst_cont(field_x, B1, B2, rrate, [e1,e2])
    
    set_pol_cplus(offset = off_cplus)
    set_file('Tb_hyst_plus')
    hyst_cont(field_x, B2, B1, rrate, [e1,e2])

    loop_count = loop_count + 1
    print "Loop # ",loop_count," finished."

close_valve()


#XMCD in remanence after -4T

#Fe XMCD
off_plus=-0.5
off_minus=2.0

Estart = 690
Eend   = 750

sample_name = 'Fe'
nr_loops=1

set_har(1.0)
set_slit(5.0)
set_range('k1', 9)
set_range('k2', 9)
open_valve()

loop_count = 0

while loop_count < nr_loops:
    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')
    loop_count = loop_count + 1   

 
close_valve()

#Tb XMCD

off_plus=-6.5
off_minus=-2.5


Estart=1215
Eend=1300

sample_name = 'Tb'

nr_loops=2

set_slit(50.0)
set_har(3.0)
set_range('k1', 10)
set_range('k2', 10)
open_valve()

loop_count = 0
while loop_count < nr_loops:
    otf(start=Estart, end=Eend, time=3, delay=10, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=3, delay=10, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=3, delay=10, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')   
    otf(start=Estart, end=Eend, time=3, delay=10, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')   
    otf(start=Estart, end=Eend, time=3, delay=10, mode = "CIRC +", offset = off_plus, name = sample_name + '_plus')
    otf(start=Estart, end=Eend, time=3, delay=10, mode = "CIRC -", offset = off_minus, name = sample_name + '_minus')
    loop_count = loop_count + 1

close_valve()