Changeset 7009 for trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

Timestamp:

05/10/05 14:41:37 (19 years ago)

Author:

meyer

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

@@ -16 +16 @@
 !
 !
+!   Author(s): Keiichi Mase, 10/2004
 !   Author(s): Markus Meyer, 02/2005 <mailto:meyer@astro.uni-wuerzburg.de>
 !   Author(s): Thomas Bretz, 04/2005 <mailto:tbretz@astro.uni-wuerzburg.de>
@@ -28 +29 @@
 // MHSingleMuon
 //
-// This class is a histogram class for a displaying muonparameters.
-// The folowing histgrams will be plotted:
-// - Radius (TH1F)
-// - ArcWidth (TH1F)
+// This class is a histogram class for displaying the radial (fHistWidth)
+// and the azimuthal (fHistPhi) intensity distribution for one muon.
+// You can retrieve the histogram (TH1F) using the function GetHistPhi()
+// or GetHistWidth().
+// From these histograms the fraction of the ring segment (ArcPhi) and the
+// Width of the muon ring (ArcWidth) is calculated.
+//
+// First, the radius and center of the ring has to be calculted by
+// MMuonSearchParCalc
+// After that the histograms has to be filled in the following way:
+//
+// MFillH fillmuon("MHSingleMuon", "", "FillMuon");
+//
+// The allowed region to estimate ArcPhi is a certain margin around the
+// radius. The default value is 0.2 deg (60mm). If the estimated radius
+// of the arc is 1.0 deg, the pixel contents in the radius range from
+// 0.8 deg to 1.2 deg are fill in the histogram.
+//
+// For ArcPhi only bins over a certain threshold are supposed to be part
+// of the ring.
+// For ArcWidth, the same algorithm is used to determine the fit region
+// for a gaussian fit to the radial intensity distribution. The ArcWidth
+// is defined as the sigma value of the gaussian fit.
+//
+// The binning of the histograms can be changed in the following way:
+//
+// MBinning bins1("BinningMuonWidth");
+// MBinning bins2("BinningArcPhi");
+// bins1.SetEdges(28, 0.3, 1.7);
+// bins2.SetEdges(20, -180,180);
+// plist.AddToList(&bins1);
+// plist.AddToList(&bins2);
+//
+// The values for the thresholds and the margin are saved in MMuonSetup.
+// They can be easily changed in star.rc.
+//
+// Please have in mind, that changes in this basic parameters will change
+// your results!!
 //
 // Inputcontainer:
 //   - MGeomCam
 //   - MMuonSearchPar
-//   - MMuonCalibPar
+//
 //
 ////////////////////////////////////////////////////////////////////////////
@@ -319 +354 @@
     f1.SetParameter(0, fHistWidth.GetMaximumBin());
     f1.SetParameter(1, fHistWidth.GetBinCenter(m));
-    f1.SetParameter(2, (endfitval-startfitval)/2);
+//    f1.SetParameter(2, (endfitval-startfitval)/2);
+    f1.SetParameter(2, 0.1);
 
     // options : N  do not store the function, do not draw
@@ -338 +374 @@
 // --------------------------------------------------------------------------
 //
-//  An impact parameter is calculated by fitting the histogram of photon
+// An impact parameter is calculated by fitting the histogram of photon
 // distribution along the circle with a theoritical model. 
 // (See G. Vacanti et. al., Astroparticle Physics 2, 1994, 1-11. 
 // The function (6) is used here.) 
 //
-//  By default this calculation is suppressed because this calculation is 
+// By default this calculation is suppressed because this calculation is
 // very time consuming. If you want to calculate an impact parameter,
 // you can call the function of EnableImpactCalc().