StPicoEpdHit.h

#ifndef StPicoEpdHit_h
#define StPicoEpdHit_h

// C++ headers
#include <cstdlib>

// ROOT headers
#include "TObject.h"

//
// \class StPicoEpdHit
// \author Mike Lisa
// \date 6 March 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
 *   that has a nonzero ADC value
 * 
 * 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 StPicoEpdHit in the StMuDST is basically the same as
 *  the StMuEpdHit object in the muDst and
 *  the StEpdHit object in StEvent
 * Except 
 *   1) it does not inherit from StObject, so can be used in "vanilla root"
 *   2) there is no "Truth ID" which is used for simulations.
 * 
 *
 * To access GEOMETRICAL information (e.g. location of a tile, etc)
 *  use the StEpdGeom class.  All you have to do is send StEpdGeom
 *  the StPicoEpdHit::id(), and you'll get all you want.
 *  (Also works for StEpdHit::id() and StMuEpdHit::id())
 *
 * - Mike Lisa March 2018
 ************************************************/

class StPicoEpdHit : public TObject {

 public:

  StPicoEpdHit();
  StPicoEpdHit(Int_t position, Int_t tile, Int_t EW, Int_t ADC, Int_t TAC, Int_t TDC,
               Bool_t hasTAC, Float_t nMIP, Bool_t statusIsGood);
  StPicoEpdHit(Short_t id, Int_t QTdata, Float_t nMIP);
  StPicoEpdHit(const StPicoEpdHit &hit);
  virtual ~StPicoEpdHit();
  virtual void Print(const Char_t *option = "") const;

  Bool_t hasTac() const             { return (mQTdata >> 29) & 0x1; }
  Int_t  adc() const                { return mQTdata & 0x0FFF; }
  Int_t  tac() const                { return (mQTdata >> 12) & 0x0FFF; }
  Int_t  tdc() const                { return (mQTdata >> 24) & 0x001F; }
  Short_t side() const              { return mId < 0 ? -1 : +1;}

  Short_t id() const                { return mId; }
  Int_t position() const            { return std::abs(mId / 100); }
  Int_t tile() const                { return std::abs(mId % 100); }
  Int_t row() const                 { return (std::abs(mId % 100))/2 + 1; }
  Int_t qtData() const              { return mQTdata; }
  Float_t nMIP() const              { return this->isGood() ? mnMIP : 0.0; }
  Bool_t isGood() const             { return (mQTdata >> 30) & 0x1; }
  Float_t TnMIP(float MAX=2.0, float thresh=0.3) const   {float nm = this->nMIP(); if (nm<thresh) return 0.0; else return nm<MAX ? nm : MAX;}

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

 protected:

  Short_t mId;

  Int_t mQTdata;

  Float_t mnMIP;

  ClassDef(StPicoEpdHit, 1)
};

#endif