Changeset 3779 for trunk/MagicSoft


Ignore:
Timestamp:
04/19/04 11:36:48 (21 years ago)
Author:
tbretz
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
2 added
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mhist/MHAlpha.cc

    r3778 r3779  
    2727// MHAlpha
    2828//
    29 // Create a false source plot. For the Binning in x,y the object MBinning
    30 // "BinningFalseSource" is taken from the paremeter list.
    31 //
    32 // The binning in alpha is currently fixed as 18bins from 0 to 90deg.
    33 //
    34 // If MTime, MObservatory and MPointingPos is available in the paremeter
    35 // list each image is derotated.
    36 //
    37 // MHAlpha fills a 3D histogram with alpha distribution for
    38 // each (derotated) x and y position.
    39 //
    40 // Only a radius of 1.2deg around the camera center is taken into account.
    41 //
    42 // The displayed histogram is the projection of alpha<fAlphaCut into
    43 // the x,y plain and the projection of alpha>90-fAlphaCut
    44 //
    45 // By using the context menu (find it in a small region between the single
    46 // pads) you can change the AlphaCut 'online'
    47 //
    48 // Here is a slightly simplified version of the algorithm:
    49 // ------------------------------------------------------
    50 //    MHillas hil; // Taken as second argument in MFillH
    51 //
    52 //    MSrcPosCam src;
    53 //    MHillasSrc hsrc;
    54 //    hsrc.SetSrcPos(&src);
    55 //
    56 //    for (int ix=0; ix<nx; ix++)
    57 //        for (int iy=0; iy<ny; iy++)
    58 //        {
    59 //            TVector2 v(cx[ix], cy[iy]); //cx center of x-bin ix
    60 //            if (rho!=0)                 //cy center of y-bin iy
    61 //                v.Rotate(-rho);         //image rotation angle
    62 //
    63 //            src.SetXY(v);               //source position in the camera
    64 //
    65 //            if (!hsrc.Calc(hil))        //calc source dependant hillas
    66 //                return;
    67 //
    68 //            //fill absolute alpha into histogram
    69 //            const Double_t alpha = hsrc.GetAlpha();
    70 //            fHist.Fill(cx[ix], cy[iy], TMath::Abs(alpha), w);
    71 //        }
    72 //    }
    73 //
    74 // Use MHFalse Source like this:
    75 // -----------------------------
    76 //    MFillH fill("MHAlpha", "MHillas");
    77 // or
    78 //    MHAlpha hist;
    79 //    hist.SetAlphaCut(12.5);  // The default value
    80 //    hist.SetBgmean(55);      // The default value
    81 //    MFillH fill(&hist, "MHillas");
    82 //
    83 // To be implemented:
    84 // ------------------
    85 //  - a switch to use alpha or |alpha|
    86 //  - taking the binning for alpha from the parlist (binning is
    87 //    currently fixed)
    88 //  - a possibility to change the fit-function (eg using a different TF1)
    89 //  - a possibility to change the fit algorithm (eg which paremeters
    90 //    are fixed at which time)
    91 //  - use a different varaible than alpha
    92 //  - a possiblity to change the algorithm which is used to calculate
    93 //    alpha (or another parameter) is missing (this could be done using
    94 //    a tasklist somewhere...)
    95 //  - a more clever (and faster) algorithm to fill the histogram, eg by
    96 //    calculating alpha once and fill the two coils around the mean
    97 //  - more drawing options...
    98 //  - Move Significance() to a more 'general' place and implement
    99 //    also different algorithms like (Li/Ma)
    100 //  - implement fit for best alpha distribution -- online (Paint)
     29// Create a single Alpha-Plot. The alpha-plot is fitted online. You can
     30// check the result when it is filles in the MStatusDisplay
     31// For more information see MHFalseSource::FitSignificance
     32//
     33// For convinience (fit) the output significance is stored in a
     34// container in the parlisrt
     35//
     36// PRELIMINARY!
    10137//
    10238//////////////////////////////////////////////////////////////////////////////
  • trunk/MagicSoft/Mars/mhist/MHAlpha.h

    r3778 r3779  
    4848    static Double_t Significance(Double_t s, Double_t b);
    4949
    50     ClassDef(MHAlpha, 1) //3D-histogram in alpha, x and y
     50    ClassDef(MHAlpha, 1) // Alpha-Plot which is fitted online
    5151};
    5252
Note: See TracChangeset for help on using the changeset viewer.