]>
| Syntax: | [phi, nph] = readMCHis(filename, mua_vec); | |
|---|---|---|
| Inputs: | filename | History file to read (without the .his extension) |
| mua_vec | Matrix of tissue absorption coefficients | |
| Outputs: | phi | Total fluence at each detector |
| nph | Total number of photons reaching each detector | |
The absorption can be adjusted after the fact because the path lengths are recorded in the history file and changing the absorption only changes the photon weightings, not its path. Changing the scattering, though, would change the actual path length distributions, and therefore to change the scattering requires re-running the simulation.