| 1 | void optimcrab()
|
|---|
| 2 | {
|
|---|
| 3 | MJOptimize opt;
|
|---|
| 4 | opt.SetDebug(2);
|
|---|
| 5 | opt.SetOptimizer(MJOptimize::kSimplex);
|
|---|
| 6 |
|
|---|
| 7 | // ---------------------------------------------------
|
|---|
| 8 |
|
|---|
| 9 | opt.FixParameter(0, -0.10011, -5, 5);
|
|---|
| 10 | opt.FixParameter(1, 0.0669669, -5, 5);
|
|---|
| 11 |
|
|---|
| 12 | opt.FixParameter(2, 0.923, -5,5);
|
|---|
| 13 | opt.FixParameter(3, -0.8900478, -5, 5);
|
|---|
| 14 |
|
|---|
| 15 | opt.FixParameter(4, 0.215533, -5, 5);
|
|---|
| 16 |
|
|---|
| 17 | // ---------------------------------------------------
|
|---|
| 18 | opt.FixParameter(5, 0.106519, -1, 1);
|
|---|
| 19 | opt.FixParameter(6, 0.5);
|
|---|
| 20 | // --------------------ALPHA--------------------------
|
|---|
| 21 |
|
|---|
| 22 | opt.FixParameter(7, 2.13579, 0, 20);
|
|---|
| 23 | opt.FixParameter(8, 4.78171, 0, 20);
|
|---|
| 24 | opt.FixParameter(9, 1.17461, 0, 20);
|
|---|
| 25 | // ---------------------AREA--------------------------
|
|---|
| 26 |
|
|---|
| 27 | opt.FixParameter(10, 0.263628, 0, 20);
|
|---|
| 28 | opt.FixParameter(11, 5.25951, 0, 20);
|
|---|
| 29 | opt.FixParameter(12, 0.111246, 0, 20);
|
|---|
| 30 | // --------------------D1/D2--------------------------
|
|---|
| 31 |
|
|---|
| 32 | opt.FixParameter(13, 0.188982, 0, 20);
|
|---|
| 33 | opt.FixParameter(14, 1.10405, 0, 20);
|
|---|
| 34 |
|
|---|
| 35 | opt.FixParameter(15, 0.0960984, 0, 10);
|
|---|
| 36 | opt.FixParameter(16, 0.0991635, 0, 10);
|
|---|
| 37 |
|
|---|
| 38 | opt.FixParameter(17, 0, -10, 10);
|
|---|
| 39 | opt.FixParameter(18, 0, -10, 10);
|
|---|
| 40 | // --------------------MD1/MD2--------------------------
|
|---|
| 41 |
|
|---|
| 42 | opt.SetParameter(20, 0.4, -10, 10);
|
|---|
| 43 | opt.FixParameter(21, 0, -10, 10);
|
|---|
| 44 | opt.SetParameter(22, 1.0, -10, 10);
|
|---|
| 45 | opt.FixParameter(23, 0, -10, 10);
|
|---|
| 46 |
|
|---|
| 47 | opt.SetParameter(24, -1, -10, 10);
|
|---|
| 48 | opt.FixParameter(25, 0, -10, 10);
|
|---|
| 49 | opt.FixParameter(26, 0, -10, 10);
|
|---|
| 50 |
|
|---|
| 51 | // S=14.7, E=485, B=285
|
|---|
| 52 | // ---------------------------------------------------
|
|---|
| 53 |
|
|---|
| 54 | MStatusDisplay *d = new MStatusDisplay;
|
|---|
| 55 | opt.SetDisplay(d);
|
|---|
| 56 |
|
|---|
| 57 | MAlphaFitter fit;
|
|---|
| 58 | fit.SetScaleMode(MAlphaFitter::kNone);
|
|---|
| 59 | fit.SetSignalIntegralMax(99);
|
|---|
| 60 | //fit.SetMinimizationStrategy(MAlphaFitter::kExcess);
|
|---|
| 61 | fit.SetMinimizationStrategy(MAlphaFitter::kSignificanceLogExcess);
|
|---|
| 62 | //fit.SetMinimizationStrategy(MAlphaFitter::kSignificance);
|
|---|
| 63 |
|
|---|
| 64 | // A88A
|
|---|
| 65 | // DIST vs. atan(width/w, length/l)
|
|---|
| 66 | // Image cleaning with AT to get better size! (IMPORTANCE: E-est!)
|
|---|
| 67 |
|
|---|
| 68 | MFMagicCuts cuts;
|
|---|
| 69 | cuts.SetAlphaCut(MFMagicCuts::kWobble);
|
|---|
| 70 | opt.RunOnOff("wobblelza-abs.root", &cuts, &fit);
|
|---|
| 71 | }
|
|---|
