Changeset 3113

Timestamp:

02/12/04 11:40:54 (21 years ago)

Author:

jlopez

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/macros/pointspreadfunction.C

@@ -37 +37 @@
 //     (It's possible not to use the calibration file and then the gain     //
 //      of the pmts are supouse to be the same for all of them.             //
+//  3- The third argument is just the number of dc measurements you want    //
+//     analize. If you put 0 it just stops after each fit and show you      //
+//     results.                                                             //
 //                                                                          //
 //--------------------------------------------------------------------------// 
 
 const Int_t numPixels = 577;
-Int_t nPixelsFitted; // 396[inners] 36[3rings] 60[4rings] 90[5rings]
+Int_t nPixelsFitted; 
 Bool_t isPixelsFitted[numPixels];
 Float_t z[numPixels],x[numPixels],y[numPixels],errorz[numPixels];
@@ -47 +50 @@
 
 //______________________________________________________________________________
+//
+// Function use for Minuit to do the fit
+//
 void fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
 {
@@ -73 +79 @@
 
 //______________________________________________________________________________
+//
+// The 2D gaussian fucntion use to fit the spot of the star
+//
 Double_t func(float x,float y,Double_t *par)
 {
@@ -102 +111 @@
 }
 
+
 void pointspreadfunction(TString fname, TString clname = "NULL", Int_t userNumLines = 1000)
 {
-
-// First run over continuos light files to have a DC calibration
 
 
@@ -131 +139 @@
     if (clname != "NULL")
     {
+
+        //
+        // First run over continuos light files to have a DC calibration
+        //
+
         MParList plist0;
         
@@ -179 +192 @@
             currrms[swpix] /= meanofcurrmean;   
             
-//      cout << "swpixel " << swpix << " dc fCalib " << currmean[swpix] << " +- " << currrms[swpix] << endl;
         }
         
         
-        //    curr.Print();
 /*     MHCamera display0(geomcam);
        display0.SetPrettyPalette();
        display0.Draw();
        
+        //    curr.Print();
        display0.SetCamContent(currbis);
        cout << "PSF>> DC mean values drawn" << endl;
@@ -198 +210 @@
     else
     {
+        // If you don't use the continuous light this is the currrms[] array 
+        // is the error associated to the currents TMinuit will use.
         for (Int_t swpix=0; swpix<numPixels; swpix++)
         {
             currmean[swpix] = 1.;
-            currrms[swpix] = 0.2;
+            currrms[swpix] = 0.2; 
         }
         
     }
 
-// Now we can run over the dc data to extract the psf rms. 
-        
+//
+// Now we can run over the dc data to extract the psf. 
+//      
     const Int_t maxNumLines = 10000;
         
@@ -241 +256 @@
     MHCamera display(geomcam);
     display.SetPrettyPalette();
-//    display.SetInvDeepBlueSeaPalette()
     display.Draw();
     gPad->SetLogy();
@@ -252 +266 @@
     Double_t xmean,xsigma,ymean,ysigma;
     Double_t xmeanerror,xsigmaerror,ymeanerror,ysigmaerror;
-
-//    TVector xellipsecenter;
-//    TVector yellipsecenter;
 
     TEllipse ellipse;
@@ -274 +285 @@
 
     TCanvas *psfcanvas = new TCanvas("psfcanvas","Point Spread Funtion 2D",200,20,900,700);
+
+
+    //
+    // 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.
+    // 
 
     tlist.Process();
@@ -362 +379 @@
         break;
 */
+
+    // Minuit initialization
+
     TMinuit *gMinuit = new TMinuit(7);  //initialize TMinuit with a maximum of 5 params
     gMinuit->SetFCN(fcn);
@@ -402 +422 @@
     gMinuit->mnparm(3, "ymean", vstart[3], step[3], lowlimit[3], uplimit[3],ierflg);
     gMinuit->mnparm(4, "ysigma", vstart[4], step[4], lowlimit[4], uplimit[4],ierflg);
-
-    numLines++;
 
     while (tlist.Process() && numLines < maxNumLines)
@@ -479 +497 @@
         }
 */
-            numLines++;
             // Remove the comments if you want to go through the file
             //event-by-event:
@@ -492 +509 @@
                     break;
             }
+            numLines++;
         }
     }
@@ -498 +516 @@
     evtloop.PostProcess();
 
-//Draw the ditributions of the sigmas the point spread function
-
-    cout<<"Number of lines "<<numLines<<endl;
+    //
+    // Draw the ditributions of the sigmas the point spread function
+    //
+
+    cout<<"PSF>> Number of lines "<<numLines<<endl;
 
     gROOT->Reset();