---

McPositionHit.h

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/GlastEvent/GlastEvent/MonteCarlo/McPositionHit.h,v 1.6 2001/04/18 19:07:25 igable Exp $
#ifndef GlastEvent_McPositionHit_H
#define GlastEvent_McPositionHit_H 1

// If you wish to introduce the namespace `GlastEvent', uncomment
// the lines commented as `NameSpace'.


// Include files
#include <iostream>
#include <math.h>
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/ContainedObject.h"
#include "GaudiKernel/SmartRef.h"
#include "GlastEvent/TopLevel/Definitions.h"
#include "CLHEP/Geometry/Point3D.h"
#include "GlastEvent/Utilities/CLHEPStreams.h"
// Include all Gaudi container types here
//   to simplify inlude statements in algorithms
#include "GaudiKernel/ObjectVector.h"
#include "GaudiKernel/ObjectList.h"
#include "GlastEvent/Utilities/VolumeID.h"


#include "GlastEvent/MonteCarlo/McParticle.h"
#include "GlastEvent/MonteCarlo/McConstants.h"



//namespace GlastEvent { // NameSpace

// Forward declarations
class McParticle;

extern const CLID& CLID_McPositionHit;

class McPositionHit : virtual public ContainedObject {
  public:

    virtual const CLID& clID() const   { return McPositionHit::classID(); }
    static const CLID& classID()       { return CLID_McPositionHit; }
    McPositionHit()
      : m_depositedEnergy(0.),
        m_timeOfFlight(0.)
    {}
    virtual ~McPositionHit() { }


    const VolumeID volumeID() const;
    void setVolumeID( VolumeID value );

    const HepPoint3D& entryPoint() const;
          HepPoint3D& entryPoint();
    void setEntryPoint( const HepPoint3D& value );

    const HepPoint3D& exitPoint() const;
          HepPoint3D& exitPoint();
    void setExitPoint( const HepPoint3D& value );

    double depositedEnergy() const;
    void setDepositedEnergy( double value );
    double particleEnergy() const;
    void setParticleEnergy( double value );

    bool primaryOrigin() const;
    void setPrimaryOrigin( bool value );
    bool caloShowerOrigin() const;
    void setCaloShowerOrigin( bool value );

    bool needDigi() const;
    void setNeedDigi( bool value );

    double directionCosine() const;

    double timeOfFlight() const;
    void setTimeOfFlight( double value );

    const McParticle* mcParticle() const;
          McParticle* mcParticle();
    void setMcParticle( McParticle* value );
    void setMcParticle( SmartRef<McParticle> value );

    const McParticle* originMcParticle() const;
          McParticle* originMcParticle();
    void setOriginMcParticle( McParticle* value );
    void setOriginMcParticle( SmartRef<McParticle> value );

    virtual StreamBuffer& serialize( StreamBuffer& s ) const;
    virtual StreamBuffer& serialize( StreamBuffer& s );
    virtual std::ostream& fillStream( std::ostream& s ) const;

  private:
    VolumeID                m_volumeID;
    HepPoint3D              m_entry;
    HepPoint3D              m_exit;
    double                  m_depositedEnergy;
    double                  m_particleEnergy;
    double                  m_timeOfFlight;
    SmartRef<McParticle>    m_mcParticle;
    SmartRef<McParticle>    m_originMcParticle;
    unsigned long           m_packedFlags;
};


typedef ObjectVector<McPositionHit> McPositionHitVector;
typedef ObjectList<McPositionHit>   McPositionHitList;

inline StreamBuffer& McPositionHit::serialize( StreamBuffer& s ) const
{
  ContainedObject::serialize(s);
  return s
    << m_volumeID
    << m_entry
    << m_exit
    << m_depositedEnergy
    << m_particleEnergy
    << m_timeOfFlight
    << m_mcParticle(this)
    << m_originMcParticle(this)
    << m_packedFlags;
}


inline StreamBuffer& McPositionHit::serialize( StreamBuffer& s )
{
  ContainedObject::serialize(s);
  return s
    >> m_volumeID
    >> m_entry
    >> m_exit
    >> m_depositedEnergy
    >> m_particleEnergy
    >> m_timeOfFlight
    >> m_mcParticle(this)
    >> m_originMcParticle(this)
    >> m_packedFlags;
}


inline std::ostream& McPositionHit::fillStream( std::ostream& s ) const
{
  return s
    << "    base class McPositionHit :"
    << "\n        Volume ID             = " << m_volumeID
    << "\n        Entry point (x, y, z) = ( "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_entry.x() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_entry.y() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_entry.z() << " )"
    << "\n        Deposited Energy      = "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_depositedEnergy
    << "\n        Particle Energy       = "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_particleEnergy
    << "\n        Time of flight        = "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_timeOfFlight
    << "\n        Exit point (x, y, z)  = ( "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_exit.x() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_exit.y() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_exit.z() << " )"
    << "\n        McParticle            = " << m_mcParticle(this)
    << "\n        ancestor McParticle   = " << m_originMcParticle(this);
}

#endif    // GlastEvent_McPositionHit_H

---