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

Timestamp:

02/09/04 22:50:52 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mcalib/MHCalibrationPixel.cc

@@ -61 +61 @@
 const Int_t   MHCalibrationPixel::fNDFLimit          = 5;
 
-const Axis_t  MHCalibrationPixel::fMaxFreq           = 30.5;
 const Axis_t  MHCalibrationPixel::fMinFreq           = 0.;
-const Int_t   MHCalibrationPixel::fPSDNbins          = 30;
+const Int_t   MHCalibrationPixel::fPSDNbins          = 50;
   
 // --------------------------------------------------------------------------
@@ -260 +259 @@
   //
   
-  // This cuts only the non-used zero's, but MFFT will cut the rest
+  // This cuts only the non-used zero's, but MFFT will later cut the rest
   CutArrayBorder(fHiGains);
   CutArrayBorder(fLoGains);  
 
-
   MFFT fourier;
 
@@ -272 +270 @@
   Int_t entries;
   TArrayF *array;
-  
+
+  const Double_t max_freq = (fChargeSigma > 0.) ? fChargeSigma*0.2 : 2.;
+  const Double_t min_freq = 0.;
+
   if (IsUseLoGain())
     {
       fHPSD = new TH1F(Form("%s%d","HPSD",fPixId),
                        Form("%s%s","Power Spectrum Density Projection ","LoGain"),
-                       fPSDNbins,fMinFreq,fMaxFreq);
+                       fPSDNbins,min_freq,max_freq);
       
       array = fPSDLoGain;
@@ -286 +287 @@
       fHPSD = new TH1F(Form("%s%d","HPSD",fPixId),
                        Form("%s%s","Power Spectrum Density Projection ","HiGain"),
-                       fPSDNbins,fMinFreq,fMaxFreq);
+                       fPSDNbins,min_freq,max_freq);
       
       array = fPSDHiGain;
@@ -305 +306 @@
   fPSDExpFit->SetParameters(area_guess,10.);
   fPSDExpFit->SetParNames("Area","slope");
-  fPSDExpFit->SetParLimits(0,0.,3.*area_guess);
-  fPSDExpFit->SetParLimits(1,0.,20.);
-  fPSDExpFit->SetRange(fMinFreq,fMaxFreq);
+  fPSDExpFit->SetParLimits(0,0.,20.*area_guess);
+  fPSDExpFit->SetParLimits(1,0.,30.);
+  fPSDExpFit->SetRange(min_freq,max_freq);
 
   fHPSD->Fit(fPSDExpFit,"RQL0");
@@ -351 +352 @@
 
   const Int_t diff = fChargeXaxis->GetSize()-n;
+
   fChargeXaxis->Set(n);
   if (diff < 0)
-    for (Int_t i=n;i<n+diff;i++)
+    for (Int_t i=n-1;i>n-diff-1;i--)
       fChargeXaxis->AddAt((Float_t)i,i);
 }
@@ -582 +584 @@
   
   TCanvas *c = MH::MakeDefCanvas(this,600,900); 
+  c->SetBit(kCanDelete);
   
   c->Divide(2,5);
@@ -714 +717 @@
   
   gStyle->SetOptStat(111111);
+  gStyle->SetOptFit();
   gPad->SetTicks();  
 
@@ -729 +733 @@
   c->Modified();
   c->Update();
+
+  c->Iconify();
   
   return;