Changeset 3334 for trunk/MagicSoft/Mars/macros

Timestamp:

02/26/04 16:42:03 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/starmcstereo.C (modified) (3 diffs)

trunk/MagicSoft/Mars/macros/starmcstereo.C

@@ -33 +33 @@
 /////////////////////////////////////////////////////////////////////////////
 
-void starmcstereo(Int_t ct1 = 2, Int_t ct2 = 3)
+//
+// User change. 
+//
+
+Float_t ctx[7] = {0., 0., 0., 0., 0., 0., 0.}; 
+Float_t cty[7] = {-70., -40., -30., 30., 50., 60., 70.}; // in meters
+//
+// FIXME: unfortunately present version of reflector was not prepared for
+// stereo configurations and keeps no track of CT position. So the positions
+// must be set above by the user, making sure that they correspond to the
+// files one is analysing.
+//
+
+void starmcstereo(Int_t ct1 = 1, Int_t ct2 = 2)
 {
-  Int_t CT[2] = {ct1, ct2};
-
+  if (ct1 > sizeof(ctx)/sizeof(ctx[0]) || 
+      ct2 > sizeof(ctx)/sizeof(ctx[0]) )
+    {
+      cout << endl << "Wrong CT id number!" << endl;
+      return;
+    }
+
+  Int_t CT[2] = {ct1, ct2};  // Only 2-telescope analysis for the moment
   Int_t NCTs = sizeof(CT)/sizeof(CT[0]);
 
@@ -47 +66 @@
   Float_t CleanLev[2] = {4., 3.}; // Tail cuts for image analysis
 
-  Int_t BinsHigh[2] = {0, 5}; // First and last FADC bin of the range to be integrated,
-  Int_t BinsLow[2]  = {0, 5}; // for high and low gain respectively.
+  Int_t BinsHigh[2] = {5, 9}; // First and last FADC bin of the range to be integrated,
+  Int_t BinsLow[2]  = {5, 9}; // for high and low gain respectively.
 
   // -------------------------------------------
@@ -204 +223 @@
   //
 
-  if (CT[0] == 1)
-    stereocalc.SetCT1coor(0.,0.);
-  else if (CT[0] == 2)
-    stereocalc.SetCT1coor(0.,0.);
-  else if (CT[0] == 3)
-    stereocalc.SetCT1coor(60.,60.);  // in meters
-  else if (CT[0] == 4)
-    stereocalc.SetCT1coor(60.,-60.);
-  else if (CT[0] == 5)
-    stereocalc.SetCT1coor(-60.,60.);
-  else if (CT[0] == 6)
-    stereocalc.SetCT1coor(-60.,-60.);
-  else
-    {
-      cout << "Unknown CT id!" << endl;
-      exit;
-    }
-
-  if (NCTs==2)
-    {
-      if (CT[1] == 1)
-        stereocalc.SetCT2coor(0.,0.);
-      else if (CT[1] == 2)
-        stereocalc.SetCT2coor(0.,0.);
-      else if (CT[1] == 3)
-        stereocalc.SetCT2coor(60.,60.);  // in meters
-      else if (CT[1] == 4)
-        stereocalc.SetCT2coor(60.,-60.);
-      else if (CT[1] == 5)
-        stereocalc.SetCT2coor(-60.,60.);
-      else if (CT[1] == 6)
-        stereocalc.SetCT2coor(-60.,-60.);
-      else
-        {
-          cout << "Unknown CT id!" << endl;
-          exit;
-        }
-
-      tlist.AddToList(&stereocalc);
-    }
+  stereocalc.SetCT1coor(ctx[CT[0]-1],cty[CT[0]-1]);
+  stereocalc.SetCT2coor(ctx[CT[1]-1],cty[CT[1]-1]);
+
+  tlist.AddToList(&stereocalc);
+
 
   MF filter1("{MMcEvt;1.fEvtNumber%2}<0.5");