Changeset 3114 for trunk/MagicSoft/Mars

Timestamp:

02/12/04 11:45:59 (21 years ago)

Author:

jlopez

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/pointspreadfunction.C (modified) (10 diffs)

trunk/MagicSoft/Mars/macros/pointspreadfunction.C

@@ -51 +51 @@
 //______________________________________________________________________________
 //
-// Function use for Minuit to do the fit
+// Function used by Minuit to do the fit
 //
 void fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
@@ -80 +80 @@
 //______________________________________________________________________________
 //
-// The 2D gaussian fucntion use to fit the spot of the star
+// The 2D gaussian fucntion used to fit the spot of the star
 //
 Double_t func(float x,float y,Double_t *par)
@@ -136 +136 @@
     
     const Float_t conv4mm2deg = geomcam.GetConvMm2Deg();
+
          
     if (clname != "NULL")
@@ -288 +289 @@
 
     //
-    // Using the first dc meausrement we search the star and we define the area to be fitted
-    // around 3 rings of neightbords around the peak of the star.
+    // Using the first dc measurement we search the pixels which contains the star and define
+    // an area to be fitted by Minuit which is 3 rings of neightbords around the peak of the star.
     // 
 
@@ -362 +363 @@
    
     for(Int_t swpixel=1; swpixel<numPixels; swpixel++)
-    {
-//              cout << "isPixelsFitted[" << swpixel << "] ";
-//              cout << (Int_t)isPixelsFitted[swpixel] << endl;
         curr[swpixel] = (Float_t)isPixelsFitted[swpixel]; 
-    }
     
 /*    MHCamera display0(geomcam);
@@ -411 +408 @@
             step[i] = TMath::Abs(vstart[i]/sqrt(2));
         else
-            step[i] = 1;
-
-        cout << "DGB vstat[" << i << "] " << vstart[i] << " step[" << i << "] " << step[i] <<endl;
+            step[i] = uplimit[i]/2;
     }
-//    cout << endl;
   
     gMinuit->mnparm(0, "max", vstart[0], step[0], lowlimit[0], uplimit[0],ierflg);
@@ -424 +418 @@
 
     while (tlist.Process() && numLines < maxNumLines)
-        {
-//        curr.Print();
+    {
 
         for (Int_t swpixel=1; swpixel<577; swpixel++)
         {
-//          cout << swpixel << '\t' << curr[swpixel] << '\t' << geomcam[swpixel].GetX() << '\t' << geomcam[swpixel].GetY() << endl;
             
             x[swpixel] = geomcam[swpixel].GetX()*conv4mm2deg;
@@ -480 +472 @@
             display.SetCamContent(curr);
             gPad->cd(1);
-//      ellipse.DrawEllipse(xmean,ymean,TMath::Abs(xsigma),TMath::Abs(ysigma),0,360,0);
             ellipse.SetX1(xmean/conv4mm2deg);
             ellipse.SetY1(ymean/conv4mm2deg);
@@ -489 +480 @@
             gPad->Update();
 
-/*      if (numLines%10==0)
-        {
-            char imagepath[60];
-            sprintf(&imagepath,"/home/Javi/mnt_users/Data/CaCo/img/ZetaTau_%04d.gif",numLines/10);
-            cout << imagepath << endl;
-            gPad->SaveAs(imagepath);
-        }
-*/
             // Remove the comments if you want to go through the file
             //event-by-event:
@@ -673 +656 @@
 
     c1->cd(5);
-
-//    TCanvas *sigmacanvas = new TCanvas("sigmacanvas","Point Spread Funtion RMS",100,10,800,600);
 
     TMath math;