Changeset 3073 for trunk/MagicSoft/Mars/mcalib/MHCalibrationPixel.cc

Timestamp:

02/09/04 17:57:31 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationPixel.cc (modified) (12 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationPixel.cc

@@ -61 +61 @@
 const Int_t   MHCalibrationPixel::fNDFLimit          = 5;
 
-const Axis_t  MHCalibrationPixel::fNyquistFreq       = 1.0;
+const Axis_t  MHCalibrationPixel::fMaxFreq           = 3.5;
 const Axis_t  MHCalibrationPixel::fMinFreq           = 0.;
 const Int_t   MHCalibrationPixel::fPSDNbins          = 30;
@@ -79 +79 @@
         fChargeXaxis(NULL),
         fPSDXaxis(NULL),
-        fCurrentSize(1024),
-        fFitLegend(NULL)
+        fFitLegend(NULL),
+        fCurrentSize(1024)
 { 
 
@@ -260 +260 @@
       fHPSD = new TH1F(Form("%s%d","HPSD",fPixId),
                        Form("%s%s","Power Spectrum Density Projection ","LoGain"),
-                       fPSDNbins,fMinFreq,fNyquistFreq);
+                       fPSDNbins,fMinFreq,fMaxFreq);
       
       array = fPSDLoGain;
@@ -269 +269 @@
       fHPSD = new TH1F(Form("%s%d","HPSD",fPixId),
                        Form("%s%s","Power Spectrum Density Projection ","HiGain"),
-                       fPSDNbins,fMinFreq,fNyquistFreq);
+                       fPSDNbins,fMinFreq,fMaxFreq);
       
       array = fPSDHiGain;
@@ -286 +286 @@
   fPSDExpFit = new TF1(Form("%s%d","PSDExpFit",fPixId),"[0]*exp(-[1]*x)",0.,1.);
 
-  fPSDExpFit->SetParameters(entries,10.);
+  fPSDExpFit->SetParameters(area_guess,10.);
   fPSDExpFit->SetParNames("Area","slope");
   fPSDExpFit->SetParLimits(0,0.,3.*area_guess);
-  fPSDExpFit->SetParLimits(1,5.,20.);
-  fPSDExpFit->SetRange(fMinFreq,fNyquistFreq);
+  fPSDExpFit->SetParLimits(1,0.,20.);
+  fPSDExpFit->SetRange(fMinFreq,fMaxFreq);
 
   fHPSD->Fit(fPSDExpFit,"RQL0");
@@ -638 +638 @@
   fHAbsTimeLoGain->Draw(opt);
 
-  CutArrayBorder(fHiGains);
+  //  CutArrayBorder(fHiGains);
   CreateChargeXaxis(fHiGains->GetSize());
 
@@ -651 +651 @@
   
   CutArrayBorder(fLoGains);  
-  CreateChargeXaxis(fHiGains->GetSize());
+  CreateChargeXaxis(fLoGains->GetSize());
 
   c->cd(8);
@@ -658 +658 @@
                            fChargeXaxis->GetArray(),
                            fLoGains->GetArray());  
-  gr2->SetTitle("Evolution of HiGain charges with event number");
+  gr2->SetTitle("Evolution of LoGain charges with event number");
   gr2->SetBit(kCanDelete);
   gr2->Draw("AL");
@@ -668 +668 @@
   
   TArrayF *array;
+  TString title = "Power Spectrum Density ";
 
   if(IsUseLoGain())
@@ -674 +675 @@
         return;
       array = fPSDLoGain;
+      title += "LoGain";
     }
   else
@@ -680 +682 @@
         return;
       array = fPSDHiGain;
+      title += "HiGain";
     }
   
@@ -686 +689 @@
 
   TGraph *gr3 = new TGraph(fPSDXaxis->GetSize(),fPSDXaxis->GetArray(),array->GetArray());
-  gr3->SetTitle("Power Spectrum Density");
+  gr3->SetTitle(title.Data());
   gr3->SetBit(kCanDelete);
   gr3->Draw("AL");