message: Adding D+ to 3-body decays analysis
diff: #endif
  -enum ESys { kAPL, kPAL, kApXz, kPrAxz, kKKpi, nSys };
  +enum ESys { kAPL, kPAL, kApXz, kPrAxz, kKKpi, kDplus, nSys };
  -const char *sysNames[nSys] = { "APL", "PAL", "ApXz", "PrAxz", "KKpi"};
  -const bool runSys[nSys] = { 1 , 1 , 1 , 1 , 0 };
  -const double distMin[nSys] = { 1.5 , 1.5 , 1.5 , 1.5 , 1.0 }; // everything in GeV here
  -const double distMax[nSys] = { 2.5 , 2.5 , 3.0 , 3.0 , 2.0 };
  -const double distBinWidth[nSys] = { 0.001, 0.001, 0.001 , 0.001 , 0.001 };
  +const char *sysNames[nSys] = { "APL", "PAL", "ApXz", "PrAxz", "D+" , "KKpi"};
  +const bool runSys[nSys] = { 1 , 1 , 0 , 0 , 1 , 0 };
  +const double distMin[nSys] = { 1.5 , 1.5 , 1.5 , 1.5 , 1.0 , 1.0 }; // everything in GeV here
  +const double distMax[nSys] = { 2.5 , 2.5 , 3.0 , 3.0 , 2.8 , 2.0 };
  +const double distBinWidth[nSys] = { 0.001, 0.001, 0.001 , 0.001 , 0.001 , 0.001 };
  -const double dalitzCutMass[nSys]= { 2.148, 2.148, 2.330 , 2.330 , 0 };
  -const double dalitzCutGamma[nSys]={ 0.134, 0.134, 0.100 , 0.100 , 0 };
  +const double dalitzCutMass[nSys]= { 1.148, 2.148, 2.330 , 2.330 ,1.86959, 0 };
  +const double dalitzCutGamma[nSys]={ 0.134, 0.134, 0.100 , 0.100 ,0.00063, 0 };
  const double dalitzMin = 0.0;
  const double dalitzMax = 6.0;
  @@ -439,6 +439,11 @@ void GetParticlesForSystem(ESys system, AliFemtoTrio::EPart &firstParticle, AliF
  secondParticle = AliFemtoTrio::kKaonMinus;
  thirdParticle = AliFemtoTrio::kPionMinus;
  }
  + if(system == kDplus){
  + firstParticle = AliFemtoTrio::kKaonMinus;
  + secondParticle = AliFemtoTrio::kPionPlus;
  + thirdParticle = AliFemtoTrio::kPionPlus;
  + }
  }
  