Changeset 5135 for trunk/MagicSoft/Mars/mhcalib/MHGausEvents.cc

Timestamp:

09/25/04 13:51:15 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhcalib/MHGausEvents.cc (modified) (9 diffs)

trunk/MagicSoft/Mars/mhcalib/MHGausEvents.cc

@@ -88 +88 @@
 #include <TCanvas.h>
 #include <TStyle.h>
+#include <TRandom.h>
 
 #include "MFFT.h"
@@ -151 +152 @@
   fHGausHist.UseCurrentStyle();
   fHGausHist.SetDirectory(NULL);
-  fHGausHist.GetYaxis()->CenterTitle();
+  //  TAxis *xaxe = fHGausHist.GetXaxis();
+  //  xaxe->Set(100,0.,100.);
+  TAxis *yaxe = fHGausHist.GetYaxis();
+  //  yaxe->SetDefaults();
+  yaxe->CenterTitle();
 }
 
@@ -173 +178 @@
 {
 
-  // delete histograms
-  if (fHPowerProbability)
-    delete fHPowerProbability;
-
   //
   // The next two lines are important for the case that 
@@ -193 +194 @@
       delete fFExpFit;
   
+  // delete histograms
+  if (fHPowerProbability)
+    delete fHPowerProbability;
+
   // delete arrays
   if (fPowerSpectrum)  
@@ -203 +208 @@
   if (fGraphPowerSpectrum)
     delete fGraphPowerSpectrum;
+
 
 }
@@ -522 +528 @@
     {
       pad->cd(cwin++);
-      DrawPowerSpectrum(*pad,cwin);
-    }
-}
-
+      DrawPowerSpectrum();
+      pad->cd(cwin);
+      DrawPowerProjection();
+    }
+}
+
+// -----------------------------------------------------------------------------
+// 
+// DrawEvents:
+//
+// Will draw the graph with the option "A", unless the option:
+// "SAME" has been chosen 
+//
 void MHGausEvents::DrawEvents(Option_t *opt)
 {
@@ -551 +566 @@
 
 
-void MHGausEvents::DrawPowerSpectrum(TVirtualPad &pad, Int_t i)
-{
-  
+// -----------------------------------------------------------------------------
+// 
+// DrawPowerSpectrum
+//
+// Will draw the fourier spectrum of the events sequence with the option "A", unless the option:
+// "SAME" has been chosen 
+//
+void MHGausEvents::DrawPowerSpectrum(Option_t *option)
+{
+
+  TString opt(option);
+  
+  if (!fPowerSpectrum)
+    CreateFourierSpectrum();
+
   if (fPowerSpectrum)
     {
       if (!fGraphPowerSpectrum)
         CreateGraphPowerSpectrum();
+
+      if (!fGraphPowerSpectrum)
+        return;
       
-      fGraphPowerSpectrum->Draw("AL");          
-      fGraphPowerSpectrum->SetBit(kCanDelete);
-    }
-  
-  pad.cd(i);
+      if (opt.Contains("same"))
+        {
+          opt.ReplaceAll("same","");      
+          fGraphPowerSpectrum->Draw(opt+"L");
+        }
+      else
+        {
+          fGraphPowerSpectrum->Draw(opt+"AL");  
+          fGraphPowerSpectrum->SetBit(kCanDelete);
+        }
+    }
+}
+
+// -----------------------------------------------------------------------------
+// 
+// DrawPowerProjection
+//
+// Will draw the projection of the fourier spectrum onto the power probability axis
+// with the possible options of TH1D
+//
+void MHGausEvents::DrawPowerProjection(Option_t *option)
+{
+  
+  TString opt(option);
+
+  if (!fHPowerProbability)
+    CreateFourierSpectrum();
 
   if (fHPowerProbability && fHPowerProbability->GetEntries() > 0)
     {
       gPad->SetLogy();
-      fHPowerProbability->Draw();
+      fHPowerProbability->Draw(opt.Data());
       if (fFExpFit)
         {
@@ -772 +824 @@
 void MHGausEvents::InitBins()
 {
+  //  const TAttAxis att(fHGausHist.GetXaxis());
   fHGausHist.SetBins(fNbins,fFirst,fLast);
+  //  att.Copy(fHGausHist.GetXaxis());
 }
 
@@ -904 +958 @@
 }
 
+// ----------------------------------------------------------------------------
+//
+// Simulates Gaussian events and fills them  into the histogram and the array
+// In order to do a fourier analysis, call CreateFourierSpectrum()
+//
+void  MHGausEvents::SimulateGausEvents(const Float_t mean, const Float_t sigma, const Int_t nevts)
+{
+
+  if (!IsEmpty())
+    *fLog << warn << "The histogram is already filled, will superimpose simulated events on it..." << endl;
+  
+  for (Int_t i=0;i<nevts;i++) {
+    const Double_t ran = gRandom->Gaus(mean,sigma);
+    FillHistAndArray(ran);
+  }
+  
+}