PndEmcDigi.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id:$
//
// Description:
//      EMC Digi.
//
//       Software developed for the BaBar Detector at the SLAC B-Factory.
// Adapted for the PANDA experiment at GSI              
//
// Author List:
//      Xiaorong Shi            Lawrence Livermore National Lab
//      Steve Playfer           University of Edinburgh
//      Stephen Gowdy           University of Edinburgh
// Dima Melnichuk - adaption for PANDA
//
// Copyright Information:
//      Copyright (C) 1994      Lawrence Livermore National Lab
//
//#pragma once
#ifndef PNDEMCDIGI_H
#define PNDEMCDIGI_H

#include "TObject.h"
#include "PndEmcTwoCoordIndex.h"
#include "PndEmcHit.h"
#include "FairTimeStamp.h"
#include "TVector3.h"
#include "TClonesArray.h"

class TVector3;
class PndEmcXtal;
class PndEmcSharedDigi;

class PndEmcDigi : public FairTimeStamp
{

        friend std::ostream& operator<< (std::ostream& out, PndEmcDigi& digi){
                out << "PndEmc Digi in crystal: " << digi.GetDetectorId() 
                        << " energy: " << digi.GetEnergy()
                        << " timestamp: "<<digi.GetTimeStamp();
                return out;
        }

        public:    
        //for time based simulation
        static TClonesArray* fDigiArrayTBD;
        static void InitDigiArrayTBD();

        PndEmcDigi();

        PndEmcDigi(Int_t trackid, Int_t id, Float_t energy, Float_t time, Int_t hitIndex=-1);
        PndEmcDigi(const PndEmcDigi&);

        virtual ~PndEmcDigi();    
        // Operators

        virtual bool operator==(const PndEmcDigi & otherDigi) const;
        virtual bool operator!=(const PndEmcDigi & otherDigi) const;
        virtual bool operator<( const PndEmcDigi & otherDigi) const;

        virtual bool equal(FairTimeStamp* data){
                PndEmcDigi* myDigi = dynamic_cast <PndEmcDigi*> (data);
                if (myDigi != 0){
                        if (fDetectorId == myDigi->GetDetectorId())
                                return true;
                }
                return false;
        }

        virtual void Print(const Option_t* opt ="") const;

        enum PositionMethod { surface, depth };

        void SetEnergy(Double32_t energy) { fEnergy     = energy ;};
        void SetTrackId(Int_t id)         { fTrackId    = id     ;};
        void SetDetectorId(Int_t id);

        virtual void AbsorbEnergy(PndEmcDigi& otherDigi);

        //check if 2 digis are neigbour or not (prelim for testing)
        bool isNeighbour(const PndEmcDigi* theDigi) const;


        virtual Double_t GetEnergy()  const;
        Int_t GetTrackId()    const { return fTrackId    ;};
        Int_t GetDetectorId() const { return fDetectorId ;};
        PndEmcTwoCoordIndex* GetTCI() const;
        Int_t GetThetaInt() const {return fThetaInd;};
        Int_t GetPhiInt() const {return fPhiInd;};
        Double_t GetTheta() const {return fTheta;};
        Double_t GetPhi() const {return fPhi;};
        Short_t GetModule()      const { return (fDetectorId/100000000);};
        Short_t GetRow()         const { return ((fDetectorId/1000000)%100);};
        Short_t GetCrystal()     const { return (fDetectorId%10000);};
        Short_t GetCopy()        const { return ((fDetectorId/10000)%100);};
        Short_t GetXPad()        const;
        Short_t GetYPad()        const;

        Int_t GetHitIndex() {return fHitIndex;}
        const TVector3 &where() const {return fWhere;};

        virtual PndEmcSharedDigi* dynamic_cast_PndEmcSharedDigi();
        virtual const PndEmcSharedDigi* dynamic_cast_PndEmcSharedDigi() const;

        static Double_t getRescaleFactor(){ return fRescaleFactor; };
        static Double_t getPositionDepthPWO(){ return fPositionDepthPWO; };
        static Double_t getPositionDepthShashlyk(){ return fPositionDepthShashlyk; };


        Int_t fEvtNo;

        protected:
        Double_t fEnergy;    // digi amplitude
        Int_t fTrackId;
        Int_t fDetectorId;
        Int_t fHitIndex; //  Index of hit which is converted to digi
        TVector3 fWhere;
        Int_t fThetaInd;
        Int_t fPhiInd;
        Double_t fTheta;
        Double_t fPhi;


        private:

        friend class PndEmcWaveformToDigi;
        friend class PndEmcWaveformToDigiFPGA;
        friend class PndEmcWaveformToCalibratedDigi;
        friend class PndEmcMultiWaveformToCalibratedDigi;
        friend class PndEmcFWEndcapDigi;
        friend class PndEmcTmpWaveformToDigi;
        friend class PndEmcMakeDigi;
        friend class PndEmcCorrBump;

        static void selectDigiPositionMethod( PositionMethod , 
                        double positionDepthPWO = 0., 
                        double positionDepthShahslyk = 0.,
                        double rescaleFactor = 1.);

        static TVector3 surfacePosition( const PndEmcXtal* xtal );
        static TVector3 depthPosition( const PndEmcXtal* xtal );

        static TVector3 ( *&algPointer() ) ( const PndEmcXtal* );

        static double fRescaleFactor;
        static double fPositionDepthPWO;
        static double fPositionDepthShashlyk;


        ClassDef(PndEmcDigi,5);
};

#endif //PndEmcDigi_H