Changeset 3779 for trunk/MagicSoft/Mars/mhist
- Timestamp:
- 04/19/04 11:36:48 (21 years ago)
- Location:
- trunk/MagicSoft/Mars/mhist
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mhist/MHAlpha.cc
r3778 r3779 27 27 // MHAlpha 28 28 // 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! 101 37 // 102 38 ////////////////////////////////////////////////////////////////////////////// -
trunk/MagicSoft/Mars/mhist/MHAlpha.h
r3778 r3779 48 48 static Double_t Significance(Double_t s, Double_t b); 49 49 50 ClassDef(MHAlpha, 1) // 3D-histogram in alpha, x and y50 ClassDef(MHAlpha, 1) // Alpha-Plot which is fitted online 51 51 }; 52 52
Note:
See TracChangeset
for help on using the changeset viewer.