StMuEpdHit.h

#ifndef StMuEpdHit_h
#define StMuEpdHit_h

#include <cstdlib>
#include "StObject.h"

//
// \class StMuEpdHit
// \author Mike Lisa
// \date 14 Jan 2018
// \brief Holds signals for tiles in STAR Event Plane Detector

/*************************************************
 * Hit Class for the STAR Event Plane Detector.
 * There is one "hit" for each (good) tile in the detector
 * 
 * Total size of this object is 10 bytes
 * This object contains only
 *  1) the id of the tile, compacted into a Short_t
 *  2) the QT information (ADC, TDC, TAC, hasTac flag, status flag)
 *     compacted into 32 bits (Int_t
 *  3) the gain-corrected energy loss, in units of the
 *     most probable value (MPV) of a single MIP in the tile,
 *     according to a Landau distribution.  Stored as Float_t
 *
 * The StMuEpdHit in the StMuDST is basically the same as
 *  the StEpdHit object in StEvent
 *
 * - Mike Lisa Jan 2018
 ************************************************/
class StEpdHit;

class StMuEpdHit : public StObject
{
public:

  StMuEpdHit();
  StMuEpdHit(Int_t position, Int_t tile, Short_t EW, Int_t ADC, Int_t TAC, Int_t TDC, bool hasTAC, Float_t nMIP, bool statusIsGood, Int_t truthId);

  StMuEpdHit(StEpdHit* epdHit);

  bool hasTac() const;
  Int_t  adc() const;
  Int_t  tac() const;
  Int_t  tdc() const;
  Short_t side() const;

  Short_t id() const;
  Int_t position() const;     // 1...12
  Int_t tile() const;         // 1...31
  Int_t qtData() {return mQTdata;}
  Float_t nMIP() {return mnMIP;}
  bool isGood() const;

  void setQTdata(Int_t packedData){mQTdata=packedData;}
  void setId(Short_t id){mId = id;}
  void SetnMIP(Float_t nMIP){mnMIP = nMIP;}

  void setIdTruth(Int_t id){mTruthId = id;}

  Int_t idTruth(){return mTruthId;}

protected:

  Short_t mId;

  Int_t mQTdata;

  Float_t mnMIP;

  Int_t mTruthId;

  ClassDef(StMuEpdHit, 1)
};

inline Short_t StMuEpdHit::side() const { return mId < 0 ? -1 : +1;}

inline Short_t StMuEpdHit::id() const { return mId; }
inline Int_t   StMuEpdHit::position() const { return std::abs(mId / 100); }
inline Int_t   StMuEpdHit::tile() const { return std::abs(mId % 100); }

inline Int_t   StMuEpdHit::adc() const { return mQTdata & 0x0FFF; }
inline Int_t   StMuEpdHit::tac() const { return (mQTdata >> 12) & 0x0FFF; }
inline Int_t   StMuEpdHit::tdc() const { return (mQTdata >> 24) & 0x001F; }
inline bool    StMuEpdHit::hasTac() const { return (mQTdata >> 29) & 0x1; }
inline bool    StMuEpdHit::isGood() const { return (mQTdata >> 30) & 0x1; }
#endif