Possibility to fix some of the parameters. New method to get the number of free param...

[u/mrichter/AliRoot.git] / STEER / AliESDcascade.cxx

diff --git a/STEER/AliESDcascade.cxx b/STEER/AliESDcascade.cxx
index 1a01046b8c7c8796115e3afcfaf600a2f62a50be..5e96cb1e066245454664d3a49a971ff98a8339ae 100644 (file)
--- a/STEER/AliESDcascade.cxx
+++ b/STEER/AliESDcascade.cxx
@@ -35,10 +35,10 @@ ClassImp(AliESDcascade)
 
 AliESDcascade::AliESDcascade() : 
   AliESDv0(),
-  fPdgCode(kXiMinus),
-  fEffMass(TDatabasePDG::Instance()->GetParticle(kXiMinus)->Mass()),
-  fChi2Xi(1.e+33),
-  fDcaXiDaughters(999),
+  fEffMassXi(TDatabasePDG::Instance()->GetParticle(kXiMinus)->Mass()),
+  fChi2Xi(1024),
+  fDcaXiDaughters(1024),
+  fPdgCodeXi(kXiMinus),
   fBachIdx(-1)
 {
   //--------------------------------------------------------------------
@@ -49,17 +49,17 @@ AliESDcascade::AliESDcascade() :
     fBachMom[j]=0.;
   }
 
-  fPosCovXi[0]=1e10;
+  fPosCovXi[0]=1024;
   fPosCovXi[1]=fPosCovXi[2]=0.;
-  fPosCovXi[3]=1e10;
+  fPosCovXi[3]=1024;
   fPosCovXi[4]=0.;
-  fPosCovXi[5]=1e10;
+  fPosCovXi[5]=1024;
 
-  fBachMomCov[0]=1e10;
+  fBachMomCov[0]=1024;
   fBachMomCov[1]=fBachMomCov[2]=0.;
-  fBachMomCov[3]=1e10;
+  fBachMomCov[3]=1024;
   fBachMomCov[4]=0.;
-  fBachMomCov[5]=1e10;
+  fBachMomCov[5]=1024;
 }
 
 AliESDcascade::~AliESDcascade() {
@@ -68,10 +68,10 @@ AliESDcascade::~AliESDcascade() {
 AliESDcascade::AliESDcascade(const AliESDv0 &v,
                             const AliExternalTrackParam &t, Int_t i) : 
   AliESDv0(v),
-  fPdgCode(kXiMinus),
-  fEffMass(-1),
-  fChi2Xi(1.e+33),
-  fDcaXiDaughters(-1),
+  fEffMassXi(TDatabasePDG::Instance()->GetParticle(kXiMinus)->Mass()),
+  fChi2Xi(1024),
+  fDcaXiDaughters(1024),
+  fPdgCodeXi(kXiMinus),
   fBachIdx(i)
 {
   //---------------------------------------------------------------------------------------------
@@ -108,7 +108,7 @@ AliESDcascade::AliESDcascade(const AliESDv0 &v,
   Double_t e1=TMath::Sqrt(0.13957*0.13957 + px1*px1 + py1*py1 + pz1*pz1);
   Double_t e2=TMath::Sqrt(1.11568*1.11568 + px2*px2 + py2*py2 + pz2*pz2);
   
-  fEffMass=TMath::Sqrt((e1+e2)*(e1+e2)-
+  fEffMassXi=TMath::Sqrt((e1+e2)*(e1+e2)-
     (px1+px2)*(px1+px2)-(py1+py2)*(py1+py2)-(pz1+pz2)*(pz1+pz2));
 
 
@@ -117,28 +117,28 @@ AliESDcascade::AliESDcascade(const AliESDv0 &v,
   fBachMom[0]=px1; fBachMom[1]=py1; fBachMom[2]=pz1; 
 
   //PH Covariance matrices: to be calculated correctly in the future
-  fPosCovXi[0]=1e10;
+  fPosCovXi[0]=1024;
   fPosCovXi[1]=fPosCovXi[2]=0.;
-  fPosCovXi[3]=1e10;
+  fPosCovXi[3]=1024;
   fPosCovXi[4]=0.;
-  fPosCovXi[5]=1e10;
+  fPosCovXi[5]=1024;
 
-  fBachMomCov[0]=1e10;
+  fBachMomCov[0]=1024;
   fBachMomCov[1]=fBachMomCov[2]=0.;
-  fBachMomCov[3]=1e10;
+  fBachMomCov[3]=1024;
   fBachMomCov[4]=0.;
-  fBachMomCov[5]=1e10;
+  fBachMomCov[5]=1024;
 
-  fChi2Xi=7.;   
+  fChi2Xi=1024.;   
 
 }
 
 AliESDcascade::AliESDcascade(const AliESDcascade& cas) :
   AliESDv0(cas),
-  fPdgCode(cas.fPdgCode),
-  fEffMass(cas.fEffMass),
+  fEffMassXi(cas.fEffMassXi),
   fChi2Xi(cas.fChi2Xi),
   fDcaXiDaughters(cas.fDcaXiDaughters),
+  fPdgCodeXi(cas.fPdgCodeXi),
   fBachIdx(cas.fBachIdx)
 {
   //copy constructor
@@ -152,6 +152,44 @@ AliESDcascade::AliESDcascade(const AliESDcascade& cas) :
   }
 }
 
+AliESDcascade& AliESDcascade::operator=(const AliESDcascade& cas){
+
+  // -------
+  // Assigmnet oeprator
+  // -----
+
+  if(this==&cas) return *this;
+  AliESDv0::operator=(cas);
+
+  fEffMassXi = cas.fEffMassXi;
+  fChi2Xi    = cas.fChi2Xi;
+  fDcaXiDaughters = cas.fDcaXiDaughters;
+  fPdgCodeXi      = cas.fPdgCodeXi;
+  fBachIdx        = cas.fBachIdx;
+  for (int i=0; i<3; i++) {
+    fPosXi[i]     = cas.fPosXi[i];
+    fBachMom[i]   = cas.fBachMom[i];
+  }
+  for (int i=0; i<6; i++) {
+    fPosCovXi[i]   = cas.fPosCovXi[i];
+    fBachMomCov[i] = cas.fBachMomCov[i];
+  }
+  return *this;
+}
+
+void AliESDcascade::Copy(TObject &obj) const {
+  
+  // this overwrites the virtual TOBject::Copy()
+  // to allow run time copying without casting
+  // in AliESDEvent
+
+  if(this==&obj)return;
+  AliESDcascade *robj = dynamic_cast<AliESDcascade*>(&obj);
+  if(!robj)return; // not an AliESDcascade
+  *robj = *this;
+
+}
+
 Double_t AliESDcascade::ChangeMassHypothesis(Double_t &v0q, Int_t code) {
   //--------------------------------------------------------------------
   // This function changes the mass hypothesis for this cascade
@@ -161,7 +199,7 @@ Double_t AliESDcascade::ChangeMassHypothesis(Double_t &v0q, Int_t code) {
   Double_t nmass=0.13957, pmass=0.93827, ps0=0.101; 
   Double_t bmass=0.13957, mass =1.3213,  ps =0.139;
 
-  fPdgCode=code;
+  fPdgCodeXi=code;
 
   switch (code) {
   case 213: 
@@ -181,12 +219,13 @@ Double_t AliESDcascade::ChangeMassHypothesis(Double_t &v0q, Int_t code) {
        break;
   default:
        AliError("Invalide PDG code !  Assuming XiMinus's...");
-       fPdgCode=kXiMinus;
+       fPdgCodeXi=kXiMinus;
     break;
   }
 
   Double_t pxn=fNmom[0], pyn=fNmom[1], pzn=fNmom[2];
   Double_t pxp=fPmom[0], pyp=fPmom[1], pzp=fPmom[2];
+
   Double_t px0=pxn+pxp, py0=pyn+pyp, pz0=pzn+pzp;
   Double_t p0=TMath::Sqrt(px0*px0 + py0*py0 + pz0*pz0);
 
@@ -210,7 +249,7 @@ Double_t AliESDcascade::ChangeMassHypothesis(Double_t &v0q, Int_t code) {
   Double_t pxl=px0+pxb, pyl=py0+pyb, pzl=pz0+pzb;
   Double_t pl=TMath::Sqrt(pxl*pxl + pyl*pyl + pzl*pzl);
   
-  fEffMass=TMath::Sqrt((e0+eb)*(e0+eb) - pl*pl);
+  fEffMassXi=TMath::Sqrt((e0+eb)*(e0+eb) - pl*pl);
 
   Double_t beta=pl/(e0+eb);
   Double_t pl0=(px0*pxl + py0*pyl + pz0*pzl)/pl;
@@ -229,7 +268,7 @@ AliESDcascade::GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const {
   //--------------------------------------------------------------------
   // This function returns the cascade momentum (global)
   //--------------------------------------------------------------------
-  px=fNmom[0]+fPmom[0]+fBachMom[0]; 
+  px=fNmom[0]+fPmom[0]+fBachMom[0];
   py=fNmom[1]+fPmom[1]+fBachMom[1]; 
   pz=fNmom[2]+fPmom[2]+fBachMom[2]; 
 }
@@ -282,3 +321,8 @@ Double_t AliESDcascade::GetCascadeCosineOfPointingAngle(Double_t& refPointX, Dou
   
   return cosinePointingAngle;
 }
+
+void AliESDcascade::GetPosCovXi(Double_t cov[6]) const {
+
+  for (Int_t i=0; i<6; ++i) cov[i] = fPosCovXi[i];
+}