Changeset 3053

Timestamp:

02/08/04 20:41:14 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mtools/MFFT.h (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -20 +20 @@
    * manalysis/MHPedestalPixel.cc
      - increased number of bins in for charges plot
+
+   * mtools/MFFT.[h,cc]
+     - possibility to have PSD directly from TArray
 
 

trunk/MagicSoft/Mars/mtools/MFFT.h

@@ -29 +29 @@
   void Swap(Double_t &a, Double_t &b)   { Double_t c = a;  a = b;  b = c;  }
 
-  void TransformF(const Int_t isign);
-  void TransformD(const Int_t isign);  
+  void TransformF(const Int_t isign, TArrayF &data);
+  void TransformD(const Int_t isign, TArrayD &data);  
   void RealFTF(const Int_t isign);
   void RealFTD(const Int_t isign);
@@ -36 +36 @@
   void CheckDim(Int_t a);
   TH1 *CheckHist(const TH1 *hist, const Int_t flag);
+
+  Float_t ApplyWindow(const Int_t j, const Float_t a, const Float_t b) const 
+    {
+
+      return 1.0-TMath::Abs((j-a)*b);      // Bartlett
+      // return 1.0;                        // Square
+      // return 1.0-(((j-a)*b)*((j-a)*b));   // Welch
+      
+    }
   
   Int_t   fDim;
   TArrayF fDataF;
   TArrayD fDataD;  
+  TArrayF fWindowF;
+  TArrayD fWindowD;  
 
 public:
@@ -58 +69 @@
   Double_t* RealFunctionIFFT(const Int_t n, const Double_t *data);  
   
+  TH1F* PowerSpectrumDensity(const TH1 *hist);
   TH1F* PowerSpectrumDensity(const TH1F *hist);
+  TH1F* PowerSpectrumDensity(const TH1I *hist);  
   TH1D* PowerSpectrumDensity(const TH1D *hist);
+
+  TArrayF* PowerSpectrumDensity(const TArrayF *array);
+  TArrayD* PowerSpectrumDensity(const TArrayD *array);
+  
+  TArrayF*  RealFunctionSpectrum(const TArrayF *data);
   
   ClassDef(MFFT,0)  // Class to perform a Fast Fourier Transform