Changeset 9308 for trunk/MagicSoft/Mars/msim

Timestamp:

02/09/09 11:30:54 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msim

Files:

• MHPhotonEvent.cc (modified) (6 diffs)
• MHPhotonEvent.h (modified) (2 diffs)
• MPhotonData.h (modified) (1 diff)
• MPhotonEvent.h (modified) (1 diff)

trunk/MagicSoft/Mars/msim/MHPhotonEvent.cc

@@ -54 +54 @@
 // Fill expects a MPhotonEvent (the second argumnet in MFillH).
 //
+//
+// Class Version 2:
+// ----------------
+//  + TH1D fHistH;
+//
+//
 /////////////////////////////////////////////////////////////////////////////
 #include "MHPhotonEvent.h"
@@ -113 +119 @@
     fHistWL.SetDirectory(NULL);
 
+    fHistH.SetName("Height");
+    fHistH.SetTitle("Production Height");
+    fHistH.SetXTitle("h [km]");
+    fHistH.SetYTitle("Counts");
+    fHistH.SetDirectory(NULL);
+
     // FIXME: Get this information from the corsika run-header
     MBinning(70, 275, 625).Apply(fHistWL);
+    MBinning(100, 0, 25).Apply(fHistH);
 }
 
@@ -259 +272 @@
         const MPhotonData &ph = (*evt)[idx];
 
+        if (ph.GetPrimary()==MMcEvtBasic::kNightSky)
+            continue;
+
         const Double_t x = ph.GetPosX();
         const Double_t y = ph.GetPosY();
@@ -265 +281 @@
         const Double_t t = ph.GetTime()-min;
         const Double_t w = ph.GetWavelength();
+        const Double_t h = ph.GetProductionHeight()/100000;
 
         //TVector3 dir = dat->GetDir3();
@@ -273 +290 @@
         fHistT.Fill(x, y, t);
         fHistWL.Fill(w);
+        fHistH.Fill(h);
     }
 
@@ -341 +359 @@
     gPad->SetGrid();
     fHistWL.Draw();
-}
+
+    pad->cd(6);
+    gPad->SetBorderMode(0);
+    gPad->SetGrid();
+    fHistH.Draw();
+}

trunk/MagicSoft/Mars/msim/MHPhotonEvent.h

@@ -23 +23 @@
     TProfile2D fHistT;
     TH1D       fHistWL;
+    TH1D       fHistH;
 
     Int_t      fType;
@@ -40 +41 @@
     void Paint(Option_t *o="");
 
-    ClassDef(MHPhotonEvent, 1) // Histogram to display the information of MPhotonEvents
+    ClassDef(MHPhotonEvent, 2) // Histogram to display the information of MPhotonEvents
 };
 

trunk/MagicSoft/Mars/msim/MPhotonData.h

@@ -83 +83 @@
     // Getter Others
     UShort_t GetWavelength() const { return fWavelength; }
+    Float_t GetProductionHeight() const { return fProductionHeight; }
     MMcEvtBasic::ParticleId_t GetPrimary() const { return fPrimary; }
 

trunk/MagicSoft/Mars/msim/MPhotonEvent.h

@@ -157 +157 @@
     Int_t GetMaxIndex() const { return fMaxIndex; }
 
-    Bool_t IsValid() const { return fTimeLast>=fTimeFirst; }
+//    Bool_t IsValid() const { return fTimeLast>=fTimeFirst; }
 
     ClassDef(MPhotonStatistics, 1)