Retrofit from Release developement
authorskowron <skowron@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Feb 2005 17:02:53 +0000 (17:02 +0000)
committerskowron <skowron@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Feb 2005 17:02:53 +0000 (17:02 +0000)
ANALYSIS/ANALYSISLinkDef.h
ANALYSIS/AliAODPair.cxx
ANALYSIS/AliAODPair.h
ANALYSIS/AliAODPairBaseCut.h
ANALYSIS/AliAODPairCut.cxx
ANALYSIS/AliAODPairCut.h

index 32844bf..accfef2 100644 (file)
@@ -70,6 +70,9 @@
 #pragma link C++ class AliAODQInvCut+;
 #pragma link C++ class AliAODKtCut+;
 #pragma link C++ class AliAODKStarCut+;
+#pragma link C++ class AliAODKStarOutCut+;
+#pragma link C++ class AliAODKStarSideCut+;
+#pragma link C++ class AliAODKStarLongCut+;
 #pragma link C++ class AliAODQSideLCMSCut+;
 #pragma link C++ class AliAODQOutLCMSCut+;
 #pragma link C++ class AliAODQLongLCMSCut+;
index e508f4c..57fac0f 100644 (file)
@@ -467,6 +467,16 @@ Double_t AliAODPair::GetGammaToLCMS()
 }
 /************************************************************************/
 
+Double_t AliAODPair::GetGammaToTransverse()
+{
+  //calculates gamma factor of the boost to LCMS
+  Double_t beta = 2.0*GetKt() / GetMt();
+  Double_t gamma = 1.0/TMath::Sqrt(1.0 - beta*beta);
+
+  return gamma;
+}
+/************************************************************************/
+
 Double_t AliAODPair::GetMt()
 {
   //Calculates transverse mass of the pair
@@ -624,3 +634,20 @@ void   AliAODPair::DeleteSecond()
   delete fPart2;
   fPart2 = 0x0;
 }
+
+void   AliAODPair::Print(Option_t* /*option*/) 
+{
+  if (fPart1) fPart1->Print();
+  if (fPart2) fPart2->Print();
+  
+  Info("Print","GetKStar() %f",GetKStar());
+  Info("Print","GetKt() %f",GetKt() );
+  Info("Print","QInv %f", GetQInv() );
+  Info("Print","GetQOutLCMS() %f",GetQOutLCMS() );
+  Info("Print","GetQSideLCMS %f",GetQSideLCMS() );
+  Info("Print","GetQLongLCMS() %f", GetQLongLCMS());
+  Info("Print","GetDeltaTheta() %f", GetDeltaTheta());
+  Info("Print","GetDeltaPhi() %f", GetDeltaPhi());
+  
+  
+}
index a928c62..69cb2d8 100644 (file)
@@ -67,6 +67,7 @@ class AliAODPair: public TObject
    virtual Double_t GetDeltaPhi();
    
    virtual Double_t GetGammaToLCMS();
+   virtual Double_t GetGammaToTransverse();
    virtual Double_t GetPIDProb() const {return fPart1->GetPidProb()*fPart2->GetPidProb();}
    
    virtual Double_t GetRStar() ;
@@ -74,6 +75,9 @@ class AliAODPair: public TObject
    
    void   MirrorSecond();
    void   DeleteSecond();
+   
+   void   Print(Option_t* /*option*/ = "" ) ;
+   
  protected:
    AliVAODParticle* fPart1;  //pointer to first particle
    AliVAODParticle* fPart2;  //pointer to second particle
index 4f7917d..6680165 100644 (file)
@@ -17,6 +17,9 @@ class AliAODPairBaseCut: public TObject
       kHbtPairCutPropQInv, //Q invariant
       kHbtPairCutPropKt,
       kHbtPairCutPropKStar,
+      kHbtPairCutPropKStarOut,
+      kHbtPairCutPropKStarSide,
+      kHbtPairCutPropKStarLong,
       kHbtPairCutPropQSideLCMS,
       kHbtPairCutPropQOutLCMS,
       kHbtPairCutPropQLongLCMS,
@@ -111,6 +114,42 @@ class AliAODKStarCut: public AliAODPairBaseCut
 };
 /******************************************************************/
 
+class AliAODKStarOutCut: public AliAODPairBaseCut
+{
+ public:
+  AliAODKStarOutCut(Double_t min = 0.0, Double_t max = 0.0):AliAODPairBaseCut(min,max,kHbtPairCutPropKStarOut){}
+  virtual ~AliAODKStarOutCut(){}
+ protected:
+  virtual Double_t  GetValue(AliAODPair* pair) const {return pair->GetKStarOut();}
+
+  ClassDef(AliAODKStarOutCut,1)
+};
+/******************************************************************/
+class AliAODKStarSideCut: public AliAODPairBaseCut
+{
+ public:
+  AliAODKStarSideCut(Double_t min = 0.0, Double_t max = 0.0):AliAODPairBaseCut(min,max,kHbtPairCutPropKStarSide){}
+  virtual ~AliAODKStarSideCut(){}
+ protected:
+  virtual Double_t  GetValue(AliAODPair* pair) const {return pair->GetKStarSide();}
+
+  ClassDef(AliAODKStarSideCut,1)
+};
+
+/******************************************************************/
+
+class AliAODKStarLongCut: public AliAODPairBaseCut
+{
+ public:
+  AliAODKStarLongCut(Double_t min = 0.0, Double_t max = 0.0):AliAODPairBaseCut(min,max,kHbtPairCutPropKStarLong){}
+  virtual ~AliAODKStarLongCut(){}
+ protected:
+  virtual Double_t  GetValue(AliAODPair* pair) const {return pair->GetKStarLong();}
+
+  ClassDef(AliAODKStarLongCut,1)
+};
+/******************************************************************/
+
 class AliAODQSideLCMSCut: public AliAODPairBaseCut
 {
  public:
index b51b787..033100a 100644 (file)
@@ -262,6 +262,33 @@ void AliAODPairCut::SetKStarRange(Double_t min, Double_t max)
 }
 /**********************************************************/
 
+void AliAODPairCut::SetKStarOutRange(Double_t min, Double_t max)
+{
+  // set range of accepted KStar (invariant pair momentum difference (fourvector))
+  AliAODKStarOutCut* cut= (AliAODKStarOutCut*)FindCut(AliAODPairBaseCut::kHbtPairCutPropKStarOut);
+  if(cut) cut->SetRange(min,max);
+  else fCuts[fNCuts++] = new AliAODKStarOutCut(min,max);
+}
+/**********************************************************/
+
+void AliAODPairCut::SetKStarSideRange(Double_t min, Double_t max)
+{
+  // set range of accepted KStar (invariant pair momentum difference (fourvector))
+  AliAODKStarSideCut* cut= (AliAODKStarSideCut*)FindCut(AliAODPairBaseCut::kHbtPairCutPropKStarSide);
+  if(cut) cut->SetRange(min,max);
+  else fCuts[fNCuts++] = new AliAODKStarSideCut(min,max);
+}
+/**********************************************************/
+
+void AliAODPairCut::SetKStarLongRange(Double_t min, Double_t max)
+{
+  // set range of accepted KStar (invariant pair momentum difference (fourvector))
+  AliAODKStarLongCut* cut= (AliAODKStarLongCut*)FindCut(AliAODPairBaseCut::kHbtPairCutPropKStarLong);
+  if(cut) cut->SetRange(min,max);
+  else fCuts[fNCuts++] = new AliAODKStarLongCut(min,max);
+}
+/**********************************************************/
+
 void AliAODPairCut::SetAvSeparationRange(Double_t min, Double_t max)
 {
   //sets avarage separation cut ->Anti-Merging cut
index cd43db4..f78a520 100644 (file)
@@ -42,6 +42,9 @@ class AliAODPairCut: public TNamed
   void SetQInvRange(Double_t min, Double_t max);
   void SetKtRange(Double_t min, Double_t max);
   void SetKStarRange(Double_t min, Double_t max);
+  void SetKStarOutRange(Double_t min, Double_t max);
+  void SetKStarSideRange(Double_t min, Double_t max);
+  void SetKStarLongRange(Double_t min, Double_t max);
   void SetQOutLCMSRange(Double_t min, Double_t max);
   void SetQSideLCMSRange(Double_t min, Double_t max);
   void SetQLongLCMSRange(Double_t min, Double_t max);