Changeset 8729 for trunk


Ignore:
Timestamp:
09/03/07 15:07:42 (17 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/macros/train
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/macros/train/traindisp.C

    r8718 r8729  
    9494     MFMagicCuts cuts;
    9595     cuts.SetHadronnessCut(MFMagicCuts::kArea);
    96      cuts.SetThetaCut(MFMagicCuts::kOn);
     96     cuts.SetThetaCut(MFMagicCuts::kNone);
    9797
    9898     TArrayD arr(13);
     
    112112     cuts.SetVariables(arr);
    113113
    114      opt.AddPreCut(&cuts);
     114     opt.AddTestCut(&cuts);
    115115
    116116     // -------------------- Other cuts ----------------------
     
    119119     opt.AddPreCut("MHillasSrc.fDCA*MGeomCam.fConvMm2Deg<0.3");
    120120     opt.AddPreCut("MPointingPos.fZd<25");
     121     opt.AddTestCut("MHillasExt.fM3Long   *sign(MHillasSrc.fCosDeltaAlpha)*MGeomCam.fConvMm2Deg>-0.07");
     122     opt.AddTestCut("MHillasExt.fSlopeLong*sign(MHillasSrc.fCosDeltaAlpha)/MGeomCam.fConvMm2Deg<(MHillasSrc.fDist*MGeomCam.fConvMm2Deg-0.5)*7.2");
    121123
    122124     // ---------------------- Histogram  --------------------
    123      MHn hist("MyHist", "Energy Residual (lg E_{est} - lg E_{mc})");
     125     MHn hist("MyHist", "Dist Residual (Disp-Dist)");
    124126
    125127     MBinning bins(50, -0.5, 0.5);
     
    166168     opt.AddPreTask(&weight);
    167169
     170     // ------- A trick for testing the parametrization -------
     171     // This will overwrite the random forest calculated disp.
     172     // Therefore you can test our disp-parametrization.
     173     MParameterCalc calc("(1.15136 +"
     174                         "0.0681437*MHillasExt.fSlopeLong*sign(MHillasSrc.fCosDeltaAlpha)/3.37e-3 +"
     175                         "0.0507821*(log10(MHillas.fSize)-1.51279)^2*(log10(MHillas.fSize)>1.51279) +"
     176                         "2.62932*MNewImagePar.fLeakage1)*"
     177                         "(1-MHillas.fWidth/MHillas.fLength)");
     178     calc.SetNameParameter("Disp");
     179     opt.AddPostTask(&calc);
     180
    168181     // ------------------------------------------------------
    169182    */
     183
     184     // To allow overwrite of existing files
     185     // opt.SetOverwrite();
    170186
    171187     // The number is the number of events read from the file
  • trunk/MagicSoft/Mars/macros/train/trainenergy.C

    r8708 r8729  
    165165    */
    166166
     167     // To allow overwrite of existing files
     168     // opt.SetOverwrite();
     169
    167170     opt.Train("rf-energy.root", set, 30000);
    168171}
Note: See TracChangeset for help on using the changeset viewer.