noatoms  = float(report['ANDOR']['noatoms'])

#SEQUENCE
s=seq.sequence(stepsize)
s=gen.initial(s)

s.digichg('odt7595',0)

#Keep ODT on
ODT = gen.bstr('ODT',report)
if ODT == True:
    s.digichg('odtttl',1)
s.wait(20.0)

#LOAD ODT
s, ENDUVMOT = uvred.run(s,'ANDOR')
s.digichg('uvaom1',1)

#Insert ODT overlap with UVMOT
overlapdt = float(report['ODT']['overlapdt'])
s.wait(-overlapdt)
s.digichg('odtttl',1)
s.digichg('odt7595',1)
s.wait(overlapdt)

#~ #Leave UVMOT on for state transfer
#~ fstatedt  = float(report['ODT']['fstatedt'])
#~ s.wait(fstatedt)
#~ s.digichg('uvaom1',0)
#~ s.wait(-fstatedt) 

#Use MOT beams for fluorescence imaging
probe = 'motswitch'

#SEQUENCE

s=seq.sequence(stepsize)
s=gen.initial(s)

s.wait(20.0)


#LoadRamps refers to everything done up to loading the optical trap
#Edit loadtrap.py to change any of this
s, duration = uvred.run(s,'BASLER')


#RELEASE FROM MOT
s.digichg('motswitch',0) 
s.digichg('uvshutter',0)
s.digichg('field',0)

s.wait(tof)

#Take fluorescence imaging shot with the MOT beams. 
#PICTURE OF ATOMS
s=basler.BaslerPicture(s,preexp,texp,postexp,probe)
#PICTURE OF BACKGROUND
s.wait(noatoms)
s=basler.BaslerPicture(s,preexp,texp,postexp,probe)
noatoms = float(report['ANDOR']['noatoms'])

#SEQUENCE
s = seq.sequence(stepsize)
s = gen.initial(s)

s.digichg('odt7595', 0)

#Keep ODT on
ODT = gen.bstr('ODT', report)
if ODT == True:
    s.digichg('odtttl', 1)
s.wait(20.0)

#LOAD ODT
s, ENDUVMOT = uvred.run(s, 'ANDOR')
s.digichg('uvaom1', 1)

#Insert ODT overlap with UVMOT
overlapdt = float(report['ODT']['overlapdt'])
s.wait(-overlapdt)
s.digichg('odtttl', 1)
s.digichg('odt7595', 1)
s.wait(overlapdt)

#~ #Leave UVMOT on for state transfer
#~ fstatedt  = float(report['ODT']['fstatedt'])
#~ s.wait(fstatedt)
#~ s.digichg('uvaom1',0)
#~ s.wait(-fstatedt)
noatoms = float(report["ANDOR"]["noatoms"])

# SEQUENCE
s = seq.sequence(stepsize)
s = gen.initial(s)

s.digichg("odt7595", 0)

# Keep ODT on
ODT = gen.bstr("ODT", report)
if ODT == True:
    s.digichg("odtttl", 1)
s.wait(20.0)

# LOAD ODT
s, ENDUVMOT = uvred.run(s, "ANDOR")
s.digichg("uvaom1", 1)

# Insert ODT overlap with UVMOT
overlapdt = float(report["ODT"]["overlapdt"])
s.wait(-overlapdt)
s.digichg("odtttl", 1)
s.digichg("odt7595", 1)
s.wait(overlapdt)

# ~ #Leave UVMOT on for state transfer
# ~ fstatedt  = float(report['ODT']['fstatedt'])
# ~ s.wait(fstatedt)
# ~ s.digichg('uvaom1',0)
# ~ s.wait(-fstatedt)