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, Float_t local_time /*check of T res*/, Int_t hitIndex=-1, Int_t pileIndex=-1, Int_t DCid=-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 != nullptr) {
      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);
    void SetDCnumber(Int_t DCid) { fDCid = DCid ;}; 

    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 GetDCnumber() const; // for mapping into virtual Data Concentrators (DC)
    Int_t GetEventNr()    const { return fEvtNo    ;}; //

    Int_t GetHitIndex() {return fHitIndex;}
        Float_t GetLocalTime() {return flocal_time;}
        Int_t GetPileIndex() {return fPileIndex;}
    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:
        Float_t flocal_time;
    Double_t fEnergy;    // digi amplitude
    Int_t fTrackId;
    Int_t fDetectorId;
    Int_t fHitIndex; //  Index of hit which is converted to digi
        Int_t fPileIndex; // Index of pile-up: 1 - single hits, 2 - come from pile-up
    TVector3 fWhere;
    Int_t fThetaInd;
    Int_t fPhiInd;
    Double_t fTheta;
    Double_t fPhi;

    Int_t fDCid;

    private:

    friend class PndEmcWaveformToDigi;
    friend class PndEmcWaveformToDigiFPGA;
    friend class PndEmcWaveformToCalibratedDigi;
    friend class PndEmcMultiWaveformToCalibratedDigi;
    friend class PndEmcFWEndcapDigi;
    friend class PndEmcBWEndcapDigi;
    friend class PndEmcTmpWaveformToDigi;
    friend class PndEmcMakeDigi;
    friend class PndEmcMakeTbDigi;
    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,8);
};

#endif // PndEmcDigi_H