Removing non-existing AliRICHReconOld class from the list

[u/mrichter/AliRoot.git] / HBTAN / AliHBTCrab.cxx

diff --git a/HBTAN/AliHBTCrab.cxx b/HBTAN/AliHBTCrab.cxx
index b91bfef0aef1e69a19a9acdb568edbf2a1845efd..d0623c4d8761a3b77c0288ea0fefa83335d1ac49 100644 (file)
--- a/HBTAN/AliHBTCrab.cxx
+++ b/HBTAN/AliHBTCrab.cxx
@@ -10,7 +10,8 @@
 /////////////////////////////////////////////////////////////////////////
  
 #include "AliHBTPair.h"
-
+#include "AliVAODParticle.h"
+#include "TDatabasePDG.h"
 #include <TMath.h>
 
 #include "volya_complex.h"
@@ -60,6 +61,10 @@ fReducedMom(kTRUE),
 fMaxMomentum(100.0)
 {
   //ctor
+  if(fgCrab)
+   {
+     Fatal("AliHBTCrab","Do not use constructor directly. Use Instance() instead.");
+   }
 }
 //===================================================================
 
@@ -72,7 +77,7 @@ fMaxMomentum(100.0)
   //ctor
 }
 //===================================================================
-const AliHBTCrab & AliHBTCrab::operator=(const AliHBTCrab& /*source*/)
+AliHBTCrab & AliHBTCrab::operator=(const AliHBTCrab& /*source*/)
 {
 //cpy constructor
  return *AliHBTCrab::Instance();
@@ -108,7 +113,7 @@ Bool_t AliHBTCrab::SetConfig(const AliHBTPair* pair)
 }
 //===================================================================
 
-Double_t AliHBTCrab::GetWeight(const AliHBTPair* partpair)
+Double_t AliHBTCrab::GetWeight(AliHBTPair* partpair)
 {
 //returns the weight
   Double_t qred, r, qdotr, mom;
@@ -141,15 +146,15 @@ void AliHBTCrab::GetComQuantities(const AliHBTPair* pair,
  static const Double_t kCmToFm = 1.e13;
 // static const Double_t cmtoOneOverGeV = kCmToFm*fgkWcons;
  
- AliHBTParticle *part1 = pair->Particle1();
- AliHBTParticle *part2 = pair->Particle2();
+ AliVAODParticle *part1 = pair->Particle1();
+ AliVAODParticle *part2 = pair->Particle2();
 
- p1[0] = part1->Energy()*1000.0;
+ p1[0] = part1->E()*1000.0;
  p1[1] = part1->Px()*1000.0;
  p1[2] = part1->Py()*1000.0;
  p1[3] = part1->Pz()*1000.0;
  
- p2[0] = part2->Energy()*1000.0;
+ p2[0] = part2->E()*1000.0;
  p2[1] = part2->Px()*1000.0;
  p2[2] = part2->Py()*1000.0;
  p2[3] = part2->Pz()*1000.0;
@@ -223,7 +228,7 @@ void AliHBTCrab::GetComQuantities(const AliHBTPair* pair,
   else //identical
   {
   //  const double  kdotp=fMass2*fMass2-fMass1*fMass1;
-   const double  kdotp = part2->GetMass()*part2->GetMass()- part1->GetMass()*part1->GetMass();
+   const double  kdotp = part2->Mass()*part2->Mass()- part1->Mass()*part1->Mass();
    *test=1;
    *mom=-(p2[0]-p1[0])*(p2[0]-p1[0]);
    ptot2=(p1[0]+p2[0])*(p1[0]+p2[0]);