**************************************************************************/
#define AliFlowAnalysisWithLYZEventPlane_cxx
+
+// AliFlowLYZEventPlane:
+//
+// Class to calculate the event plane and event weight from the LYZ method
+// It needs input from the standard LYZ first and second run
+// author: N. van der Kolk (kolk@nikhef.nl)
#include "Riostream.h"
#include "TProfile.h"
#include "TFile.h"
#include "TComplex.h"
+#include "AliFlowVector.h"
#include "AliFlowLYZConstants.h"
#include "AliFlowEventSimple.h"
#include "AliFlowLYZEventPlane.h"
class AliFlowTrackSimple;
-// AliFlowLYZEventPlane:
-//
-// Class to calculate the event plane and event weight from the LYZ method
-//
-// author: N. van der Kolk (kolk@nikhef.nl)
-
ClassImp(AliFlowLYZEventPlane)
//-----------------------------------------------------------------------
fFirstr0theta(0)
{
// Constructor.
- fQ.Set(0.,0.); // flow vector
+
}
//-----------------------------------------------------------------------
fFirstr0theta = (TProfile*)fFirstRunFile->Get("First_FlowPro_r0theta_LYZ");
}
-
-
}
-
//-----------------------------------------------------------------------
-void AliFlowLYZEventPlane::CalculateRPandW(TVector2 fQ)
+void AliFlowLYZEventPlane::CalculateRPandW(AliFlowVector fQ)
{
//declare variables
Int_t fNtheta = AliFlowLYZConstants::kTheta;
fPsi = 0.5*TMath::ATan2(sinTerm,cosTerm); //takes care of the signs correctly!
if (fPsi < 0.) { fPsi += TMath::Pi(); } //to shift distribution from (-pi/2 to pi/2) to (0 to pi)
+}
-}
-//-----------------------------------------------------------------------
-TVector2 AliFlowLYZEventPlane::GetQ(AliFlowEventSimple* fEvent)
-{
- //get the Q vector
- TVector2 fQ = fEvent->GetQ();
-
- return fQ;
-
-}
#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 "TVector2.h"
+#include "AliFlowVector.h"
class AliFlowEventSimple;
class TProfile;
class TFile;
-// AliFlowLYZEventPlane:
-// Class to calculate the event plane and event weight from the LYZ method
-// author: N. van der Kolk (kolk@nikhef.nl)
-
class AliFlowLYZEventPlane {
public:
AliFlowLYZEventPlane();
virtual ~AliFlowLYZEventPlane();
void Init();
- void CalculateRPandW(TVector2 fQ);
+ void CalculateRPandW(AliFlowVector fQ);
Double_t GetWR() const {return this->fWR; }
Double_t GetPsi() const {return this->fPsi; }
private:
- TVector2 GetQ(AliFlowEventSimple* fEvent);
-
TFile* fFirstRunFile ; //! pointer to file from first run
TFile* fSecondRunFile ; //! pointer to file from second run
TString fFirstRunFileName; //!
TString fSecondRunFileName; //!
-
- TVector2 fQ; // flow vector
Double_t fWR; // event weight
Double_t fPsi; // reaction plane
};
#endif
+
+