DigiPhoton.h

Go to the documentation of this file.

// $Id: DigiPhoton.h,v 1.3 2005/02/01 19:09:19 tagg Exp $
//
// A single photon in the scintillator system. The coordinate system
// should always be the local coordinate frame. For example, in a
// scintillator strip x should be the distance along the width of the
// strip, y should be the the distance along the thickness of the
// strip and z should be the distance along the length.
//
// messier@indiana.edu
//
// Fundamental redesign by n.tagg1@physics.ox.ac.uk (I AM my email address.)
#ifndef DIGIPHOTON_H
#define DIGIPHOTON_H
#include "TObject.h"
#include "TRef.h"
#ifndef PLEXSTRIPENDID_H
# include "Plex/PlexStripEndId.h"
#endif
#ifndef MUNITS_H
# include "Conventions/Munits.h"
#endif

#include "TVector3.h"
#include "TMath.h"

class DigiScintHit;

class DigiPhoton : public TObject
{
public:
  DigiPhoton();
  DigiPhoton( const DigiScintHit* hit,
             double wavelength, double t, 
             double x, double y, double z,
             double u, double v, double w);

  DigiPhoton(const DigiScintHit* hit,
             double wavelength, double t,
             const TVector3& pos,
             const TVector3& dir );

  // Getters:
  const DigiScintHit*  ParentHit()     const { return fHit; };

  Double_t             WaveLength()  const { return fWavelength; };
  Double_t             E()           const;
  Double_t             T()           const { return fT;   }

  const TVector3&      GetPos()      const { return fPos; };
  Double_t             X()           const { return fPos.X();   }
  Double_t             Y()           const { return fPos.Y();   }
  Double_t             Z()           const { return fPos.Z();   }

  const TVector3&      GetDir()      const { return fDir; };
  Double_t             GetCosX()     const { return fDir.X(); };
  Double_t             GetCosY()     const { return fDir.Y(); };
  Double_t             GetCosZ()     const { return fDir.Z(); };
  Double_t             U()           const { return fDir.X(); };
  Double_t             V()           const { return fDir.Y(); };
  Double_t             W()           const { return fDir.Z(); };
  Double_t             DxDz()        const { return fDir.Z()>0?(fDir.X()/fDir.Z()):999.;}
  Double_t             DyDz()        const { return fDir.Z()>0?(fDir.Y()/fDir.Z()):999.;}

  // Setters:
  void  SetParentHit( const DigiScintHit* hit )          { fHit = hit; };
  void  SetWavelength( double l )                        { fWavelength = l; };
  void  SetTime( double t )                              { fT = t; };
  void  SetPosition(const TVector3& in)                  { fPos = in; };
  void  SetPosition(Double_t x, Double_t y, Double_t z)  { fPos.SetXYZ(x,y,z); };
  void  SetDirection(const TVector3& in)                 { fDir = in; fDir.SetMag(1.0); };
  void  SetDirection(Double_t u, Double_t v, Double_t w) { fDir.SetXYZ(u,v,w); };

  virtual void Print( Option_t* option="") const;

public:
  const DigiScintHit* fHit;  // The parent energy deposition
  
  Double_t  fWavelength; // in Munits.
  Double_t  fT;    // Time at position
  TVector3  fPos;  // Position in strip.
  TVector3  fDir;  // Direction.
    
  ClassDef(DigiPhoton,2) // A photon in the scintillator system
};


inline Double_t DigiPhoton::E() const
{
  // E = hc/lambda
  // H() in Js, C() in m/s, fWavelength in meters, convert J to energy standard..
  return  (TMath::H() * TMath::C() / fWavelength) * TMath::Qe() / Munits::eV;  
}

#endif // DIGIPHTON

---