StEEmcPi0Analysis.h

#ifndef __StEEmcPi0Analysis_h__
#define __StEEmcPi0Analysis_h__

#include "StMaker.h"

#include "StMuDSTMaker/COMMON/StMuDstMaker.h"
#include "StMuDSTMaker/COMMON/StMuDst.h"
#include "StMuDSTMaker/COMMON/StMuEvent.h"
#include "StMuDSTMaker/COMMON/StMuTriggerIdCollection.h"
#include "StMuDSTMaker/COMMON/StMuPrimaryVertex.h" 
#include "StEvent/StTriggerId.h" 

#include "StSpinPool/StSpinDbMaker/StSpinDbMaker.h"

#include "StEEmcMixEvent.h" 

class StEEmcMixMaker;
class StEEmcPointMaker;
class StRFEmcTrigMaker;
class TH1F;
class TH2F;
class TFile; 
class SpinCuts; 
class EEmcGeomSimple;
class StEEmcA2EMaker; 
class TTree; 

#include <vector>
#include "StEEmcPair.h"
#include "SpinHistos.h"
#include "TString.h" 

class StEEmcPi0Analysis : public StMaker
{

 public:
  StEEmcPi0Analysis(const Char_t *name);
  ~StEEmcPi0Analysis(){ /* nada */ };

  Int_t Init();
  Int_t Make();
  Int_t InitRun(Int_t);
  void  Clear(Option_t *opts="");
  void mixer(const Char_t *name);
  void points(const Char_t *name);
  void mudst(const Char_t *name); 
  void analysis(const Char_t *name); 
  void spin(const Char_t *name);

  void trigger(Int_t trig){ mTriggerList.push_back(trig); }

  void minbias(Int_t trig){ mMinBias = trig; }

  SpinCuts *cuts(){ return mCuts; } 

  void filename(const Char_t *fname){ mFilename=fname; } 

  void triggerSim( const Char_t *name, Int_t t ) { mTrigSim=(StRFEmcTrigMaker*)GetMaker(name); mTrigSimThreshold=t; }

 private:
 protected:

  Bool_t twoBodyCut( StEEmcPair &p );

  StEEmcMixMaker   *mEEmixer;  
  StEEmcPointMaker *mEEpoints; 
  StMuDstMaker *mMuDst; 
  StEEmcA2EMaker *mEEanalysis; 
  StSpinDbMaker  *mSpinDb;

  std::vector<Int_t> mTriggerList;   
  Int_t mMinBias;                    
  enum SpinStates { kAny=0, kUpUp, kUpDn, kDnUp, kDnDn };

  SpinHistos *mHistograms[100];
  SpinHistos *mBackgrounds[100];
  SpinCuts   *mCuts; 
  TH1F       *hFillPatternI;    
  TH1F       *hSpin4;           
  TH1F       *hBx7;             
  TH1F       *hBx48;            
  TH2F       *hBx7diffBx48;     
  TH1F       *hBxStar;          
  TH1F       *hBxStarPi0;       
  TH2F       *hMassBx;          
  TH2F       *hZvertBx;         
  TH2F       *hZggBx;           
  TH2F       *hEtaBx;           
  TH1F       *hSpin4mb;         
  enum EventCutTypes { 
      kEvent=0,
      kMinBias, 
      kTrig,
      kSoftTrig, 
      kSoftTrig36, 
      kVertex,
      kCandidate 
  }; 

  enum PairCutTypes {
      kPair=1, 
      kTwoBody,
      kAdcTow,
      kEtTow,
      kKinematics 
  }; 
      
  TH1F *hEventCounter; 

  TH1F *hPairCounter; 

  Bool_t accept(StMuDst *mu); 
  Bool_t accept(StEEmcPair pair, Bool_t fill=true ); 

  Int_t getSpinState(StMuDst *mu, Int_t &bxs);

  TString mFilename;
  TFile *mFile; 

  StRFEmcTrigMaker *mTrigSim;
  Int_t mTrigSimThreshold;

  Int_t mRunNumber;
  Bool_t mSpinSort;

  EEmcGeomSimple *mEEgeom;

  StEEmcMixEvent *mRealEvent;
  StEEmcMixEvent *mMixedEvent; 

  TTree *mRealTree;
  TTree *mMixedTree; 

  ClassDef(StEEmcPi0Analysis,1);  

};

#endif