Changeset 9348 for trunk/MagicSoft/Mars/msim

Timestamp:

02/19/09 15:38:48 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msim

Files:

• MHPhotonEvent.cc (modified) (1 diff)
• MPhotonData.cc (modified) (5 diffs)
• MPhotonData.h (modified) (4 diffs)
• MSimAbsorption.cc (modified) (1 diff)

trunk/MagicSoft/Mars/msim/MHPhotonEvent.cc

@@ -258 +258 @@
     */
 
-    const Int_t f = fPermanentReset ? 2 : 1;
-
-    Double_t xmin = h->GetWavelengthMin();
-    Double_t xmax = h->GetWavelengthMax();
-    Int_t    num  = TMath::CeilNint((xmax-xmin)/(5*f));
-
-    MH::FindGoodLimits(num, num, xmin, xmax, kFALSE);
-
-    MBinning(abs(num), xmin, xmax).Apply(fHistWL);
+    const Int_t f = fPermanentReset ? 10 : 2;
+
+    Double_t xmin = h->GetWavelengthMin()-20;
+    Double_t xmax = h->GetWavelengthMax()+20;
+    Int_t    num  = TMath::CeilNint((xmax-xmin)/f);
+
+    MH::FindGoodLimits(num, num, xmin, xmax, kTRUE);
+
+    MBinning(abs(num), xmin-.5, xmax-.5).Apply(fHistWL);
 
     return kTRUE;

trunk/MagicSoft/Mars/msim/MPhotonData.cc

@@ -30 +30 @@
 //  Storage container to store Corsika events
 //
+// For details on the coordinate systems see our Wiki.
+//
 //   Version 1:
 //   ----------
@@ -97 +99 @@
 // --------------------------------------------------------------------------
 //
+// Return the square cosine of the Theta-angle == 1-CosU^2-CosV^2
+//
+Double_t MPhotonData::GetCosW2() const
+{
+    return 1 - GetSinW2();
+}
+
+// --------------------------------------------------------------------------
+//
+// Return the square sine of the Theta-angle == CosU^2+CosV^2
+//
+Double_t MPhotonData::GetSinW2() const
+{
+    return fCosU*fCosU + fCosV*fCosV;
+}
+
+// --------------------------------------------------------------------------
+//
 // return the cosine of the Theta-angle == sqrt(1-CosU^2-CosV^2)
 //
 Double_t MPhotonData::GetCosW() const
 {
-    return TMath::Sqrt(1 - fCosU*fCosU - fCosV*fCosV);
+    return TMath::Sqrt(GetCosW2());
+}
+
+// --------------------------------------------------------------------------
+//
+// return the sine of the Theta-angle == sqrt(CosU^2+CosV^2)
+//
+Double_t MPhotonData::GetSinW() const
+{
+    return TMath::Sqrt(GetSinW2());
 }
 
@@ -110 +139 @@
 Double_t MPhotonData::GetTheta() const
 {
-    return TMath::ACos(GetCosW());
+    return TMath::ASin(GetSinW());
 }
 
@@ -186 +215 @@
     {
         // FIXME: Could be done in MPhotonEvent::ReadCorsikaEvent
-        gLog << err << "ERROR - MPhotonData::FillCorsika: fNummPhotons not 1, but " << fNumPhotons << endl;
+        gLog << err << "ERROR - MPhotonData::FillCorsika: fNumPhotons not 1, but " << fNumPhotons << endl;
         gLog << "        This is not yet supported." << endl;
         return kERROR;
@@ -196 +225 @@
     //fCosU = f[3];  // cos to x
     //fCosV = f[4];  // cos to y
-    // x=east, y=north
+    // x=west, y=south
     fPosX =  f[2];  // [cm]
     fPosY = -f[1];  // [cm]

trunk/MagicSoft/Mars/msim/MPhotonData.h

@@ -49 +49 @@
     Float_t fWeight; //! A weight for external use
 
-protected:
-    virtual Int_t FillCorsika(Float_t f[7]);
-    virtual Int_t FillRfl(Float_t f[8]);
-
 public:
     MPhotonData(/*const char *name=NULL, const char *title=NULL*/);
@@ -64 +60 @@
     Float_t  GetCosV()  const { return fCosV; }
     Double_t GetCosW()  const;
+    Double_t GetSinW()  const;
+    Double_t GetCosW2()  const;
+    Double_t GetSinW2()  const;
     Double_t GetTheta() const;
 
@@ -111 +110 @@
     Int_t GetTag() const { return fTag; }
 
+    // TObject
     //void Clear(Option_t * = NULL);
     void Print(Option_t * = NULL) const;
@@ -124 +124 @@
         return 0;
     }
-/*
-    void FillRad(TH1 &hist, Float_t scale=1) const;
-    void FillRad(TH2 &hist, Double_t x, Float_t scale=1) const;
-    void Fill(TH2 &hist, Float_t scale=1) const;
-    void Fill(TH3 &hist, Double_t z, Float_t scale=1) const;
- */
+
+    // I/O
     Int_t ReadCorsikaEvt(istream &fin);
     Int_t ReadRflEvt(istream &fin);
 
+    Int_t FillCorsika(Float_t f[7]);
+    Int_t FillRfl(Float_t f[8]);
+
     ClassDef(MPhotonData, 1) //Container to store a cherenkov photon bunch from a CORSUKA file
 };
-/*
-class MPhotonDataWeighted : public MPhotonData
-{
-private:
-    Float_t fWeight; // A weight for external use
 
-public:
-    Float_t GetWeight() const { return fWeight; }
-
-    ClassDef(MPhotonData, 1)
-};
-*/
 #endif

trunk/MagicSoft/Mars/msim/MSimAbsorption.cc

@@ -264 +264 @@
 Bool_t MSimAbsorption::ReInit(MParList *pList)
 {
+    if (fUseTheta)
+        return kTRUE;
+
     MCorsikaRunHeader *h = (MCorsikaRunHeader*)pList->FindObject("MCorsikaRunHeader");
     if (!h)