| 1 | void optimwobble()
|
|---|
| 2 | {
|
|---|
| 3 | MJOptimizeCuts opt("MHThetaSq");
|
|---|
| 4 | opt.SetDebug(2);
|
|---|
| 5 | opt.SetOptimizer(MJOptimize::kSimplex);
|
|---|
| 6 |
|
|---|
| 7 | // ------------------- Xi -----------------------
|
|---|
| 8 | opt.FixParameter( 0, 1.15136);
|
|---|
| 9 | opt.FixParameter( 8, 0.0681437);
|
|---|
| 10 | opt.FixParameter( 9, 2.62932);
|
|---|
| 11 | opt.FixParameter(10, 1.51279);
|
|---|
| 12 | opt.FixParameter(11, 0.0507821);
|
|---|
| 13 |
|
|---|
| 14 | // --------------- Theta Sq ---------------------
|
|---|
| 15 | opt.FixParameter( 1, 0.215);
|
|---|
| 16 |
|
|---|
| 17 | // ----------------- Area -----------------------
|
|---|
| 18 | opt.SetParameter( 2, 0.21);
|
|---|
| 19 | opt.SetParameter( 3, 5.6);
|
|---|
| 20 | opt.SetParameter( 4, 0.083);
|
|---|
| 21 |
|
|---|
| 22 | // ----------------- M3L ------------------------
|
|---|
| 23 | opt.FixParameter( 5, -0.07);
|
|---|
| 24 |
|
|---|
| 25 | // ---------------- Slope -----------------------
|
|---|
| 26 | opt.FixParameter( 6, 7.2);
|
|---|
| 27 | opt.FixParameter( 7, 0.5);
|
|---|
| 28 |
|
|---|
| 29 | // ----------------- Fit ------------------------
|
|---|
| 30 |
|
|---|
| 31 | MAlphaFitter fit;
|
|---|
| 32 | fit.SetScaleMode(MAlphaFitter::kNone);
|
|---|
| 33 | fit.SetSignalFunction(MAlphaFitter::kThetaSq);
|
|---|
| 34 | fit.SetMinimizationStrategy(MAlphaFitter::kSignificanceLogExcess);
|
|---|
| 35 |
|
|---|
| 36 | // For wobble mode with three off-regions
|
|---|
| 37 | // (sets also ScaleMode to kUser)
|
|---|
| 38 | fit.SetScaleUser(1./3);
|
|---|
| 39 |
|
|---|
| 40 | // ----------------- Cuts -----------------------
|
|---|
| 41 |
|
|---|
| 42 | MFMagicCuts cuts;
|
|---|
| 43 | cuts.SetThetaCut(MFMagicCuts::kOn);
|
|---|
| 44 |
|
|---|
| 45 | /* You can also try (switch on scaling):
|
|---|
| 46 |
|
|---|
| 47 | fit.SetScaleMode(MAlphaFitter::kBackground);
|
|---|
| 48 | fit.SetBackgroundFitMin(0.12);
|
|---|
| 49 | fit.SetBackgroundFitMax(0.5);
|
|---|
| 50 |
|
|---|
| 51 | cuts.SetThetaCut(MFMagicCuts::kOff);
|
|---|
| 52 | */
|
|---|
| 53 |
|
|---|
| 54 | // -----------------------------------------------------------
|
|---|
| 55 |
|
|---|
| 56 | MStatusDisplay *d = new MStatusDisplay;
|
|---|
| 57 | opt.SetDisplay(d);
|
|---|
| 58 |
|
|---|
| 59 | opt.RunOnOff("ganymed00000001-summary.root", &cuts, &fit);
|
|---|
| 60 | }
|
|---|
