/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
/* $Id$ */

#ifndef AliFlowLYZEventPlane_H
#define AliFlowLYZEventPlane_H

// AliFlowLYZEventPlane:
// Class to calculate the event plane and event weight from the LYZ method
// author: N. van der Kolk (kolk@nikhef.nl)

#include "TString.h"
#include "AliFlowVector.h"

class AliFlowEventSimple;
class TProfile;
class TFile;

class AliFlowLYZEventPlane {
 public:
  AliFlowLYZEventPlane();
  virtual ~AliFlowLYZEventPlane();

  void Init();
  void CalculateRPandW(AliFlowVector aQ);

  Double_t GetWR() const  {return this->fWR; }
  Double_t GetPsi() const {return this->fPsi; }
  
  // input file
  void	   SetSecondRunFileName(TString name) 	
    { this->fSecondRunFileName = name ; }     // Sets input file name
  TString  GetSecondRunFileName() const		
    { return this->fSecondRunFileName ; }     // Gets output file name
  void     SetSecondRunFile(TFile* file)         
    { this->fSecondRunFile = file ; }         // Sets second run file


 private:
  
  AliFlowLYZEventPlane(const AliFlowLYZEventPlane& aAnalysis);             // copy constructor
  AliFlowLYZEventPlane& operator=(const AliFlowLYZEventPlane& aAnalysis);  // assignment operator
  
  TFile*   fSecondRunFile ;         //! pointer to file from second run
  TString  fSecondRunFileName;      //!

  Double_t fWR;            // event weight
  Double_t fPsi;           // reaction plane

  TProfile* fSecondReDtheta; // holds Re of Dtheta
  TProfile* fSecondImDtheta; // holds Im of Dtheta
  TProfile* fFirstr0theta;   // holds r0(theta)

  ClassDef(AliFlowLYZEventPlane, 0);          
};

#endif