Some more effc++ warnings.

[u/mrichter/AliRoot.git] / PWG3 / AliBtoJPSItoEle.cxx

diff --git a/PWG3/AliBtoJPSItoEle.cxx b/PWG3/AliBtoJPSItoEle.cxx
index 842e154bcfb6e28aea35895514272a6ce7d0f23d..d31e26c1d92337f182dd27246b2ac33f1cd17e1d 100644 (file)
--- a/PWG3/AliBtoJPSItoEle.cxx
+++ b/PWG3/AliBtoJPSItoEle.cxx
@@ -29,6 +29,7 @@
 #include <TCanvas.h>
 #include <TPaveLabel.h>
 #include <TVector3.h>
+#include <TString.h>
 
 #include "AliBtoJPSItoEle.h"
 
@@ -50,6 +51,7 @@ fWgtJPsi(0.)
 {
   // Default constructor
   
+
   fTrkNum[0] = 0;
   fTrkNum[1] = 0;
 
@@ -82,7 +84,7 @@ fWgtJPsi(0.)
 AliBtoJPSItoEle::AliBtoJPSItoEle(Int_t ev,Int_t trkNum[2],
                       Double_t v1[3],Double_t v2[3], 
                       Double_t dca,
-                      Double_t mom[6],Double_t d0[2]):
+                      Double_t mom[6],Double_t d0[2]) :
 fSignal(kFALSE),
 fJpsiPrimary(kFALSE),
 fEvent(ev),
@@ -100,6 +102,13 @@ fWgtJPsi(0.)
   fTrkNum[0] = trkNum[0];
   fTrkNum[1] = trkNum[1];
 
+  fV1x = v1[0];
+  fV1y = v1[1];
+  fV1z = v1[2];
+  fV2x = v2[0];
+  fV2y = v2[1];
+  fV2z = v2[2];
+
   fPx[0] = mom[0];
   fPy[0] = mom[1];
   fPz[0] = mom[2];
@@ -123,6 +132,7 @@ fWgtJPsi(0.)
   fTagKa[1]  = 0.;
   fTagNid[0] = 0.;
   fTagNid[1] = 0.;
+
 }
 //----------------------------------------------------------------------------
 AliBtoJPSItoEle::~AliBtoJPSItoEle() {}
@@ -380,6 +390,14 @@ void AliBtoJPSItoEle::SetPIDresponse(Double_t resp0[5],Double_t resp1[5]) {
   return;
 } 
 //-----------------------------------------------------------------------------
+Double_t AliBtoJPSItoEle::GetTagEl(Int_t child) const {
+  // Get tag probabilities for electrons
+                                                                                                                          
+  if(child!=0 && child !=1) return -1;
+
+  return fTagEl[child];
+}
+//-----------------------------------------------------------------------------
 void AliBtoJPSItoEle::GetPIDresponse(Double_t resp0[5],Double_t resp1[5]) const {
   // Get combined PID detector response probabilities
                                                                                                                           
@@ -412,10 +430,13 @@ Double_t AliBtoJPSItoEle::TRDTPCCombinedPIDParametrization(Double_t p) const {
   Double_t p3=0.007;
   value=p2+p3*(1.-exp(-TMath::Power(p/p1,4.)));
 
-// Better estimation based on TRD test beam (as presented by Andrea at Munster)
   value/=0.01; // here remove from TPC+TRD the TRD contribution estimated to be 0.01
-  if (p<10.) value*=(1.32-0.18*p+0.076*p*p-0.0037*p*p*p)/100.;
-  if (p>10.) value*=(0.48+0.287*p)/100.;
+// Better estimation based on TRD test beam (as presented by Andrea at Munster)
+//  if (p<10.) value*=(1.32-0.18*p+0.076*p*p-0.0037*p*p*p)/100.;
+//  if (p>10.) value*=(0.48+0.287*p)/100.;
+// From Silvia MASCIOCCHI (Darmstadt) at PWG3 AliceWeek October 2007
+  if (p<10.) value*=(0.44-0.06*p+0.031*p*p-0.0008*p*p*p)/100.;
+  if (p>10.) value*=(-0.67+0.28*p)/100.;
 
   return value;
 }
@@ -424,4 +445,3 @@ Double_t AliBtoJPSItoEle::TRDTPCCombinedPIDParametrization(Double_t p) const {
 
 
 
-