DeltaE added and appropriate Pair Base Cut

diff --git a/ANALYSIS/AliAODPair.h b/ANALYSIS/AliAODPair.h
index f8de3c9cbed7ae00def663c9d19eee50ff96d8bd..8952966daade22873f78b379ae3582e5437c33b9 100644 (file)
--- a/ANALYSIS/AliAODPair.h
+++ b/ANALYSIS/AliAODPair.h
@@ -51,7 +51,9 @@ class AliAODPair: public TObject
    
    virtual Double_t GetAvarageDistance();//returns avarage distnace between two tracks
    
    
    virtual Double_t GetAvarageDistance();//returns avarage distnace between two tracks
    

-   virtual Double_t GetDeltaP(); //return difference of momenta
+   virtual Double_t GetDeltaE(); //return difference of Energies
+   virtual Double_t GetDeltaP(); //return difference of momenta (scalar difference)
+   virtual Double_t GetDeltaPvector(); //return legth of difference vector of momenta

    virtual Double_t GetDeltaPt();
    virtual Double_t GetDeltaPx();
    virtual Double_t GetDeltaPy();
    virtual Double_t GetDeltaPt();
    virtual Double_t GetDeltaPx();
    virtual Double_t GetDeltaPy();


@@ -244,10 +246,27 @@ void AliAODPair::CalculateDiffs()
  }
 
 /****************************************************************/
  }
 
 /****************************************************************/

+
+inline 
+Double_t AliAODPair::GetDeltaE() 
+{
+ //returns difference of energies
+  return fPart1->E() - fPart2->E();
+}
+/****************************************************************/
+
+inline 
+Double_t AliAODPair::GetDeltaP() 
+{
+ //returns difference of momenta (scalars)
+ return fPart1->P() - fPart2->P();
+}
+/****************************************************************/
+

 inline 
 inline 

-Double_t AliAODPair::GetDeltaP() //return difference of momenta
+Double_t AliAODPair::GetDeltaPvector() //return difference of momenta

 {
 {

- //returns difference of momenta (length of vector)
+ //returns legth of the momenta difference vector 

  CalculateDiffs();
  return TMath::Sqrt(fPxDiff*fPxDiff + fPyDiff*fPyDiff + fPzDiff*fPzDiff);
 }
  CalculateDiffs();
  return TMath::Sqrt(fPxDiff*fPxDiff + fPyDiff*fPyDiff + fPzDiff*fPzDiff);
 }

diff --git a/ANALYSIS/AliAODPairBaseCut.cxx b/ANALYSIS/AliAODPairBaseCut.cxx
index 1a9cc037a769588e693341f95e906d539d571386..8642d016afa2a1fb701ebc4c8fca0a7bc0628092 100644 (file)
--- a/ANALYSIS/AliAODPairBaseCut.cxx
+++ b/ANALYSIS/AliAODPairBaseCut.cxx
@@ -10,6 +10,11 @@ ClassImp(AliAODKtCut)
 ClassImp(AliAODQSideLCMSCut)
 ClassImp(AliAODQOutLCMSCut)
 ClassImp(AliAODQLongLCMSCut)
 ClassImp(AliAODQSideLCMSCut)
 ClassImp(AliAODQOutLCMSCut)
 ClassImp(AliAODQLongLCMSCut)

+ClassImp(AliAODDeltaECut)
+ClassImp(AliAODDeltaPCut)
+ClassImp(AliAODDeltaPvectorCut)
+ClassImp(AliAODDeltaPhiCut)
+ClassImp(AliAODDeltaThetaCut)

 
 /******************************************************************/
 ClassImp(AliAODAvSeparationCut)
 
 /******************************************************************/
 ClassImp(AliAODAvSeparationCut)

diff --git a/ANALYSIS/AliAODPairBaseCut.h b/ANALYSIS/AliAODPairBaseCut.h
index 2b1f40404ac7504fae21ea3194e8b8d266fbf854..4f7917d156f9a7ec214458370687c5136b5f9a6c 100644 (file)
--- a/ANALYSIS/AliAODPairBaseCut.h
+++ b/ANALYSIS/AliAODPairBaseCut.h
@@ -22,7 +22,9 @@ class AliAODPairBaseCut: public TObject
       kHbtPairCutPropQLongLCMS,
       kHbtPairCutPropDeltaPhi,
       kHbtPairCutPropDeltaTheta,
       kHbtPairCutPropQLongLCMS,
       kHbtPairCutPropDeltaPhi,
       kHbtPairCutPropDeltaTheta,

-      kHbtPairCutPropDeltaP,
+      kHbtPairCutPropDeltaE,
+      kHbtPairCutPropDeltaP,//scalar difference
+      kHbtPairCutPropDeltaPvector,//legth of the momenta difference vector

       kHbtPairCutPropDeltaPt,
       kHbtPairCutPropAvSepar,
       kHbtPairCutPropSepar,
       kHbtPairCutPropDeltaPt,
       kHbtPairCutPropAvSepar,
       kHbtPairCutPropSepar,


@@ -152,6 +154,48 @@ class AliAODQLongLCMSCut: public AliAODPairBaseCut
 };
 /******************************************************************/
 
 };
 /******************************************************************/
 

+class AliAODDeltaECut: public AliAODPairBaseCut
+{
+ public:
+  AliAODDeltaECut(Double_t min = 0.0, Double_t max = 0.0):
+    AliAODPairBaseCut(min,max,kHbtPairCutPropDeltaE){}
+  virtual ~AliAODDeltaECut(){}
+ protected:
+  virtual Double_t  GetValue(AliAODPair* pair) const 
+    {return pair->GetDeltaE();}
+
+  ClassDef(AliAODDeltaECut,1)
+};
+/******************************************************************/
+
+class AliAODDeltaPCut: public AliAODPairBaseCut
+{
+ public:
+  AliAODDeltaPCut(Double_t min = 0.0, Double_t max = 0.0):
+    AliAODPairBaseCut(min,max,kHbtPairCutPropDeltaP){}
+  virtual ~AliAODDeltaPCut(){}
+ protected:
+  virtual Double_t  GetValue(AliAODPair* pair) const 
+    {return pair->GetDeltaP();}
+
+  ClassDef(AliAODDeltaPCut,1)
+};
+/******************************************************************/
+
+class AliAODDeltaPvectorCut: public AliAODPairBaseCut
+{
+ public:
+  AliAODDeltaPvectorCut(Double_t min = 0.0, Double_t max = 0.0):
+    AliAODPairBaseCut(min,max,kHbtPairCutPropDeltaPvector){}
+  virtual ~AliAODDeltaPvectorCut(){}
+ protected:
+  virtual Double_t  GetValue(AliAODPair* pair) const 
+    {return pair->GetDeltaPvector();}
+
+  ClassDef(AliAODDeltaPvectorCut,1)
+};
+/******************************************************************/
+

 class AliAODDeltaPhiCut: public AliAODPairBaseCut
 {
  public:
 class AliAODDeltaPhiCut: public AliAODPairBaseCut
 {
  public: