---

SubMCEvent.h

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/GlastEvent/GlastEvent/TopLevel/SubMCEvent.h,v 1.2 2001/04/18 19:07:26 igable Exp $
#ifndef LHCBEVENT_SUBMCEVENT_H
#define LHCBEVENT_SUBMCEVENT_H 1


// Include files
#include <iostream>
#include <vector>
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/StreamBuffer.h"
#include "GlastEvent/Utilities/RandomNumberSeed.h"
#include "GlastEvent/Utilities/ProcessingVersion.h"
#include "CLHEP/Geometry/Point3D.h"
#include "GlastEvent/Utilities/CLHEPStreams.h"
#include "GlastEvent/TopLevel/Definitions.h"


//------------------------------------------------------------------------------
//
// ClassName:   SubMCEvent
//  
// Description: Essential information of the sub- Monte Carlo event
//              (the Monte Carlo event can be assembled from several
//              of these sub-events)
//
// The class SubMCEvent uses the Class Library for HEP
// <A HREF="http://wwwinfo.cern.ch/asd/lhc++/clhep/manual/RefGuide/index.html">CLHEP</A>
//              
//              It contains:
//              - detector characteristics
//              - primary vertex
//              - Monte Carlo event weight
//              - random number seed
//              - processing version
//
// Author:      Pavel Binko
// Changes:     P.Binko 19/10/1999 : Formating of ASCII output
//
//------------------------------------------------------------------------------


class SubMCEvent                                                               {

public:
  SubMCEvent()                                                               { }
  ~SubMCEvent()                                                              { }

  SubMCEvent& operator= ( const SubMCEvent clone )                             {
    m_detectorCharacteristics = clone.detectorCharacteristics();
    m_processingVersion = clone.processingVersion();
    m_randomNumberSeed = clone.randomNumberSeed();
    m_primaryVertex = clone.primaryVertex();
    m_weight = clone.weight();
    return *this;
  }

  bool operator== ( const SubMCEvent& t ) const                                {
    return 
      detectorCharacteristics() == t.detectorCharacteristics()
      && processingVersion() == t.processingVersion()
      && randomNumberSeed() == t.randomNumberSeed()
      && primaryVertex() == t.primaryVertex()
      && weight() == t.weight();
  }

  long detectorCharacteristics () const                                        {
    return m_detectorCharacteristics;
  }
  void setDetectorCharacteristics (long value)                                 {
    m_detectorCharacteristics = value;
  }

  const HepPoint3D& primaryVertex () const                                     {
    return m_primaryVertex;
  }
  HepPoint3D& primaryVertex ()                                                 {
    return m_primaryVertex;
  }
  void setPrimaryVertex (const HepPoint3D& value)                              {
    m_primaryVertex = value;
  }

  double weight () const                                                       {
    return m_weight;
  }
  void setWeight (double value)                                                {
    m_weight = value;
  }

  const RandomNumberSeed& randomNumberSeed () const                            {
    return m_randomNumberSeed; 
  }
  void setRandomNumberSeed (const RandomNumberSeed& value)                     {
    m_randomNumberSeed = value;
  }

  const ProcessingVersion& processingVersion () const                          {
    return m_processingVersion; 
  }
  void setProcessingVersion (const ProcessingVersion& value)                   {
    m_processingVersion = value; 
  }

  friend StreamBuffer& operator<< ( StreamBuffer& s, const SubMCEvent& obj )   {
    return s
      << obj.m_detectorCharacteristics
      << obj.m_primaryVertex
      << obj.m_weight
      << obj.m_randomNumberSeed
      << obj.m_processingVersion;
  }
  friend StreamBuffer& operator>> ( StreamBuffer& s, SubMCEvent& obj )         {
    return s
      >> obj.m_detectorCharacteristics
      >> obj.m_primaryVertex
      >> obj.m_weight
      >> obj.m_randomNumberSeed
      >> obj.m_processingVersion;
  }


  friend std::ostream& operator<< ( std::ostream& s, const SubMCEvent& obj )   {
    return obj.fillStream(s);
  }  
  std::ostream& fillStream( std::ostream& s ) const;

private: 
  long                m_detectorCharacteristics;
  HepPoint3D          m_primaryVertex;
  double              m_weight;
  RandomNumberSeed    m_randomNumberSeed;
  ProcessingVersion   m_processingVersion;

};


//
// Inline code must be outside the class definition
//


inline std::ostream& SubMCEvent::fillStream( std::ostream& s ) const           {
  return s
    << "class SubMCEvent :\n"
    << "\n    Detector characteristics = "
    << GlastEventField( GlastEvent::field12 )
    << m_detectorCharacteristics
    << "\n    Primary vertex location  = ( "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_primaryVertex.x() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_primaryVertex.y() << ", "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_primaryVertex.z() << " )"
    << "\n    Monte Carlo event weight = "
    << GlastEventFloatFormat( GlastEvent::width, GlastEvent::precision )
    << m_weight
    << "\n    Random number seed       = " << m_randomNumberSeed
    << "\n    Processing version       = " << m_processingVersion;
}


#endif    // LHCBEVENT_SUBMCEVENT_H

---