Changeset 9567 for trunk/MagicSoft/Mars/msim/MHPhotonEvent.cc

Timestamp:

04/13/10 11:22:26 (15 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msim/MHPhotonEvent.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/msim/MHPhotonEvent.cc

@@ -38 +38 @@
 // Type 1:
 //   The maximum in x and y is determined from MCorsikaRunHeader
-//   (not yet implemented. Fixed to 25000)
+//   (not yet implemented. Fixed to 25000, 100 bins)
 //
 // Type 2:
 //   The maximum in x and y is determined from MReflector->GetMaxR();
+//   (100 bins)
 //
 // Type 3:
 //   The maximum in x and y is determined from MGeomCam->GetMaxR();
+//   (roughly 10 bins per pixel)
 //
 // Type 4:
-//   As type 3 but divided by 10.
-//
-// The binning is optimized using MH::FindGoodLimits. The number of bins
-// in 100 in the default case and 50 for type 3-4.
+//   Two times the traversal size of pixel #0 ( 4*MGeomCam[0].GetT() )
+//   (101 bins, units cm)
+//
+// Type 5:
+//   As type 3 but in cm instead of mm
+//
+// The binning is optimized using MH::FindGoodLimits.
 //
 // Fill expects a MPhotonEvent (the second argumnet in MFillH).
@@ -59 +64 @@
 //  + TH1D fHistH;
 //
+// Class Version 3:
+// ----------------
+//  + replaced TH1D by TH1F
+//  + replaced TH2D by TH2F
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -173 +182 @@
 
     MH::FindGoodLimits(num, num, min, max, kFALSE);
-
+    MH::FindGoodLimits(num, num, min, max, kFALSE);
     MBinning binsd, binsa, binsz;
     binsd.SetEdges(num, min, max);
@@ -213 +222 @@
             break;
         }
-    case 3: // The maximum radius
+    case 3: // The maximum radius (mm)
+        fHistXY.SetXTitle("X [mm]");
+        fHistXY.SetYTitle("Y [mm]");
+        fHistT.SetXTitle("X [mm]");
+        fHistT.SetYTitle("Y [mm]");
+        // *Fallthrough*
+
     case 4: // Two times the pixel-0 traversal size
+    case 5: // The maximum radius (cm)
         {
             MGeomCam *c = (MGeomCam*)pList->FindObject("MGeomCam");
@@ -222 +238 @@
                 return kFALSE;
             }
-            xmax = fType==3 ? c->GetMaxRadius()/10 : 2*(*c)[0].GetT()/10;
-            num  = 50;
+            // Type 3: Define ~10 bins per pixel
+            xmax = fType==3 ? c->GetMaxRadius() : 4*(*c)[0].GetT()/10;
+            num  = fType==3 ? TMath::Nint(10*(2*c->GetMaxRadius())/(2*(*c)[0].GetT())) : 101;
+
+            if (fType==5)
+                xmax /= 10;
 
             break;