source: Tutorials/PythonScript/data/sim.py @ 3209

'''Sample script to demonstrate use of GSASIIscriptable to simulate a powder pattern

This script is described in this tutorial: 
https://subversion.xray.aps.anl.gov/pyGSAS/Tutorials/PythonScript/Scripting.htm
'''

import os,sys
sys.path.insert(0,'/Users/toby/software/G2/GSASII')
import GSASIIscriptable as G2sc

workdir = "/Users/toby/Scratch/PythonScript"
datadir = "/Users/toby/software/G2/Tutorials/PythonScript/data"

gpx = G2sc.G2Project(filename='PbSO4sim.gpx') # create a project

# setup step 1: add a phase to the project
phase0 = gpx.add_phase(os.path.join(datadir,"PbSO4-Wyckoff.cif"),
                      phasename="PbSO4",fmthint='CIF')

# setup step 2: add a simulated histogram and link it to the previous phase(s)
hist1 = gpx.add_simulated_powder_histogram("PbSO4 simulation",
                          os.path.join(datadir,"inst_d1a.prm"),
                          5.,120.,0.01,
                          phases=gpx.phases())
gpx.data['Controls']['data']['max cyc'] = 0 # refinement not needed
gpx.do_refinements([{}])
gpx.save()