Changeset 2178 for trunk/MagicSoft/Mars/macros

Timestamp:

06/16/03 13:45:01 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/macros

Files:

• readCT1.C (modified) (2 diffs)
• readMagic.C (modified) (2 diffs)
• readcurrents.C (modified) (2 diffs)
• sumcurrents.C (modified) (2 diffs)

trunk/MagicSoft/Mars/macros/readCT1.C

@@ -89 +89 @@
             continue;
 
-        display.FillPhotNum(*evt);
+        display.Fill(*evt);
         display.Update();
 
@@ -100 +100 @@
         newimgpar.Print();
 
-        display.FillPhotNum(*(MCerPhotEvt*)plist.FindObject("MCerPhotEvt"));
+        display.Fill(*(MCerPhotEvt*)plist.FindObject("MCerPhotEvt"));
         display.Update();
 

trunk/MagicSoft/Mars/macros/readMagic.C

@@ -104 +104 @@
         cout << "Event #" << read.GetNumEntry() ":" << endl;
 
-        display.FillPhotNum(*(MCerPhotEvt*)clone.GetClone());
+        display.Fill(*(MCerPhotEvt*)clone.GetClone());
         gPad->Modified();
         gPad->Update();
@@ -116 +116 @@
         newimgpar.Print();
 
-        display.FillPhotNum(*(MCerPhotEvt*)plist.FindObject("MCerPhotEvt"));
+        display.Fill(*(MCerPhotEvt*)plist.FindObject("MCerPhotEvt"));
         gPad->Modified();
         gPad->Update();

trunk/MagicSoft/Mars/macros/readcurrents.C

@@ -45 +45 @@
 }
 
-void readcurrents(const char *fname="../currents/dcs_vega.dat")
+void readcurrents(const char *fname="../currents/dcs_arcturus.dat")
 {
     MParList plist;
@@ -86 +86 @@
     {
         // cur.Print();
-        display.FillCurrents(cur);
+        display.Fill(cur);
         gPad->Modified();
         gPad->Update();

trunk/MagicSoft/Mars/macros/sumcurrents.C

@@ -29 +29 @@
 //  This macro shows how to fill and display a histogram using Mars
 //
-void sumcurrents(const char *fname="../currents/dcs_arcturus.dat")
+void sumcurrents(const char *fname="../currents/spica_center.txt")
 {
     //
@@ -83 +83 @@
     MHCurrents &h = *(MHCurrents*)plist->FindObject("MHCurrents");
 
-    MCamDisplay *disp = new MCamDisplay(&geom);
-    disp->FillCurrents(h.GetSum());
-    disp->Draw();
+    TCanvas *c = MH::MakeDefCanvas();
+    c->Divide(3, 2);
+
+    MCamDisplay *disp1=new MCamDisplay(&geom);
+    MCamDisplay *disp2=new MCamDisplay(&geom);
+    MCamDisplay *disp3=new MCamDisplay(&geom);
+    disp1->Fill(h.GetSum());
+    disp2->Fill(h.GetRms());
+
+    TArrayF arr(577);
+    for (int i=0;i<577;i++)
+    {
+        TArrayF &r = h.GetRms();
+        TArrayF &v = h.GetSum();
+        arr[i] = (r[i]==0 ? 0 : v[i]/r[i]);
+    }
+    disp3->Fill(arr);
+
+    c->cd(1);
+    gPad->SetBorderMode(0);
+    gPad->Divide(1,1);
+    gPad->cd(1);
+    gPad->SetLogz();
+    disp1->Draw();
+    disp1->SetBit(kCanDelete);
+    c->cd(2);
+    gPad->SetBorderMode(0);
+    gPad->Divide(1,1);
+    gPad->cd(1);
+    gPad->SetLogz();
+    disp2->Draw();
+    disp2->SetBit(kCanDelete);
+    c->cd(3);
+    gPad->SetBorderMode(0);
+    gPad->Divide(1,1);
+    gPad->cd(1);
+    gPad->SetLogz();
+    disp3->Draw();
+    disp3->SetBit(kCanDelete);
+    c->cd(4);
+    gPad->SetBorderMode(0);
+    h.GetHist().DrawClone();
+    c->cd(5);
+    gPad->SetBorderMode(0);
+    TH1F h1("currents;rms", "Currents Rms", 577, -0.5, 576.5);
+    for (int i=1;i<=577; i++)
+        h1.SetBinContent(i, h.GetHist().GetBinError(i));
+    h1.DrawCopy();
+    c->cd(6);
+    gPad->SetBorderMode(0);
+    TH1F h2("currents;relerr", "Currents rel. Error [%]", 577, -0.5, 576.5);
+    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.DrawCopy();
+
+    c->SaveAs("spica_center-avg.ps");
 }