]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG2/FEMTOSCOPY/AliFemto/AliFemtoV0.cxx
Making the directory structure of AliFemto flat. All files go into one common directory
[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemto / AliFemtoV0.cxx
similarity index 88%
rename from PWG2/FEMTOSCOPY/AliFemto/Infrastructure/AliFemtoV0.cxx
rename to PWG2/FEMTOSCOPY/AliFemto/AliFemtoV0.cxx
index b178dc5eb7692ece3eea80205dbf5b47f78daf7d..75b33d78e65e1743120c06a18378c6c200e2a08e 100644 (file)
@@ -1,3 +1,12 @@
+///////////////////////////////////////////////////////////////////////////
+//                                                                       //
+// AliFemtoV0: special type of particle desling with the specifics       //
+// of the V0 type of particle                                            //
+// It stores the information both about the V0 itself and about it's     //
+// daughters, so that the caut betwen the daughter characteristics is    //
+// possible.
+//                                                                       //
+///////////////////////////////////////////////////////////////////////////
 #include "AliFemtoV0.h"
 #include "phys_constants.h"
 
@@ -26,7 +35,10 @@ AliFemtoV0::AliFemtoV0():
   fPtNeg(0),  fPtotNeg(0),
   fKeyNeg(0),   fKeyPos(0),
   fHiddenInfo(0)  /***/
-{ /* no-op */}
+{ 
+  // Default empty constructor
+  /* no-op */
+}
 // -----------------------------------------------------------------------
 AliFemtoV0::AliFemtoV0(const AliFemtoV0& v) :
   fDecayLengthV0(0), fDecayVertexV0(0), fPrimaryVertex(0),
@@ -52,7 +64,8 @@ AliFemtoV0::AliFemtoV0(const AliFemtoV0& v) :
   fPtNeg(0),  fPtotNeg(0),
   fKeyNeg(0),   fKeyPos(0),
   fHiddenInfo(0)  /***/
-{ // copy constructor
+{ 
+  // copy constructor
   fDecayLengthV0 = v.fDecayLengthV0;
   fDecayVertexV0 = v.fDecayVertexV0;
   fDcaV0Daughters = v.fDcaV0Daughters;
@@ -91,11 +104,12 @@ AliFemtoV0::AliFemtoV0(const AliFemtoV0& v) :
 
   fHelixPos = v.fHelixPos;// Gael 12 Sept
   fHelixNeg = v.fHelixNeg;// Gael 12 Sept
-  fHiddenInfo = v.fHiddenInfo? v.fHiddenInfo->clone() : 0;// GR 11 DEC 02
+  fHiddenInfo = v.fHiddenInfo? v.fHiddenInfo->Clone() : 0;// GR 11 DEC 02
   UpdateV0();
 }
 AliFemtoV0& AliFemtoV0::operator=(const AliFemtoV0& aV0)
 {
+  // assignment operator
   if (this == &aV0)
     return *this;
   fDecayLengthV0 = aV0.fDecayLengthV0;
@@ -137,7 +151,7 @@ AliFemtoV0& AliFemtoV0::operator=(const AliFemtoV0& aV0)
   fHelixPos = aV0.fHelixPos;// Gael 12 Sept
   fHelixNeg = aV0.fHelixNeg;// Gael 12 Sept
   if (fHiddenInfo) delete fHiddenInfo;
-  fHiddenInfo = aV0.fHiddenInfo? aV0.fHiddenInfo->clone() : 0;// GR 11 DEC 02
+  fHiddenInfo = aV0.fHiddenInfo? aV0.fHiddenInfo->Clone() : 0;// GR 11 DEC 02
   UpdateV0();
   
   return *this;
@@ -146,7 +160,7 @@ AliFemtoV0& AliFemtoV0::operator=(const AliFemtoV0& aV0)
 // -----------------------------------------------------------------------
 void AliFemtoV0::UpdateV0(){
   //Calc. derived memebers of the v0 class
-  float MomNegAlongV0, MomPosAlongV0;
+  float tMomNegAlongV0, tMomPosAlongV0;
 
    fMomV0  = fMomPos + fMomNeg;
    fPtV0   = fMomV0.perp();
@@ -162,11 +176,11 @@ void AliFemtoV0::UpdateV0(){
    fEPosPion = ::sqrt(fPtotPos*fPtotPos+M_PION_PLUS*M_PION_PLUS);
    fENegPion = ::sqrt(fPtotNeg*fPtotNeg+M_PION_MINUS*M_PION_MINUS);
   
-   MomNegAlongV0 =  fMomNeg*fMomV0 / ::sqrt(::pow(fPtotV0,2));
-   MomPosAlongV0 =  fMomPos*fMomV0 / ::sqrt(::pow(fPtotV0,2));
+   tMomNegAlongV0 =  fMomNeg*fMomV0 / ::sqrt(::pow(fPtotV0,2));
+   tMomPosAlongV0 =  fMomPos*fMomV0 / ::sqrt(::pow(fPtotV0,2));
 
-   fAlphaV0 = (MomPosAlongV0-MomNegAlongV0)/(MomPosAlongV0+MomNegAlongV0);
-   fPtArmV0 =  ::sqrt(fPtotPos*fPtotPos - MomPosAlongV0*MomPosAlongV0);
+   fAlphaV0 = (tMomPosAlongV0-tMomNegAlongV0)/(tMomPosAlongV0+tMomNegAlongV0);
+   fPtArmV0 =  ::sqrt(fPtotPos*fPtotPos - tMomPosAlongV0*tMomPosAlongV0);
    fMassLambda = ::sqrt(::pow(fEPosProton+fENegPion,2)-::pow(fPtotV0,2));
    fMassAntiLambda = ::sqrt(::pow(fENegProton+fEPosPion,2)-::pow(fPtotV0,2));
    fMassK0Short = ::sqrt(::pow(fENegPion+fEPosPion,2)-::pow(fPtotV0,2));
@@ -268,5 +282,5 @@ const AliFmPhysicalHelixD& AliFemtoV0::HelixNeg() const {return fHelixNeg;}// Ga
 
 void AliFemtoV0::SetHiddenInfo(AliFemtoHiddenInfo* aHiddenInfo) {fHiddenInfo=aHiddenInfo;}
 bool AliFemtoV0::ValidHiddenInfo() const { if (fHiddenInfo) return true; else return false; }
-AliFemtoHiddenInfo* AliFemtoV0::getHiddenInfo() const {return fHiddenInfo;}
+AliFemtoHiddenInfo* AliFemtoV0::GetHiddenInfo() const {return fHiddenInfo;}