Changeset 5827 for trunk/MagicSoft/Mars/mtemp/mifae/programs/makeHillas.cc

Timestamp:

01/13/05 19:28:53 (20 years ago)

Author:

rico

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/programs/makeHillas.cc (modified) (12 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/programs/makeHillas.cc

@@ -61 +61 @@
 #include "MContinue.h"
 #include "MReportDrive.h"
+#include "MTopology.h"
+#include "MTopologyCalc.h"
 
 #include "TApplication.h"
@@ -93 +95 @@
 Float_t  ltail     = 1.5;
 Int_t    islflag   = 0;
+Int_t    topflag   = 0;
 Float_t  lnew      = 40;
 Int_t    kmethod   = 1;
@@ -207 +210 @@
           pedloop2.SetNoStorage();
           pedloop2.SetExtractor(extractor);
-          pedloop2.SetExtractorResolution();
           pedloop2.SetPedestals(pedloop1.GetPedestalCam());  
           pedloop2.SetInput(&pedcaliter);
@@ -374 +376 @@
   MIslands      isl2;
   MIslands      isl3;
+  MTopology     topology;
   
   isl.SetName("MIslands");  
@@ -401 +404 @@
   plist5.AddToList(&runhead);
   plist5.AddToList(&reportdrive);
+  plist5.AddToList(&topology);
   
   // cuts
@@ -435 +439 @@
   MHillasCalc       hcalc;
   MHillasSrcCalc    csrc1;
+
+  MTopologyCalc     topcalc;
   
   MContinue applycut(&cut);
@@ -453 +459 @@
   write.AddContainer("MIslands"       , "Parameters");
   write.AddContainer("MReportDrive"   , "Parameters");
+  write.AddContainer("MTopology",       "Parameters");
   
   if (islflag == 2) 
@@ -494 +501 @@
   if(filter.Length())
     tlist5.AddToList(&applycut);
+  if(topflag>0)
+    tlist5.AddToList(&topcalc);
   tlist5.AddToList(&write);
   if(display)
@@ -503 +512 @@
       tlist5.AddToList(disphillas);
     }
+
   
   // Create and setup the eventloop
@@ -635 +645 @@
         ifun >> sext >> hifirst >> hilast >> lofirst >> lolast >> wsize;
       
+      if(strcmp(word.Data(),"TOPFLAG")==0)
+        ifun >> topflag;
+
       if(strcmp(word.Data(),"ISLFLAG")==0)
         {
@@ -674 +687 @@
   cout << "Maximum number of events: " << nmaxevents << endl;
   if(filter.Length())
-    cout << "Applying rejection cut: " << filter << endl;
+    cout << "Applying selection cut: " << filter << endl;
   cout << "Output file name: " << outname << endl;
   if(display)
@@ -692 +705 @@
   if (islflag == 1 || islflag == 2)
     cout << "Island calcultation..." << "using algorithm #" << kalgorithm <<endl;
+  if (topflag > 0)
+    cout << "Topology parameters will be computed" << endl;
   if (islflag == 2)
     {