/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Event/McTruth/McTruth-00-02-19/McTruth/McParticle.h

Go to the documentation of this file.

#ifndef Event_McParticle_H
#define Event_McParticle_H 

#include <iostream>
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/ContainedObject.h"
#include "GaudiKernel/SmartRef.h"
#include "GaudiKernel/SmartRefVector.h"
#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Geometry/Point3D.h"
#include "GaudiKernel/ObjectVector.h"
#include "GaudiKernel/ObjectList.h"
#include "EventModel/EventModel.h"
using namespace CLHEP;
extern const CLID& CLID_McParticle;

namespace Event {  // NameSpace


class McParticle  : virtual public ContainedObject  {
  public:
    typedef int  StdHepId;
    

    enum StatusBits{  
        DECAYED =1 ,    
        DECAYFLT=1<<1,  
        PRIMARY =1<<2,  
        LEAF = 1<<3,    
        ERROR = 1<<4    
        //MISSED=  1<<2,  //! Does not hit detector
        //NOINTER =1<<3,  //! Traverses detector w/o interacting 
        //STOPPED =1<<4,  //! Energy below cut; other bits may say why 
        //INTERACT=1<<5,  //! Interacted, no further decision to be made
        //INTSHDEP=1<<6,  //! Interacted, further decision depends on ! selection of shower deposition  
        //PRIMARY =1<<7,  //! primary particle 
        //SWERROR =1<<8,  //! Error occurred in swimming the track 
        //BCKSPL=1<<9,    //! The particle is backsplashed from the CAL back in the TKR region
        //POSHIT =1<<10,     //! The particle released a hit a PositionHit detector
        //NOTTRACK=1<<11, //! Not tracked by user request 
        //Swum =   1<<12,  //! this particle was produced by the swimmer
        //LEAF = 1<<13     //! this particle is a leaf in the particle tree 
    };
    

    virtual const CLID& clID() const   { return McParticle::classID(); }
    static const CLID& classID()       { return CLID_McParticle; }
    McParticle() :
      m_statusFlags(0)
    {}
    virtual ~McParticle() {}

    void init( McParticle* mother, 
        StdHepId id, 
        unsigned int statusBits,
        const HepLorentzVector& initialMomentum,
        const HepLorentzVector& finalMomentum,
        const HepLorentzVector& initialPosition,
        const HepLorentzVector& finalPosition,
        const std::string process = "");

    void initialize(
        StdHepId id, 
        unsigned int statusBits,
        const HepLorentzVector& initialMomentum,
        const HepLorentzVector& initialPosition,
        const std::string process = "");


    void initialize( McParticle* mother, 
        StdHepId id, 
        unsigned int statusBits,
        const HepLorentzVector& initialMomentum,
        const HepLorentzVector& initialPosition,
        const std::string process = "");

    //void finalize( const HepLorentzVector& finalMomentum,const HepLorentzVector& finalPosition);

    void finalize( const HepLorentzVector& finalPosition);

    StdHepId particleProperty() const;

    unsigned int statusFlags()const;

    void addStatusFlag(unsigned int flag){m_statusFlags |= flag;};

    bool primaryParticle() const;

    bool leafParticle() const;

    bool decayFromGenerator() const;

    bool decayInFlight() const;

    void setVertexIndex0(int index0) { m_vertexIndex0 = index0; }

    int  vertexIndex0() const {return m_vertexIndex0;}

    void setVertexIndex1(int index1) { m_vertexIndex1 = index1; }

    int  vertexIndex1() const {return m_vertexIndex1;}
    
    // Set track index
    void setTrackIndex(int trackIndex) { m_trackIndex = trackIndex; }
 
    // Get track index
    int trackIndex() const {return m_trackIndex;}

    const HepLorentzVector& initialPosition() const;
    const HepLorentzVector& finalPosition() const;

    const HepLorentzVector&  initialFourMomentum()const;
    //const HepLorentzVector&  finalFourMomentum()const;

    const McParticle& mother()const; 

    void setMother(const SmartRef<McParticle> m);

    void addDaughter(const SmartRef<McParticle> d){m_daughters.push_back(d);};

    void removeDaughter(const SmartRef<McParticle> mcPart);

    //const std::string getProcess()const{return m_process;};

    // void setInitialId(idents::VolumeIdentifier id){m_initialId = id;};
    // void setFinalId(idents::VolumeIdentifier id){m_finalId = id;};

    // idents::VolumeIdentifier getInitialId(){return m_initialId;};
    // idents::VolumeIdentifier getFinalId(){return m_finalId;};

    const SmartRefVector<McParticle>& daughterList()const{return m_daughters;};


  private:

    StdHepId                  m_particleID;

    int m_trackIndex;

    int m_vertexIndex0;
    
    int m_vertexIndex1;

    unsigned long             m_statusFlags;
    HepLorentzVector                 m_initialPosition;
    HepLorentzVector                 m_finalPosition;
    
    HepLorentzVector           m_initialFourMomentum;
    HepLorentzVector           m_finalFourMomentum;
    SmartRef<McParticle>       m_mother;
    SmartRefVector<McParticle> m_daughters;
    std::string                m_process;

    // idents::VolumeIdentifier   m_initialId;

    // idents::VolumeIdentifier   m_finalId;  
};

typedef ObjectList<McParticle>       McParticleCol;

} // NameSpace Event


#endif    // Event_McParticle_H

---