Changeset 2190 for trunk/MagicSoft/Mars/macros/sumcurrents.C

Timestamp:

06/18/03 09:29:45 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/sumcurrents.C (modified) (5 diffs)

trunk/MagicSoft/Mars/macros/sumcurrents.C

@@ -23 +23 @@
 \* ======================================================================== */
 
-// -------------------------------------------------------------------------
-//
-//  plot.C
-//
-//  This macro shows how to fill and display a histogram using Mars
-//
-void sumcurrents(const char *fname="../currents/spica_center.txt")
+void ProcessFile(TString fname)
 {
     //
@@ -51 +45 @@
     // (created with the star.C macro)
     MReadCurrents read(fname);
-    /*
-     read.AddFile("../currents/dcs_arcturus2.dat");
-     read.AddFile("../currents/dcs_arcturus3.dat");
-     read.AddFile("../currents/dcs_arcturus4.dat");
-     read.AddFile("../currents/dcs_arcturus5.dat");
-     read.AddFile("../currents/dcs_arcturus6.dat");
-     read.AddFile("../currents/dcs_arcturus7.dat");
-     read.AddFile("../currents/dcs_arcturus8.dat");
-     read.AddFile("../currents/dcs_arcturus9.dat");
-     read.AddFile("../currents/dcs_arcturus10.dat");
-     */
     tlist.AddToList(&read);
 
@@ -97 +80 @@
         TArrayF &r = h.GetRms();
         TArrayF &v = h.GetSum();
-        arr[i] = (r[i]==0 ? 0 : v[i]/r[i]);
+        arr[i] = (r[i]==0 ? 0 : r[i]/v[i]);
     }
     disp3->Fill(arr);
@@ -124 +107 @@
     c->cd(4);
     gPad->SetBorderMode(0);
-    h.GetHist().DrawClone();
+    h.GetHist().SetStats(kFALSE);
+    h.GetHist().DrawCopy();
     c->cd(5);
     gPad->SetBorderMode(0);
+    gPad->SetLogy();
     TH1F h1("currents;rms", "Currents Rms", 577, -0.5, 576.5);
+    h1.SetXTitle("Pixel Index");
+    h1.SetStats(kFALSE);
     for (int i=1;i<=577; i++)
         h1.SetBinContent(i, h.GetHist().GetBinError(i));
@@ -133 +120 @@
     c->cd(6);
     gPad->SetBorderMode(0);
+    gPad->SetLogy();
     TH1F h2("currents;relerr", "Currents rel. Error [%]", 577, -0.5, 576.5);
+    h2.SetXTitle("Pixel Index");
+    h2.SetStats(kFALSE);
     for (int i=1;i<=577; i++)
-    {
-        Float_t val = h.GetHist().GetBinContent(i);
-        Float_t rms = h.GetHist().GetBinError(i);
-        h2.SetBinContent(i, rms==0?0:val/rms);
-    }
+        h2.SetBinContent(i, arr[i]*100);
     h2.DrawCopy();
 
-    c->SaveAs("spica_center-avg.ps");
+    c->SaveAs(fname(0, fname.Last('.')+1) + "ps");
+    c->SaveAs(fname(0, fname.Last('.')+1) + "root");
 }
+
+// -------------------------------------------------------------------------
+//
+//  plot.C
+//
+//  This macro shows how to fill and display a histogram using Mars
+//
+void sumcurrents(const char *dirname="/home/MAGIC/online_data/ccdata/")
+{
+    MDirIter Next;
+    Next.AddDirectory(dirname, "dc_*.txt", -1);
+
+    TString fname;
+    while (1)
+    {
+        fname = Next();
+        if (fname.IsNull())
+            break;
+
+        ProcessFile(fname);
+    }
+}