StPxlDigmapsSim.h

/*
 * Author: M. Mustafa
 *
 */
#ifndef STAR_StPxlDigmapsSim
#define STAR_StPxlDigmapsSim

#include "StPxlUtil/StPxlConstants.h"
#include "StPxlISim.h"
#include "TVector3.h"
class StPxlDb;
class StMcPxlHit;
class StMcTrack;
class TObject;
class TRandom3;
class TF1;
class DIGPlane;
class DIGADC;
class DIGTransport;
class DIGEvent;

class StPxlDigmapsSim: public StPxlISim
{
  public:

    StPxlDigmapsSim(const Char_t *name = "pxlDigmapsSim");
    virtual ~StPxlDigmapsSim();

    virtual int initRun(TDataSet const& calib_db, TObjectSet const* pxlDbDataSet, Int_t run);
    
    void setSeed(unsigned int seed) { mOwnRndSeed = seed; }

    virtual const char *GetCVS() const
    {
      static const char cvs[] = "Tag $Name:  $ $Id: StPxlDigmapsSim.h,v 1.3 2018/03/15 21:37:42 dongx Exp $ built " __DATE__ " " __TIME__ ;
      return cvs;
    }

  private:
    void fillDigmapsEvent(int, StMcPxlHit const*, DIGEvent&) const;
    void calculateIncidencePositions(int, StMcPxlHit const*, TVector3&, TVector3&) const;
    float calculateDepositedEnergy(float totalLength, float betagamma) const;
    float betaGamma(StMcTrack const*) const;
    double dEdxvsBGNorm(double* x, double* par);
    bool goodPixel(int sec, int lad, int sen, int ix, int iy) const;
    bool goodSensor(int sec, int lad, int sen) const;
    virtual int addPxlRawHits(StMcPxlHitCollection const& in, StPxlRawHitCollection& out);

    TRandom3* mRndGen;
    unsigned int mOwnRndSeed;
    
    DIGPlane* mDigPlane;
    DIGADC* mDigAdc;
    DIGTransport* mDigTransport;

    StPxlDb* mPxlDb;

    short mSensorGoodStatusMin;
    short mSensorGoodStatusMax;

    short mRowColumnGoodStatus;

    float mEnergyLandauMean;
    float mEnergyLandauSigma;
    double mScalePar[10]; // parameters to function mdEdxvsBGNorm
    
    float mResAddX;  // additional resolution contribution from mechanical vibration and calibration precision
    float mResAddZ;  // additional resolution contribution from mechanical vibration and calibration precision
    // These two offset values are initialized per run according to the resolution above, fixed for each sensor in each run
    float mOffsetX[kNumberOfPxlSectors][kNumberOfPxlLaddersPerSector][kNumberOfPxlSensorsPerLadder];
    float mOffsetZ[kNumberOfPxlSectors][kNumberOfPxlLaddersPerSector][kNumberOfPxlSensorsPerLadder]; 
    TF1* mdEdxvsBGNorm;  // dEdx vs. beta*gamma function normalized to MIP - for momentum and particle species depenence
    //
    short mHitEffMode;
    float mMomCut;      // momentum cut for momentum depenent efficiency
    float mHitEffInner; // single hit efficiency - 0.97 (best knowledge seen from ZF cosmic ray test)
    float mHitEffOuter; // single hit efficiency - 0.97 (best knowledge seen from ZF cosmic ray test)
};
#endif