McIntegratingHit.h

Go to the documentation of this file.

#ifndef MCINTEGRATINGHIT_H
#define MCINTEGRATINGHIT_H

#include "TObject.h"
#include "TVector.h"
#include "TRefArray.h"

#include "commonRootData/idents/VolumeIdentifier.h"
#include "McParticle.h"

#include <vector>
#ifndef R__GLOBALSTL
#ifndef WIN32
using std::vector;
#else
using namespace std;
#endif
#endif


class McIntegratingHit: public TObject {
    
public:
    
    enum Particle{
        PRIMARY = 0,
        ELECTRON = 1,
        POSITRON = 2
    };

    McIntegratingHit();
    
    virtual ~McIntegratingHit();

    void Clear(Option_t * ="");
    void Fake( Int_t ievent, UInt_t rank, Float_t randNum ) ; // for tests
    Bool_t CompareInRange( const McIntegratingHit &, const std::string & name = "" ) const ; // for tests
    void Print(Option_t * ="") const;

    void initialize(const VolumeIdentifier& id);
    
    const VolumeIdentifier getVolumeId() const { return m_volumeId; };
    
    Double_t getTotalEnergy() const { return m_totalEnergy; };
    
    void addEnergyItem( const Double_t& energy, McParticle* t, const TVector3& pos );

    void setEnergyItems( const Double_t &totE, const Double_t *energyArr, 
        const TVector3& moment1, const TVector3& moment2);

    const McParticle* itemizedEnergyNext (Double_t &energy) const ;
    void itemizedEnergyReset() const { m_mapPtr = 0; };
    UInt_t itemizedEnergySize() const { return m_mcPartArr.GetEntries(); };

    Double_t getMcParticleEnergy(Particle p) const;

    const TVector3 getMoment1 () const;
    const TVector3 getMoment2 () const;

    void* operator new(size_t size);

    void* operator new(size_t size, void* vp);

    McIntegratingHit& operator=(const McIntegratingHit& rhs);
    
private:
    Double_t m_totalEnergy;
    Double_t m_energyArray[3];
    UInt_t m_packedFlags;
    VolumeIdentifier m_volumeId;
    TVector3 m_moment1Seed;
    TVector3 m_moment2Seed;
    
    TRefArray m_mcPartArr;
    vector<Double_t> m_energyPtrArr;
   
    //std::vector<Int_t> m_mcIdArr;

    //energyDepositMap m_energyItem; //!
    
    mutable UInt_t m_mapPtr; 

    ClassDef(McIntegratingHit,2)  // Monte Carlo Integrating Hit Class
};

#endif

---