StFgtSlowSimuMaker.h

// $Id: StFgtSlowSimuMaker.h,v 1.3 2014/08/06 11:43:13 jeromel Exp $


#ifndef STAR_StFgtSlowSimuMaker
#define STAR_StFgtSlowSimuMaker

#ifndef StMaker_H
#include "StMaker.h"
#endif

#include "StFgtUtil/StFgtConsts.h"
#include <TVector3.h>

class St_g2t_fgt_hit;
class g2t_fgt_hit_st;
class TH2F;
class TH1F;

#include <TRandom3.h>
class StFgtDisc;
class StFgtDb;

class StFgtStripCollection;


//#define __FGT_QA_HISTO__
/* enable this switch to generate many QA histos. 
Note, you need to pass pointer to the valid  TObjArray top get them out of this maker.
In bfc.C you should have:
  StFgtSlowSimuMaker    *myMkSM =  get valid pointer....
   HList=new  TObjArray;
   myMkSM->setHList(HList); 
.....
after event loop save them:
  TString &fnameOut1=chain->GetFileOut();  TString fnameOut2=fnameOut1;
  fnameOut2.ReplaceAll(".root",".fgt");
  myMkSM->saveHisto(fnameOut2);// jan , save one maker is enough, all histos are merged
*/


struct fgt_g2t_auxil {
  TVector3 Rlab;
  TVector3 Rloc, Dloc; 
  g2t_fgt_hit_st *hitPtr; 
  Int_t iQuad; 
};


class StFgtSlowSimuMaker : public StMaker {   
 protected:

 static const Float_t pulseShape[20];

 private:


  enum {mxH=32};
  TH1 *hA[mxH];
  TObjArray *HList;

  
  StFgtDb *fgtDb;

  TH2F *quadDigitizationXY; 
  TH1F *quadDigitizationPhi; 
  TH1F *quadDigitizationRad; 

  TH2F *digXYAll;  
  TH2F *digPadcAll, *digRadcAll;
  TH1F *digPAll, *digRAll; 

  Int_t mInpEve;
  Int_t mEventId;


  Double_t   par_stripThreshAdc; 
  Double_t   par_XYamplSigma; 
 
  Int_t      par_cutoffOfBichel; 
  Double_t   par_pairsPerCm ; 
  Double_t   par_trackTOFcutoff; 
  Double_t   par_trackPcutoff; 

  Double_t   par_transDiffusionPerPath; 
  Int_t      par_binStep; 
  

  Double_t   par_2DampCutoffScale; 
  Double_t   par_overalGain; 
  Double_t   par_PplaneChargeFraction; 

  Int_t      par_badSetup;
  Int_t      switch_addPeds; 

  TRandom3* mRnd;

  void   addHit(TVector3 rLocal, Double_t  ampl=1.) ;
  void   responseMipModel(TVector3 Rloc, TVector3 Dloc);// detailed, w/ fluct
  vector<fgt_g2t_auxil> mG2tHitList[kFgtNumDiscs][kFgtNumQuads];

  void InitHisto1();
  void InitHisto2();
  void CloseHisto();

  void unpack_g2t_hits(St_g2t_fgt_hit *);
  void projectQuad2strips( Int_t idisc, Int_t iquad);
  void exportStripPlane2StEvent(TH1F *hInp, Int_t stripIdOffset, StFgtStripCollection *out);

 
  /* model transverse response in 1D
     Frank: However, from test beam data I calculate the RMS of each cluster, 
     and there I get something like 0.5 - 0.6 strips, which is around 350 um.
     all clusters (both in my diploma thesis and in my email) are 1D. A
     cluster is defined as something on ONE projection of the detector. So
     the FWHM and the sigma is all 1D. So you should not divide by Sqrt(2) to
     get the 1 D value.
  */

  Double_t  amplFunc(Double_t del_cm) {
    Double_t relD=del_cm/par_XYamplSigma;
    return exp(-relD*relD/2.);
  }


 public:
  StFgtSlowSimuMaker(const char *name="FgtSlowSimu");
  virtual       ~StFgtSlowSimuMaker();
  virtual Int_t  Init();
  virtual Int_t  InitRun(Int_t runNumber);
  virtual Int_t  Finish(); 
  virtual Int_t  Make();
  virtual void  Clear(Option_t *option="");
  void  setHList(TObjArray * x){HList=x;}
  void  setFgtDb(StFgtDb *x) { fgtDb=x;}
  void  enableAddPeds(){switch_addPeds=1;}
  void  setRndSeed(int x) { mRnd->SetSeed(x); } // force reproducible rnd sequence
  void  saveHisto(TString fname);

 
  virtual const char *GetCVS() const {
    static const char cvs[]="Tag $Name:  $ $Id: StFgtSlowSimuMaker.h,v 1.3 2014/08/06 11:43:13 jeromel Exp $ built " __DATE__ " " __TIME__ ; 
    return cvs;
  }
 private:
  
  ClassDef(StFgtSlowSimuMaker,0)   
};
    
#endif



// $Log: StFgtSlowSimuMaker.h,v $
// Revision 1.3  2014/08/06 11:43:13  jeromel
// Suffix on literals need to be space (later gcc compiler makes it an error) - first wave of fixes
//
// Revision 1.2  2012/06/20 18:32:40  avossen
// setting elec ids for strips now, implemented pulse shape over 7 timebins
//
// Revision 1.1  2012/06/06 20:35:09  jeromel
// Code  review closed (requested Anselm/Jan; reviewed Jonathan/Jason)
//