---

McIntegratingHit.h

Go to the documentation of this file.

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

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


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



//namespace GlastEvent {  // NameSpace

// Posible rehash Ian
// Forward declarations
//class McParticle;
#include "GlastEvent/MonteCarlo/McParticle.h"

extern const CLID& CLID_McIntegratingHit;


class McIntegratingHit : virtual public ContainedObject {
  public:

    virtual const CLID& clID() const   { return McIntegratingHit::classID(); }
    static const CLID& classID()       { return CLID_McIntegratingHit; }
    // Here I needed to takeout SmartRef<McParticle> because the 
    // template library was not giving a compile error.
    typedef std::map<McParticle*,double>  energyDepositMap;

    McIntegratingHit() : m_packedFlags(0)
    {}
    ~McIntegratingHit(){}


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

    double totalEnergy() const;
    const HepPoint3D moment1 () const;
          HepPoint3D moment1 ();
    const HepPoint3D moment2 () const;
          HepPoint3D moment2 ();

    const energyDepositMap& itemizedEnergy() const;
          energyDepositMap& itemizedEnergy();
    void setEnergyItems( const energyDepositMap& value );
    void clearEnergyItems();
    void addEnergyItem( const double& energy, McParticle* t, const HepPoint3D& position );
    void addEnergyItem( const double& energy, SmartRef<McParticle> t, const HepPoint3D& position );

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

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

    virtual StreamBuffer& serialize( StreamBuffer& s );
    virtual StreamBuffer& serialize( StreamBuffer& s ) const;
    virtual std::ostream& fillStream( std::ostream& s ) const;    
       
  private:
      VolumeID                      m_volumeID;
      energyDepositMap              m_energyItem;
      double                        m_totalEnergy;
      HepPoint3D                    m_moment1seed;
      HepPoint3D                    m_moment2seed;
      unsigned long                 m_packedFlags;
};

  
// Definition of all container types of McIntegratingHit
template <class TYPE> class ObjectVector;
typedef ObjectVector<McIntegratingHit>     McIntegratingHitVector;
template <class TYPE> class ObjectList;
typedef ObjectList<McIntegratingHit>       McIntegratingHitList;

inline StreamBuffer& McIntegratingHit::serialize( StreamBuffer& s ) const
{
    ContainedObject::serialize(s);
    s
      << m_volumeID
      << m_totalEnergy
      << m_moment1seed
      << m_moment2seed
//    << m_energyItem(this)     // The operator<< has not implemented yet. FIX ME!!
      << m_energyItem.size();
    energyDepositMap::const_iterator it;
    for (it = m_energyItem.begin(); it != m_energyItem.end(); it++){
        s << it->first//(this)
          << it->second;
    }
    return s
      << m_packedFlags;
}


inline StreamBuffer& McIntegratingHit::serialize( StreamBuffer& s )
{
    energyDepositMap::size_type m_energyItem_size;
    ContainedObject::serialize(s);
    s
      >> m_volumeID
      >> m_totalEnergy
      >> m_moment1seed
      >> m_moment2seed
//    >> m_energyItem(this)     // The operator<< has not implemented yet. FIX ME!!
      >> m_energyItem_size;
    m_energyItem.clear();
    energyDepositMap::size_type i;
    for (i = 0; i < m_energyItem_size; i++){//for (i = 0; i < m_energyItem_size; i++){ //Nasté
        SmartRef<McParticle> first;
        double               second;
        s >> first(this)
          >> second;
        m_energyItem[first] = second;
    }
        return s
      >> m_packedFlags;
}


inline std::ostream& McIntegratingHit::fillStream( std::ostream& s ) const
{
    s << "class McCaloHitBase :"
      << "\n    Deposited Energy        = "
      << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
      << m_totalEnergy
      << "\n    First moment (x, y, z)  = "
      << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
      << m_moment1seed.x() / m_totalEnergy << ", "
      << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
      << m_moment1seed.y() / m_totalEnergy << ", "
      << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
      << m_moment1seed.z() / m_totalEnergy << " )"
      << "\n    Second moment (x, y, z) = "
      << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
      << m_moment2seed.x() / m_totalEnergy << ", "
      << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
      << m_moment2seed.y() / m_totalEnergy << ", "
      << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
      << m_moment2seed.z() / m_totalEnergy << " )"
      << "\n    Volume ID               = " << m_volumeID;
    return s;
}

#endif // GlastEvent_McIntegratingHit_H

---