StFmsFilterMaker.h

// $Id: StFmsFilterMaker.h,v 1.1 2018/05/22 19:55:42 akio Exp $

#ifndef STAR_StFmsFilterMaker
#define STAR_StFmsFilterMaker

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

class TH1F;
class StFmsDbMaker;
class StFmsCollection;

//class FmsGeomSimple;
//class StFmsTower;

class StFmsFilterMaker : public StMaker {
 private:
  enum {kUnknown=1, kUseFixedVertex, kUseRecoVertex};
  int myMode;
  int nInpEve, nRecVert,nZverOK, nAccEve;
  TH1F * mH0;
  //  FmsGeomSimple*  mGeomE;
  //  Float_t transverseNRG(Float_t vertexPosZ, StFmsTower * tower);
  Int_t triggerCondition();// Float_t vertexPosZ, StFmsTower *highTow,  Float_t &patchEt);
  float par_E_thres, par_Z0_vert, par_delZ_vert;


  StFmsDbMaker* mFmsDbMaker;
  StFmsCollection* mFmsColl;

  
 public: 
  void setEThres(float x) { par_E_thres=x;}  
  void setZvertCut(float z0, float dz) {  par_Z0_vert=z0; par_delZ_vert=dz;myMode=kUseRecoVertex;}
  void setFixedVertex(float z0) {  par_Z0_vert=z0; myMode=kUseFixedVertex;}

  StFmsFilterMaker(const char *name="FmsFilter");
  virtual       ~StFmsFilterMaker();
  virtual Int_t Init();
  virtual Int_t  Make();

  // virtual Int_t InitRun  (int runumber){return 0;}; // Overload empty StMaker::InitRun 
  virtual Int_t FinishRun(int runumber);

  virtual const char *GetCVS() const {
    static const char cvs[]="Tag $Name:  $ $Id: StFmsFilterMaker.h,v 1.1 2018/05/22 19:55:42 akio Exp $ built " __DATE__ " " __TIME__ ; 
    return cvs;
  }

  ClassDef(StFmsFilterMaker,0)  
};

#endif


/* ===================================

#if 1 
  gSystem->Load("StFmsA2EMaker");
  gSystem->Load("StFmsPoolFmsFilterMaker");

  StFmsA2EMaker *EEa2eMK=new StFmsA2EMaker("EE_A2E");
  EEa2eMK->database("eeDb");   // sets db connection
  EEa2eMK->source("StEventMaker",2);
  EEa2eMK->scale(1.0);      // scale reco Endcap energy by a factor

  //  this filter-1 _assumes_ fixed Zvertex,
  //    needs only Endcap slow simu 
  StFmsFilterMaker *eeFlt1Mk=new StFmsFilterMaker("FmsFilterPreTPC"); 
  eeFlt1Mk->setEtThres(14.);// (GeV), event-eta used 
  eeFlt1Mk->setFixedVertex(-60.); // (cm), Z0

  //  this filter-2 requires reco vertex
  StFmsFilterMaker *eeFlt2Mk=new StFmsFilterMaker("FmsFilterPostTPC"); 
  eeFlt2Mk->setEtThres(16.);// (GeV), event-eta used 
  eeFlt2Mk->setZvertCut(-60.,15.); // (cm), Z0, delatZ

  // position makers at the right location in the chain
  // WARN, order of operations below is important
  StMaker *tpcMk = chain->GetMaker("tpcChain");
  chain->AddAfter("eefs",tpcMk); // now TPC-chain is after EMC-makers
 
  chain->AddAfter("eefs",eeFlt1Mk);
  chain->AddAfter("eefs",EEa2eMK); // filter 1 is before TPC-chain
  chain->AddBefore("MuDst",eeFlt2Mk); // filter 2 is after vertex finder
  chain->ls(3);
 #endif

 
  ========================== */
// $Log: StFmsFilterMaker.h,v $
// Revision 1.1  2018/05/22 19:55:42  akio
// Initial version of after GEANT filter for FMS by Mriganka
//
// Revision 1.5  2014/08/06 11:42:57  jeromel
// Suffix on literals need to be space (later gcc compiler makes it an error) - first wave of fixes
//
// Revision 1.4  2008/05/12 15:33:01  balewski
// *** empty log message ***
//
// Revision 1.3  2008/05/11 18:49:18  balewski
// merged 2 makers, now one can us it before and after TPC tracking, run 2 independent copies
//
// Revision 1.2  2008/05/09 22:14:35  balewski
// new there are 2 filters
//
// Revision 1.1  2008/04/21 15:47:09  balewski
// star
//