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

Legend:

Unmodified
Added
Removed

  • 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.