Changeset 6993 for trunk/MagicSoft/Mars

Timestamp:

05/03/05 18:02:37 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mbadpixels/MBadPixelsCam.cc (modified) (1 diff)
• mjobs/MJCut.cc (modified) (4 diffs)
• mjobs/MJSpectrum.cc (modified) (6 diffs)
• mjobs/MJStar.cc (modified) (4 diffs)
• mjobs/MJStar.h (modified) (1 diff)
• mjobs/Makefile (modified) (1 diff)
• star.cc (modified) (3 diffs)
• star.rc (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -29 +29 @@
    * mpointing/MSrcPosCalc.cc:
      - removed dependancy on MTime for MC
+
+   * star.cc:
+     - addde possibility to switch off  muon analysis
+
+   * star.rc:
+     - added parameters for muon analysis
+
+   * mbadpixels/MBadPixelsCam.cc:
+     - fixed a typo in a comment
+
+   * mjobs/MJStar.[h,cc]:
+     - added possibility to switch off muon analysis
+
+   * mjobs/Makefile:
+     - added -I../mmuon
+
+   * mjobs/MJCut.cc, mjobs/MJSpectrum.cc:
+     - renamed MHAlpha to Hist (to have a common naming for all options)
+     - renamed MHAlphaOff to HistOff
 
 

trunk/MagicSoft/Mars/mbadpixels/MBadPixelsCam.cc

@@ -299 +299 @@
 // --------------------------------------------------------------------------
 //
-// Returns the size of the biggest cluster with the given USuitableType
+// Returns the size of the biggest cluster with the given UnsuitableType
 // type and the given area index.
 //

trunk/MagicSoft/Mars/mjobs/MJCut.cc

@@ -258 +258 @@
 
     // Save also the result, not only the setup
-    const MHAlpha *halpha = (MHAlpha*)plist.FindObject(fNameHist, "MHAlpha");
+    const MHAlpha *halpha = (MHAlpha*)plist.FindObject("Hist", "MHAlpha");
     if (halpha)
         arr.Add((TObject*)(&halpha->GetAlphaFitter()));
@@ -365 +365 @@
     }
 
-    const TString objname(Form("%s%s", fNameHist.Data(), name));
+    const TString objname(Form("%s%s", "Hist"/*fNameHist.Data()*/, name));
     MHAlpha *h = (MHAlpha*)plist.FindCreateObj(fNameHist, objname);
     if (!h)
@@ -446 +446 @@
 
     MHAlpha *halphaoff = CreateNewHist(plist, "Off");
-    MFillH falpha(halphaoff, "MHillasSrc", "FillAlpha");
+    MFillH falpha(halphaoff, "", "FillHist");
     MFillH ffs("MHFalseSourceOff [MHFalseSource]", "MHillas", "FillFS");
 
@@ -675 +675 @@
     */
     MHAlpha *halphaon=CreateNewHist(plist);
-    MFillH falpha2(halphaon, "MHillasSrc", "FillAlpha");
+    MFillH falpha2(halphaon, "", "FillHist");
     MFillH ffs2("MHFalseSource", "MHillas", "FillFS");
 

trunk/MagicSoft/Mars/mjobs/MJSpectrum.cc

@@ -66 +66 @@
 // Tasks/Filter
 #include "MReadMarsFile.h"
+#include "MReadMarsFile.h"
 #include "MFEventSelector2.h"
 #include "MFDataMember.h"
@@ -173 +174 @@
 
     TH1D *vstime = (TH1D*)arr.FindObjectInCanvas("Theta",  "TH1D", "OnTime");
-    TH1D *size   = (TH1D*)arr.FindObjectInCanvas("Excess", "TH1D", "MHAlpha");
+    TH1D *size   = (TH1D*)arr.FindObjectInCanvas("Excess", "TH1D", "Hist");
     if (!vstime || !size)
         return -1;
@@ -183 +184 @@
 
     if (fDisplay)
-        arr.DisplayIn(*fDisplay, "MHAlpha");
+        arr.DisplayIn(*fDisplay, "Hist");
 
     if (!ReadTask(fCut0, "Cut0"))
@@ -343 +344 @@
     taskenv1.SetDefault(fEstimateEnergy ? fEstimateEnergy : &est);
 
-    MFillH fill1("MHAlphaOff [MHAlpha]", "MHillasSrc", "FillAlphaOff");
-    MFillH fill2("MHAlphaOn  [MHAlpha]", "MHillasSrc", "FillAlphaOn");
+    MFillH fill1("HistOff [MHAlpha]", "MHillasSrc", "FillHistOff");
+    MFillH fill2("HistOn  [MHAlpha]", "MHillasSrc", "FillHistOn");
 
     MFDataMember f0("DataType.fVal", '<', 0.5, "FilterOffData");
@@ -381 +382 @@
     }
 
-    const MHAlpha *halpha = (MHAlpha *)plist.FindObject("MHAlphaOn");
+    const MHAlpha *halpha = (MHAlpha *)plist.FindObject("HistOn");
     if (!halpha)
     {
-        *fLog << err << GetDescriptor() << ": MHAlphaOn [MHAlpha] not found... abort." << endl;
+        *fLog << err << GetDescriptor() << ": HistOn [MHAlpha] not found... abort." << endl;
         return kFALSE;
     }
@@ -615 +616 @@
     cout << "Searching..." << endl;
 
-    TH1 *excess = (TH1D*)arr.FindObjectInCanvas("Excess", "TH1D", "MHAlpha");
+    TH1 *excess = (TH1D*)arr.FindObjectInCanvas("Excess", "TH1D", "Hist");
     if (!excess)
         return kFALSE;

trunk/MagicSoft/Mars/mjobs/MJStar.cc

@@ -62 +62 @@
 #include "MSrcPosCalc.h"
 #include "MHillasCalc.h"
-#include "../mmuon/MMuonSearchParCalc.h"
-#include "../mmuon/MMuonCalibParCalc.h"
-#include "../mmuon/MMuonSetup.h"
+#include "MMuonSearchParCalc.h"
+#include "MMuonCalibParCalc.h"
 #include "MFillH.h"
 #include "MWriteRootFile.h"
 
+#include "MMuonSetup.h"
 #include "MObservatory.h"
 #include "MPointingPosCalc.h"
@@ -81 +81 @@
 // Sets fRuns to 0, fExtractor to NULL, fDataCheck to kFALSE
 //
-MJStar::MJStar(const char *name, const char *title) 
+MJStar::MJStar(const char *name, const char *title) : fMuonAnalysis(kTRUE)
 {
     fName  = name  ? name  : "MJStar";
     fTitle = title ? title : "Standard analysis and reconstruction";
+}
+
+Bool_t MJStar::CheckEnvLocal()
+{
+    DisableMuonAnalysis(!GetEnv("MuonAnalysis", fMuonAnalysis));
+    return kTRUE;
 }
 
@@ -357 +363 @@
     tlist2.AddToList(&fill5);
 
-    // *********************************************************
-    MMuonSetup setup;
-    plist.AddToList(&setup);
-
+    // ----------------------- Muon Analysis ----------------------
     MF fmuon1("MHillas.fSize>150", "MuonPreCut");
-
-    MF fmuon2("(MMuonSearchPar.fRadius>180) &&"
-              "(MMuonSearchPar.fRadius<400) &&"
+    MF fmuon2("(MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) &&"
               "(MMuonSearchPar.fDeviation<45)", "MuonSearchCut");
-    fmuon2.SetFilter(&fmuon1);
-
-    MF fmuon3("(MMuonCalibPar.fArcPhi>175)  &&"
-              "(MMuonSearchPar.fDeviation<37)", "MuonFinalCut");
+    MF fmuon3("(MMuonCalibPar.fArcPhi>175)  && (MMuonSearchPar.fDeviation<37)",
+              "MuonFinalCut");
 
     MMuonSearchParCalc muscalc;
@@ -377 +376 @@
     mcalc.SetFilter(&fmuon2);
 
-    MFillH fill11("MHSingleMuon", "", "FillMuon");
-    MFillH fill10("MHMuonPar",    "", "FillMuonPar");
-    fill11.SetFilter(&fmuon2);
-    fill10.SetFilter(&fmuon3);
-
-    tlist2.AddToList(&fmuon1);
-    tlist2.AddToList(&muscalc);
-    tlist2.AddToList(&fmuon2);
-    tlist2.AddToList(&fill11);
-    tlist2.AddToList(&mcalc);
-    tlist2.AddToList(&fmuon3);
-    tlist2.AddToList(&fill10);
-
-    // *********************************************************
+    MFillH fillmuon("MHSingleMuon", "", "FillMuon");
+    MFillH fillmpar("MHMuonPar",    "", "FillMuonPar");
+    fillmuon.SetFilter(&fmuon2);
+    fillmpar.SetFilter(&fmuon3);
+
+    if (fMuonAnalysis)
+    {
+        tlist2.AddToList(&fmuon1);
+        tlist2.AddToList(&muscalc);
+        tlist2.AddToList(&fmuon2);
+        tlist2.AddToList(&fillmuon);
+        tlist2.AddToList(&mcalc);
+        tlist2.AddToList(&fmuon3);
+        tlist2.AddToList(&fillmpar);
+    }
+    // ------------------------------------------------------------
 
     //tlist2.AddToList(&fill6);

trunk/MagicSoft/Mars/mjobs/MJStar.h

@@ -9 +9 @@
 {
 private:
+    Bool_t fMuonAnalysis;
+
     Bool_t WriteResult();
+    Bool_t CheckEnvLocal();
 
 public:
     MJStar(const char *name=NULL, const char *title=NULL);
 
+    // Setup
+    void DisableMuonAnalysis(Bool_t b=kTRUE) { fMuonAnalysis = !b; }
+
+    // Process
     Bool_t Process(Bool_t ismc=kFALSE);
 

trunk/MagicSoft/Mars/mjobs/Makefile

@@ -23 +23 @@
            -I../mbadpixels -I../msignal -I../mraw -I../mpedestal -I../mtools \
            -I../mimage -I../mpointing -I../mastro -I../mfbase -I../mhvstime \
-           -I../mtrigger
+           -I../mtrigger -I../mmuon
 
 SRCFILES = MSequence.cc \

trunk/MagicSoft/Mars/star.cc

@@ -66 +66 @@
     gLog << "   --out=path                Path to write the all results to [def=local path]" << endl;
     gLog << "                             (overwrites --outc and --outy)" << endl;
+    gLog << "   --no-muons                Switch off Muon analysis (for fast tests)" << endl;
     gLog << "   --print-seq               Print Sequence information" << endl;
     gLog << "   --print-files             Print Files taken from Sequence" << endl;
@@ -122 +123 @@
     //const Bool_t  kForceExec  = arg.HasOnlyAndRemove("-ff");
     const Bool_t  kIsMC       = arg.HasOnlyAndRemove("-mc");
+    const Bool_t  kNoMuons    = arg.HasOnlyAndRemove("--no-muons");
 
     const TString kInpath     = arg.GetStringAndRemove("--ind=", "");
@@ -224 +226 @@
         job.SetPathData(kInpath);
         // job.SetPathIn(kInpath); // not yet needed
+        if (kNoMuons)
+            job.DisableMuonAnalysis();
 
         if (!job.Process(kIsMC))

trunk/MagicSoft/Mars/star.rc

@@ -43 +43 @@
 #MJStar.MaxEvents: 10000
 #MJStar.Overwrite: yes,no
-#MJStar.
 
 # -------------------------------------------------------------------------
-# Setup the muon analysis here
+# Setup or switch off the muon analysis here
 # -------------------------------------------------------------------------
-# MMuonSetup.Margin:            0.2
-# MMuonSetup.ThresholdArcPhi:   30
-# MMuonSetup.ThresholdArcWidth: 30
-# BinningRadius.Raw:            20  0.5  1.5
-# BinningArcWidth.Raw:          60  0.0  0.3
-# BinningRingBroadening.Raw:    20  0.5  1.5
-# BinningSizeVsArcRadius.Raw:   20  0.5  1.5
-# BinningsArcPhi.Raw:           21 -180  180
-# BinningsMuonWidth.Raw:        28  0.3  1.7
+
+#MJStar.MuonAnalysis: On
+
+#MMuonSetup.Margin:            0.2
+#MMuonSetup.ThresholdArcPhi:   30
+#MMuonSetup.ThresholdArcWidth: 30
+
+#BinningRadius.Raw:            20  0.5  1.5
+#BinningArcWidth.Raw:          60  0.0  0.3
+#BinningRingBroadening.Raw:    20  0.5  1.5
+#BinningSizeVsArcRadius.Raw:   20  0.5  1.5
+#BinningsArcPhi.Raw:           21 -180  180
+#BinningsMuonWidth.Raw:        28  0.3  1.7
+
+#MuonPreCut.Condition:   MHillas.fSize>150
+#MuonSearcCut.Condition: (MMuonSearchPar.fRadius>180) && (MMuonSearchPar.fRadius<400) && (MMuonSearchPar.fDeviation<45)
+#MuonFinalCut.Condition: (MMuonCalibPar.fArcPhi>175)  && (MMuonSearchPar.fDeviation<37)
 
 # -------------------------------------------------------------------------
 # -------------------------------------------------------------------------
-MJStar.MImgCleanStd.CleanLevel1: 4.5
-MJStar.MImgCleanStd.CleanLevel2: 3.0
-MJStar.MImgCleanStd.CleanMethod: Scaled
+MImgCleanStd.CleanLevel1: 4.5
+MImgCleanStd.CleanLevel2: 3.0
+MImgCleanStd.CleanMethod: Scaled
 
 # -------------------------------------------------------------------------