Changeset 6907


Ignore:
Timestamp:
04/04/05 10:18:01 (20 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/Changelog

    r6906 r6907  
    3939   * mjobs/MJStar.cc:
    4040    - added MPointingPosCalc to tasklist for "Events" in case of MC
     41
     42   * ganymed.cc, mjobs/MJCut.[h,cc]:
     43     - implemented an option to enhance the display
     44     - do not display FS plot in wobble mode
     45     - made sure files are not written when they should not be written
     46
     47   * mhflux/MHAlpha.[h,cc]:
     48     - let FitEnergyBins return the mean of the significance of
     49       all Energy Bins
     50
     51   * mimage/MHHillasExt.cc, mimage/MHHillasSrc.cc,
     52     mimage/MHImagePar.cc, mimage/MHNewImagePar.cc:
     53     - implemented or fixed implementation of "same" Draw-option
     54     - slightly modified binning for Dist in MHHillasSrc
     55
     56   * mjobs/MJCalibrateSignal.cc:
     57     - do not plot pulse time twice
     58
     59   * mjobs/MJOptimize.cc:
     60     - print all parameters in minuit after they have been
     61       setup
    4162
    4263

  • trunk/MagicSoft/Mars/ganymed.cc

    r6890 r6907  
    6969    gLog << "   --print-ds                Print Dataset information" << endl;
    7070    gLog << "   --print-files             Print Files taken from Sequences ('+' found, '-' missing)" << endl;
     71    gLog << "   --full-display            Show as many plots as possible" << endl;
    7172    gLog << "   --config=ganymed.rc       Resource file [default=ganymed.rc]" << endl;
    7273    gLog << endl;
    … …  
    110111    const Bool_t  kDebugMem     = arg.HasOnlyAndRemove("--debug-mem");
    111112    const Bool_t  kWriteOnly    = arg.HasOnlyAndRemove("--write-only");
     113    const Bool_t  kFullDisplay  = arg.HasOnlyAndRemove("--full-display");
    112114    Int_t  kDebugEnv = arg.HasOnlyAndRemove("--debug-env") ? 1 : 0;
    113115    kDebugEnv = arg.GetIntAndRemove("--debug-env=", kDebugEnv);
    … …  
    235237    job.EnableWriteOnly(kWriteOnly);
    236238    job.EnableWobbleMode(iswobble);
     239    job.EnableFullDisplay(kFullDisplay);
    237240    if (kWriteSummary) // Don't change flag set in SetNameSummaryFile
    238241        job.EnableStorageOfSummary();

  • trunk/MagicSoft/Mars/mhflux/MHAlpha.cc

    r6906 r6907  
    184184}
    185185
    186 void MHAlpha::FitEnergyBins(Bool_t paint)
     186Float_t MHAlpha::FitEnergyBins(Bool_t paint)
    187187{
    188188    // Do not store the 'final' result in fFit
    … …  
    193193    fHEnergy.SetEntries(0);
    194194
     195    Float_t mean = 0;
     196    Int_t   num  = 0;
    195197    for (int i=1; i<=n; i++)
    196198    {
    … …  
    202204            else
    203205                fHEnergy.SetBinError(i, fit.GetEventsExcess());
    204         }
    205     }
     206
     207            if (fit.GetSignificance()>1)
     208            {
     209                mean += fit.GetSignificance();
     210                num++;
     211            }
     212        }
     213    }
     214    return mean/num;
    206215}
    207216
    … …  
    412421    fHAlphaTime.SetEntries(0);
    413422
    414     // Get number of bins
    415     const Int_t n = fHTime.GetNbinsX();
    416 
    417423    //
    418424    // Prepare Histogram
    … …  
    420426    if (final)
    421427        time->Plus1ns();
     428
     429    // Get number of bins
     430    const Int_t n = fHTime.GetNbinsX();
    422431
    423432    // Enhance binning

  • trunk/MagicSoft/Mars/mhflux/MHAlpha.h

    r6240 r6907  
    117117
    118118    void FitEnergySpec(Bool_t paint=kFALSE);
    119     void FitEnergyBins(Bool_t paint=kFALSE);
     119    Float_t FitEnergyBins(Bool_t paint=kFALSE);
    120120    void FitThetaBins(Bool_t paint=kFALSE);
    121121

  • trunk/MagicSoft/Mars/mimage/MHHillasExt.cc

    r5142 r6907  
    244244// together with the canvas.
    245245//
    246 void MHHillasExt::Draw(Option_t *)
     246void MHHillasExt::Draw(Option_t *o)
    247247{
    248248    TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this);
    … …  
    251251    AppendPad("");
    252252
    253     pad->Divide(2, 2);
     253    // FIXME: If same-option given make two independant y-axis!
     254    const TString opt(o);
     255    const Bool_t same = opt.Contains("same");
     256
     257    if (!same)
     258        pad->Divide(2,2);
     259    else
     260    {
     261        fHM3Long.SetLineColor(kMagenta);
     262        fHM3Trans.SetLineColor(kCyan);
     263
     264        fHAsym.SetLineColor(kBlue);
     265        fHMaxDist.SetLineColor(kBlue);
     266    }
    254267
    255268    pad->cd(1);
    256269    gPad->SetBorderMode(0);
    257     MH::DrawSame(fHM3Long, fHM3Trans, "3^{rd} Moments");
     270    MH::DrawSame(fHM3Long, fHM3Trans, "3^{rd} Moments", same);
    258271
    259272    pad->cd(3);
    260273    gPad->SetBorderMode(0);
    261     fHAsym.Draw();
     274    fHAsym.Draw(same?"same":"");
    262275
    263276    pad->cd(2);
    264277    gPad->SetBorderMode(0);
    265     fHMaxDist.Draw();
     278    fHMaxDist.Draw(same?"same":"");
    266279
    267280    delete pad->GetPad(4);

  • trunk/MagicSoft/Mars/mimage/MHHillasSrc.cc

    r5142 r6907  
    7171    //
    7272    fAlpha    = new TH1F("Alpha",    "Alpha of Ellipse",                 90,   -90,  90);
    73     fDist     = new TH1F("Dist",     "Dist of Ellipse",                 100,     0, 445);
     73    fDist     = new TH1F("Dist",     "Dist of Ellipse",                  70,     0, 623);
    7474    fCosDA    = new TH1F("CosDA",    "cos(Delta,Alpha) of Ellipse",     101,    -1,   1);
    7575    fDCA      = new TH1F("DCA",      "Distance of closest aproach",     101,  -500, 500);
    … …  
    214214// together with the canvas.
    215215//
    216 void MHHillasSrc::Draw(Option_t *)
     216void MHHillasSrc::Draw(Option_t *o)
    217217{
    218218    TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this);
    … …  
    223223    // FIXME: Display Source position
    224224
    225     pad->Divide(2, 2);
     225    // FIXME: If same-option given make two independant y-axis!
     226    const TString opt(o);
     227    const Bool_t same = opt.Contains("same");
     228
     229    if (!same)
     230        pad->Divide(2,2);
     231    else
     232    {
     233        fAlpha->SetLineColor(kGreen);
     234        fDist->SetLineColor(kGreen);
     235        fDCA->SetLineColor(kGreen);
     236        fCosDA->SetLineColor(kGreen);
     237        fDCADelta->SetLineColor(kGreen);
     238    }
    226239
    227240    pad->cd(1);
    228241    gPad->SetBorderMode(0);
    229     fAlpha->Draw();
     242    fAlpha->Draw(same?"same":"");
    230243
    231244    pad->cd(2);
    232245    gPad->SetBorderMode(0);
    233     fDist->Draw();
     246    fDist->Draw(same?"same":"");
    234247
    235248    pad->cd(3);
    236249    gPad->SetBorderMode(0);
    237     fDCA->Draw();
     250    fDCA->Draw(same?"same":"");
    238251
    239252    pad->cd(4);
    … …  
    241254
    242255    TVirtualPad *p = gPad;
    243     p->Divide(1,2);
     256    if (!same)
     257        p->Divide(1,2);
    244258    p->cd(1);
    245259    gPad->SetBorderMode(0);
    246     fCosDA->Draw();
     260    fCosDA->Draw(same?"same":"");
    247261
    248262    p->cd(2);
    249263    gPad->SetBorderMode(0);
    250     fDCADelta->Draw();
     264    fDCADelta->Draw(same?"same":"");
    251265}
    252266

  • trunk/MagicSoft/Mars/mimage/MHImagePar.cc

    r6869 r6907  
    231231        pad->Divide(2,2);
    232232    else
     233    {
     234        fHistSatHi.SetLineColor(kMagenta);
     235        fHistSatLo.SetLineColor(kCyan);
     236
     237        fHistSizeMainIsland.SetLineColor(kMagenta);
     238        fHistSizeSubIslands.SetLineColor(kCyan);
     239
    233240        fHistIslands.SetLineColor(kGreen);
    234 
     241        fHistNumSP.SetLineColor(kGreen);
     242        fHistSizeSP.SetLineColor(kGreen);
     243    }
     244
     245    pad->cd(1);
     246    gPad->SetBorderMode(0);
     247    MH::DrawSame(fHistSatHi, fHistSatLo, "Saturating Pixels", same);
     248    fHistSatHi.SetMinimum();    // switch off to allow log-scale
     249    fHistSatLo.SetMinimum();    // switch off to allow log-scale
     250    fHistSatLo.SetMaximum(0.1); // dummy value to allow log-scale
     251
     252    pad->cd(4);
     253    gPad->SetBorderMode(0);
     254    MH::DrawSame(fHistSizeMainIsland, fHistSizeSubIslands, "Sizes...", same);
     255    fHistSizeMainIsland.SetMinimum();    // switch off to allow log-scale
     256    fHistSizeSubIslands.SetMinimum();    // switch off to allow log-scale
     257    fHistSizeSubIslands.SetMaximum(0.1); // dummy value to allow log-scale
     258
     259    pad->cd(2);
     260    gPad->SetBorderMode(0);
    235261    if (!same)
     262        pad->GetPad(2)->Divide(1,2,1e-10,1e-10);
     263    if (pad->GetPad(2))
    236264    {
    237         pad->cd(1);
     265        pad->GetPad(2)->cd(1);
    238266        gPad->SetBorderMode(0);
    239         MH::DrawSame(fHistSatHi, fHistSatLo, "Saturating Pixels");
    240         fHistSatHi.SetMinimum();    // switch off to allow log-scale
    241         fHistSatLo.SetMinimum();    // switch off to allow log-scale
    242         fHistSatLo.SetMaximum(0.1); // dummy value to allow log-scale
    243 
    244         pad->cd(4);
     267        fHistNumSP.Draw(same?"same":"");
     268        pad->GetPad(2)->cd(2);
    245269        gPad->SetBorderMode(0);
    246         MH::DrawSame(fHistSizeMainIsland, fHistSizeSubIslands, "Sizes...");
    247         fHistSizeMainIsland.SetMinimum();    // switch off to allow log-scale
    248         fHistSizeSubIslands.SetMinimum();    // switch off to allow log-scale
    249         fHistSizeSubIslands.SetMaximum(0.1); // dummy value to allow log-scale
     270        fHistSizeSP.Draw(same?"same":"");
    250271    }
    251 
    252     pad->cd(2);
    253     gPad->SetBorderMode(0);
    254     pad->GetPad(2)->Divide(1,2,1e-10,1e-10);
    255     pad->GetPad(2)->cd(1);
    256     gPad->SetBorderMode(0);
    257     fHistNumSP.Draw(same?"same":"");
    258     pad->GetPad(2)->cd(2);
    259     gPad->SetBorderMode(0);
    260     fHistSizeSP.Draw(same?"same":"");
    261272
    262273    pad->cd(3);

  • trunk/MagicSoft/Mars/mimage/MHNewImagePar.cc

    r6869 r6907  
    284284void MHNewImagePar::Paint(Option_t *o)
    285285{
    286     if (TString(o)==(TString)"log" && fHistLeakage1.GetMaximum()>0)
    287         gPad->SetLogy();
     286    if (fHistLeakage1.GetMaximum()>0 && gPad->GetPad(1))
     287        gPad->GetPad(1)->SetLogy();
    288288}
    289289
    … …  
    294294// together with the canvas.
    295295//
    296 void MHNewImagePar::Draw(Option_t *)
     296void MHNewImagePar::Draw(Option_t *o)
    297297{
    298298    TVirtualPad *pad = gPad ? gPad : MakeDefCanvas(this);
    … …  
    301301    AppendPad("");
    302302
    303     pad->Divide(2,2);
     303    // FIXME: If same-option given make two independant y-axis!
     304    const TString opt(o);
     305    const Bool_t same = opt.Contains("same");
     306
     307    if (!same)
     308        pad->Divide(2,2);
     309    else
     310    {
     311        fHistLeakage1.SetLineColor(kMagenta);
     312        fHistLeakage1.SetLineColor(kCyan);
     313
     314        fHistCorePix.SetLineColor(kMagenta);
     315        fHistUsedPix.SetLineColor(kCyan);
     316
     317        fHistConc1.SetLineColor(kMagenta);
     318        fHistConc.SetLineColor(kCyan);
     319
     320        fHistCoreArea.SetLineColor(kMagenta);
     321        fHistUsedArea.SetLineColor(kCyan);
     322    }
    304323
    305324    pad->cd(1);
    … …  
    307326    TAxis &x = *fHistLeakage1.GetXaxis();
    308327    x.SetRangeUser(0.0, x.GetXmax());
    309     MH::DrawSame(fHistLeakage1, fHistLeakage2, "Leakage1 and Leakage2");
     328    MH::DrawSame(fHistLeakage1, fHistLeakage2, "Leakage1 and Leakage2", same);
    310329    fHistLeakage1.SetMinimum();
    311330    fHistLeakage2.SetMinimum();
    312331    fHistLeakage2.SetMaximum(0.1);   // dummy value to allow log-scale
    313     AppendPad("log");
    314332
    315333    pad->cd(2);
    316334    gPad->SetBorderMode(0);
    317     MH::DrawSame(fHistCorePix, fHistUsedPix, "Number of core/used Pixels");
     335    MH::DrawSame(fHistCorePix, fHistUsedPix, "Number of core/used Pixels", same);
    318336
    319337    pad->cd(3);
    320338    gPad->SetBorderMode(0);
    321     MH::DrawSame(fHistConc1, fHistConc, "Concentrations");
     339    MH::DrawSame(fHistConc1, fHistConc, "Concentrations", same);
    322340
    323341    pad->cd(4);
    324342    gPad->SetBorderMode(0);
    325     MH::DrawSame(fHistCoreArea, fHistUsedArea, "Area of core/used Pixels");
     343    MH::DrawSame(fHistCoreArea, fHistUsedArea, "Area of core/used Pixels", same);
    326344}
    327345

  • trunk/MagicSoft/Mars/mjobs/MJCalibrateSignal.cc

    r6906 r6907  
    866866
    867867    MGeomCam &geomcam = *(MGeomCam*)plist.FindObject("MGeomCam");
    868 
     868/*
    869869    //
    870870    // Get container from list
    … …  
    894894        disp21.CamDraw(c0, 2, 2, 1);
    895895      }
    896 
     896 */
    897897    if (!fIsHiLoCalibration)
    898898      return;

  • trunk/MagicSoft/Mars/mjobs/MJCut.cc

    r6892 r6907  
    5656#include "../mhflux/MAlphaFitter.h"
    5757#include "MBinning.h"
     58#include "MWeight.h"
    5859#include "MDataSet.h"
    5960#include "MParameters.h"
    … …  
    7172MJCut::MJCut(const char *name, const char *title)
    7273    : fStoreSummary(kFALSE), fStoreResult(kFALSE), fWriteOnly(kFALSE),
    73     fIsWobble(kFALSE), fEstimateEnergy(0), fCalcHadronness(0)
     74    fIsWobble(kFALSE), fFullDisplay(kFALSE), fEstimateEnergy(0),
     75    fCalcHadronness(0)
    7476{
    7577    fName  = name  ? name  : "MJCut";
    76     fTitle = title ? title : "Standard program to perform g/h-seperation cuts";
     78    fTitle = title ? title : "Standard program to perform g/h-separation cuts";
    7779}
    7880
    … …  
    260262//   "DataType"     to  "Events"
    261263//
    262 void MJCut::SetupWriter(MWriteRootFile &write, const char *name) const
    263 {
    264     write.SetName(name);
    265     write.AddContainer("MHillas",        "Events");
    266     write.AddContainer("MHillasSrc",     "Events");
    267     write.AddContainer("MHillasSrcAnti", "Events", kFALSE);
    268     write.AddContainer("MHadronness",    "Events", kFALSE);
    269     write.AddContainer("MEnergyEst",     "Events", kFALSE);
    270     write.AddContainer("DataType",       "Events");
     264void MJCut::SetupWriter(MWriteRootFile *write, const char *name) const
     265{
     266    if (!write)
     267        return;
     268
     269    write->SetName(name);
     270    write->AddContainer("MHillas",        "Events");
     271    write->AddContainer("MHillasSrc",     "Events");
     272    write->AddContainer("MHillasSrcAnti", "Events", kFALSE);
     273    write->AddContainer("MNewImagePar",   "Events", kFALSE);
     274    write->AddContainer("MNewImagePar2",  "Events", kFALSE);
     275    write->AddContainer("MHadronness",    "Events", kFALSE);
     276    write->AddContainer("MEnergyEst",     "Events", kFALSE);
     277    write->AddContainer("DataType",       "Events");
    271278
    272279    // Should not be the default: Either as option, or as
    … …  
    326333    MBinning bins3(50, 0,  60,   "BinningTheta",     "cos");
    327334    MBinning bins4("BinningFalseSource");
     335    //MBinning bins5("BinningWidth");
     336    //MBinning bins6("BinningLength");
     337    //MBinning bins7("BinningDist");
    328338    plist.AddToList(&bins1);
    329339    plist.AddToList(&bins2);
    330340    plist.AddToList(&bins3);
    331341    plist.AddToList(&bins4);
     342    //plist.AddToList(&bins5);
     343    //plist.AddToList(&bins6);
     344    //plist.AddToList(&bins7);
    332345
    333346    // --------------------------------------------------------------------------------
    … …  
    368381    fname1 += fNameResult.IsNull()  ?  (TString) Form("ganymed%08d-result.root",  set.GetNumAnalysis()) : fNameResult;
    369382
    370     MWriteRootFile write0(fPathOut.IsNull()?0:fname0.Data(), fOverwrite?"RECREATE":"NEW");
    371     MWriteRootFile write1(fPathOut.IsNull()?0:fname1.Data(), fOverwrite?"RECREATE":"NEW");
    372     if (CanStoreSummary())
    373         SetupWriter(write0, "WriteAfterCut0");
    374     if (CanStoreSummary())
    375         SetupWriter(write1, "WriteAfterCut3");
     383    MWriteRootFile *write0 = CanStoreSummary() ? new MWriteRootFile(fPathOut.IsNull()?0:fname0.Data(), fOverwrite?"RECREATE":"NEW") : 0;
     384    MWriteRootFile *write1 = CanStoreResult()  ? new MWriteRootFile(fPathOut.IsNull()?0:fname1.Data(), fOverwrite?"RECREATE":"NEW") : 0;
     385    SetupWriter(write0, "WriteAfterCut0");
     386    SetupWriter(write1, "WriteAfterCut3");
    376387
    377388
    … …  
    384395    taskenv2.SetDefault(fCalcHadronness);
    385396
    386     MFillH fill1a("MHHillasOffPre  [MHHillas]", "MHillas",    "FillHillasPre");
    387     MFillH fill2a("MHHillasOffPost [MHHillas]", "MHillas",    "FillHillasPost");
    388     MFillH fill3a("MHVsSizeOffPost [MHVsSize]", "MHillasSrc", "FillVsSizePost");
     397    MFillH fill1a("MHHillasOffPre  [MHHillas]",      "MHillas",      "FillHillasPre");
     398    MFillH fill2a("MHHillasOffPost [MHHillas]",      "MHillas",      "FillHillasPost");
     399    MFillH fill3a("MHVsSizeOffPost [MHVsSize]",      "MHillasSrc",   "FillVsSizePost");
     400    MFillH fill4a("MHHilExtOffPost [MHHillasExt]",   "MHillasSrc",   "FillHilExtPost");
     401    MFillH fill5a("MHHilSrcOffPost [MHHillasSrc]",   "MHillasSrc",   "FillHilSrcPost");
     402    MFillH fill6a("MHImgParOffPost [MHImagePar]",    "MImagePar",    "FillImgParPost");
     403    MFillH fill7a("MHNewParOffPost [MHNewImagePar]", "MNewImagePar", "FillNewParPost");
    389404    fill1a.SetNameTab("PreCut");
    390405    fill2a.SetNameTab("PostCut");
    391406    fill3a.SetNameTab("VsSize");
     407    fill4a.SetNameTab("HilExt");
     408    fill5a.SetNameTab("HilSrc");
     409    fill6a.SetNameTab("ImgPar");
     410    fill7a.SetNameTab("NewPar");
    392411
    393412    MPrint print2("MEffectiveOnTime");
    … …  
    412431    tlist2.AddToList(&taskenv2);
    413432    tlist2.AddToList(&cont0);
    414     if (CanStoreSummary())
    415         tlist2.AddToList(&write0);
     433    if (write0)
     434        tlist2.AddToList(write0);
    416435    if (!fWriteOnly)
    417436        tlist2.AddToList(&fill1a);
    418437    tlist2.AddToList(&cont1);
    419     if (!fWriteOnly)
     438    if (!fWriteOnly && !fIsWobble)
    420439        tlist2.AddToList(&ffs);
    421440    tlist2.AddToList(&cont2);
    … …  
    423442    {
    424443        tlist2.AddToList(&fill2a);
    425         tlist2.AddToList(&fill3a);
     444        if (fFullDisplay)
     445        {
     446            tlist2.AddToList(&fill3a);
     447            tlist2.AddToList(&fill4a);
     448            tlist2.AddToList(&fill5a);
     449            tlist2.AddToList(&fill6a);
     450            tlist2.AddToList(&fill7a);
     451        }
    426452    }
    427453    if (!fWriteOnly)
    428454        tlist2.AddToList(&falpha);
    429455    tlist2.AddToList(&cont3);
    430     if (CanStoreResult())
    431         tlist2.AddToList(&write1);
     456    if (write1)
     457        tlist2.AddToList(write1);
    432458
    433459    tlist.AddToList(&readoff);
    … …  
    461487            return kFALSE;
    462488        }
     489
     490        //plist.FindObject("MTimeEffectiveOnTime")->Clear();
    463491    }
    464492
    … …  
    473501        scalc.SetWobbleMode(kFALSE); /********************/
    474502
    475     MFillH fill1b("MHHillasOnPre  [MHHillas]", "MHillas",    "FillHillasPre");
    476     MFillH fill2b("MHHillasOnPost [MHHillas]", "MHillas",    "FillHillasPost");
    477     MFillH fill3b("MHVsSizeOnPost [MHVsSize]", "MHillasSrc", "FillVsSizePost");
     503    MFillH fill1b("MHHillasOnPre  [MHHillas]",      "MHillas",      "FillHillasPre");
     504    MFillH fill2b("MHHillasOnPost [MHHillas]",      "MHillas",      "FillHillasPost");
     505    MFillH fill3b("MHVsSizeOnPost [MHVsSize]",      "MHillasSrc",   "FillVsSizePost");
     506    MFillH fill4b("MHHilExtOnPost [MHHillasExt]",   "MHillasSrc",   "FillHilExtPost");
     507    MFillH fill5b("MHHilSrcOnPost [MHHillasSrc]",   "MHillasSrc",   "FillHilSrcPost");
     508    MFillH fill6b("MHImgParOnPost [MHImagePar]",    "MImagePar",    "FillImgParPost");
     509    MFillH fill7b("MHNewParOnPost [MHNewImagePar]", "MNewImagePar", "FillNewParPost");
    478510    fill1b.SetNameTab("PreCut");
    479511    fill2b.SetNameTab("PostCut");
    480512    fill3b.SetNameTab("VsSize");
     513    fill4b.SetNameTab("HilExt");
     514    fill5b.SetNameTab("HilSrc");
     515    fill6b.SetNameTab("ImgPar");
     516    fill7b.SetNameTab("NewPar");
    481517    fill1b.SetDrawOption(set.HasOffSequences()||fIsWobble?"same":"");
    482518    fill2b.SetDrawOption(set.HasOffSequences()||fIsWobble?"same":"");
    483519    fill3b.SetDrawOption(set.HasOffSequences()||fIsWobble?"same":"");
    484 
     520    fill4b.SetDrawOption(set.HasOffSequences()||fIsWobble?"same":"");
     521    fill5b.SetDrawOption(set.HasOffSequences()||fIsWobble?"same":"");
     522    fill6b.SetDrawOption(set.HasOffSequences()||fIsWobble?"same":"");
     523    fill7b.SetDrawOption(set.HasOffSequences()||fIsWobble?"same":"");
     524/*
     525    MWeight weight;
     526    weight.SetWeight(-1);
     527    fill2a.SetWeight(&weight);
     528    fill3a.SetWeight(&weight);
     529  */
    485530    MFillH falpha2("MHAlpha", "MHillasSrc", "FillAlpha");
    486531    MFillH ffs2("MHFalseSource", "MHillas", "FillFS");
    … …  
    490535    {
    491536        tlist2.Replace(&fill1b);
     537/*        if (fIsWobble)
     538        {
     539            tlist2.AddToListAfter(&fill2b, &fill2a);
     540            tlist2.AddToListAfter(&fill3b, &fill3a);
     541        }
     542        else
     543        */
    492544        tlist2.Replace(&fill2b);
    493         tlist2.Replace(&fill3b);
     545        if (fFullDisplay)
     546        {
     547            tlist2.Replace(&fill3b);
     548            tlist2.Replace(&fill4b);
     549            tlist2.Replace(&fill5b);
     550            tlist2.Replace(&fill6b);
     551            tlist2.Replace(&fill7b);
     552        }
    494553        tlist2.Replace(&falpha2);
    495         tlist2.Replace(&ffs2);
     554        if (!fIsWobble)
     555            tlist2.Replace(&ffs2);
    496556    }
    497557
    … …  
    518578    }
    519579
     580    if (write0)
     581        delete write0;
     582    if (write1)
     583        delete write1;
     584
    520585    tlist.PrintStatistics();
    521586

  • trunk/MagicSoft/Mars/mjobs/MJCut.h

    r6874 r6907  
    1717    Bool_t  fWriteOnly;
    1818    Bool_t  fIsWobble;
     19    Bool_t  fFullDisplay;
    1920
    2021    TString fNameSummary;
    … …  
    2728    TString GetOutputFile(UInt_t num) const;
    2829    Bool_t  CheckEnvLocal();
    29     void    SetupWriter(MWriteRootFile &write, const char *name) const;
     30    void    SetupWriter(MWriteRootFile *write, const char *name) const;
    3031    Bool_t  WriteTasks(UInt_t num, TObjArray &cont) const;
    3132    Bool_t  WriteResult(UInt_t num) const;
    … …  
    4445    void EnableWriteOnly(Bool_t b=kTRUE)         { fWriteOnly    = b; }
    4546    void EnableWobbleMode(Bool_t b=kTRUE)        { fIsWobble     = b; }
     47    void EnableFullDisplay(Bool_t b=kTRUE)       { fFullDisplay  = b; }
    4648
    4749    void SetNameSummaryFile(const char *name="");
    … …  
    5254    void SetHadronnessCalculator(const MTask *task=0);
    5355
    54     ClassDef(MJCut, 0) // Standard program to perform g/h-seperation cuts
     56    ClassDef(MJCut, 0) // Standard program to perform g/h-separation cuts
    5557};
    5658

  • trunk/MagicSoft/Mars/mjobs/MJOptimize.cc

    r6900 r6907  
    497497
    498498    gMinuit = new TMinuit(fParameters.GetSize());
    499     //gMinuit->SetPrintLevel(-1);
    500499
    501500    gMinuit->SetFCN(fcn);
    502501    gMinuit->SetObjectFit(this);
     502    gMinuit->SetPrintLevel(-1); // Don't print when DefineParameter
    503503
    504504    //
    … …  
    526526            gMinuit->FixParameter(i);
    527527    }
     528
     529    gMinuit->SetPrintLevel(1); // Switch on pritning again
     530    gMinuit->mnprin(1,0);      // Print all parameters
    528531
    529532    fEvtLoop = &evtloop;
Note: See TracChangeset for help on using the changeset viewer.