Changeset 4707 for trunk/MagicSoft/Mars

Timestamp:

08/23/04 13:44:36 (20 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/findTelAxisFromStars.C (modified) (8 diffs)

trunk/MagicSoft/Mars/mtemp/findTelAxisFromStars.C

@@ -25 +25 @@
   MGeomCamMagic geomcam;
   MCameraDC     dccam;
-  MStarLocalCam starcam;
-
-  //$$$$$$$$$$$$$$$$ ww
+  MStarCam      starcam;
+
+  //$$$$$$$$$$$$$$$$ ww
+  MStarCam      sourcecam;
+  sourcecam.SetName("MSourceCam");
+
   MHTelAxisFromStars htelaxis;
   plist.AddToList(&htelaxis);
@@ -36 +39 @@
   plist.AddToList(&dccam);
   plist.AddToList(&starcam);
+  plist.AddToList(&sourcecam);
 
   // Reads the trees of the root file and the analysed branches
@@ -82 +86 @@
   MSourceDirections sdirs;
   sdirs.SetGeometryFile(geometryfile);
-  const Double_t ra  = MAstro::Hms2Rad(11, 4, 26);
+
+  const Double_t ra  = MAstro::Hms2Rad(11,  4, 26);
   const Double_t dec = MAstro::Dms2Rad(38, 12, 36);
   sdirs.SetRaDec(ra,dec);
   sdirs.AddDirection(ra,dec,1,"Mkn 421");
-  const Double_t ra  = MAstro::Hms2Rad(11, 4, 31);
+
+  const Double_t ra  = MAstro::Hms2Rad(11,  0, 50);
+  const Double_t dec = MAstro::Dms2Rad(39, 12, 44);
+  sdirs.AddDirection(ra,dec,1,"My_UMa 49");
+
+  const Double_t ra  = MAstro::Hms2Rad(11,  4, 31);
   const Double_t dec = MAstro::Dms2Rad(38, 14, 29);
   sdirs.AddDirection(ra,dec,1,"My_UMa 51");
+
+
   sdirs.SetRadiusFOV(3);
   //$$$$$$$$$$$$$$$$ ww
@@ -109 +121 @@
 
   //$$$$$$$$$$$$$$$$ ww
-  //Int_t InputType = 0;
+  Int_t InputType = 1;
   //findstars.SetUseCorrelatedGauss(kFALSE);
 
@@ -115 +127 @@
   //telaxis.FixRotationAngleAt(-1.0);
   //telaxis.FixScaleFactorAt(-1.0);
-  //telaxis.SetInputType(InputType);
+  telaxis.SetInputType(InputType);
 
   MFillH fillhisto("MHTelAxisFromStars[MHTelAxisFromStars]","");
-  //htelaxis.SetInputType(InputType);
+  htelaxis.SetInputType(InputType);
   //$$$$$$$$$$$$$$$$ ww
 
@@ -167 +179 @@
   while (tlist.Process())
     {
+      gLog << "---------------------------------------------" << endl;
+      gLog << "Macro : MStarPos content = " << endl;
       starcam.Print("namepossizchierr");
+      gLog << "---------------------------------------------" << endl;
+      gLog << "Macro : MSourcePos content = " << endl;
+      sourcecam.Print("namepossizchierr");
+      gLog << "---------------------------------------------" << endl;
 
       numevents++;
@@ -182 +200 @@
         // This is how to iterate over stars found:
         //        TIter Next(starlist);
-        //        MStarLocalPos* star;
+        //        MStarPos* star;
         //        UInt_t starnum = 0;
         //        cout << filename << " ";
         //        cout << "Iterating over list" << endl;
-        //        while ((star=(MStarLocalPos*)Next())) 
+        //        while ((star=(MStarPos*)Next())) 
         //          {
         //            cout << "star[" << starnum << "] ";
@@ -217 +235 @@
   tlist.PrintStatistics();
     
+  plist.Print();
+
   //$$$$$$$$$$$$$$$$ ww
   
   gLog << "Event loop finished; call DrawClone of MHTelAxisFromStars" << endl;
+  TObject *srccam = plist.FindObject("MSourceCam");
+
+  gLog << "srccam = " << srccam << endl;
+
   TObject *obj = plist.FindObject("MHTelAxisFromStars");
   if (obj)