StFlowLeeYangZerosMaker.h

//
// $Id: StFlowLeeYangZerosMaker.h,v 1.5 2014/08/06 11:43:14 jeromel Exp $
//
// Authors: Markus Oldenberg and Art Poskanzer, LBNL
//
//
// Description:  Maker to analyze Flow by the LeeYangZeros method
//

#ifndef StFlowLeeYangZerosMaker_H
#define StFlowLeeYangZerosMaker_H
#include <Stiostream.h>
#include "StMaker.h"
#include "StFlowMaker/StFlowConstants.h"
#include "TVector2.h"
#include "TComplex.h"
#include "TString.h"
//#include "StTimer.hh"
class StFlowEvent;
class StFlowSelection;
class TH1F;
class TH1D;
class TProfile;
class TProfile2D;

class StFlowLeeYangZerosMaker : public StMaker {

public:

  StFlowLeeYangZerosMaker(const Char_t* name="FlowLeeYangZeros");

  StFlowLeeYangZerosMaker(const Char_t* name,
                          const StFlowSelection& pFlowSelect);
  StFlowLeeYangZerosMaker(const StFlowLeeYangZerosMaker &from) {};

  virtual  ~StFlowLeeYangZerosMaker();

  Int_t    Init();
  Int_t    Make();
  Int_t    Finish();
  void     SetHistoRanges(Bool_t ftpc_included = kFALSE);
  void     SetPtRange_for_vEta(Float_t lo, Float_t hi);
  void     SetEtaRange_for_vPt(Float_t lo, Float_t hi);
  static   void SetV1Mixed(Bool_t);
  virtual  const char *GetCVS() const {static const char cvs[]=
    "Tag $Name:  $  StFlowLeeYangZerosMaker.h, " __DATE__ " " __TIME__ ;
    return cvs;}

private:

  Bool_t   FillFromFlowEvent();
  void     FillParticleHistograms();
  Bool_t   mZeroPass;
  Bool_t   mFirstPass;
  static   Bool_t   mV1Mixed;        // flag for v1 mixed harmonic
#ifndef __CINT__
  TVector2 mQ[Flow::nSels][Flow::nHars];                     
  TVector2 mqTPC[Flow::nSels][Flow::nHars];                  
  TVector2 mqTPCE[Flow::nSels][Flow::nHars];                 
  TVector2 mqTPCW[Flow::nSels][Flow::nHars];                 
  Double_t mQ2[Flow::nSels][Flow::nHars];                    
  Float_t  mQtheta[Flow::nSels][Flow::nHars][Flow::nTheta];  
  Float_t  mr0theta[Flow::nSels][Flow::nHars][Flow::nTheta]; 
  TComplex mGr0theta[Flow::nSels][Flow::nHars][Flow::nTheta];
  TComplex mV1neum[Flow::nTheta];                            
  TComplex mGV1r0theta[Flow::nTheta1][Flow::nTheta];         
#endif /*__CINT__*/
  Int_t            mNEvents;    
  Int_t            mMult;       
  Int_t            mPtBinsPart; 
  TString          xLabel;      
  StFlowEvent*     pFlowEvent;  
  StFlowSelection* pFlowSelect; 

  // for single histograms
  TH1D*     mHistYieldPartPt;          
  TH1D*     mHistYieldPartEta;         
  TH1F*     mHistMult;                 

  // for each harmonic, each selection, and each theta
  struct histThetas {
    TH1D*       mHistGtheta;
    TProfile*   mHistProReGtheta;
    TProfile*   mHistProImGtheta;    
    TH1D*       mHistReGtheta;
    TH1D*       mHistImGtheta;    
    TProfile2D* mHistReNumer2D;
    TProfile*   mHistReNumerEta;
    TProfile2D* mHistImNumer2D;
    TProfile*   mHistImNumerEta;
  };

  // for each harmonic and each selection
  struct histFullHars {
    TProfile*   mHistPro_vEta;
    TProfile*   mHistPro_vPt;
    TProfile*   mHistPro_Vtheta;
    TProfile*   mHistPro_r0theta;
    TH1D*       mHist_vEta;
    TH1D*       mHist_vPt;
    TH1D*       mHist_Vtheta;
    TH1D*       mHist_r0theta;
    TProfile*   mHistReDenom;
    TProfile*   mHistImDenom;
    TProfile*   mHistCentX;
    TProfile*   mHistCentY;
    TProfile*   mHistQCent;
    struct histThetas histTheta[Flow::nTheta]; 
  };

  // for each selection
  struct histFulls;     
  friend struct histFulls;
  struct histFulls {
    TProfile*   mHistPro_V;
    TProfile*   mHistPro_vr0;
    TH1D*       mHist_V;
    TH1D*       mHist_vr0;
    TH1F*       mHist_v;
    struct histFullHars histFullHar[Flow::nHars];
  };
  struct histFulls histFull[Flow::nSels]; 

  TString      MakerName;

  Float_t mEtaMin;
  Float_t mEtaMax;
    Int_t mNEtaBins;

  Float_t mPtRange_for_vEta[2];
  Float_t mEtaRange_for_vPt[2];

/*   StTimer timeInit; */
/*   StTimer timeEvent; */
/*   StTimer timePart; */
/*   StTimer timeFinish; */

  ClassDef(StFlowLeeYangZerosMaker,0)              // macro for rootcint
};

inline void StFlowLeeYangZerosMaker::SetV1Mixed(Bool_t V1Mixed) { mV1Mixed = V1Mixed; }

#endif