Changeset 6723 for trunk/MagicSoft/Mars/mtemp/mifae

Timestamp:

03/03/05 17:29:38 (20 years ago)

Author:

rico

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mifae

Files:

• Changelog (modified) (1 diff)
• programs/makeHillas.cc (modified) (8 diffs)
• programs/makehillas.datacard (modified) (2 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/Changelog

@@ -18 +18 @@
 
                                                  -*-*- END OF LINE -*-*-
+ 2005/03/03 Javier Rico
+    * programs/makeHillas.cc programs/makehillas.datacard
+      - read initial list of bad pixels from ascii file
+      - include MImagePar as output container
 
  2005/02/15 Javier Rico

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

@@ -98 +98 @@
 TString  filter;
 TString  psfilename("makehillas.ps");
+TString  bpfilename;
 MRunIter pedcaliter;
 MRunIter caliter;
@@ -202 +203 @@
   /* FIRST LOOP: PEDESTAL COMPUTATION FOR CALIBRATION */
   /****************************************************/
-
+  
   // If you want to exclude pixels from the beginning, read 
   // an ascii-file with the corr. pixel numbers (see MBadPixelsCam)  
   //badcam.AsciiRead("badpixels.dat");
-
+  MBadPixelsCam badcam;
+  MGeomCamMagic geomcam;
+  MGeomApply      geomapl;
+  badcam.InitSize(geomcam.GetNumPixels());
+
+  if(bpfilename.Length()>0)
+    {
+      ifstream fin("/home/emma/Mars/mtemp/mifae/programs/badpixels.dat");
+      badcam.AsciiRead((istream&)fin);
+      badcam.Print();
+    }
+  
+  
   MJPedestal pedloop1;
   MJPedestal pedloop2;
@@ -215 +228 @@
       pedloop1.SetInput(&pedcaliter);
       pedloop1.SetExtractor(extractor);
-      pedloop1.SetNoStorage();
+      pedloop1.SetBadPixels(badcam);
+     pedloop1.SetNoStorage();
       if (!pedloop1.Process())
         return;
@@ -226 +240 @@
           /***********************************************************/
           pedloop2.SetNoStorage();
-          pedloop2.SetExtractor(extractor);
-          pedloop2.SetPedestals(pedloop1.GetPedestalCam());  
+          pedloop2.SetExtractor(extractor); 
+          pedloop2.SetBadPixels(pedloop1.GetBadPixels());
+          pedloop2.SetPedestals(pedloop1.GetPedestalCam());  
           pedloop2.SetInput(&pedcaliter);
           
@@ -265 +280 @@
   MParList  plist3;
   MTaskList tlist3;
-  
-  // containers
-  MGeomCamMagic geomcam;
-  MGeomApply      geomapl;
-  
+    
   plist3.AddToList(&tlist3);
   plist3.AddToList(&geomcam);
@@ -547 +558 @@
   write.AddContainer("MTopology"      , "Parameters");
   write.AddContainer("MPointingPos"   , "Parameters");
+  write.AddContainer("MImagePar"      , "Parameters");
   if(rffilename.Length())
     write.AddContainer("MHadronness",     "Parameters");
@@ -714 +726 @@
       if(strcmp(word.Data(),"PSFILENAME")==0)
         ifun >> psfilename;
+
+      // excluded pixels file name
+      if(strcmp(word.Data(),"BADPIXELFILE")==0)
+        ifun >> bpfilename;
 
       // calibration flag
@@ -804 +820 @@
   if(display)
     cout << "Generating PS file: " << psfilename << endl;
+  if(bpfilename.Length())
+    cout << "Bad pixels will be read from " << bpfilename << endl;
   cout << "Signal Extractor: " << chext[sext] << " with bounds ("<<hifirst<<","<<hilast<<","<<lofirst<<","<<lolast<<"), window size " << wsize << endl;
   cout << "Calibration: ";

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

@@ -11 +11 @@
 
 // Maximun number of (data) events to be processed)
-NEVENTS 10000
+NEVENTS 9999999
 
 // data file directory
@@ -41 +41 @@
 // PS file name
 PSFILENAME makehillas.ps
+
+// If want to include a file with a list of bad pixels, add next word
+// followed by the name of the file (provide the whole path)
+BADPIXELFILE /home/emma/Mars/mtemp/mifae/programs/badpixels.dat
 
 // Signal extractor, higain first, higainlast, logainfirst, logainlast, window size