]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGDQ/dielectron/AliDielectronPair.h
-iarsene: FMD includes commented out
[u/mrichter/AliRoot.git] / PWGDQ / dielectron / AliDielectronPair.h
index 9c3f916bb2ed896f2292eb29a75a8dbc973e8855..541716f014bb12cca0e114ed215620b47b314e43 100644 (file)
@@ -26,6 +26,7 @@
 
 #include <AliKFParticle.h>
 #include <AliVParticle.h>
+#include <AliVEvent.h>
 
 class AliVVertex;
 class AliVTrack;
@@ -52,6 +53,10 @@ public:
   void SetTracks(AliVTrack * const particle1, Int_t pid1,
                  AliVTrack * const particle2, Int_t pid2);
 
+  void SetGammaTracks(AliVTrack * const particle1, Int_t pid1,
+                     AliVTrack * const particle2, Int_t pid2);
+
+
   void SetTracks(const AliKFParticle * const particle1,
                  const AliKFParticle * const particle2,
                  AliVTrack * const refParticle1,
@@ -97,6 +102,7 @@ public:
 
   UChar_t GetType() const { return fType; }
   void SetType(Char_t type) { fType=type; }
+  static void SetBeamEnergy(AliVEvent *ev, Double_t beamEbyHand=-1.);
 
   // MC information
   void SetLabel(Int_t label) {fLabel=label;}
@@ -120,13 +126,25 @@ public:
   // calculate cos(theta*) and phi* in HE and CS pictures
   void GetThetaPhiCM(Double_t &thetaHE, Double_t &phiHE, Double_t &thetaCS, Double_t &phiCS) const;
   
+
   Double_t ThetaPhiCM(Bool_t isHE, Bool_t isTheta) const;
   static Double_t ThetaPhiCM(const AliVParticle* d1, const AliVParticle* d2, 
                                               const Bool_t isHE, const Bool_t isTheta);
 
   Double_t PsiPair(Double_t MagField)const; //Angle cut w.r.t. to magnetic field
   Double_t PhivPair(Double_t MagField)const; //Angle of ee plane w.r.t. to magnetic field
+  //  Double_t PairPlanev0rpH2Angle(Double_t kv0CrpH2)const;
+  // Double_t PairPlaneMagAngle(Double_t kv0CrpH2)const;
+  Double_t PairPlaneAngle(Double_t kv0CrpH2)const;
+
   Double_t GetCosPointingAngle(const AliVVertex *primVtx) const;
+  Double_t GetArmAlpha() const;
+  Double_t GetArmPt()    const;
+
+  // calculate rotation of p1 p2
+  void GetRotPair(Double_t &RotPairx, Double_t &RotPairy, Double_t &RotPairz) const;
+  
+
   // internal KF particle
   const AliKFParticle& GetKFParticle()       const { return fPair; }
   const AliKFParticle& GetKFFirstDaughter()  const { return fD1;   }
@@ -148,6 +166,7 @@ private:
   Char_t   fType;         // type of the pair e.g. like sign SE, unlike sign SE, ... see AliDielectron
   Int_t    fLabel;        // MC label
   Int_t    fPdgCode;      // pdg code in case it is a MC particle
+  static Double_t fBeamEnergy; //!beam energy
   
   AliKFParticle fPair;   // KF particle internally used for pair calculation
   AliKFParticle fD1;     // KF particle first daughter
@@ -158,7 +177,8 @@ private:
 
   Bool_t fKFUsage;       // Use KF for vertexing
   
-  ClassDef(AliDielectronPair,4)
+  ClassDef(AliDielectronPair,5)
 };
 
 #endif
+