project
SUM_PARTNER_MAX_BACKUP = '_sum_partner_previous_max'
module-attribute
USER_CONSTRAINT_FLAG = '_easyreflectometry_user_constraint'
module-attribute
logger = logging.getLogger(__name__)
module-attribute
Q_MIN = 0.001
module-attribute
Q_MAX = 0.3
module-attribute
Q_RESOLUTION = 500
module-attribute
GUIDE_FIELD_ANGLE = 270.0
module-attribute
SPIN_ASYMMETRY_SIGNIFICANCE = 3.0
module-attribute
SPIN_ASYMMETRY_CANCELLATION_FRACTION = 0.001
module-attribute
MAGNETIC_MOMENT_FLOOR_FRACTION = 0.01
module-attribute
DEFAULT_MINIMIZER = AvailableMinimizers.LMFit_leastsq
module-attribute
_PATH_SKIPPED_PROPERTIES = frozenset({'interface', 'parent', 'front_layer', 'back_layer', 'head_layer', 'tail_layer'})
module-attribute
CalculatorFactory
Bases: InterfaceFactoryTemplate
Source code in src/easyreflectometry/calculators/factory.py
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 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 | |
fit_func
property
Fit func.
__init__()
Init function.
Source code in src/easyreflectometry/calculators/factory.py
15 16 17 | |
reset_storage()
Reset storage.
Source code in src/easyreflectometry/calculators/factory.py
19 20 21 | |
sld_profile(model_id)
Sld profile.
Source code in src/easyreflectometry/calculators/factory.py
23 24 25 | |
polarized_reflectivity_profiles(x_array, model_id)
Reflectivity profiles of all four spin channels ('pp', 'pm', 'mp', 'mm').
Source code in src/easyreflectometry/calculators/factory.py
27 28 29 | |
reflectivity_profile_channel(x_array, model_id, channel)
Reflectivity profile of one explicit spin channel ('pp', 'pm', 'mp' or 'mm').
Source code in src/easyreflectometry/calculators/factory.py
31 32 33 | |
magnetic_sld_profile(model_id)
Nuclear and magnetic sld profiles: z, sld(z), rhoM(z) and thetaM(z).
Source code in src/easyreflectometry/calculators/factory.py
35 36 37 | |
fit_func_for_channel(channel)
A fit function evaluating one explicit spin channel.
Used for simultaneous multi-channel fitting: each channel dataset gets its own fit function while all of them share the same model (and hence the same parameters).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
channel
|
PolarizationChannel | str
|
One of 'pp', 'pm', 'mp', 'mm' (or the corresponding enum member). |
required |
Returns:
| Type | Description |
|---|---|
Callable
|
Function of (x_array, model_id) returning the channel reflectivity. |
Source code in src/easyreflectometry/calculators/factory.py
61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 | |
PolarizationChannel
Bases: str, Enum
Spin cross-section channels for polarized neutron reflectometry.
The accepted spellings are exactly the enum values ('pp', 'pm', 'mp', 'mm'); uppercase strings are rejected.
Source code in src/easyreflectometry/calculators/polarization.py
7 8 9 10 11 12 13 14 15 16 17 | |
CalculatorBase
Bases: SerializerComponent
This class is a template and defines all properties that a calculator should have.
Source code in src/easyreflectometry/calculators/calculator_base.py
24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 | |
supports_magnetism
property
Whether this calculator backend can model magnetic samples.
include_magnetism
property
writable
Include magnetism.
polarization_channel
property
writable
The spin channel ('pp', 'pm', 'mp' or 'mm') used by reflectity_profile when magnetism is enabled.
Note: this state belongs to the currently-active calculator instance; switching calculators via the factory constructs a fresh instance and resets it.
__init_subclass__(is_abstract=False, **kwargs)
Initialise all subclasses so that they can be created in the factory.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
cls
|
|
required | |
is_abstract
|
bool
|
Is this a subclass which shouldn't be dded. By default, False. |
False
|
**kwargs
|
Key word arguments. |
{}
|
Source code in src/easyreflectometry/calculators/calculator_base.py
33 34 35 36 37 38 39 40 41 42 43 44 45 46 | |
__init__()
Init function.
Source code in src/easyreflectometry/calculators/calculator_base.py
48 49 50 51 | |
reset_storage()
Reset the storage area of the calculator.
Source code in src/easyreflectometry/calculators/calculator_base.py
53 54 55 | |
create(model)
Creation function.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
model
|
Material | Layer | Multilayer | Model
|
Object to be created. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 | |
assign_material_to_layer(material_id, layer_id)
Assign a material to a layer.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
material_id
|
str
|
The material name. |
required |
layer_id
|
str
|
The layer name. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
132 133 134 135 136 137 138 139 140 141 142 | |
add_layer_to_item(layer_id, item_id)
Add a layer to the item stack.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
item_id
|
str
|
The item id. |
required |
layer_id
|
str
|
The layer id. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
144 145 146 147 148 149 150 151 152 153 154 | |
remove_layer_from_item(layer_id, item_id)
Remove a layer from an item stack.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
item_id
|
str
|
The item id. |
required |
layer_id
|
str
|
The layer id. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
156 157 158 159 160 161 162 163 164 165 166 | |
add_item_to_model(item_id, model_id)
Add a layer to the item stack.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
item_id
|
str
|
The item id. |
required |
model_id
|
str
|
The model id. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
168 169 170 171 172 173 174 175 176 177 178 | |
remove_item_from_model(item_id, model_id)
Remove an item from the model.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
item_id
|
str
|
The item id. |
required |
model_id
|
str
|
The model id. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
180 181 182 183 184 185 186 187 188 189 190 | |
reflectity_profile(x_array, model_id)
Determines the reflectivity profile for the given range and model.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
x_array
|
np.ndarray
|
Points to be calculated at. |
required |
model_id
|
str
|
The model id. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
192 193 194 195 196 197 198 199 200 201 202 | |
reflectivity_profile_channel(x_array, model_id, channel)
Determine the reflectivity profile of one explicit spin channel.
Unlike polarization_channel (global calculator state), the channel is an
argument, so several channels can be evaluated against the same model —
one per dataset in a simultaneous multi-channel fit.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
x_array
|
np.ndarray
|
Points to be calculated at. |
required |
model_id
|
str
|
The model id. |
required |
channel
|
PolarizationChannel | str
|
One of 'pp', 'pm', 'mp', 'mm' (or the corresponding enum member). |
required |
Returns:
| Type | Description |
|---|---|
np.ndarray
|
Reflectivity of the requested channel at q. |
Source code in src/easyreflectometry/calculators/calculator_base.py
204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 | |
polarized_reflectivity_profiles(x_array, model_id)
Determines the reflectivity profiles of all four spin channels for the given range and model.
Requires include_magnetism to be enabled and a calculator that supports it (refl1d).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
x_array
|
np.ndarray
|
Points to be calculated at. |
required |
model_id
|
str
|
The model id. |
required |
Returns:
| Type | Description |
|---|---|
dict[str, np.ndarray]
|
Reflectivity per spin channel, keyed 'pp', 'pm', 'mp', 'mm' (in that order). |
Source code in src/easyreflectometry/calculators/calculator_base.py
227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 | |
sld_profile(model_id)
Return the scattering length density profile.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
model_id
|
str
|
The model id. |
required |
Returns:
| Type | Description |
|---|---|
tuple[np.ndarray, np.ndarray]
|
z and sld(z). |
Source code in src/easyreflectometry/calculators/calculator_base.py
246 247 248 249 250 251 252 253 254 255 256 257 258 259 | |
magnetic_sld_profile(model_id)
Return the nuclear and magnetic scattering length density profiles.
Requires include_magnetism to be enabled and a calculator that supports it (refl1d).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
model_id
|
str
|
The model id. |
required |
Returns:
| Type | Description |
|---|---|
tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray]
|
z, sld(z), magnetic sld rhoM(z) and magnetic angle thetaM(z). |
Source code in src/easyreflectometry/calculators/calculator_base.py
261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 | |
set_resolution_function(resolution_function)
Set resolution function.
Source code in src/easyreflectometry/calculators/calculator_base.py
278 279 280 | |
remove_layer_magnetism(layer_id)
Remove the magnetic state of one layer; disables include_magnetism
when no magnetic layer is left.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
layer_id
|
str
|
The layer id. |
required |
Source code in src/easyreflectometry/calculators/calculator_base.py
287 288 289 290 291 292 293 294 295 296 | |
DataSet1D
Bases: SerializerComponent
Source code in src/easyreflectometry/data/data_store.py
15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 | |
model
property
writable
Model function.
is_experiment
property
Is experiment.
is_simulation
property
Is simulation.
__init__(name='Series', x=None, y=None, ye=None, xe=None, model=None, x_label='x', y_label='y', auto_background=True)
Init function.
Source code in src/easyreflectometry/data/data_store.py
16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 | |
data_points()
Data points.
Source code in src/easyreflectometry/data/data_store.py
85 86 87 | |
__repr__()
Repr function.
Source code in src/easyreflectometry/data/data_store.py
89 90 91 | |
PolarizedDataSet
A polarized experiment: one dataset per measured spin channel, one shared model.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
name
|
str
|
Name of the experiment. By default, 'PolarizedSeries'. |
'PolarizedSeries'
|
channels
|
dict[PolarizationChannel | str, DataSet1D]
|
The measured channels; at least one. Keys are 'pp', 'pm', 'mp', 'mm' (or the corresponding enum members). Two-channel NSF-only experiments simply provide 'pp' and 'mm'. |
None
|
model
|
Model
|
The model shared by all channels. By default, None. |
None
|
Source code in src/easyreflectometry/data/polarized.py
88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 | |
model
property
writable
The model shared by all channels.
channels
property
The measured channels, in canonical order (pp, pm, mp, mm).
Read-only view: use :meth:set_channel / :meth:remove_channel to
modify the channel set, so validation and model propagation apply.
available_channels
property
The measured channels, in canonical order (pp, pm, mp, mm).
is_experiment
property
Is experiment.
is_simulation
property
Is simulation.
set_channel(channel, dataset)
Add or replace one channel dataset.
The dataset is validated, adopts the shared model, and the canonical channel order is preserved.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
channel
|
Union[PolarizationChannel, str]
|
One of 'pp', 'pm', 'mp', 'mm' (or the corresponding enum member). |
required |
dataset
|
DataSet1D
|
The channel data. |
required |
Source code in src/easyreflectometry/data/polarized.py
157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 | |
remove_channel(channel)
Remove one channel dataset; the last remaining channel cannot be removed.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
channel
|
Union[PolarizationChannel, str]
|
One of 'pp', 'pm', 'mp', 'mm' (or the corresponding enum member). |
required |
Source code in src/easyreflectometry/data/polarized.py
177 178 179 180 181 182 183 184 185 186 187 188 189 190 | |
MultiFitter
Source code in src/easyreflectometry/fitting.py
294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 | |
sampler
property
The Sampler behind the most recent :meth:mcmc_sample call, or None.
Holds the live BUMPS chain state, so the sampling run can be continued
with fitter.sampler.extend(additional_samples=...) instead of
starting a fresh chain.
chi2
property
Total chi-squared across all fitted datasets, or None if no fit has been performed.
reduced_chi
property
Reduced chi-squared from the most recent fit, or None if no fit has been performed.
classical_chi2
property
Classical chi-squared using only points with positive variances.
classical_reduced_chi
property
Reduced classical chi-squared using only points with positive variances.
objective_chi2
property
Objective-space chi-squared returned by the minimizer.
objective_reduced_chi
property
Objective-space reduced chi-squared returned by the minimizer.
__init__(*args, objective='hybrid')
A convenience class for the class:
easyscience.Fitting.Fitting
which will populate the class:
sc.DataGroup appropriately
after the fitting is performed.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
*args
|
Model
|
Reflectometry model(s). |
()
|
objective
|
str
|
Zero-variance handling strategy. One of
|
'hybrid'
|
Source code in src/easyreflectometry/fitting.py
295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 | |
for_experiments(experiments, objective='hybrid', constraints_factory_provider=None)
classmethod
Build a fitter for a mixed list of unpolarized and polarized experiments.
Every experiment contributes one fit function per dataset: an ordinary
experiment one, a polarized experiment one per measured spin channel,
each evaluating that channel's spin cross-section against the single
model the channels share. Structural parameters are therefore common to
all channels and the magnetic ones are constrained by all of them at
once, exactly as in :meth:fit_polarized — but here several experiments
(and several models) can be fitted together, which is what an
application's "fit everything that is loaded" action needs.
The resulting fitter is not run: the caller supplies the data arrays
to easy_science_multi_fitter.fit(...) in the order given by
:attr:fit_datasets, which lets a GUI drive it from a worker thread.
That call resolves :attr:constraints_factory_provider at call time,
so inequality constraints are applied on this path too — pass
constraints_factory_provider (e.g.
project.build_constraints_factory) or set the attribute before
fitting; with none set, no inequality constraints are enforced.
Note
Built via cls(*models, objective=objective) and then overwrites
_fit_func / easy_science_multi_fitter with the per-channel
versions — the __init__-built pair is briefly constructed and
discarded. Unlike :meth:fit_polarized, the minimizer selection and
its tolerance / max_evaluations are not carried over: the
returned fitter starts from easy_science_multi_fitter's defaults,
so a caller that needs a specific minimizer must set it explicitly
before calling .fit(...).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
experiments
|
list[DataSet1D | PolarizedDataSet]
|
The loaded experiments, in the order they should be fitted. |
required |
objective
|
str
|
Zero-variance handling strategy, see :meth: |
'hybrid'
|
constraints_factory_provider
|
Callable[[], Callable | None] | None
|
Zero-argument callable returning the |
None
|
Returns:
| Type | Description |
|---|---|
MultiFitter
|
Fitter whose |
Source code in src/easyreflectometry/fitting.py
399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 | |
fit(data, id=0, objective=None, constraints_factory=None)
Perform the fitting and populate the DataGroups with the result.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
sc.DataGroup
|
DataGroup to be fitted to and populated. |
required |
id
|
int
|
Unused parameter kept for backward compatibility. By default, 0. |
0
|
objective
|
str | None
|
Per-call override for the zero-variance objective.
If |
None
|
constraints_factory
|
Callable | None
|
Inequality constraints to enforce (BUMPS engines only); see
:mod: |
None
|
Returns:
| Type | Description |
|---|---|
sc.DataGroup
|
A new DataGroup with fitted model curves, SLD profiles, and fit statistics. |
Source code in src/easyreflectometry/fitting.py
494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 | |
fit_single_data_set_1d(data, objective=None, constraints_factory=None)
Perform fitting on a single 1D dataset.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
DataSet1D
|
The 1D dataset to fit. Note that |
required |
objective
|
str | None
|
Per-call override for the zero-variance objective.
If |
None
|
constraints_factory
|
Callable | None
|
Inequality constraints to enforce (BUMPS engines only). Defaults
to what :attr: |
None
|
Returns:
| Type | Description |
|---|---|
FitResults
|
Fit results from the minimizer. |
Source code in src/easyreflectometry/fitting.py
573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 | |
fit_polarized(data, objective=None, constraints_factory=None)
Fit all measured spin channels of a polarized experiment simultaneously.
Each channel dataset gets its own fit function evaluating the
corresponding spin cross-section, while every channel shares the single
model — so structural parameters (thickness, roughness, nuclear SLD,
scale, background) are common, and the magnetic parameters
(rho_m/theta_m) are constrained by all channels at once. The refl1d
backend computes all four cross-sections in one kernel evaluation and
caches them per iteration, so fitting N channels costs about as much as
fitting one.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data
|
PolarizedDataSet
|
The polarized experiment (its model must be the model this fitter
was constructed with). Note that per-channel |
required |
objective
|
str | None
|
Per-call override for the zero-variance objective.
If |
None
|
constraints_factory
|
Callable | None
|
Inequality constraints to enforce (BUMPS engines only); see
:mod: |
None
|
Returns:
| Type | Description |
|---|---|
dict[str, FitResults]
|
Fit results per channel, keyed 'pp', 'pm', 'mp', 'mm' (measured channels only, in that order). |
Note
Unlike :meth:for_experiments, this method does not populate
:attr:fit_datasets / :attr:fit_channels — those are set only by the
caller-driven, for_experiments-built flow.
Source code in src/easyreflectometry/fitting.py
619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 | |
mcmc_sample(data, samples=10000, burn=2000, thin=10, population=None, objective=None, initializer=None, progress_callback=None, abort_test=None, constraints_factory=None)
Run Bayesian MCMC sampling on reflectometry data using the DREAM sampler.
Requires that the minimizer is a BUMPS instance (i.e. the minimizer was
switched to AvailableMinimizers.Bumps).
:param data: DataGroup with reflectivity data.
:param samples: Number of retained DREAM samples requested from BUMPS.
:param burn: Burn-in steps.
:param thin: Thinning interval.
:param population: BUMPS DREAM population count for advanced users.
:param objective: Zero-variance handling strategy.
:param initializer: DREAM population initializer. One of 'eps',
'cov', 'lhs', or 'random'. By default, None (BUMPS
uses 'eps').
:param progress_callback: Optional callback for progress updates during
sampling. Forwarded to the core MultiFitter.
:param abort_test: Optional callable returning True to abort sampling.
:param constraints_factory: Inequality constraints to enforce (see
:mod:easyreflectometry.inequality_constraints); the posterior is
penalised in the infeasible region. Defaults to what
:attr:constraints_factory_provider returns.
:return: Dictionary with keys 'draws', 'param_names', 'state',
and 'logp'.
:raises RuntimeError: If the current minimizer is not a BUMPS instance.
The underlying :class:~easyscience.fitting.Sampler is retained on
:attr:sampler, so the chain can be continued without re-running the
burn-in::
fitter.mcmc_sample(data, samples=2000, burn=500, thin=10)
extended = fitter.sampler.extend(additional_samples=8000, thin=10)
Source code in src/easyreflectometry/fitting.py
721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 | |
record_fit_results(results)
Adopt fit results produced elsewhere, so this fitter reports on them.
An application that drives easy_science_multi_fitter.fit(...) itself
— to run it in a worker thread, for instance — leaves the MultiFitter
that owns the goodness-of-fit properties none the wiser. Handing the
results back here makes :attr:chi2 / :attr:reduced_chi describe that
fit instead of reporting "no fit performed".
Only the minimizer-reported metrics are restored: the classical ones
need the original data arrays, which are not part of FitResults, so
:attr:classical_chi2 and :attr:classical_reduced_chi stay None.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
results
|
list[FitResults] | None
|
Results of the fit, one per fitted dataset. None clears them. |
required |
Source code in src/easyreflectometry/fitting.py
903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 | |
switch_minimizer(minimizer)
Switch the minimizer for the fitting.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
minimizer
|
AvailableMinimizers
|
Minimizer to be switched to. |
required |
Source code in src/easyreflectometry/fitting.py
924 925 926 927 928 929 930 931 932 | |
InequalityEvaluation
dataclass
Result of evaluating a spec against the current parameter values.
Source code in src/easyreflectometry/inequality_constraints.py
242 243 244 245 246 247 248 249 250 251 252 253 | |
InequalitySpec
dataclass
Declarative cross-parameter inequality, e.g. t_head < t_tail.
Attributes:
| Name | Type | Description |
|---|---|---|
lhs_expression |
str
|
Expression over the aliases in |
op |
str
|
One of |
rhs_expression |
str
|
Expression over the aliases in |
lhs_paths, rhs_paths |
dict[str, str]
|
Alias to structural parameter path ( |
name |
str
|
Optional user label. |
enabled |
bool
|
Disabled specs are kept but not applied to fits. |
Source code in src/easyreflectometry/inequality_constraints.py
136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 | |
paths
property
All alias → path pairs of both sides.
validate_syntax()
Check that both sides parse and every identifier has an alias mapping.
Source code in src/easyreflectometry/inequality_constraints.py
176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 | |
Model
Bases: BaseCore
Model is the class that represents the experiment.
It is used to store the information about the experiment and to perform the calculations.
Source code in src/easyreflectometry/model/model.py
60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 | |
total_thickness
property
Read-only parameter: the summed thickness of the film layers.
The superphase (first layer of the sample) and the subphase (last layer) are semi-infinite and excluded. The parameter is dependent: it tracks the layer thicknesses, never enters a fit and cannot be set; it can be referenced from constraint expressions and inequality constraints like any other parameter.
has_magnetism
property
Whether any layer in the sample carries magnetic properties.
is_default
property
writable
Whether this model was created as a default placeholder.
resolution_function
property
writable
Return the resolution function.
__init__(sample=None, scale=None, background=None, resolution_function=None, name='Model', color=COLORS[0], unique_name=None, interface=None)
Constructor.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
unique_name
|
Optional[str]
|
By default, None. |
None
|
color
|
str
|
By default, COLORS[0]. |
COLORS[0]
|
sample
|
Union[Sample, None]
|
The sample being modelled. By default, None. |
None
|
scale
|
Union[Parameter, Number, None]
|
Scaling factor of profile. By default, None. |
None
|
background
|
Union[Parameter, Number, None]
|
Linear background magnitude. By default, None. |
None
|
name
|
str
|
Name of the model. By default, 'Model'. |
'Model'
|
resolution_function
|
Union[ResolutionFunction, None]
|
Resolution function. By default, None. |
None
|
interface
|
Calculator interface. By default, None. |
None
|
Source code in src/easyreflectometry/model/model.py
66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 | |
interface(new_interface)
Set the interface; runs generate_bindings and then refreshes the
calculator's resolution function.
Source code in src/easyreflectometry/model/model.py
294 295 296 297 298 299 300 301 302 | |
refresh_derived()
Re-derive the model's computed parameters from the current layer structure.
Cheap (a few string operations) and idempotent: the dependency is
only rebuilt when the set of contributing layers changed. Called on
every access to :attr:total_thickness, so layers added or removed
through any path (Model.add_assemblies, assembly.layers.append,
...) are picked up without explicit notification.
Source code in src/easyreflectometry/model/model.py
170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 | |
add_assemblies(*assemblies)
Add assemblies to the model sample.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
*assemblies
|
list[BaseAssembly]
|
Assemblies to add to model sample. |
()
|
Source code in src/easyreflectometry/model/model.py
215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 | |
duplicate_assembly(index)
Duplicate a given item or layer in a sample.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
|
required |
idx
|
Index of the item or layer to duplicate. |
required |
Source code in src/easyreflectometry/model/model.py
236 237 238 239 240 241 242 243 244 245 246 247 | |
remove_assembly(index)
Remove an assembly from the model.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
|
required |
idx
|
Index of the item to remove. |
required |
Source code in src/easyreflectometry/model/model.py
249 250 251 252 253 254 255 256 257 258 259 260 261 | |
to_dict(skip=None)
Serialize the model, encoding the resolution function and interface name.
Source code in src/easyreflectometry/model/model.py
327 328 329 330 331 332 333 334 335 336 337 338 339 340 | |
as_dict(skip=None)
Compatibility alias for :meth:to_dict.
Source code in src/easyreflectometry/model/model.py
342 343 344 | |
as_orso()
The sample as an ORSO simple-model (sample.model) dictionary.
Slab representation: lengths in angstrom, SLDs in 1/angstrom^2,
repeating multilayers via the inline N ( ... ) stack syntax.
Returns:
| Type | Description |
|---|---|
dict
|
The ORSO model-language dictionary (the content of an .ort file's
|
Source code in src/easyreflectometry/model/model.py
346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 | |
from_dict(passed_dict)
classmethod
Create a Model from a dictionary.
Source code in src/easyreflectometry/model/model.py
365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 | |
ModelCollection
Bases: BaseCollection
Source code in src/easyreflectometry/model/model_collection.py
21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 | |
next_color_index
property
Index of the next colour to assign — kept around so it round-trips.
__init__(*models, name='Models', interface=None, unique_name=None, populate_if_none=True, next_color_index=None, **kwargs)
Init function.
Source code in src/easyreflectometry/model/model_collection.py
22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 | |
next_color()
Colour the next appended model should get.
Appending advances the cycle, so callers that build a Model
themselves can keep the per-model colours distinct::
model = Model(sample=sample, color=collection.next_color())
collection.add_model(model)
Source code in src/easyreflectometry/model/model_collection.py
66 67 68 69 70 71 72 73 74 75 | |
add_model(model=None)
Add a model to the collection.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
model
|
Optional[Model]
|
Model to add. By default, None (a new model is created with
the collection's next colour; a supplied model keeps its own
colour — use :meth: |
None
|
Source code in src/easyreflectometry/model/model_collection.py
77 78 79 80 81 82 83 84 85 86 87 88 89 | |
duplicate_model(index)
Duplicate a model in the collection.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Model to duplicate. |
required |
Source code in src/easyreflectometry/model/model_collection.py
91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 | |
from_dict(this_dict)
classmethod
Create an instance of a collection from a dictionary.
Source code in src/easyreflectometry/model/model_collection.py
107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 | |
append(model)
Append function.
Source code in src/easyreflectometry/model/model_collection.py
136 137 138 | |
PercentageFwhm
Bases: ResolutionFunction
Source code in src/easyreflectometry/model/resolution_functions.py
59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 | |
__init__(constant=None)
Init function.
Source code in src/easyreflectometry/model/resolution_functions.py
60 61 62 63 64 | |
smearing(q)
Return per-point sigma values from the constant FWHM percentage.
constant is a FWHM percentage of q; it is converted to an
absolute sigma so the smearing() contract is sigma for all types.
Source code in src/easyreflectometry/model/resolution_functions.py
66 67 68 69 70 71 72 73 74 | |
as_dict(skip=None)
As dict.
Source code in src/easyreflectometry/model/resolution_functions.py
76 77 78 79 80 | |
Pointwise
Bases: ResolutionFunction
Pointwise resolution defined by a per-point resolution provided with the data.
The resolution is supplied as the variance of the Qz values (sQz) at the
measured Qz data points, which is the form produced by data reduction (e.g.
Qz_0.variances). The resolution width at each point is sqrt(sQz).
For a requested q the width is obtained by linearly interpolating onto
q, exactly as :class:LinearSpline does for explicitly provided widths.
This is a convenience wrapper around :class:LinearSpline that derives the
widths from the [Qz, R, sQz] triple loaded from a data file; the returned
widths are consumed by the calculators (refnx x_err / refl1d dq),
which perform the actual convolution against the model.
Serialization contract: as_dict/from_dict store sQz_data_points
as variances (sigma squared). This is deliberately unchanged by the
ORSO value_is: FWHM support — FWHM columns are converted to sigma at
load time, so stored values are always sigma squared and saved projects
round-trip without migration.
Source code in src/easyreflectometry/model/resolution_functions.py
109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 | |
__init__(q_data_points)
Init function.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
q_data_points
|
List[np.ndarray]
|
|
required |
Source code in src/easyreflectometry/model/resolution_functions.py
130 131 132 133 134 135 136 137 138 139 140 | |
smearing(q=None)
Return the resolution sigma interpolated onto q.
sQz is the variance of Qz, so the sigma at each data point is
sqrt(sQz); values are linearly interpolated onto the requested
q. This already satisfies the sigma smearing() contract, so no
FWHM conversion is applied. When q is None the sigma values
are returned at the stored data points.
Source code in src/easyreflectometry/model/resolution_functions.py
142 143 144 145 146 147 148 149 150 151 152 153 154 155 | |
as_dict(skip=None)
As dict.
Source code in src/easyreflectometry/model/resolution_functions.py
157 158 159 160 161 162 163 164 165 166 | |
Layer
Bases: BaseCore
Source code in src/easyreflectometry/sample/elements/layers/layer.py
40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 | |
__init__(material=None, thickness=None, roughness=None, name='EasyLayer', unique_name=None, interface=None, magnetism=None)
Constructor.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
unique_name
|
Optional[str]
|
By default, None. |
None
|
material
|
Union[Material, None]
|
The material for the layer. By default, None. |
None
|
thickness
|
Union[Parameter, float, None]
|
Layer thickness in Angstrom. By default, None. |
None
|
roughness
|
Union[Parameter, float, None]
|
Upper roughness on the layer in Angstrom. By default, None. |
None
|
magnetism
|
Union[LayerMagnetism, None]
|
Magnetic properties of the layer; None for a non-magnetic layer. By default, None. |
None
|
name
|
str
|
Name of the layer. By default, 'EasyLayer'. |
'EasyLayer'
|
interface
|
Interface object. By default, None. |
None
|
Source code in src/easyreflectometry/sample/elements/layers/layer.py
41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 | |
interface(new_interface)
Set the interface; runs generate_bindings and, for a magnetic layer,
enables magnetism on the calculator (raising if it does not support it).
Source code in src/easyreflectometry/sample/elements/layers/layer.py
106 107 108 109 110 111 112 113 | |
assign_material(material)
Assign a material to the layer interface.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
material
|
Material
|
The material to assign to the layer. |
required |
Source code in src/easyreflectometry/sample/elements/layers/layer.py
175 176 177 178 179 180 181 182 183 184 185 | |
Material
Bases: BaseCore
Source code in src/easyreflectometry/sample/elements/materials/material.py
39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 | |
__init__(sld=None, isld=None, name='EasyMaterial', unique_name=None, interface=None)
Constructor.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
unique_name
|
Optional[str]
|
By default, None. |
None
|
sld
|
Union[Parameter, float, None]
|
Real scattering length density. By default, None. |
None
|
isld
|
Union[Parameter, float, None]
|
Imaginary scattering length density. By default, None. |
None
|
name
|
str
|
Name of the material. By default, 'EasyMaterial'. |
'EasyMaterial'
|
interface
|
Calculator interface. By default, None. |
None
|
Source code in src/easyreflectometry/sample/elements/materials/material.py
40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 | |
MaterialCollection
Bases: BaseCollection
Source code in src/easyreflectometry/sample/collections/material_collection.py
22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 | |
__init__(*materials, name='EasyMaterials', interface=None, unique_name=None, populate_if_none=True, **kwargs)
Init function.
Source code in src/easyreflectometry/sample/collections/material_collection.py
23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 | |
add_material(material=None)
Add a material to the collection.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
material
|
Optional[Material]
|
Material to add. By default, None. |
None
|
Source code in src/easyreflectometry/sample/collections/material_collection.py
48 49 50 51 52 53 54 55 56 57 58 59 | |
duplicate_material(index)
Duplicate a material in the collection.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
|
required |
material
|
Assembly to add. |
required |
Source code in src/easyreflectometry/sample/collections/material_collection.py
61 62 63 64 65 66 67 68 69 70 71 72 73 | |
Multilayer
Bases: BaseAssembly
A multi layer is build from a single or a list of Layer or LayerCollection.
The multi layer will arrange the layers as slabs, allowing the reflectometry to be determined from them. The front layer is where the neutron beam starts in, it has an index of 0.
More information about the usage of this assembly is available in the
multilayer documentation_
.. _multilayer documentation: ../sample/assemblies_library.html#multilayer
Source code in src/easyreflectometry/sample/assemblies/multilayer.py
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 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 | |
__init__(layers=None, name='EasyMultilayer', unique_name=None, interface=None, type='Multi-layer', populate_if_none=True, conformal_thickness=False, conformal_roughness=False)
Constructor.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
populate_if_none
|
Optional[bool]
|
By default, True. |
True
|
unique_name
|
Optional[str]
|
By default, None. |
None
|
layers
|
Union[Layer, list[Layer], LayerCollection, None]
|
The layers that make up the multi-layer. By default, None. |
None
|
name
|
str
|
Name for multi layer. By default, 'EasyMultilayer'. |
'EasyMultilayer'
|
interface
|
Calculator interface. By default, None. |
None
|
|
type
|
str
|
Type of the constructed instance. By default, 'Multi-layer'. |
'Multi-layer'
|
conformal_thickness
|
bool
|
Tie every layer's thickness to the front layer's. Serialization
reads the current graph state through the matching property, so
the ties are rebuilt on |
False
|
conformal_roughness
|
bool
|
Tie every layer's roughness to the front layer's, likewise persistent. By default, False. |
False
|
Source code in src/easyreflectometry/sample/assemblies/multilayer.py
26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 | |
add_layer(*layers)
Add a layer to the multi layer.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
*layers
|
tuple[Layer]
|
Layers to add to the multi layer. |
()
|
Source code in src/easyreflectometry/sample/assemblies/multilayer.py
86 87 88 89 90 91 92 93 94 95 96 97 98 | |
duplicate_layer(idx)
Duplicate a given layer.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
idx
|
int
|
Index of layer to duplicate. |
required |
Source code in src/easyreflectometry/sample/assemblies/multilayer.py
100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 | |
remove_layer(idx)
Remove a layer from the item.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
idx
|
int
|
Index of layer to remove. |
required |
Source code in src/easyreflectometry/sample/assemblies/multilayer.py
117 118 119 120 121 122 123 124 125 126 127 | |
from_dict(data)
classmethod
Create a Multilayer from a dictionary.
Source code in src/easyreflectometry/sample/assemblies/multilayer.py
135 136 137 138 139 | |
Sample
Bases: BaseCollection
A sample is a collection of assemblies that represent the structure for which experimental measurements exist.
Source code in src/easyreflectometry/sample/collections/sample.py
26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 | |
superphase
property
The superphase of the sample.
subphase
property
The subphase of the sample.
__init__(*assemblies, name='EasySample', interface=None, unique_name=None, populate_if_none=True, **kwargs)
Constructor.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
**kwargs
|
|
{}
|
|
populate_if_none
|
bool
|
By default, True. |
True
|
unique_name
|
Optional[str]
|
By default, None. |
None
|
*assemblies
|
Optional[List[BaseAssembly]]
|
|
()
|
args
|
The assemblies in the sample. |
required | |
name
|
str
|
Name of the sample. By default, 'EasySample'. |
'EasySample'
|
interface
|
Calculator interface. By default, None. |
None
|
Source code in src/easyreflectometry/sample/collections/sample.py
29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 | |
add_assembly(assembly=None)
Add an assembly to the sample.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
assembly
|
Optional[BaseAssembly]
|
Assembly to add. By default, None. |
None
|
Source code in src/easyreflectometry/sample/collections/sample.py
78 79 80 81 82 83 84 85 86 87 88 89 90 91 | |
duplicate_assembly(index)
Add an assembly to the sample.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
|
required |
assembly
|
Assembly to add. |
required |
Source code in src/easyreflectometry/sample/collections/sample.py
93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 | |
move_up(index)
Move the assembly at the given index up in the sample.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the assembly to move up. |
required |
Source code in src/easyreflectometry/sample/collections/sample.py
118 119 120 121 122 123 124 125 126 | |
move_down(index)
Move the assembly at the given index down in the sample.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the assembly to move down. |
required |
Source code in src/easyreflectometry/sample/collections/sample.py
128 129 130 131 132 133 134 135 136 | |
remove_assembly(index)
Remove the assembly at the given index from the sample.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the assembly to remove. |
required |
Source code in src/easyreflectometry/sample/collections/sample.py
138 139 140 141 142 143 144 145 146 | |
BaseCollection
Bases: EasyList
Local base for sample-tree collections (Material/Layer/Assembly/Model collections).
Built on top of easyscience.base_classes.EasyList (the replacement for
the deprecated CollectionBase). On top of EasyList this class adds:
- a
nameproperty - an
interfaceproperty whose setter propagates the calculator interface to every contained item - propagation of the current
interfaceto newly inserted items - a
populate_if_noneflag preserved for serialization round-trip - convenience helpers
names,move_up,move_down,remove_at - a yaml-formatted
__repr__driven by_dict_repr - an
as_dictalias forto_dictwithskip=support andinterfaceexcluded by default
Subclasses (LayerCollection, MaterialCollection, Sample,
ModelCollection) keep their existing constructor shape — they pass
name and interface positionally to this class, items as *args, and
additional configuration as kwargs.
Source code in src/easyreflectometry/sample/collections/base_collection.py
17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 | |
names
property
List of item names.
data
property
Read-only view of the underlying item list.
Provided for compatibility with code (and tests) written against the
legacy CollectionBase shape, which exposed items via .data.
insert(index, value)
Insert and (if an interface is set) propagate it to the new item.
Legacy CollectionBase.insert did value.interface = self.interface
after registering the item; we replicate the same behaviour here so
downstream calculator state stays in sync when items are appended
after the collection's interface was already set.
The type check from EasyList.insert is bypassed: each subclass
accepts a single item type (Layer / Material / BaseAssembly / Model)
and enforces it elsewhere, while the EasyList check would require
every item to be a NewBase subclass — which was historically not
guaranteed and forces an extra coupling we don't need.
Source code in src/easyreflectometry/sample/collections/base_collection.py
124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 | |
move_up(index)
Move the element at the given index up in the collection.
Source code in src/easyreflectometry/sample/collections/base_collection.py
167 168 169 170 171 | |
move_down(index)
Move the element at the given index down in the collection.
Source code in src/easyreflectometry/sample/collections/base_collection.py
173 174 175 176 177 | |
remove_at(index)
Remove the item at index from the collection.
Renamed from the legacy BaseCollection.remove(index) which shadowed
MutableSequence.remove(value) (remove-by-value, inherited from
EasyList). Callers that meant "remove by index" should use this; the
standard remove(value) is still available with its usual semantics.
Source code in src/easyreflectometry/sample/collections/base_collection.py
179 180 181 182 183 184 185 186 187 | |
get_parameters()
Compatibility alias for legacy callers; prefer get_all_parameters.
Source code in src/easyreflectometry/sample/collections/base_collection.py
191 192 193 | |
get_all_variables()
Flat list of every Parameter/Descriptor across all items.
Walks each item in the collection and unions whatever each item
exposes via its own get_all_variables (for ModelBase /
BaseCore children) or, for objects that lack that hook,
unions any direct DescriptorBase attributes.
Source code in src/easyreflectometry/sample/collections/base_collection.py
195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 | |
get_fit_parameters()
Alias kept for the minimizer; matches ModelBase.get_fit_parameters.
Source code in src/easyreflectometry/sample/collections/base_collection.py
225 226 227 | |
to_dict(skip=None)
Serialize with skip support; interface excluded by default.
EasyList.to_dict doesn't accept a skip argument and is hard-wired
to dump data plus the parent's _arg_spec view. We reimplement here
so existing callers (Project, Model.as_dict, etc.) can keep
passing skip=['unique_name'] or similar.
NewBase.to_dict mutates the skip list in-place (e.g. it
appends 'unique_name' when the collection's own unique_name is
default-generated). We therefore pass a copy to it, otherwise the
mutation would leak into the per-item serialization below and force
every item Parameter dict to drop its unique_name — breaking the
from_dict round-trip.
Source code in src/easyreflectometry/sample/collections/base_collection.py
279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 | |
as_dict(skip=None)
Compatibility alias for :meth:to_dict.
Source code in src/easyreflectometry/sample/collections/base_collection.py
319 320 321 | |
__deepcopy__(memo)
Round-trip via dict-skip-unique to get a fresh copy.
NewBase.__copy__ already does this; the override is kept (rather
than deleted) to mirror the legacy BaseCollection.__deepcopy__
semantics — callers that relied on copy.deepcopy(collection) still
get a clone built from from_dict(as_dict(skip=['unique_name'])).
Source code in src/easyreflectometry/sample/collections/base_collection.py
323 324 325 326 327 328 329 330 331 | |
Project
Source code in src/easyreflectometry/project.py
113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1867 1868 1869 1870 1871 1872 1873 1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945 1946 1947 1948 1949 1950 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124 2125 2126 2127 2128 2129 2130 | |
parameters
property
Get all unique parameters from all models in the project.
Parameters shared across multiple models (e.g. material SLD) are only included once to avoid double-counting.
q_min
property
writable
Q min.
q_max
property
writable
Q max.
q_resolution
property
writable
Q resolution.
current_material_index
property
writable
Current material index.
current_model_index
property
writable
Current model index.
current_assembly_index
property
writable
Current assembly index.
current_layer_index
property
writable
Current layer index.
current_experiment_index
property
writable
Current experiment index.
created
property
Created function.
path
property
Path function.
models
property
writable
Models function.
fitter
property
Fitter function.
inequality_constraints
property
The project's inequality constraints (a copy of the list).
calculator
property
writable
Calculator function.
calculator_supports_magnetism
property
Whether the active calculator can model magnetic samples.
The GUI uses this to gate magnetism-related controls (e.g. when the refnx or bornagain backend is selected).
calculators_supporting_magnetism
property
Names of the available calculators that can model magnetic samples.
Lets an application offer the switch a magnetic sample needs ("this requires refl1d — change to it?") instead of only reporting that the current calculator cannot do it. The active calculator is not touched.
models_have_magnetism
property
Whether any model in the project carries a magnetic layer.
A calculator without magnetism cannot be selected while this holds — the binding it would have to build does not exist.
minimizer
property
writable
Minimizer function.
experiments
property
writable
Experiments function.
path_json
property
Path json.
__init__()
Init function.
Source code in src/easyreflectometry/project.py
114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 | |
reset()
Reset function.
Source code in src/easyreflectometry/project.py
141 142 143 144 145 146 147 | |
set_path_project_parent(path)
Set path project parent.
Source code in src/easyreflectometry/project.py
335 336 337 | |
parameter_path(parameter)
Structural path of parameter within this project, e.g.
models/0/sample/1/layers/0/thickness.
Paths are stable across save/load (unlike unique names, which are
regenerated) and are the way inequality constraints reference
parameters. Returns None when the parameter is not reachable
from the project's models.
Source code in src/easyreflectometry/project.py
424 425 426 427 428 429 430 431 432 433 | |
resolve_parameter_path(path)
Return the parameter at a structural path (see :meth:parameter_path).
Source code in src/easyreflectometry/project.py
435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 | |
add_inequality_constraint(spec, validate=True)
Register an inequality constraint; returns it.
With validate the paths are resolved and the units of both sides
compared, raising KeyError/ValueError on problems.
Source code in src/easyreflectometry/project.py
462 463 464 465 466 467 468 469 470 471 472 473 | |
remove_inequality_constraint(which)
Remove a constraint by index, by name or by identity.
A name removes every spec carrying that name; use the index or the spec object to remove a single one when names are shared.
Source code in src/easyreflectometry/project.py
475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 | |
evaluate_inequality_constraints()
Evaluate every constraint (enabled or not) at the current values.
Source code in src/easyreflectometry/project.py
496 497 498 | |
violated_inequality_constraints()
Enabled constraints that the current parameter values violate.
A fit started from an infeasible point begins on the BUMPS penalty plateau; callers should refuse or warn before launching.
Source code in src/easyreflectometry/project.py
500 501 502 503 504 505 506 507 508 509 510 | |
build_constraints_factory()
constraints_factory hook for the enabled inequality constraints, or None.
Source code in src/easyreflectometry/project.py
512 513 514 | |
get_index_air()
Get index air.
Source code in src/easyreflectometry/project.py
759 760 761 762 763 | |
get_index_si()
Get index si.
Source code in src/easyreflectometry/project.py
765 766 767 768 769 | |
get_index_sio2()
Get index sio2.
Source code in src/easyreflectometry/project.py
771 772 773 774 775 | |
get_index_d2o()
Get index d2o.
Source code in src/easyreflectometry/project.py
777 778 779 780 781 | |
load_orso_file(path)
Load an ORSO file, creating a model from its sample.model (when present) and an experiment.
.. deprecated::
Use :meth:load_new_experiment (data) together with
:meth:set_sample_from_orso /
:func:easyreflectometry.orso_utils.load_orso_model (model)
instead. This wrapper now routes through the same
DataSet1D + title + resolution path as every other importer.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
Union[Path, str]
|
Path to the ORSO file. |
required |
Source code in src/easyreflectometry/project.py
783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 | |
set_sample_from_orso(sample)
Replace the current project model collection with a single model built from an ORSO-parsed sample.
This is a convenience helper for the ORSO import pipeline where a complete
:class:~easyreflectometry.sample.Sample is constructed elsewhere.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
sample
|
Sample
|
Sample to set as the project's (single) model. |
required |
Returns:
| Type | Description |
|---|---|
None
|
|
Source code in src/easyreflectometry/project.py
815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 | |
add_sample_from_orso(sample)
Add a new model with the given sample to the existing model collection.
The created model is appended to :attr:models, its calculator interface is
set to the project's current calculator, and any materials referenced in the
sample are added to the project's material collection.
After adding the model, :attr:current_model_index is updated to point to
the newly added model.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
sample
|
Sample
|
Sample to add as a new model. |
required |
Returns:
| Type | Description |
|---|---|
None
|
|
Source code in src/easyreflectometry/project.py
834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 | |
replace_models_from_orso(sample)
Replace all models and materials with a single model from an ORSO sample.
All existing models and their associated materials are removed. A new model is created from sample, assigned to the project's calculator, and the material collection is rebuilt from the new model only.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
sample
|
Sample
|
Sample to set as the project's only model. |
required |
Returns:
| Type | Description |
|---|---|
None
|
|
Source code in src/easyreflectometry/project.py
869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 | |
load_new_experiment(path, data_group=None)
Load new experiment.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
Union[Path, str]
|
Path to the experiment data file. |
required |
data_group
|
Pre-loaded scipp DataGroup for path (avoids re-parsing the file). By default, None. |
None
|
Source code in src/easyreflectometry/project.py
975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 | |
count_datasets_in_file(path)
Return the number of datasets contained in the file at path.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
Union[Path, str]
|
Path to the data file. |
required |
Returns:
| Type | Description |
|---|---|
int
|
Number of datasets found; 1 if a non-ORSO file cannot be introspected. A corrupt ORSO file (banner present but unparsable) raises instead of being silently miscounted. |
Source code in src/easyreflectometry/project.py
1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 | |
load_all_experiments_from_file(path)
Load all datasets from a file as separate experiments sharing the current model.
For a multi-dataset ORSO file (e.g. a multi-angle measurement), each dataset is
registered as an independent experiment. All experiments share the model that is
currently selected. Single-dataset files go through
:meth:load_new_experiment. A corrupt ORSO file raises (no silent
plain-text fallback).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
Union[Path, str]
|
Path to the data file. |
required |
Returns:
| Type | Description |
|---|---|
int
|
Number of experiments that were added. |
Source code in src/easyreflectometry/project.py
1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 | |
suggest_polarized_channel_assignment(paths)
Suggest a spin-channel assignment for a set of data files.
For each file, the ORSO header polarization is used when present, falling back to filename heuristics ('_uu'/'_up' → pp, '_dd'/'_down' → mm, ...). Files that cannot be identified map to None; a GUI should present the result as an editable file → channel table.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
paths
|
List[Union[Path, str]]
|
Paths of the per-channel data files. |
required |
Returns:
| Type | Description |
|---|---|
Dict[str, Optional[PolarizationChannel]]
|
Detected channel (or None) per path. |
Source code in src/easyreflectometry/project.py
1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 | |
load_polarized_experiment(paths, model_index=None)
Load a polarized experiment from one data file per spin channel.
The channel datasets form a single :class:PolarizedDataSet experiment
sharing one model. Two-channel NSF-only experiments simply pass 'pp' and
'mm' entries; spin-flip channels are optional.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
paths
|
Dict[Union[PolarizationChannel, str], Union[Path, str]]
|
Explicit channel → file assignment (e.g. from the GUI import dialog,
pre-filled via :meth: |
required |
model_index
|
Optional[int]
|
Index of the model the experiment belongs to. By default, the current model. |
None
|
Returns:
| Type | Description |
|---|---|
int
|
Index of the newly loaded experiment, so a caller can make it current. |
Source code in src/easyreflectometry/project.py
1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 | |
load_polarized_experiment_from_file(path, model_index=None)
Load a polarized experiment from a single multi-dataset ORSO file.
Multi-dataset packing is the ORSO format's intended way to store spin
states: each data_set: block is classified by its own header
(instrument_settings.polarization), so per-dataset overrides are
honoured. Only the fully-analysed cross-sections pp/pm/mp/mm are
mapped; po/mo/op/om/unpolarized/vector are not coerced.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
Union[Path, str]
|
Path to the multi-dataset ORSO file. |
required |
model_index
|
Optional[int]
|
Index of the model the experiment belongs to. By default, the current model. |
None
|
Returns:
| Type | Description |
|---|---|
int
|
Index of the newly loaded experiment. |
Raises:
| Type | Description |
|---|---|
ValueError :
|
If any dataset lacks a mappable polarization header, or two
datasets declare the same channel. Such files cannot be classified
unambiguously; use :meth: |
Source code in src/easyreflectometry/project.py
1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 | |
save_experiment_as_orso(path, index=None)
Write an experiment to an ORSO file (.ort text, or .orb binary).
The experiment's model (when set) is exported as the ORSO
sample.model (slab representation); data columns are written as
sigma. A polarized experiment becomes one file with one data_set:
block per spin channel.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
Union[Path, str]
|
Destination path; a |
required |
index
|
Optional[int]
|
Index of the experiment to save. By default, the current one. |
None
|
Source code in src/easyreflectometry/project.py
1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 | |
load_experiment_for_model_at_index(path, index=0)
Load experiment for model at index.
Source code in src/easyreflectometry/project.py
1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 | |
sld_data_for_model_at_index(index=0)
Sld data for model at index.
Source code in src/easyreflectometry/project.py
1300 1301 1302 1303 1304 1305 1306 1307 1308 | |
model_has_magnetism_at_index(index=0)
Whether the model at index carries magnetism (False when there is no such model).
Source code in src/easyreflectometry/project.py
1310 1311 1312 1313 1314 1315 | |
magnetic_sld_data_for_model_at_index(index=0)
Nuclear and magnetic depth profiles of a magnetic model.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the model. |
0
|
Returns:
| Type | Description |
|---|---|
Dict[str, DataSet1D]
|
Profiles versus depth z, keyed:
|
Raises:
| Type | Description |
|---|---|
ValueError
|
The model has no magnetic layer, so there is no magnetic profile. |
NotImplementedError
|
The active calculator cannot model magnetism. |
Source code in src/easyreflectometry/project.py
1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 | |
sample_data_for_model_at_index(index=0, q_range=None)
Sample data for model at index.
Source code in src/easyreflectometry/project.py
1378 1379 1380 1381 1382 1383 1384 1385 | |
model_data_for_model_at_index(index=0, q_range=None, channel=None)
Model data for model at index.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the model. |
0
|
q_range
|
Optional[np.array]
|
Points to calculate at; the project q range by default. |
None
|
channel
|
Optional[Union[PolarizationChannel, str]]
|
Spin cross-section to calculate. |
None
|
Source code in src/easyreflectometry/project.py
1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 | |
experimental_data_for_model_at_index(index=0, channel=None)
Experimental data for model at index.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the experiment. |
0
|
channel
|
Optional[Union[PolarizationChannel, str]]
|
Spin channel of a polarized experiment. |
None
|
Returns:
| Type | Description |
|---|---|
Union[DataSet1D, PolarizedDataSet]
|
The experiment, or the |
Raises:
| Type | Description |
|---|---|
IndexError
|
No experiment is loaded at |
KeyError
|
|
ValueError
|
|
Source code in src/easyreflectometry/project.py
1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 | |
experiment_is_polarized_at_index(index=0)
Whether the experiment at index holds per-channel (polarized) data.
Returns False when no experiment is loaded at index, so consumers can
use it as a plain predicate.
Source code in src/easyreflectometry/project.py
1480 1481 1482 1483 1484 1485 1486 | |
experiment_channels_at_index(index=0)
Measured spin channels of the experiment at index ([] when unpolarized).
Source code in src/easyreflectometry/project.py
1488 1489 1490 1491 1492 1493 | |
experiment_supports_spin_asymmetry_at_index(index=0)
Whether SA can be formed for the experiment at index.
Spin asymmetry needs both non-spin-flip channels; an NSF-incomplete or spin-flip-only experiment has no SA. The two channel datasets must also be structurally usable (non-empty, matching lengths, a strictly ordered q grid) — an application asks this before offering a spin-asymmetry view, so a dataset that cannot be turned into one must not be advertised and then fail.
Source code in src/easyreflectometry/project.py
1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 | |
spin_asymmetry_for_experiment_at_index(index=0)
Spin asymmetry SA = (R⁺⁺ − R⁻⁻) / (R⁺⁺ + R⁻⁻) of a polarized experiment.
The measured SA is formed on the q grid of the pp channel. A mm channel
measured on a different grid is linearly interpolated onto it — values
with the usual weights, variances with the squared weights, which is
the propagation rule for independent endpoints (the covariance this
introduces between neighbouring SA points is not representable in
DataSet1D and is discarded). Interpolation happens only inside the q
range both channels cover: np.interp would otherwise clamp to the
edge value and turn extrapolated points into fabricated measurements.
Points are dropped when the denominator cannot carry a meaningful asymmetry:
- it is not positive, or is lost to cancellation between the two
channels (see :data:
SPIN_ASYMMETRY_CANCELLATION_FRACTION) — this guard needs no uncertainties and is what keeps a background-subtracted tail from throwing ±10³ values onto the axis; - it is not above :data:
SPIN_ASYMMETRY_SIGNIFICANCEtimes its own uncertainty, when the channels carry usable uncertainties; - the point itself is not usable (non-finite reflectivity, or a negative or non-finite variance).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the experiment. |
0
|
Returns:
| Type | Description |
|---|---|
Dict[str, object]
|
|
Raises:
| Type | Description |
|---|---|
IndexError
|
No experiment is loaded at |
ValueError
|
The experiment does not carry both pp and mm channels, or their datasets are not structurally usable. |
Source code in src/easyreflectometry/project.py
1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 | |
default_model()
Default model.
Source code in src/easyreflectometry/project.py
1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 | |
is_default_model(index)
Check if the model at the given index is a default model.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the model to check. |
required |
Returns:
| Type | Description |
|---|---|
bool
|
True if the model was created as a default placeholder. |
Source code in src/easyreflectometry/project.py
1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 | |
remove_model_at_index(index)
Remove the model at the given index.
Removes the model from the model collection, removes the experiment at the same index (if any), and reindexes experiments above the removed index so model/experiment indices stay aligned.
Adjusts the current model index if necessary.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
index
|
int
|
Index of the model to remove. |
required |
Raises:
| Type | Description |
|---|---|
IndexError :
|
If the index is out of range. |
ValueError :
|
If trying to remove the last remaining model. |
Source code in src/easyreflectometry/project.py
1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 | |
add_material(material)
Add a material to the project material collection.
:param material: the material to add. :raises ValueError: if the material is already in the collection. The collection is left untouched in that case, so a caller cannot mistake a refused add for a successful one.
Source code in src/easyreflectometry/project.py
1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 | |
remove_material(index)
Remove the material at index from the project material collection.
:param index: position of the material in the collection. :raises IndexError: if there is no material at index. :raises ValueError: if the material is used by one of the models. Removing it would leave the model referring to a material the project no longer holds, so the collection is left untouched.
Source code in src/easyreflectometry/project.py
1770 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 | |
create()
Create the project directory tree on disk.
:raises FileExistsError: if the project directory already exists. Nothing is written in
that case and created stays False, so a caller cannot mistake a collision for a new
project.
Source code in src/easyreflectometry/project.py
1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 | |
save_as_json(overwrite=False)
Save the project as json.
The write is atomic: the serialized project is written to a temporary file in the
destination directory and then moved into place with os.replace. Any failure while
serializing or writing therefore leaves a previously saved project file untouched.
Failures are raised rather than reported on stdout, so that callers (notably the GUI) can tell a successful save from a failed one.
:param overwrite: whether an existing project file may be replaced.
:raises FileExistsError: if the project file exists and overwrite is False.
:raises ValueError: if the project cannot be serialized, e.g. when a constraint depends
on a parameter that is not reachable from the models.
:raises OSError: if the file cannot be written or moved into place.
Source code in src/easyreflectometry/project.py
1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 | |
load_from_json(path=None)
Load a project file, replacing the current project state.
:param path: the project file; defaults to this project's own path_json.
:raises FileNotFoundError: if there is no file at path. The current project is left
untouched.
:raises ValueError: if the file is not valid JSON or predates the supported file format.
Source code in src/easyreflectometry/project.py
1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890 1891 1892 | |
as_dict(include_materials_not_in_model=False)
As dict.
Source code in src/easyreflectometry/project.py
1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 1930 1931 | |
from_dict(project_dict)
From dict.
Source code in src/easyreflectometry/project.py
1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 | |
constrain(parameter, expression, **parameters)
Tie parameter to an arbitrary expression of other parameters.
Placeholders in expression are supplied as keyword arguments:
.. code-block:: python
constrain(layer_b.thickness, '2 * t', t=layer_a.thickness)
parameter becomes dependent; its value, unit, variance, min and max
are replaced by the evaluated expression and its fixed flag is
cleared. Placeholder names must be valid Python identifiers and not
Python keywords. An unmapped name in the expression that matches a
mathematical builtin (e, pi, sin, ...) evaluates silently
instead of raising NameError, so prefer descriptive placeholder
names.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
parameter
|
Parameter
|
Parameter to make dependent (the follower). |
required |
expression
|
str
|
Mathematical expression to evaluate, e.g. |
required |
parameters
|
DescriptorNumber
|
Placeholder-name to parameter mapping for |
{}
|
Source code in src/easyreflectometry/constraints.py
72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 | |
detect_polarization_channel(path)
Detect the spin channel of a data file.
Tries the ORSO header (data_source.measurement.instrument_settings.polarization)
first, then filename heuristics ('_uu'/'_pp'/'_up' → pp, '_dd'/'_mm'/'_down' → mm,
'_ud'/'_pm' → pm, '_du'/'_mp' → mp, ...). Returns None when neither yields a
channel; a GUI should then ask the user to assign one.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
path
|
str
|
Path to the data file. |
required |
Returns:
| Type | Description |
|---|---|
Optional[PolarizationChannel]
|
The detected channel, or None. |
Source code in src/easyreflectometry/data/polarized.py
219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 | |
detect_polarization_channels_per_dataset(orso_data)
Classify every dataset of a parsed ORSO file by its own header.
Unlike :func:detect_polarization_channel, which reads only the first
dataset, this honours per-dataset polarization: overrides in
multi-dataset files.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
orso_data
|
list
|
Parsed ORSO dataset list (as returned by |
required |
Returns:
| Type | Description |
|---|---|
list[tuple[bool, Optional[PolarizationChannel]]]
|
Per dataset: (header declares a polarization, mapped channel or None). |
Source code in src/easyreflectometry/data/polarized.py
275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 | |
load_as_dataset(fname, data_group=None)
Load data from an ORSO .ort file as a DataSet1D.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
fname
|
Union[TextIO, str]
|
The file to be read. |
required |
data_group
|
Optional[sc.DataGroup]
|
Pre-loaded DataGroup for fname (avoids re-parsing the file). By default, None. |
None
|
Returns:
| Type | Description |
|---|---|
DataSet1D
|
The (first) dataset in the file. |
Source code in src/easyreflectometry/data/measurement.py
86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 | |
dataset_from_datagroup(data_group, data_key=None)
Build a DataSet1D from one dataset of a loaded DataGroup.
The ORSO header (when present) is attached to the returned dataset as the
orso_header attribute (a plain dict), so exporters can reuse the
original data_source/reduction provenance.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
data_group
|
sc.DataGroup
|
A DataGroup as returned by :func: |
required |
data_key
|
Optional[str]
|
The data entry to use (e.g. |
None
|
Returns:
| Type | Description |
|---|---|
DataSet1D
|
The dataset. |
Source code in src/easyreflectometry/data/measurement.py
46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 | |
extract_orso_title(data_group, data_name)
Extract orso title.
Source code in src/easyreflectometry/data/measurement.py
110 111 112 113 114 115 116 117 118 119 120 | |
load_measurement_file(fname)
Load data from an ORSO file (.ort/.orb) or a plain text file.
The discriminator is the ORSO banner line (or the HDF5 magic for binary files), not the file extension: a file carrying the banner that fails to parse raises instead of being silently re-read as plain text (which would drop the entire header, including polarization).
For a plain text file the columns are read by position in the ORSO order
Qz, R, sR, sQz -- further numeric columns are ignored -- and the error
columns must hold standard deviations, not variances. See
:func:_load_txt.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
fname
|
Union[TextIO, str]
|
The file to be read. |
required |
Returns:
| Type | Description |
|---|---|
sc.DataGroup
|
The loaded data. |
Source code in src/easyreflectometry/data/measurement.py
18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 | |
load_data_from_orso_file(fname)
Load data from an ORSO file (.ort text or .orb binary).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
fname
|
str
|
Path to the file. |
required |
Returns:
| Type | Description |
|---|---|
sc.DataGroup
|
A scipp DataGroup with data, coords, and attrs. |
Raises:
| Type | Description |
|---|---|
ValueError :
|
If the file cannot be parsed as ORSO. Parse failures are not swallowed here; falling back to plain-text loading is the caller's decision and only valid when the file carries no ORSO banner. |
Source code in src/easyreflectometry/orso_utils.py
177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 | |
build_constraints_factory(specs, resolve)
Build the constraints_factory hook for the given specs.
The returned callable is what easyscience's BUMPS engine invokes
with the {prefixed unique name: bumps.Parameter} mapping of a
freshly built problem; it returns one penalty object per enabled spec.
Paths are resolved and dependent parameters expanded at that moment, so
the factory always reflects the model as it is when the fit starts.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
specs
|
Iterable[InequalitySpec]
|
Specs to apply; disabled ones are skipped. |
required |
resolve
|
Callable[[str], DescriptorNumber]
|
Structural-path resolver, typically |
required |
Returns:
| Type | Description |
|---|---|
Callable | None
|
The factory, or |
Source code in src/easyreflectometry/inequality_constraints.py
370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 | |
check_units(spec, resolve)
Raise :class:UnitError when the two sides of spec have incompatible units.
Each side is evaluated with the unit-carrying DescriptorNumber
objects themselves (the same arithmetic the equality constraints use),
so t1 + sld is rejected by EasyScience and t_head < t_tail
passes. A plain numeric side is accepted against any unit: it is read in
the unit of the other side. The same applies to a side mixing literals
with parameters (90 - b): it is checked numerically and its literals
are read in the unit of the other side.
Source code in src/easyreflectometry/inequality_constraints.py
429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 | |
evaluate_spec(spec, resolve)
Evaluate a spec against the current parameter values.
Used for the start-point feasibility check before a fit is launched and for displaying the constraint state; dependent parameters contribute their current (already propagated) value.
Source code in src/easyreflectometry/inequality_constraints.py
414 415 416 417 418 419 420 421 422 423 424 425 426 | |
apply_default_limits(parameter, kind)
Apply default min/max to a parameter if current bounds are infinite.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
parameter
|
Parameter
|
The parameter to adjust. |
required |
kind
|
str
|
One of 'thickness', 'roughness', 'sld', 'isld', 'rho_m', 'scale'. |
required |
Source code in src/easyreflectometry/limits.py
12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | |
is_orso_file(fname)
Whether fname is an ORSO file: text with the ORSO banner line, or HDF5 (.orb).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
fname
|
str
|
Path to the file. |
required |
Returns:
| Type | Description |
|---|---|
bool
|
True when the file starts with the ORSO banner or the HDF5 magic bytes. |
Source code in src/easyreflectometry/orso_utils.py
95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 | |
save_orso_experiment(experiment, fname, model=None)
Write an experiment to an ORSO file (.ort text, or .orb binary).
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
experiment
|
A DataSet1D or PolarizedDataSet (the latter is written as one file with
one |
required | |
fname
|
str
|
Destination path; a |
required |
model
|
optional
|
The model whose sample is exported as the ORSO |
None
|
Source code in src/easyreflectometry/orso_utils.py
871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 | |
_weak_constraints_provider(project)
Source code in src/easyreflectometry/project.py
103 104 105 106 107 108 109 110 | |
_unwrapped_angle(angle, mask)
Make an angle profile continuous along z within each magnetic region.
The moment angle is periodic, so a profile that turns smoothly from 359° to
1° comes back from atan2 as a jump of nearly 360°. Drawn as a line that is
a full sweep across the chart where the moment barely moves. Each contiguous
run of mask (a magnetic region) is therefore unwrapped on its own and then
shifted by whole turns so it sits as close to the conventional [0, 360)
range as possible — neighbouring regions stay independent, since the angle
between them is not defined.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
angle
|
np.ndarray
|
Wrapped angles in degrees. |
required |
mask
|
np.ndarray
|
Which points carry a moment. |
required |
Returns:
| Type | Description |
|---|---|
np.ndarray
|
Angles in degrees, continuous within each masked region. Points outside the mask are returned unchanged (callers drop them). |
Source code in src/easyreflectometry/project.py
2133 2134 2135 2136 2137 2138 2139 2140 2141 2142 2143 2144 2145 2146 2147 2148 2149 2150 2151 2152 2153 2154 2155 2156 2157 2158 2159 2160 2161 2162 2163 2164 2165 2166 2167 2168 2169 | |
_has_both_nsf_channels(channels)
Whether both non-spin-flip channels (pp, mm) — the pair spin asymmetry needs — are present.
Source code in src/easyreflectometry/project.py
2172 2173 2174 | |
_ordered_channel_arrays(dataset)
A channel's (q, reflectivity, variance) arrays, ordered and checked.
Spin asymmetry pairs two channels point by point and interpolates one onto
the other, both of which assume well-formed, strictly increasing q. Rather
than trusting that, the arrays are checked here and sorted when needed —
np.interp silently returns nonsense for a descending grid.
Raises:
| Type | Description |
|---|---|
ValueError
|
The dataset is empty, its arrays disagree in length or shape, its q values are not finite, or it visits the same q twice. |
Source code in src/easyreflectometry/project.py
2177 2178 2179 2180 2181 2182 2183 2184 2185 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2197 2198 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210 2211 2212 2213 2214 2215 2216 2217 2218 2219 2220 2221 2222 2223 2224 2225 2226 | |
_interpolate_with_variance(q, q_source, values, variances)
Linear interpolation of values and their variances onto q.
Values use the linear weights (1−t, t); variances use their squares,
which is the propagation rule for independent endpoints — interpolating a
variance linearly (as np.interp would) overestimates it, by a factor 2 at
the midpoint of two equal variances.
q must lie inside q_source; the caller restricts it to the overlap.
Source code in src/easyreflectometry/project.py
2229 2230 2231 2232 2233 2234 2235 2236 2237 2238 2239 2240 2241 2242 2243 2244 2245 2246 | |
_spin_asymmetry(r_pp, r_mm, var_pp=None, var_mm=None)
Spin asymmetry, its variance, which points to keep, and why not.
Parameters:
| Name | Type | Description | Default |
|---|---|---|---|
r_pp
|
np.ndarray
|
Non-spin-flip reflectivities on a common q grid. |
required |
r_mm
|
np.ndarray
|
Non-spin-flip reflectivities on a common q grid. |
required |
var_pp
|
Optional[np.ndarray]
|
Their variances ( |
None
|
var_mm
|
Optional[np.ndarray]
|
Their variances ( |
None
|
Returns:
| Type | Description |
|---|---|
tuple
|
SA, its variance (zeros without input variances), a boolean mask of the points to keep, and a dict counting the dropped ones by reason. |
Source code in src/easyreflectometry/project.py
2249 2250 2251 2252 2253 2254 2255 2256 2257 2258 2259 2260 2261 2262 2263 2264 2265 2266 2267 2268 2269 2270 2271 2272 2273 2274 2275 2276 2277 2278 2279 2280 2281 2282 2283 2284 2285 2286 2287 2288 2289 2290 2291 2292 2293 2294 2295 2296 2297 2298 2299 2300 2301 2302 2303 2304 2305 2306 2307 2308 2309 2310 2311 2312 2313 2314 2315 2316 2317 | |