Changeset 3081 for trunk/MagicSoft/Mars/mhbase/MH.cc

Timestamp:

02/10/04 16:10:52 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhbase/MH.cc (modified) (1 diff)

trunk/MagicSoft/Mars/mhbase/MH.cc

@@ -1081 +1081 @@
     return NULL;
 }
+
+TH1I* MH::ProjectArray(const TArrayF *array, Int_t nbins, const char* name, const char* title)
+{
+
+  const Int_t size  = array->GetSize();
+  TH1I *h1=0;
+
+  //check if histogram with identical name exist
+  TObject *h1obj = gROOT->FindObject(name);
+  if (h1obj && h1obj->InheritsFrom("TH1I")) {
+    h1 = (TH1I*)h1obj;
+    h1->Reset();
+  }
+
+  Double_t max = 0.;
+  Double_t min = 0.;
+  Int_t newbins = 0;
+
+  // first loop over array to find the maximum:
+  for (Int_t i=0; i<size;i++)
+    if (array->At(i) > max)
+      max = array->At(i);
+
+  FindGoodLimits(nbins, newbins, min, max, kFALSE);
+
+  if (!h1)
+    {
+      h1 = new TH1I(name, title, newbins, min, max);
+      h1->SetXTitle("");
+      h1->SetYTitle("Counts");
+      h1->SetDirectory(NULL);
+    }
+  
+    // Second loop to fill the histogram
+    for (Int_t i=0;i<size;i++) 
+      h1->Fill(array->At(i));
+
+    h1->SetEntries(size);
+    
+    return h1;
+}
+
+TH1I* MH::ProjectArray(const TArrayD *array, Int_t nbins, const char* name, const char* title)
+{
+
+  const Int_t size  = array->GetSize();
+  
+  TH1I *h1=0;
+
+  //check if histogram with identical name exist
+  TObject *h1obj = gROOT->FindObject(name);
+  if (h1obj && h1obj->InheritsFrom("TH1I")) {
+    h1 = (TH1I*)h1obj;
+    h1->Reset();
+  }
+
+  Double_t max = 0.;
+  Double_t min = 0.;
+  Int_t newbins = 0;
+
+  // first loop over array to find the maximum:
+  for (Int_t i=0; i<size;i++)
+    if (array->At(i) > max)
+      max = array->At(i);
+
+  FindGoodLimits( nbins, newbins, min, max, kFALSE);
+
+  if (!h1)
+    {
+      h1 = new TH1I(name, title, newbins, min, max);
+      h1->SetXTitle("");
+      h1->SetYTitle("Counts");
+      h1->SetDirectory(NULL);
+    }
+  
+    // Second loop to fill the histogram
+    for (Int_t i=0;i<size;i++) 
+      h1->Fill(array->At(i));
+
+    h1->SetEntries(size);
+    
+    return h1;
+}
+