Changeset 5827 for trunk/MagicSoft/Mars/mtemp

Timestamp:

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

Author:

rico

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mifae/programs

Files:

• makeHillas.cc (modified) (12 diffs)
• makehillas.datacard (modified) (4 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)
     {

trunk/MagicSoft/Mars/mtemp/mifae/programs/makehillas.datacard

@@ -11 +11 @@
 
 // Maximun number of (data) events to be processed)
-NEVENTS 5000
+NEVENTS 999999
 
 // data file directory
@@ -30 +30 @@
 // Enclose all conditions between brakets, like: "(x<y) && (z<5)"
 // (see MF description class for more details)
-// FILTER (MHillas.fLength<100) && (MHillas.fLength>50)
+FILTER (MNewImagePar.fNumUsedPixels<400)
 
 // Display flag 
@@ -61 +61 @@
 
 // Cleaning level (tail cut, boundary cut, number of rings, cleaning method - 1=standard, 2=democratic)
-CLEANLEVEL 2.5 2.0 1 1
+CLEANLEVEL 4.0 3.5 1 1
 
 //Island calculations
@@ -74 +74 @@
 ISLANDCLEAN 1 40
 
+// topology flag
+// 0: do not calculate topology parameters
+// 1: calculate topology parameters
+TOPFLAG 1
+