Changeset 5876 for trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.cc

Timestamp:

01/17/05 16:39:19 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCam.cc

@@ -691 +691 @@
   TH1D *h = convcam.ProjectionS(TArrayI(),TArrayI(1,&aidx),"_py",750);
   h->SetDirectory(NULL);
-  h->Fit("gaus","Q");
-  TF1 *fit = h->GetFunction("gaus");
-  
+
+  TF1 *fit = NULL;
+
+  if (geom.InheritsFrom("MGeomCamMagic"))
+    {
+
+      fit = new TF1("fit","gaus",0.4,5.);
+      
+      // Fix the ranges, as found by Nadia
+      if(aidx == 0)
+        {h->Fit("fit","REQ", "",0.4,1.5);}
+      else
+        {h->Fit("fit","REQ", "",1.,5.);}
+    }
+  else
+    {
+      h->Fit("gaus","Q");
+      fit = h->GetFunction("gaus");
+    }
+
   Float_t ci2   = fit->GetChisquare();
   Float_t sigma = fit->GetParameter(2);
@@ -699 +716 @@
   if (ci2 > 500. || sigma > sg)
     {
-      h->Fit("gaus","QM");
-      fit = h->GetFunction("gaus");
+      if (geom.InheritsFrom("MGeomCamMagic"))
+        {
+          // Fix the ranges, as found by Nadia
+          if(aidx == 0)
+            {h->Fit("fit","REQ", "",0.4,1.5);}
+          else
+            {h->Fit("fit","REQ", "",1.,5.);}
+        }
+      else
+        {
+          h->Fit("gaus","MREQ");
+          fit = h->GetFunction("gaus");
+        }
       
       ci2   = fit->GetChisquare();
@@ -720 +748 @@
         << "  Chisquare: " << Form("%4.3f",fit->GetChisquare()) << "  NDF  : " << ndf << endl;          
   
+  delete fit;
   delete h;
   gROOT->GetListOfFunctions()->Remove(fit);