LMU/musrfit
LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 3 Activity

shift the calculated field by the dead layer

Browse Source
This commit is contained in:
suter_a 2012-05-15 07:57:06 +00:00
parent f3214c3020
commit 7605aab46f
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/tests/nonlocal/PPippard.cpp
View File

@ -420,12 +420,18 @@ void PPippard::SaveField()
fprintf(fp, "%% Data --------------------------------------\n");
fprintf(fp, "%% z (nm), B/B_0 \n");
if (fParams.specular) {
for (Int_t i=0; i<(Int_t)(fParams.deadLayer/f_dz); i++) {
fprintf(fp, "%lf, %lf\n", f_dz*(Double_t)i, 1.0);
}
for (Int_t i=0; i<PippardFourierPoints/2; i++) {
fprintf(fp, "%lf, %lf\n", f_dz*(Double_t)i, fFieldB[i+fShift][1]);
fprintf(fp, "%lf, %lf\n", fParams.deadLayer + f_dz*(Double_t)i, fFieldB[i+fShift][1]);
}
} else {
for (Int_t i=0; i<(Int_t)(fParams.deadLayer/(f_dz * XiP_T(fParams.t))); i++) {
fprintf(fp, "%lf, %lf\n", f_dz * XiP_T(fParams.t) * (Double_t)i, 1.0);
}
for (Int_t i=0; i<PippardDiffusePoints; i++) {
fprintf(fp, "%lf, %lf\n", f_dz * XiP_T(fParams.t) * (Double_t)i, (*fFieldDiffuse)(i));
fprintf(fp, "%lf, %lf\n", fParams.deadLayer + f_dz * XiP_T(fParams.t) * (Double_t)i, (*fFieldDiffuse)(i));
}
}