StFmsPointMaker.h

Go to the documentation of this file.

// $Id: StFmsPointMaker.h,v 1.10 2019/06/26 16:49:43 akio Exp $
//
// $Log: StFmsPointMaker.h,v $
// Revision 1.10  2019/06/26 16:49:43  akio
// shower shape scaling for all shapes
//
// Revision 1.9  2018/03/08 20:07:27  akio
// initialization of mObjectCount and mMaxEnergySum was missing
//
// Revision 1.8  2018/03/02 20:26:44  akio
// Big update from Zhanwen Zhu with new shower shape and six z slices
//
// Revision 1.7  2015/11/05 17:53:09  akio
// Adding setScaleShowerShape() option for scaling up shower shape function for large cell
//
// Revision 1.6  2015/11/02 22:40:04  akio
// adding option for new cluster categorization
//
// Revision 1.5  2015/10/21 15:49:12  akio
// Adding 3 options to control how reconstruction works:
//   setGlobalRefit(int v=1)
//   setMergeSmallToLarge(int v=1)
//   setTry1PhotonFit(int v=1)
//
// Revision 1.4  2015/09/18 18:46:47  akio
// Move energy sum check for killing LED tail event to whole FMS, not each module
// Also make it not dependent on beam energy, so that it runs on simulation as well.
//
// Revision 1.3  2015/09/02 14:52:15  akio
// Adding readMuDst() to give options when reading back from mudst
//
#ifndef STROOT_STFMSPOINTMAKER_STFMSPOINTMAKER_H_
#define STROOT_STFMSPOINTMAKER_STFMSPOINTMAKER_H_

#include <map>
#include <vector>

#include "StMaker.h"
//#include "StRoot/StFmsUtil/StFmsGeometry.h"

class StFmsCollection;
class StFmsDbMaker;
class StFmsPoint;
class StMuDst;

namespace FMSCluster {
class StFmsFittedPhoton;
class StFmsTower;
class StFmsTowerCluster;
}  // namespace FMSCluster

class StFmsPointMaker : public StMaker {
 public:
  explicit StFmsPointMaker(const char* name = "StFmsPointMaker");
  ~StFmsPointMaker();
  Int_t InitRun(Int_t runNumber);
  Int_t Make();
  void Clear(Option_t* option = "");

  void setMaxEnergySum(Float_t v) {mMaxEnergySum=v;}

  void SetReadMuDst(int v=1) {mReadMuDst=v;} 
  
  /* Set this to perform gloal refit of all photons in a detector */
  void setGlobalRefit(int v=1) {mGlobalRefit=v;}

  /* Set this to cluster small and large cell together... experimental for now */
  void setMergeSmallToLarge(int v=1) {mMergeSmallToLarge=v;}

  /* set this to perform 1 photon fit if 2 photon fit failed */
  void setTry1PhotonFit(int v=1) {mTry1PhotonFitWhen2PhotonFitFailed=v;}

  /* set cluster categorization algo */
  void setCategorizationAlgo(int v=1) {mCategorizationAlgo=v;}

  /* set to scale shower shape for large cell, recomended values are*/
  //              large    small
  //old(1slice)    1.6      1.0
  //Zhanwen's      0.8      0.6
  //Yuxi's         1.2      0.8
  void setScaleShowerShape(float large=0.8, float small=0.6) {mScaleShowerShapeLarge=large; mScaleShowerShapeSmall=small;} 

  /* 0=original shower shape, 1=new shower shape with 6 z slices, 2=Yuxi's 6 slices */
  void setShowerShapeWithAngle(int v=1) {mShowerShapeWithAngle=v;} 

  /* 0=no vertex correction, 1=take vertex from Mudst BBC based on run11 calibration */
  void setVertexZ(int v=1) {mVertexZ=v;} 

 private:
  typedef std::vector<FMSCluster::StFmsTower> TowerList;
  typedef std::map<int, TowerList> TowerMap;
  StFmsPointMaker(const StFmsPointMaker&);
  StFmsPointMaker& operator=(const StFmsPointMaker&);
  StFmsCollection* getFmsCollection();
  int clusterEvent();
  int clusterDetector(TowerList* towers, int detectorId);

  //bool validateTowerEnergySum(const TowerList& towers) const;

  bool processTowerCluster(FMSCluster::StFmsTowerCluster* towerCluster,
                           int detectorId);
  StFmsPoint* makeFmsPoint(const FMSCluster::StFmsFittedPhoton& photon,
                           int detectorId);
  bool populateTowerLists();
  bool isValidChannel(int detector, int row, int col);
  StFmsDbMaker* mFmsDbMaker;  
  //FMSCluster::StFmsGeometry mGeometry;  //!< Access to current FMS geometry
  StFmsCollection* mFmsCollection; 
  TowerMap mTowers;  
  int mObjectCount=0;  

  Float_t mMaxEnergySum=255.0; 
   
  Int_t readMuDst();
  Int_t mReadMuDst=0;  

  Int_t mGlobalRefit=0;       
  Int_t mMergeSmallToLarge=1; 
  Int_t mTry1PhotonFitWhen2PhotonFitFailed=1; 
  Int_t mCategorizationAlgo=1;    
  Float_t mScaleShowerShapeLarge=0.8;  
  Float_t mScaleShowerShapeSmall=0.6;  
  Int_t mShowerShapeWithAngle=1;  
  Int_t mVertexZ=0;               
  Double_t vertexz=0.0;           

  StMuDst* muDst; 

  virtual const Char_t *GetCVS() const {static const Char_t cvs[]="Tag " __DATE__ " " __TIME__ ; return cvs;}
  ClassDef(StFmsPointMaker, 0)
};
#endif  // STROOT_STFMSPOINTMAKER_STFMSPOINTMAKER_H_