Changeset 3531 for trunk/MagicSoft/Mars/manalysis

Timestamp:

03/17/04 13:19:12 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MSourcePosfromStarPos.cc (modified) (6 diffs)
• MSourcePosfromStarPos.h (modified) (1 diff)

trunk/MagicSoft/Mars/manalysis/MSourcePosfromStarPos.cc

@@ -99 +99 @@
     fRaSource  = 0.0;
 
+    // these are the default values when starting the excution;
+    // later these locations contain the values of the last event
     fxSourceold  =  25.0;
     fySourceold  = -40.0;
@@ -292 +294 @@
     if (fDecSource == 0.0  ||  fRaSource == 0.0  ||  fStars == 0) 
     {
-      *fLog << err << "MSourcePosfromStarPos::PreProcess; there are no sky coordinates defined for the source or from stars; fStars = " << fStars 
+      *fLog << warn << "MSourcePosfromStarPos::PreProcess; there are no sky coordinates defined for the source or from stars; fStars = " << fStars 
             << endl;
     }
@@ -531 +533 @@
         }
 
-        MSourcePosfromStarPos::SourcefromStar( fDistCameraReflector,
-                      fDecStar, fRaStar, fDecSource, fRaSource,
-                      fThetaTel[i],   fPhiTel[i],   
-                      xStar,          yStar,
-                      dxStar,         dyStar,
-                      fxSource,       fySource,
-                      fdxSource,      fdySource);
+        SourcefromStar( fDistCameraReflector,
+                        fDecStar, fRaStar, fDecSource, fRaSource,
+                        fThetaTel[i],   fPhiTel[i],   
+                        xStar,          yStar,
+                        dxStar,         dyStar,
+                        fxSource,       fySource,
+                        fdxSource,      fdySource);
       
         fSrcPos->SetXY(fxSource, fySource);
@@ -637 +639 @@
 
     Int_t fRows = fxStar.GetNrows();
+    if (fRows <= 0)
+      fRows = 1;
     fxStar.ResizeTo(fRows, fSize);
     fyStar.ResizeTo(fRows, fSize);
@@ -665 +669 @@
     
     Int_t fRows = fxStar.GetNrows();
+    if (fRows <= 0)
+      fRows = 1;
     fxStar.ResizeTo(fRows, fSize);
     fyStar.ResizeTo(fRows, fSize);
@@ -690 +696 @@
 
 
+  while(1)
+  {
+    *fIn >> ival;
+
+    // run number must be greater than 10000
+    if (fabs(ival) < 10000)
+    {
+      *fLog << err << "===========> Error when reading file with theta and phi <========="
+            << "             ival = " << ival << endl;
+    }          
+    else
+      break;
+  }  
+
   fRuns += 1;
-
-  *fIn >> ival;
-
-
-  
   //*fLog << fRuns <<"th run : " << ival << endl;  
 

trunk/MagicSoft/Mars/manalysis/MSourcePosfromStarPos.h

@@ -55 +55 @@
 
 
-    Double_t fDecSource;
+    Double_t fDecSource;  // sky coordinates of source
     Double_t fRaSource;
-    Double_t fxSource;
+    Double_t fxSource;    // estimated position of source in camera
     Double_t fySource;
     Double_t fdxSource;
     Double_t fdySource;
 
-    TArrayI fRunNr;
-    TArrayD fThetaTel;
+    TArrayI fRunNr;       // list of run numbers
+    TArrayD fThetaTel;    // for each run : theta, phi, ...
     TArrayD fPhiTel;
     TArrayD fdThetaTel;
     TArrayD fdPhiTel;
 
-    TArrayD  fDecStar;
+    TArrayD  fDecStar;   // sky coordinates of stars
     TArrayD  fRaStar;
-    TMatrixD fxStar;
+    TMatrixD fxStar;     // for each run : positions of stars in camera (from DC currents)
     TMatrixD fyStar;
     TMatrixD fdxStar;
     TMatrixD fdyStar;
 
-    Double_t fxSourceold;
+    Double_t fxSourceold;  
     Double_t fySourceold;
     Double_t fThetaradold;