Fixed warnings and coverity + added option to control filling of THnSparse

[u/mrichter/AliRoot.git] / PWGLF / RESONANCES / AliRsnMother.cxx

diff --git a/PWGLF/RESONANCES/AliRsnMother.cxx b/PWGLF/RESONANCES/AliRsnMother.cxx
index 7be7ab05acaa1378c8413d8abefd1945d962dce9..8f365f368140c47fec46eac11b2652b6471e3260 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnMother.cxx
+++ b/PWGLF/RESONANCES/AliRsnMother.cxx
@@ -196,19 +196,19 @@ Double_t AliRsnMother::CosThetaStar(Bool_t first, Bool_t useMC)
 
    TLorentzVector &mother    = Sum(useMC);
    TLorentzVector &daughter0 = (first ? fDaughter[0]->P(useMC) : fDaughter[1]->P(useMC));
-   TLorentzVector &daughter1 = (first ? fDaughter[1]->P(useMC) : fDaughter[0]->P(useMC));
+//    TLorentzVector &daughter1 = (first ? fDaughter[1]->P(useMC) : fDaughter[0]->P(useMC));
    TVector3 momentumM(mother.Vect());
    TVector3 normal(mother.Y() / momentumM.Mag(), -mother.X() / momentumM.Mag(), 0.0);
 
-   // Computes first the invariant mass of the mother
-   Double_t mass0      = fDaughter[0]->P(useMC).M();
-   Double_t mass1      = fDaughter[1]->P(useMC).M();
-   Double_t p0         = daughter0.Vect().Mag();
-   Double_t p1         = daughter1.Vect().Mag();
-   Double_t E0         = TMath::Sqrt(mass0 * mass0 + p0 * p0);
-   Double_t E1         = TMath::Sqrt(mass1 * mass1 + p1 * p1);
-   Double_t MotherMass = TMath::Sqrt((E0 + E1) * (E0 + E1) - (p0 * p0 + 2.0 * daughter0.Vect().Dot(daughter1.Vect()) + p1 * p1));
-   MotherMass = fSum.M();
+//    // Computes first the invariant mass of the mother
+//    Double_t mass0      = fDaughter[0]->P(useMC).M();
+//    Double_t mass1      = fDaughter[1]->P(useMC).M();
+//    Double_t p0         = daughter0.Vect().Mag();
+//    Double_t p1         = daughter1.Vect().Mag();
+//    Double_t E0         = TMath::Sqrt(mass0 * mass0 + p0 * p0);
+//    Double_t E1         = TMath::Sqrt(mass1 * mass1 + p1 * p1);
+//    Double_t MotherMass = TMath::Sqrt((E0 + E1) * (E0 + E1) - (p0 * p0 + 2.0 * daughter0.Vect().Dot(daughter1.Vect()) + p1 * p1));
+//    MotherMass = fSum.M();
 
    // Computes components of beta
    Double_t betaX = -mother.X() / mother.E();