StEEmcSlowMaker.h

// $Id: StEEmcSlowMaker.h,v 2.10 2014/08/06 11:43:04 jeromel Exp $

#ifndef STAR_StEEmcSlowMaker
#define STAR_StEEmcSlowMaker

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

class TObjArray;
class StEEmcDb;
class StMuEmcCollection;
class StEmcCollection;

#include "StEEmcUtil/EEmcGeom/EEmcGeomDefs.h"
#include "StEEmcUtil/EEfeeRaw/EEdims.h"   
#include "StEEmcUtil/database/EEmcDbItem.h"                                                                                            

class StEEmcSlowMaker : public StMaker {

private: 
  // parameters, fixed :                                                                                                                                     
  Float_t mMip2ene; // conversion from mips to  energy in GeV                                                                                                   
  Float_t mSig1pe; // width of the single photoelectron peak (in p.e.)                                                                                         
  Float_t mMip2pe[MaxSmdStrips]; // conversion from mip to p.e. from ANL cosmic ray measurements for SMD strips                                                                                      
  Float_t mPmip2ene[3]; // as above for pre- and post-shower elements                                                                                              
  Float_t mPmip2pe; // as above for pre- and post-shower elements                                                                                               

  Float_t avgNumPePerMip(Int_t stripID); // avg # p.e. per mip                                                                      

  enum Source_t {kMuDst, kStEvent};
  Source_t mSource;

  StEEmcDb *mEeDb;

  Int_t mNInpEve; // private event counter

  enum {maxHist = 32};
  TH1 *mHist[maxHist]; // some global (test) histograms  
  void InitHisto();

  Int_t MakeTower(StMuEmcCollection *emc);

  Int_t MakePrePost(StMuEmcCollection *emc);

  Int_t MakeSMD(StMuEmcCollection *emc);

  Int_t MakeTower(StEmcCollection* emc);

  Int_t MakePrePost(StEmcCollection* emc);

  Int_t MakeSMD(StEmcCollection* emc);

  Bool_t checkDBped(const EEmcDbItem *x);

  Float_t getPedSmear(Float_t sigPed);

  void setZeroAdc(StMuEmcCollection *emc);

  void setZeroAdc(StEmcCollection* emc);
  
  Bool_t mEnableSMD;
  Bool_t mEnablePrePost;
  Bool_t mEnableTower;

  // Offset by pedestal (default true)           
  Bool_t mAddPed;
  // Smear the pedestals (default true)         
  Bool_t mSmearPed;
  // Drop bad channels (default false)
  Bool_t mDropBad;
  // Overwrite muDst values(default true)
  Bool_t mOverwrite;

  Bool_t mIsEmbeddingMode;  
  // set separate defaults for use in BFC      
  Bool_t mIsBFC;                                

  enum {kPre1=0, kPre2, kPost, kNumberPrepost};

  Float_t mRelativeLightYield[kNumberPrepost]; /* <N p.e.>/MIP for pre1, pre2 and post   */
  Float_t mSamplingFraction;                   /* sampling fraction from the fast simu   */
  Float_t mSamplingFractionUser;               /* user-specified sampling fraction       */
  Float_t mTowerGains[kEEmcNumSectors];        /* tower gains from the fast simu         */
  Float_t mPrepostGains;                       /* pre/post ..                            */
  Float_t mSmdGains;                           /* smd ..                                 */
  Int_t   mMaxAdc;                             /* max ADC in fast simu                   */
  Bool_t  mDoLightYield;                       /* set false to disable tower pre/post relative light yield correction */
  Float_t mTowerGainFact[kEEmcNumSectors][kEEmcNumSubSectors][kEEmcNumEtas];
  Float_t mSmdGainFact[kEEmcNumSectors][kEEmcNumSmdUVs][kEEmcNumStrips];

  Int_t mTruncatePedSmear; 

  // Left unimplemented
  StEEmcSlowMaker(const StEEmcSlowMaker &);
  StEEmcSlowMaker &operator=(const StEEmcSlowMaker &);
  
public: 

  StEEmcSlowMaker(const Char_t *name = "EEmcSlowSim", const Char_t *muDstMakerName = 0);
  virtual ~StEEmcSlowMaker();

  virtual Int_t Init();
  virtual Int_t Make();
  void setEmbeddingMode(Bool_t e = true) {mIsEmbeddingMode = e;}

  void disableTower() {mEnableTower = false;}
  void disablePrePost() {mEnablePrePost = false;}
  void disableSMD() {mEnableSMD = false;}

  void setAddPed(Bool_t a = true) {mAddPed = a;}
  void setSmearPed(Bool_t s = true) {mSmearPed = s;}
  void setDropBad(Bool_t d = true) {mDropBad = true;}

  void setOverwrite(Bool_t o = true) {mOverwrite = o;}

  void setTruncatePedSmear(Float_t nSigma) {mTruncatePedSmear = nSigma;}
  
  void setMipElossSmd(Float_t e) {mMip2ene = e;}
  void setNpePerMipSmd(Int_t strip, Float_t npe) {mMip2pe[strip] = npe;}
  void setNpePerMipSmd(Float_t npe) {for (Int_t i = 0;i < MaxSmdStrips;i++) mMip2pe[i] = npe;}

  void setMipElossPre(Int_t layer, Float_t e) {mPmip2ene[layer] = e;}
  void setNpePerMipPre(Float_t npe) {mPmip2pe = npe;}

  void setSinglePeResolution(Float_t r) {mSig1pe = r;}

  void setSource(const Char_t* name);

  void setRelativeLightYield(Float_t pre1, Float_t pre2, Float_t post) {mRelativeLightYield[kPre1] = pre1; mRelativeLightYield[kPre2] = pre2; mRelativeLightYield[kPost] = post;}

  void setDoLightYield(Bool_t ly) {mDoLightYield = ly;}

  void setSamplingFraction(Float_t f) {mSamplingFractionUser = f;}

  void setTowerGainSpread(Float_t s, Float_t mean = 1.0);

  void setSmdGainSpread(Float_t s, Int_t sector, Int_t uv, Int_t strip_index);
  void setSmdGainSpread(Float_t s, Int_t strip_index);
  void setSmdGainSpread(Float_t s);

  static Float_t getMipdEdx(); 

  virtual const char *GetCVS() const {
    static const char cvs[]="Tag $Name:  $ $Id: StEEmcSlowMaker.h,v 2.10 2014/08/06 11:43:04 jeromel Exp $ built " __DATE__ " " __TIME__ ; 
    return cvs;
  }

  ClassDef(StEEmcSlowMaker, 3)
};

#endif

// $Log: StEEmcSlowMaker.h,v $
// Revision 2.10  2014/08/06 11:43:04  jeromel
// Suffix on literals need to be space (later gcc compiler makes it an error) - first wave of fixes
//
// Revision 2.9  2010/09/07 22:24:52  stevens4
// give access to MIP dE/dx to other makers
//
// Revision 2.8  2010/08/05 21:23:45  stevens4
// Update sampling fraction to 4.8%
//
// Revision 2.7  2010/08/03 02:20:40  stevens4
// final update from peer review
//
// Revision 2.6  2010/07/29 16:12:03  ogrebeny
// Update after the peer review
//
// Revision 2.5  2010/05/03 20:47:23  ogrebeny
// Some code cleanup
//
// Revision 2.4  2010/02/12 23:02:38  ogrebeny
// By the request of the photon group, added an option to shift EEMC gains in the slow simulator.
//
// Revision 2.3  2009/02/05 20:06:53  ogrebeny
// Changed StEEmcDbMaker -> StEEmcDb
//
// Revision 2.2  2008/04/11 14:37:17  jwebb
// Added options to disable operation of individual slow simulaor subsystems.
//
// Revision 2.1  2007/11/28 16:17:34  jwebb
// Added the following features:
//
// 1. User may specify the sampling fraction.
//
// 2. Tower and SMD gain spreads.
//
// Revision 2.0  2007/01/13 00:03:03  jwebb
// Upgrade of the slow simulator.  The following changes have been made:
//
// 1. Towers will always be masked out when a "fail" bit is set in the database.
//    Previously this only happened if pedestals were being added, smeared.
//
// 2. Tower, preshower and postshower ADC values will be simulated using the
//    GEANT energy loss stored in StEmcHit and StMuEmcHit.  Previously, ADC
//    values from the fast simulator were used and energy loss recovered
//    using gains, resulting in roundoff errors.   Note that towers still use
//    the old path for MuDst-based analysis.
//
// 3. Tower simulation now accounts for the different light yields provided
//    by the brighter scintillator and two-fiber readout in the preshower
//    and postshower layers.  Previously, only the difference in thickness
//    was accounted for by the GEANT simulation.
//
// Revision 1.4  2006/12/12 20:29:14  balewski
// added hooks for Endcap embedding
//
// Revision 1.3  2006/08/07 18:50:11  balewski
// added capabilty to run on StEvent, use se-method, see macros/ for example
//
// Revision 1.2  2005/09/23 01:30:11  jwebb
// Tower peds now added  if option is set.
//
// Revision 1.1  2004/12/15 17:02:56  balewski
// try 2
//