ParisEvents.h

Go to the documentation of this file.

#ifndef _ParisEvents
#define _ParisEvents 1

#include "TObject.h"
#include "TClonesArray.h"

#include <list>

class PHit : public TObject
{       
public:
        enum EParticle { kUnknown = -1, kGamma, kGDR, kNeutron };

public:
        Double32_t fE; // energy
        Int_t fID;              // an ID
        Double32_t fX;  // X position
        Double32_t fY;  // Y position
        Double32_t fZ;  // Z position
        Double32_t fT;  // time of flight
        Int_t fFlag;    // a flag
        Int_t fUID;             // a universal ID

private:
        void Clear(Option_t *opt);

public:
        PHit();
        PHit(const PHit &);
        virtual ~PHit();        
                
        ClassDef( PHit , 1 ) // Paris Hit
};

class PEvent : public TObject
{
private:
        TClonesArray *fHits; //-> Collection of hits
        
        Double32_t fH;  // sum energy
        Int_t fK;               // fold

public:
        PEvent();
        virtual ~PEvent();      
        
        void CopyTo(std::vector <PHit *> &ordlist, Option_t *opt = "");

        // number of hits in that event
        Int_t GetNbHits() const 
                { return fHits->GetEntries(); }
        
        PHit *IsUID(Int_t uid, Int_t *which = 0x0)
        {
                TClonesArray &ar = *fHits;

                PHit *ishit = 0x0; Int_t entries = fHits->GetEntries();
                for (Int_t j = 0; j < entries; j++) {
                        PHit *hit = (PHit *)ar[j];
                        if ( hit->fUID == uid ) {
                                ishit = hit; 
                                if ( which )
                                        (*which) = j;
                                break;
                        }
                }
                return ishit;
        }
        
        PHit *GetHit(Int_t);
        PHit *AddHit();
        
        void AddHK(Double_t h, Int_t k = 1)
                { fH += h; fK += k; }
        
        void SetHK(Double_t h, Int_t k)
                { fH = h; fK = k; }
        
        Double_t GetH() const
                { return fH;}
        Int_t GetK() const
                { return fK; }
        
        void Clear(Option_t *opt);
        
        ClassDef( PEvent , 1 ) // Paris Event
};

class POpticalHit : public TObject
{       
public:
        Double32_t fX;  // X position of the last point of the track
        Double32_t fY;  // Y position of the last point of the track
        Double32_t fZ;  // Z position of the last point of the track
        Double32_t fTA; // arrival time for the last point of the track since event
        Double32_t fTL; // arrival time for the last point of the track since creation of the photon
        
        Double32_t fLength; // length of the trace
        
        Int_t fPrimaryID;       // Primary ID, in case more that one primary particle is giving scintillation
        Int_t fSecondaryID;     // detector ID where the optical photon is 
        Int_t fNbSteps;         // number of steps between start and stop
        
private:
        void Clear(Option_t *opt);
        
public:
        POpticalHit();
        POpticalHit(const POpticalHit &);
        virtual ~POpticalHit(); 
        
        ClassDef( POpticalHit , 1 ) // Paris Optical Hit
};

class POpticalEvent : public TObject
{
private:
        TClonesArray *fHits; //-> Collection of optical hits
        
public:
        POpticalEvent();
        virtual ~POpticalEvent();       

        // number of hits in that event
        Int_t GetNbHits() const 
                { return fHits->GetEntries(); }
        
        POpticalHit *GetHit(Int_t);
        POpticalHit *AddHit();
        
        void Clear(Option_t *opt);

        ClassDef( POpticalEvent , 1 ) // Paris Optical Event
        
};

#endif