StEEmcMixQAMaker.h

#ifndef __StEEmcMixQAMaker_h__
#define __StEEmcMixQAMaker_h__

#include "StMaker.h"

class StEEmcMixMaker;
class StEEmcPointMaker;
class TH1F;
class TH2F;

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

class StEEmcMixQAMaker : public StMaker
{

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

  Int_t Init();
  Int_t Make();
  void  Clear(Option_t *opts=""){ /* nada */ };
  void mixer(const Char_t *name, Float_t min=0., Float_t max=999.);
  void points(const Char_t *name);

  void background(){ mBackground=true; }
 
  
  Int_t maxPerSector;  
  Int_t maxPerEvent;   
  Int_t maxPerCluster; 
  Float_t minTowerFrac; 
  Float_t zVertexMin;
  Float_t zVertexMax;
  Float_t minZgg; 
  Float_t maxZgg; 
  void addBin(Float_t b){ assert(b>mBins.back()); mBins.push_back(b); }

 private:
 protected:

  Float_t mMin; 
  Float_t mMax; 
  Bool_t mBackground; 
  TH1F *hNcandidates; 
  TH1F *hNcandidatesR; 
  std::vector<TH2F *> hYXpair;  
  std::vector<TH2F *> hYXhigh;  
  std::vector<TH2F *> hYXlow;   
  std::vector<TH2F *> hE1E2;    
  
  std::vector<Float_t>               mBins;     

  //
  // The following are 2D vectors of histograms.  We 
  // bin the data in pT and sector.  Sector 13 integrates
  // over all 12 eemc sectors.  One final pT bin integrates
  // over all pT

  std::vector< std::vector<TH1F *> > hMassR;    
  std::vector< std::vector<TH1F *> > hZggR;     
  std::vector< std::vector<TH1F *> > hPhiggR;   
  std::vector< std::vector<TH1F *> > hEnergyR;  
  std::vector< std::vector<TH1F *> > hZvertexR; 

  TH1F *hMassRall;
  TH1F *hZvertexRall;

  Int_t ptbin( StEEmcPair p );

  Bool_t twoBodyCut( StEEmcPair p );

  StEEmcMixMaker   *mEEmixer;  
  StEEmcPointMaker *mEEpoints; 

  ClassDef(StEEmcMixQAMaker,1);  

};

#endif