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

added bkg

Browse Source
This commit is contained in:
nemu 2008-02-29 13:58:17 +00:00
parent dd2cbac068
commit 12c1aef4cd
1 changed files with 25 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
src/tests/skewedGaussianTest/skewedGaussian.C
View File

@ -20,6 +20,18 @@ Double_t PofB(Double_t B, Double_t B0, Double_t sm, Double_t sp)
val = TMath::Exp(-0.5*TMath::Power((B-B0)/sp, 2.0)); val = TMath::Exp(-0.5*TMath::Power((B-B0)/sp, 2.0));
} }
val *= TMath::Sqrt(2.0/TMath::Pi())/(sm+sp);
return val;
}
//-------------------------------------------------------------------------
Double_t PofBext(Double_t B, Double_t B0, Double_t s)
{
Double_t val = 0.0;
val = 1.0/TMath::Sqrt(2.0*TMath::Pi())/s*TMath::Exp(-0.5*TMath::Power((B-B0)/s, 2.0));
return val; return val;
} }
@ -29,11 +41,16 @@ void skewedGaussian()
FILE *fp; FILE *fp;
char fln[256]; char fln[256];
const Double_t B0 = 100.0; // (G) const Double_t w = 1.0; // weight of the skewed Gaussian
const Double_t sm = 4.0; // (G) const Double_t B0 = 100.0; // skewed Gaussian B0 (G)
const Double_t sp = 5.0; // (G) const Double_t sm = 4.5; // skewed Gaussian sigma- (G)
const Double_t sp = 4.5; // skewed Gaussian sigma+ (G)
sprintf(fln, "skewedGauss-B%0.2lf-sm%0.2lf-sp%0.2lf.dat", B0, sm, sp); const Double_t B0ext = 110.0; // external field Gaussian B0 (G)
const Double_t sext = 1.2; // external field Gaussian sigma (G)
sprintf(fln, "skewedGauss-B%0.2lf-sm%0.2lf-sp%0.2lf-w%0.1lf-Bext%0.2lf-sext%0.2lf.dat",
B0, sm, sp, w, B0ext, sext);
const Int_t noOfPoints = 1000; const Int_t noOfPoints = 1000;
@ -41,14 +58,15 @@ void skewedGaussian()
fp = fopen(fln, "w"); fp = fopen(fln, "w");
fprintf(fp, "# title: skewed Gaussian\n"); fprintf(fp, "# title: skewed Gaussian + background\n");
fprintf(fp, "# parameter: B0 = %lf, sm = %lf, sp = %lf\n", B0, sm, sp); fprintf(fp, "# parameter: B0 = %lf, sm = %lf, sp = %lf, w = %lf\n", B0, sm, sp, w);
fprintf(fp, "# B0ext = %lf, sext = %lf, \n", B0ext, sext);
fprintf(fp, "# -----------------------------------------------------\n"); fprintf(fp, "# -----------------------------------------------------\n");
fprintf(fp, "# B, pB\n"); fprintf(fp, "# B, pB\n");
Double_t b, pb; Double_t b, pb;
for (Int_t i=0; i<noOfPoints; i++) { for (Int_t i=0; i<noOfPoints; i++) {
b = B0-10.0*sm + res * i; b = B0-10.0*sm + res * i;
pb = PofB(b, B0, sm, sp); pb = w*PofB(b, B0, sm, sp) + (1.0-w)*PofBext(b, B0ext, sext);
fprintf(fp, "%lf, %lf\n", b, pb); fprintf(fp, "%lf, %lf\n", b, pb);
} }
fprintf(fp, "# end of file\n"); fprintf(fp, "# end of file\n");