| 1 | void optime()
|
|---|
| 2 | {
|
|---|
| 3 | MJOptimize opt;
|
|---|
| 4 | opt.SetDebug(2);
|
|---|
| 5 | opt.SetOptimizer(MJOptimize::kMigrad);
|
|---|
| 6 | opt.EnableTestTrain();
|
|---|
| 7 |
|
|---|
| 8 | opt.AddParameter("MHillas.fSize");
|
|---|
| 9 | opt.AddParameter("MPointingPos.fZd");
|
|---|
| 10 |
|
|---|
| 11 | opt.SetParameter(0, 0.6, 0, 1000);
|
|---|
| 12 | opt.SetParameter(1, 0.0005, 0, 1000);
|
|---|
| 13 | opt.SetParameter(2, 0.89, 0, 1000);
|
|---|
| 14 | opt.FixParameter(3, 0, -1000, 1000);
|
|---|
| 15 | opt.FixParameter(4, 0, -1000, 1000);
|
|---|
| 16 | opt.FixParameter(5, 0, -1000, 1000);
|
|---|
| 17 | opt.FixParameter(6, 0, -1000, 1000);
|
|---|
| 18 | //opt.SetParameter(4, 0, -1000, 1000);
|
|---|
| 19 |
|
|---|
| 20 | MStatusDisplay *d = new MStatusDisplay;
|
|---|
| 21 | opt.SetDisplay(d);
|
|---|
| 22 |
|
|---|
| 23 | char *r = "([0]+(M[1]^2*[1]))*pow(M[0] + [3]*M[0]^2, [2])";
|
|---|
| 24 |
|
|---|
| 25 | opt.RunEnergy("ganymed00000999-result.root", r);
|
|---|
| 26 | //opt.RunEnergy("~/Software/mcwobble/img-abs/*_I_MonteCarlo[0123456789]_E.root", r);
|
|---|
| 27 | }
|
|---|
