GEMC_histogram_analysis/plot_trajectory.py at master · dejac001/GEMC_histogram_analysis · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
from histogram.density_histogram import parse_arguments
from histogram.trajectory_histogram import Trajectory


def example_butane():
    args = parse_arguments()
    args['files'] = [
        'example_data/butane_400K/1/fort12.prod-%i' % i for i in range(1, 4)
    ]
    args['temperature'] = 400.
    args['plot'] = 'Yes'
    instance = Trajectory(**args)
    instance.main()
    return instance


def example_pentanal():
    args = parse_arguments()
    args['files'] = [
        'example_data/pentanal_510K/par1/fort12.prod%i' % i for i in range(1, 20)
    ]
    args['temperature'] = 510.
    args['plot'] = 'Yes'
    instance = Trajectory(**args)
    instance.main()
    return instance


def run_examples():
    example_pentanal()
    example_butane()


if __name__ == '__main__':
    # args = parse_arguments()
    # instance = Trajectory(**args)
    # instance.main()
    run_examples()