Index: /trunk/MagicSoft/Mars/Changelog
===================================================================
--- /trunk/MagicSoft/Mars/Changelog	(revision 7404)
+++ /trunk/MagicSoft/Mars/Changelog	(revision 7405)
@@ -18,4 +18,11 @@
 
                                                  -*-*- END OF LINE -*-*-
+ 2005/11/15 Daniela Dorner
+
+   * datacenter/macros/fillstar.C, setupdb.C:
+     - added new variable fInhomogeneity
+
+
+
  2005/11/15 Thomas Bretz
 
Index: /trunk/MagicSoft/Mars/datacenter/macros/fillstar.C
===================================================================
--- /trunk/MagicSoft/Mars/datacenter/macros/fillstar.C	(revision 7404)
+++ /trunk/MagicSoft/Mars/datacenter/macros/fillstar.C	(revision 7405)
@@ -45,4 +45,5 @@
 
 #include <TH1.h>
+#include <TH2.h>
 #include <TGraph.h>
 #include <TProfile.h>
@@ -106,4 +107,11 @@
     }
 
+    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,4 +120,18 @@
         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,4 +199,5 @@
 
     cout << "Sequence #" << seq << endl;
+    cout << "  Inhomogeneity         " << inhomogen << endl;
     cout << "  PSF [mm]              " << PSF       << endl;
     cout << "  Island Quality        " << islands   << endl;
@@ -199,9 +222,11 @@
                      " 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,10 +240,12 @@
                      " 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);
     }
 
Index: /trunk/MagicSoft/Mars/datacenter/macros/setupdb.C
===================================================================
--- /trunk/MagicSoft/Mars/datacenter/macros/setupdb.C	(revision 7404)
+++ /trunk/MagicSoft/Mars/datacenter/macros/setupdb.C	(revision 7405)
@@ -565,4 +565,5 @@
          "  fDataRate           INT        UNSIGNED   NOT NULL,"
          "  fMaxHumidity        FLOAT(6,1)            NOT NULL,"
+         "  fInhomogeneity      FLOAT(5,1)            NOT NULL,"
          "  fLastUpdate         TIMESTAMP"
          ")"));
