Changeset 4073 for trunk/MagicSoft/Mars/mtemp/mifae/macros

Timestamp:

05/14/04 16:45:04 (21 years ago)

Author:

jlopez

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mifae/macros

Files:

• findstars.C (modified) (4 diffs)
• psffit.C (modified) (2 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/macros/findstars.C

@@ -47 +47 @@
 
 
-void findstars(const TString filename="dc_2004_03_18_20_37_19_20663_CrabNebula.root", const TString directory="/nfs/magic/CaCodata/rootdata/CrabNebula/Period015/2004_03_19/", const UInt_t numEvents = 11)
+void findstars(const TString filename="dc_2004_03_19_00_36_50_20781_Mrk421.root", const TString directory="/nfs/magic/CaCodata/rootdata/Mrk421/Period015/2004_03_19/", const UInt_t numEvents = 0)
 {
 
@@ -80 +80 @@
 
   MGeomApply geomapl;
-
-//   const Int_t numrings = 3;
-  const Int_t numblind = 5;
-  const Short_t x[numblind] = {  8,  27, 224, 279, 339};
+  TString continuoslightfile = 
+    "/nfs/magic/CaCodata/rootdata/Miscellaneous/Period016/2004_04_16/dc_2004_04_16_04_46_18_22368_Off3c279-2CL100.root";
+  MCalibrateDC dccal;
+  dccal.SetFileName(continuoslightfile);
+  
+  const Int_t numblind = 12;
+  const Short_t x[numblind] = {  8,  27,
+                               507, 508, 509, 510, 511,
+                               543, 559, 560, 561, 567};
   const TArrayS blindpixels(numblind,(Short_t*)x);
   MFindStars findstars;
-//   findstars.SetNumRings(numrings);
-//   findstars.SetBlindPixels(blindpixels);
+  findstars.SetBlindPixels(blindpixels);
+  Float_t ringinterest = 100; //[mm]
+  Float_t tailcut = 3.5;
+  UInt_t integratedevents = 10;
+  findstars.SetRingInterest(ringinterest);
+  findstars.SetDCTailCut(tailcut);
+  findstars.SetNumIntegratedEvents(integratedevents);
+  findstars.SetMinuitPrintOutLevel(0);
 
   // prints
@@ -95 +106 @@
   tlist.AddToList(&geomapl);
   tlist.AddToList(&read);
+  tlist.AddToList(&dccal);
   //  tlist.AddToList(&pdc, "Currents");
   tlist.AddToList(&findstars, "Currents");
@@ -119 +131 @@
           return;
       
+       MHCamera display0(geomcam);
+       display0.SetPrettyPalette();
+       display0.Draw();
+//        display0.SetCamContent(dccal.GetDisplay());
+       display0.SetCamContent(dccal.GetDCCalibrationFactor());
+       gPad->Modified();
+       gPad->Update();
+
+       // Remove the comments if you want to go through the file
+       // event-by-event:
+       if (!HandleInput())
+         break;
+
        MHCamera display(geomcam);
        display.SetPrettyPalette();
        display.Draw();
        gPad->cd(1);
-//       starscam.Draw();
+       starcam.Draw();
       
-      while (tlist.Process())
-      {
-          display.SetCamContent(dccam);
-          gPad->Modified();
-          gPad->Update();
-          starcam.Print();
-          // Remove the comments if you want to go through the file
-          // event-by-event:
-          if (!HandleInput())
-              break;
+       UInt_t numevents=0;
+       
+       while (tlist.Process())
+         {
+           numevents++;
+           if (numevents%integratedevents==0)
+             {
+               //               display.SetCamContent(dccam);
+               display.SetCamContent(findstars.GetDisplay());
+               gPad->Modified();
+               gPad->Update();
+               starcam.Print("maxpossizechi");
+               // Remove the comments if you want to go through the file
+               // event-by-event:
+               if (!HandleInput())
+                 break;
+             }
       } 
 

trunk/MagicSoft/Mars/mtemp/mifae/macros/psffit.C

@@ -47 +47 @@
 
 
-void psffit(const TString filename="dc_2004_02_15_01_51_22_17286_Mrk421.root", const TString directory="/nfs/magic/CaCodata/2004_02_15/", const UInt_t numEvents = 100000)
+void psffit(const TString filename="dc_2004_03_19_00_36_50_20781_Mrk421.root", const TString directory="/nfs/magic/CaCodata/rootdata/Mrk421/Period015/2004_03_19/", const UInt_t numEvents = 0)
 {
 
@@ -130 +130 @@
           gPad->Modified();
           gPad->Update();
+          psffit.Print();
           // Remove the comments if you want to go through the file
           // event-by-event: