Changeset 4097 for trunk/MagicSoft/Mars/mtemp/mmpi/macros

Timestamp:

05/18/04 13:39:53 (21 years ago)

Author:

mazin

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mmpi/macros

Files:

• callalphaplot.C (added)
• cleaninghillas.C (modified) (2 diffs)
• falsesourcemethod.C (modified) (1 diff)
• hillasfileread.C (modified) (2 diffs)
• skyplot.h (modified) (9 diffs)

trunk/MagicSoft/Mars/mtemp/mmpi/macros/cleaninghillas.C

@@ -75 +75 @@
 
     MObservatory observ;
+    pliston.AddToList(&observ);
+
     MPointingPos mpos;
     pliston.AddToList(&mpos);
@@ -105 +107 @@
     MSourcePosfromStarPos sourcefromstar;
 
-//    if (typeInput == "ON")
-//    {
-            sourcefromstar.AddFile("~mazin/mars/MarsCVS/Mars/data/positions2103.txt", 0);
-//    }
-//    else if (typeInput == "OFF")
-//    {
-//        sourcefromstar.AddFile("~mazin/data/dataMrk421_15feb/positionsFeb04_0", 0);
-//    }
+    if (typeInput == "ON")
+    {
+            sourcefromstar.AddFile("~mazin/mars/MarsCVS/Mars/positionsON", 0);
+    }
+    else if (typeInput == "OFF")
+    {
+        sourcefromstar.AddFile("~mazin/mars/MarsCVS/Mars/positionsOFF", 0);
+    }
 
 

trunk/MagicSoft/Mars/mtemp/mmpi/macros/falsesourcemethod.C

@@ -130 +130 @@
   if (TYPEOPTION == kTRUE)  // output of OFF data sample
 //      DATPARAMIN = "data/ressignOffMrk3004.dat";  // 3000
-      DATPARAMIN = "data/ressignOffCrab1502_1005.dat";   // 2000
+      DATPARAMIN = "paramOffCrab2701.dat";   // 2000
 
 
    HILLASFILENAME = "~/data/Crab/2004_01_27/CrabHillasOFF.root";
-   DATPARAMOUT = "data/paramOffCrab2701.dat";
+   DATPARAMOUT = "paramOffCrab2701.dat";
    ROOTPLOTNAME = "/.magic/magicserv01/scratch/Daniel/plots/RootPlots/Crab/2004_01_27/CrabOffSkyPlots.root";
    hillasfileread();

trunk/MagicSoft/Mars/mtemp/mmpi/macros/hillasfileread.C

@@ -90 +90 @@
     plist.AddToList(&mcevt); 
 
+    MPointingPos mpoint;
+    plist.AddToList(&mpoint); 
+
     MObservatory observ;
     plist.AddToList(&observ);
 
-//    printf("observ: %f\n",observ);
+    MHStarMap starmap;
+    plist.AddToList(&starmap);
 
-    MHStarMap starmap;
-//    plist.AddToList(&starmap); 
-    starmap.ReInit(&plist);
 //    MRawRunHeader header;
 //    plist.AddToList(&header);
@@ -131 +132 @@
      if (ROTOPTION == kTRUE) 
      {
-        theta = mcevt->GetTelescopeTheta();
-        phi = mcevt->GetTelescopePhi();
+//        theta = mcevt->GetTelescopeTheta();
+        theta = mpoint->GetZd();
+//        phi = mcevt->GetTelescopePhi();
+        phi = mpoint->GetAz();
 //   printf("theta: %5.3f, phi: %5.3f\n", theta*180./4.1415, phi*180./4.1415);
         observ.RotationAngle(theta, phi, singam(imnum), cosgam(imnum)); 

trunk/MagicSoft/Mars/mtemp/mmpi/macros/skyplot.h

@@ -67 +67 @@
 
 
-TH2F histskyplot("Sky Plot 1","Sky Plot of significance (2nd order fit)", 
+TH2F histskyplot("SPSignif2ndOrder","Sky Plot of significance (2nd order fit)", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -76 +76 @@
 histskyplot.UseCurrentStyle();
 
-TH2F histNexOnOff("Sky Plot 2","Sky Plot of Number of excess events (using ON and OFF data)", 
+TH2F histNexOnOff("SPNexOnOff","Sky Plot of Number of excess events (using ON and OFF data)", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -85 +85 @@
 histNexOnOff.UseCurrentStyle();
 
-TH2F histChi2Off("Sky Plot 3","Sky Plot of Chi2/ndf (2nd order fit, OFF data)", 
+TH2F histChi2Off("SPChi2PerNdf2ndOrderOff","Sky Plot of Chi2/ndf (2nd order fit, OFF data)", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -94 +94 @@
 histChi2Off.UseCurrentStyle();
 
-TH2F histskyplotOnOff("Sky Plot 4","Sky Plot of significance (using ON and OFF data)", 
+TH2F histskyplotOnOff("SPSignifOnOff","Sky Plot of significance (using ON and OFF data)", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -103 +103 @@
 histskyplotOnOff.UseCurrentStyle();
 
-TH2F histNex("Sky Plot 5","Sky Plot of Number of excess events (2nd order fit)",
+TH2F histNex("SPNex2ndOrder","Sky Plot of Number of excess events (2nd order fit)",
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -112 +112 @@
 histNex.UseCurrentStyle();
 
-TH2F histChi2("Sky Plot 6","Sky Plot of Chi2/ndf (2nd order fit)", 
+TH2F histChi2("SPChi2PerNdf2ndOrderOn","Sky Plot of Chi2/ndf (2nd order fit)", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -121 +121 @@
 histChi2.UseCurrentStyle();
 
-TH2F histskyplot4("Sky Plot 7","Sky Plot of significance (4th order fit)", 
+TH2F histskyplot4("SPSignif4thOrder","Sky Plot of significance (4th order fit)", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -130 +130 @@
 histskyplot4.UseCurrentStyle();
 
-TH2F histskyLiMa("Sky Plot 9","Sky Plot of significance (ON - OFF LiMa17)", 
+TH2F histskyLiMa("SPSignifLiMaOnOff","Sky Plot of significance (ON - OFF LiMa17)", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);
@@ -139 +139 @@
 histskyLiMa.UseCurrentStyle();
 
-TH2F histBerlin("Sky Plot Berlin","Sky Plot of Events below alpha 8 degrees", 
+TH2F histBerlin("SPBerlinMethod","Sky Plot of Events below alpha 8 degrees", 
                  NUMSTEPS, MINXGRID - 0.5*STEPGRID, MAXXGRID + 0.5*STEPGRID,
                  NUMSTEPS, MINYGRID - 0.5*STEPGRID, MAXYGRID + 0.5*STEPGRID);