Added DptDpt QA Task

[u/mrichter/AliRoot.git] / TEvtGen / AliDecayerEvtGen.cxx

diff --git a/TEvtGen/AliDecayerEvtGen.cxx b/TEvtGen/AliDecayerEvtGen.cxx
index 2b4dbbf41560f2219dd5c74256c58b889933bdf3..72968d274f36d6c12be897c4535b02bb26b052ae 100644 (file)
--- a/TEvtGen/AliDecayerEvtGen.cxx
+++ b/TEvtGen/AliDecayerEvtGen.cxx
@@ -158,11 +158,12 @@ Int_t AliDecayerEvtGen::ImportParticles(TClonesArray *particles)
   py=p4.get(2);
   pz=p4.get(3);
   e=p4.get(0);
-
-  x=x4.get(1);//[mm]
-  y=x4.get(2);//[mm]
-  z=x4.get(3);//[mm]
-  t=x4.get(0);//[mm]
+  const Float_t kconvT=0.001/2.999792458e8; // mm/c to seconds conversion
+  const Float_t kconvL=1./10; // mm to cm conversion
+  x=x4.get(1)*kconvL;//[cm]
+  y=x4.get(2)*kconvL;//[cm]
+  z=x4.get(3)*kconvL;//[cm]
+  t=x4.get(0)*kconvT;//[s]
 
   AliDebug(1,Form("partnum = %d istat = %d primaMadre = %d ultimaMadre = %d primaF = %d ultimaF=%d x=%f y=%f z=%f t=%f e=%f px=%f \n",partnum,istat,jmotherfirst,jmotherlast,jdaugfirst,jdauglast,x,y,z,t,e,px));
 
@@ -206,7 +207,10 @@ void AliDecayerEvtGen::ForceDecay()
      break;
      case kBSemiElectronic:
      SetDecayTablePath(gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/DecayTable/BTOELE.DEC"));
-     break;  
+     break;
+     case kHadronicD:
+      SetDecayTablePath(gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/DecayTable/BTOD.DEC"));
+      break;
      case kHardMuons:
      case kChiToJpsiGammaToMuonMuon:
      case kChiToJpsiGammaToElectronElectron:
@@ -224,7 +228,6 @@ void AliDecayerEvtGen::ForceDecay()
      case kWToCharmToMuon:
      case kZDiMuon:
      case kZDiElectron:
-     case kHadronicD:
      case kHadronicDWithout4Bodies:
      case kPhiKK:
      case kOmega:
@@ -234,7 +237,14 @@ void AliDecayerEvtGen::ForceDecay()
      case kNeutralPion:
      case kBPsiPrimeDiElectron:
      case kBeautyUpgrade:
-     AliWarning(Form("Warning: case %s not implemented for this class!",(Char_t)decay));
+     case kBJpsiUndecayed: 
+     case kDiElectronEM:
+     case kElectronEM:
+     case kGammaEM:
+     case kJpsiDiMuon:
+     case kNoDecayBeauty:
+     case kPsiPrimeJpsiDiElectron:
+      AliWarning(Form("Warning: case %d not implemented for this class!",(int)decay));
      break;
      }
      ReadDecayTable();