Changeset 7405 for trunk/MagicSoft/Mars/datacenter/macros

Timestamp:

11/16/05 16:01:44 (19 years ago)

Author:

Daniela Dorner

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/datacenter/macros

Files:

• fillstar.C (modified) (6 diffs)
• setupdb.C (modified) (1 diff)

trunk/MagicSoft/Mars/datacenter/macros/fillstar.C

@@ -45 +45 @@
 
 #include <TH1.h>
+#include <TH2.h>
 #include <TGraph.h>
 #include <TProfile.h>
@@ -106 +107 @@
     }
 
+    TH2F *hcog = (TH2F*)arr.FindObjectInCanvas("Center", "TH2F", "MHHillas");
+    if (!hcog)
+    {
+        cout << "WARNING - Reading of MHHillas failed." << endl;
+        return 0;
+    }
+
     MHMuonPar *hmuon = (MHMuonPar*)arr.FindObjectInCanvas("MHMuonPar", "MHMuonPar", "MHMuonPar");
     if (!hmuon)
@@ -112 +120 @@
         return 0;
     }
+
+    Double_t val[6];
+    for (int x=1; x<hcog->GetNbinsX(); x++)
+        for (int y=1; y<hcog->GetNbinsY(); y++)
+        {
+            Stat_t px = hcog->GetXaxis()->GetBinCenter(x);
+            Stat_t py = hcog->GetYaxis()->GetBinCenter(y);
+            Int_t  i  = (TMath::Nint(3*TMath::ATan2(px,py)/TMath::Pi())+6)%6;
+            val[i] += hcog->GetBinContent(x, y);
+        }
+
+    Double_t inhom = TMath::RMS(6, val)*6/hcog->GetEntries()*100;
+    inhom = TMath::Nint(inhom*10)/10.;
+    TString inhomogen = Form("%5.1f", inhom);
 
     Float_t mw    = hmuon->GetMeanWidth();
@@ -177 +199 @@
 
     cout << "Sequence #" << seq << endl;
+    cout << "  Inhomogeneity         " << inhomogen << endl;
     cout << "  PSF [mm]              " << PSF       << endl;
     cout << "  Island Quality        " << islands   << endl;
@@ -199 +222 @@
                      " fPSF=%s, "
                      " fDataRate=%d, "
-                     " fMaxHumidity=%s ",
+                     " fMaxHumidity=%s ,"
+                     " fInhomogeneity=%s ",
                      seq, islands.Data(), ratio.Data(),
                      num, effontime,
                      muonrate.Data(), PSF.Data(),
-                     datarate, maxhum.Data());
+                     datarate, maxhum.Data(),
+                     inhomogen.Data());
     }
     else
@@ -215 +240 @@
                      " fPSF=%s, "
                      " fDataRate=%d, "
-                     " fMaxHumidity=%s "
+                     " fMaxHumidity=%s, "
+                     " fInhomogeneity=%s "
                      " WHERE fSequenceFirst=%d ",
                      islands.Data(), ratio.Data(),
                      num, effontime,
                      muonrate.Data(), PSF.Data(),
-                     datarate, maxhum.Data(), seq);
+                     datarate, maxhum.Data(),
+                     inhomogen.Data(), seq);
     }
 

trunk/MagicSoft/Mars/datacenter/macros/setupdb.C

@@ -565 +565 @@
          "  fDataRate           INT        UNSIGNED   NOT NULL,"
          "  fMaxHumidity        FLOAT(6,1)            NOT NULL,"
+         "  fInhomogeneity      FLOAT(5,1)            NOT NULL,"
          "  fLastUpdate         TIMESTAMP"
          ")"));