Changeset 7225

Timestamp:

07/27/05 17:41:28 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• NEWS (modified) (1 diff)
• mmuon/MHSingleMuon.cc (modified) (2 diffs)
• mmuon/MHSingleMuon.h (modified) (3 diffs)
• mmuon/MMuonSetup.cc (modified) (2 diffs)
• star.rc (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -20 +20 @@
 
                                                  -*-*- END OF LINE -*-*-
+
+ 2005/07/27 Markus Meyer, Thomas Bretz
+
+   * mmuon/MHSingleMuon.[h,cc]:
+     - added third (constant) parameter to arc width gaus fit
+     - replaced nonsense GetMaximumBin as start value for fit
+       by GetMaximum
+     - replaced TH1F for fHistWidth by TProfile
+
+   * mmuon/MMuonSetup.cc:
+     - set new threshold for arc width plot to 2
+
+   * star.rc:
+     - replaced default for arc width threshold
+
+
 
  2005/07/27 Thomas Bretz

trunk/MagicSoft/Mars/NEWS

@@ -124 +124 @@
    - star: Output now contains MReportStarguider
 
-   - ganymed: corrects now (correctly?) for the misspointing in May/April 05
+   - star: Improved muon analysis. Due to the absense of a constant
+           for the gaus fit for arc width and the fact that the arcwidth
+           plot was a histogram instead of a profile the resulting psf
+           became pedestal rms dependant
+
+   - ganymed: corrects now (correctly?) for the misspointing in 
+              May/April 05
 
    - ganymed: Two new default templates for ganymed.rc are avilable:

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

@@ -244 +244 @@
 // below the threshold is found, kFALSE is returned.
 //
-Bool_t MHSingleMuon::FindRangeAboveThreshold(const TH1F &h, Float_t thres, Int_t &first, Int_t &last) const
+Bool_t MHSingleMuon::FindRangeAboveThreshold(const TH1 &h, Float_t thres, Int_t &first, Int_t &last) const
 {
     const Int_t n      = h.GetNbinsX();
@@ -348 +348 @@
 
     // Setup the function and perform the fit
-    TF1 f1("f1", "gaus", startfitval, endfitval);
+    TF1 f1("f1", "gaus + [3]", startfitval, endfitval);
     f1.SetLineColor(kBlue);
 
     // Choose starting values as accurate as possible
-    f1.SetParameter(0, fHistWidth.GetMaximumBin());
+    f1.SetParameter(0, fHistWidth.GetMaximum());
     f1.SetParameter(1, fHistWidth.GetBinCenter(m));
 //    f1.SetParameter(2, (endfitval-startfitval)/2);
     f1.SetParameter(2, 0.1);
+    f1.SetParameter(3, 1.8);
 
     // options : N  do not store the function, do not draw

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.h

@@ -5 +5 @@
 #include "MH.h"
 #endif
-#ifndef ROOT_TH1
-#include <TH1.h>
+#ifndef ROOT_TProfile
+#include <TProfile.h>
 #endif
 
@@ -24 +24 @@
     Double_t fMargin;               //!
 
-    TH1F fHistPhi;   // Histogram of photon distribution along the arc.
-    TH1F fHistWidth; // Histogram of radial photon distribution of the arc.
+    TH1F     fHistPhi;   // Histogram of photon distribution along the arc.
+    TProfile fHistWidth; // Histogram of radial photon distribution of the arc.
 
-    Bool_t FindRangeAboveThreshold(const TH1F &h, Float_t thres, Int_t &first, Int_t &last) const;
+    Bool_t FindRangeAboveThreshold(const TH1 &h, Float_t thres, Int_t &first, Int_t &last) const;
 
 public:
@@ -39 +39 @@
 
     const TH1F &GetHistPhi() const { return fHistPhi; }
-    const TH1F &GetHistWidth() const { return fHistWidth; }
+    const TProfile &GetHistWidth() const { return fHistWidth; }
 
     void Draw(Option_t *o="");

trunk/MagicSoft/Mars/mmuon/MMuonSetup.cc

@@ -49 +49 @@
 //
 MMuonSetup::MMuonSetup(const char *name, const char *title)
-    : fMargin(0.2), fThresholdArcPhi(30), fThresholdArcWidth(30)
+    : fMargin(0.2), fThresholdArcPhi(30), fThresholdArcWidth(2)
 {
     fName  = name  ? name  : "MMuonSetup";
@@ -60 +60 @@
 //   MMuonSetup.Margin:            0.2
 //   MMuonSetup.ThresholdArcPhi:   30
-//   MMuonSetup.ThresholdArcWidth: 30
+//   MMuonSetup.ThresholdArcWidth: 2
 //
 Int_t MMuonSetup::ReadEnv(const TEnv &env, TString prefix, Bool_t print)

trunk/MagicSoft/Mars/star.rc

@@ -58 +58 @@
 #MMuonSetup.Margin:            0.2
 #MMuonSetup.ThresholdArcPhi:   30
-#MMuonSetup.ThresholdArcWidth: 30
+#MMuonSetup.ThresholdArcWidth: 2
 
 #BinningRadius.Raw:            20  0.5  1.5