Fixing new line ending
authoragrigora <alina.gabriela.grigoras@cern.ch>
Fri, 18 Jul 2014 09:17:47 +0000 (11:17 +0200)
committeragrigora <alina.gabriela.grigoras@cern.ch>
Fri, 18 Jul 2014 09:17:47 +0000 (11:17 +0200)
PWGCF/Correlations/macros/pidpid/AddAnalysisTaskPidPidCorrelations.C [changed mode: 0755->0644]
PWGCF/FEMTOSCOPY/macros/Train/KchKch_MC_ESD/ConfigFemtoAnalysis.C
PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskHdibaryonLPpi.cxx

old mode 100755 (executable)
new mode 100644 (file)
index ee98b78..f3f1b19
@@ -1,99 +1,99 @@
-// -*- c++ -*-\r\r
-// $Id: AddAnalysisTaskPidPidCorrelations.C\r\r
-\r\r
-const Int_t nBinCent = 1;\r\r
-Double_t centLimits[nBinCent+1] = {0.,5.};\r\r
-// Double_t centLimits[nBinCent+1] = {0.,1.,2.,3.,4.,5.,10.,20.,30.,40.,50.,60.,70.,80.,90.,100.1};\r\r
-const Int_t nBinZvtx = 10;\r\r
-// Double_t zvtxLimits[nBinZvtx+1] = {-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.};\r\r
-Double_t zvtxLimits[nBinZvtx+1] = {-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.};\r\r
-const Int_t nBinPt = 8;\r\r
-Double_t ptLimits[nBinPt+1] = {0.0, 0.5, 1.0, 2.0, 3.0, 4.0,, 5.0, 6.0, 7.0};\r\r
-// Double_t ptLimits[nBinPt+1] = {0.0, 0.1, 0.2, 0.3, 0.35, 0.4, 0.45, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0};\r\r
-const Int_t nBinEta = 20;\r\r
-Double_t etaLimits[nBinEta+1] = {-1.0, -0.9, -0.8, -0.7, -0.6, -0.5, -0.4, -0.3, -0.2, -0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0};\r\r
-\r\r
-\r\r
-AliAnalysisTaskPidPidCorrelations*\r\r
-AddAnalysisTaskPidPidCorrelations(Bool_t mc = kFALSE\r\r
-                               , TString centralityEstimator = "V0M"\r\r
-                               , Bool_t eventMixing = kTRUE\r\r
-                               , Double_t ptMin = 0.2\r\r
-                               , Double_t ptMax = 6.0\r\r
-                               , Double_t etaMin = -0.8\r\r
-                               , Double_t etaMax = 0.8\r\r
-                               , Double_t centrMin = 0.\r\r
-                               , Double_t centrMax = 5.\r\r
-                               , Double_t vertexZ = 10.\r\r
-                               , Int_t triggerPID = 1\r\r
-                               , Int_t assocPID = 1\r\r
-                               )\r\r
-{\r\r
-  // Get the pointer to the existing analysis manager via the static access method.\r\r
-  //===========================================================================\r\r
-  AliAnalysisManager* mgr = AliAnalysisManager::GetAnalysisManager();\r\r
-  if (mgr == NULL) { ::Error("AddAnalysisTaskPidPidCorrelations", "No analysis manager to connect to."); return NULL; }\r\r
-  \r\r
-  // Check the analysis type using the event handlers connected to the analysis manager.\r\r
-  //===========================================================================\r\r
-  if (mgr->GetInputEventHandler() == NULL) { ::Error("AddAnalysisTaskPidPidCorrelations", "This task requires an input event handler"); return NULL; }\r\r
-  \r\r
-  // Create the task, add it to manager and configure it.\r\r
-  //===========================================================================\r\r
-  AliAnalysisTaskPidPidCorrelations *taskPIDCorr = new AliAnalysisTaskPidPidCorrelations(Form("AliAnalysisTaskPidPidCorrelations__%.0f-%.0f_%s",centrMin,centrMax,centralityEstimator.Data()));\r\r
-\r\r
-  taskPIDCorr -> SetTriggerMask(AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral);\r\r
-//     taskPIDCorr->SelectCollisionCandidates(AliVEvent::kMB);//MB  //now inside task\r\r
-  \r\r
-  //______ DATA or MC\r\r
-  taskPIDCorr -> SetMC(mc);\r\r
-\r\r
-  //______ kinematics cut\r\r
-  taskPIDCorr -> SetKinematicsCutsAOD(ptMin,ptMax,etaMin,etaMax);\r\r
-//   taskPIDCorr -> SetEtaRange(-0.9, 0.9);\r\r
-  //_____ vertex\r\r
-  taskPIDCorr -> SetVertexDiamond(.3,.3,vertexZ);  \r\r
-  //____ centrality\r\r
-  taskPIDCorr -> SetCentralityEstimator(centralityEstimator);\r\r
-  taskPIDCorr -> SetCentralityRange(centrMin,centrMax);\r\r
-  taskPIDCorr -> SetEventMixing(eventMixing);\r\r
-//     taskPIDCorr -> SetTriggerRestrictEta(0.5);\r\r
-  taskPIDCorr -> SetEtaOrdering(kFALSE);\r\r
-  taskPIDCorr -> SetPairCuts(kFALSE, kFALSE);\r\r
-  taskPIDCorr -> SetTwoTrackEfficiencyCut(0.02,0.8);\r\r
-  taskPIDCorr -> SetFillpT(kFALSE);\r\r
-  taskPIDCorr -> SetMixingTracks(50000,1000);  \r\r
-  taskPIDCorr -> SetWeightPerEvent(kFALSE);\r\r
-  taskPIDCorr -> SetRejectResonanceDaughters(0);\r\r
-  taskPIDCorr -> SetSelectCharge(0);\r\r
-  taskPIDCorr -> SetSelectTriggerCharge(0);\r\r
-  taskPIDCorr -> SetSelectAssociatedCharge(0);\r\r
-  taskPIDCorr -> SetOnlyOneEtaSide(0);\r\r
-  taskPIDCorr -> SetPtOrder(kFALSE);\r\r
-  taskPIDCorr -> UseMomentumDifferenceCut(kFALSE,0.01);\r\r
-  taskPIDCorr -> SetCentBinning(nBinCent, centLimits);\r\r
-  taskPIDCorr -> SetZvtxBinning(nBinZvtx, zvtxLimits);\r\r
-  taskPIDCorr -> SetPtBinning(nBinPt, ptLimits);\r\r
-  taskPIDCorr -> SetEtaBinning(nBinEta, etaLimits);\r\r
-  taskPIDCorr -> SetPIDsToCorrelate(triggerPID,assocPID);\r\r
-  \r\r
-  // ADD the task\r\r
-  //===========================================================================\r\r
-  mgr -> AddTask(taskPIDCorr);  \r\r
-  \r\r
-  // Create ONLY the output containers for the data produced by the task.\r\r
-  // Get and connect other common input/output containers via the manager as below\r\r
-  //===========================================================================\r\r
-\r\r
-  TString outputFileName = AliAnalysisManager::GetCommonFileName();\r\r
-  outputFileName += ":PWGCFEbyE.outputPIDPIDCorrelations";\r\r
-\r\r
-  AliAnalysisDataContainer *listPIDCorr1 = mgr -> CreateContainer( Form("listPIDPIDCorr_%.0f-%.0f_%s",centrMin,centrMax,centralityEstimator.Data()), TList::Class(), AliAnalysisManager::kOutputContainer, outputFileName.Data());\r\r
-  AliAnalysisDataContainer *listPIDCorr2 = mgr -> CreateContainer( Form("OutputCFCont_%.0f-%.0f_%s",centrMin,centrMax,centralityEstimator.Data()), AliCFContainer::Class(), AliAnalysisManager::kOutputContainer, outputFileName.Data());\r\r
-  \r\r
-  mgr -> ConnectInput( taskPIDCorr, 0, mgr -> GetCommonInputContainer());\r\r
-  mgr -> ConnectOutput( taskPIDCorr, 1, listPIDCorr1);\r\r
-  mgr -> ConnectOutput( taskPIDCorr, 2, listPIDCorr2);\r\r
-  \r\r
-  return taskPIDCorr;\r\r
-}\r\r
+// -*- c++ -*-
+// $Id: AddAnalysisTaskPidPidCorrelations.C
+
+const Int_t nBinCent = 1;
+Double_t centLimits[nBinCent+1] = {0.,5.};
+// Double_t centLimits[nBinCent+1] = {0.,1.,2.,3.,4.,5.,10.,20.,30.,40.,50.,60.,70.,80.,90.,100.1};
+const Int_t nBinZvtx = 10;
+// Double_t zvtxLimits[nBinZvtx+1] = {-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.};
+Double_t zvtxLimits[nBinZvtx+1] = {-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.};
+const Int_t nBinPt = 8;
+Double_t ptLimits[nBinPt+1] = {0.0, 0.5, 1.0, 2.0, 3.0, 4.0,, 5.0, 6.0, 7.0};
+// Double_t ptLimits[nBinPt+1] = {0.0, 0.1, 0.2, 0.3, 0.35, 0.4, 0.45, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0};
+const Int_t nBinEta = 20;
+Double_t etaLimits[nBinEta+1] = {-1.0, -0.9, -0.8, -0.7, -0.6, -0.5, -0.4, -0.3, -0.2, -0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0};
+
+
+AliAnalysisTaskPidPidCorrelations*
+AddAnalysisTaskPidPidCorrelations(Bool_t mc = kFALSE
+                               , TString centralityEstimator = "V0M"
+                               , Bool_t eventMixing = kTRUE
+                               , Double_t ptMin = 0.2
+                               , Double_t ptMax = 6.0
+                               , Double_t etaMin = -0.8
+                               , Double_t etaMax = 0.8
+                               , Double_t centrMin = 0.
+                               , Double_t centrMax = 5.
+                               , Double_t vertexZ = 10.
+                               , Int_t triggerPID = 1
+                               , Int_t assocPID = 1
+                               )
+{
+  // Get the pointer to the existing analysis manager via the static access method.
+  //===========================================================================
+  AliAnalysisManager* mgr = AliAnalysisManager::GetAnalysisManager();
+  if (mgr == NULL) { ::Error("AddAnalysisTaskPidPidCorrelations", "No analysis manager to connect to."); return NULL; }
+  
+  // Check the analysis type using the event handlers connected to the analysis manager.
+  //===========================================================================
+  if (mgr->GetInputEventHandler() == NULL) { ::Error("AddAnalysisTaskPidPidCorrelations", "This task requires an input event handler"); return NULL; }
+  
+  // Create the task, add it to manager and configure it.
+  //===========================================================================
+  AliAnalysisTaskPidPidCorrelations *taskPIDCorr = new AliAnalysisTaskPidPidCorrelations(Form("AliAnalysisTaskPidPidCorrelations__%.0f-%.0f_%s",centrMin,centrMax,centralityEstimator.Data()));
+
+  taskPIDCorr -> SetTriggerMask(AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral);
+//     taskPIDCorr->SelectCollisionCandidates(AliVEvent::kMB);//MB  //now inside task
+  
+  //______ DATA or MC
+  taskPIDCorr -> SetMC(mc);
+
+  //______ kinematics cut
+  taskPIDCorr -> SetKinematicsCutsAOD(ptMin,ptMax,etaMin,etaMax);
+//   taskPIDCorr -> SetEtaRange(-0.9, 0.9);
+  //_____ vertex
+  taskPIDCorr -> SetVertexDiamond(.3,.3,vertexZ);  
+  //____ centrality
+  taskPIDCorr -> SetCentralityEstimator(centralityEstimator);
+  taskPIDCorr -> SetCentralityRange(centrMin,centrMax);
+  taskPIDCorr -> SetEventMixing(eventMixing);
+//     taskPIDCorr -> SetTriggerRestrictEta(0.5);
+  taskPIDCorr -> SetEtaOrdering(kFALSE);
+  taskPIDCorr -> SetPairCuts(kFALSE, kFALSE);
+  taskPIDCorr -> SetTwoTrackEfficiencyCut(0.02,0.8);
+  taskPIDCorr -> SetFillpT(kFALSE);
+  taskPIDCorr -> SetMixingTracks(50000,1000);  
+  taskPIDCorr -> SetWeightPerEvent(kFALSE);
+  taskPIDCorr -> SetRejectResonanceDaughters(0);
+  taskPIDCorr -> SetSelectCharge(0);
+  taskPIDCorr -> SetSelectTriggerCharge(0);
+  taskPIDCorr -> SetSelectAssociatedCharge(0);
+  taskPIDCorr -> SetOnlyOneEtaSide(0);
+  taskPIDCorr -> SetPtOrder(kFALSE);
+  taskPIDCorr -> UseMomentumDifferenceCut(kFALSE,0.01);
+  taskPIDCorr -> SetCentBinning(nBinCent, centLimits);
+  taskPIDCorr -> SetZvtxBinning(nBinZvtx, zvtxLimits);
+  taskPIDCorr -> SetPtBinning(nBinPt, ptLimits);
+  taskPIDCorr -> SetEtaBinning(nBinEta, etaLimits);
+  taskPIDCorr -> SetPIDsToCorrelate(triggerPID,assocPID);
+  
+  // ADD the task
+  //===========================================================================
+  mgr -> AddTask(taskPIDCorr);  
+  
+  // Create ONLY the output containers for the data produced by the task.
+  // Get and connect other common input/output containers via the manager as below
+  //===========================================================================
+
+  TString outputFileName = AliAnalysisManager::GetCommonFileName();
+  outputFileName += ":PWGCFEbyE.outputPIDPIDCorrelations";
+
+  AliAnalysisDataContainer *listPIDCorr1 = mgr -> CreateContainer( Form("listPIDPIDCorr_%.0f-%.0f_%s",centrMin,centrMax,centralityEstimator.Data()), TList::Class(), AliAnalysisManager::kOutputContainer, outputFileName.Data());
+  AliAnalysisDataContainer *listPIDCorr2 = mgr -> CreateContainer( Form("OutputCFCont_%.0f-%.0f_%s",centrMin,centrMax,centralityEstimator.Data()), AliCFContainer::Class(), AliAnalysisManager::kOutputContainer, outputFileName.Data());
+  
+  mgr -> ConnectInput( taskPIDCorr, 0, mgr -> GetCommonInputContainer());
+  mgr -> ConnectOutput( taskPIDCorr, 1, listPIDCorr1);
+  mgr -> ConnectOutput( taskPIDCorr, 2, listPIDCorr2);
+  
+  return taskPIDCorr;
+}
index b72795a..8359b66 100644 (file)
-\r
-/*********************************************************************\r
- *                                                                   *\r
- * ConfigFemtoAnalysis.C - configuration macro for the femtoscopic   *\r
- * analysis, meant as a QA process for two-particle effects          *\r
- *                                                                   *\r
- * Author: Adam Kisiel (Adam.Kisiel@cern.ch)                         *\r
- *                                                                   *\r
- *********************************************************************/\r
-\r
-#if !defined(__CINT__) || defined(__MAKECINT_)\r
-#include "AliFemtoManager.h"\r
-#include "AliFemtoEventReaderESDChain.h"\r
-#include "AliFemtoEventReaderESDChainKine.h"\r
-#include "AliFemtoEventReaderAODChain.h"\r
-#include "AliFemtoSimpleAnalysis.h"\r
-#include "AliFemtoBasicEventCut.h"\r
-#include "AliFemtoESDTrackCut.h"\r
-#include "AliFemtoCorrFctn.h"\r
-#include "AliFemtoCutMonitorParticleYPt.h"\r
-#include "AliFemtoCutMonitorParticleVertPos.h"\r
-#include "AliFemtoCutMonitorParticleMomRes.h"\r
-#include "AliFemtoCutMonitorParticlePID.h"\r
-#include "AliFemtoCutMonitorEventMult.h"\r
-#include "AliFemtoCutMonitorEventVertex.h"\r
-#include "AliFemtoShareQualityTPCEntranceSepPairCut.h"\r
-#include "AliFemtoPairCutAntiGamma.h"\r
-#include "AliFemtoPairCutRadialDistance.h"\r
-#include "AliFemtoQinvCorrFctn.h"\r
-#include "AliFemtoShareQualityCorrFctn.h"\r
-#include "AliFemtoTPCInnerCorrFctn.h"\r
-#include "AliFemtoVertexMultAnalysis.h"\r
-#include "AliFemtoCorrFctn3DSpherical.h"\r
-#include "AliFemtoChi2CorrFctn.h"\r
-#include "AliFemtoCorrFctnTPCNcls.h"\r
-#include "AliFemtoBPLCMS3DCorrFctn.h"\r
-#include "AliFemtoCorrFctn3DLCMSSym.h"\r
-#include "AliFemtoModelBPLCMSCorrFctn.h"\r
-#include "AliFemtoModelCorrFctn3DSpherical.h"\r
-#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"\r
-#include "AliFemtoModelGausRinvFreezeOutGenerator.h"\r
-#include "AliFemtoModelManager.h"\r
-#include "AliFemtoModelWeightGeneratorBasic.h"\r
-#include "AliFemtoModelWeightGeneratorLednicky.h"\r
-#include "AliFemtoCorrFctnDirectYlm.h"\r
-#include "AliFemtoModelCorrFctnDirectYlm.h"\r
-#include "AliFemtoModelCorrFctnSource.h"\r
-#include "AliFemtoModelCorrFctn.h"\r
-#include "AliFemtoCutMonitorParticlePtPDG.h"\r
-#include "AliFemtoKTPairCut.h"\r
-#include "AliFemtoCutMonitorCollections.h"\r
-#endif\r
-\r 
-//_  6 fm  pure kaons _(0.2-0.4)______________________________________________________________________\r
-AliFemtoManager* ConfigFemtoAnalysis() {\r
-\r
-  double PionMass = 0.13956995;\r
-  double KaonMass = 0.493677;\r
-  double ProtonMass = 0.938272013;\r
-\r
-  //multiplicity bins\r
-  int runmults[10] = {1, 0, 0, 0, 0, 0, 0, 0, 0, 0};\r
-  int multbins[11] = {0.01, 200000, 400, 600, 900, 950, 500, 600, 700, 800, 900};\r
-\r
-  int runch[3] = {1, 1, 0};\r
-  const char *chrgs[3] = { "kk", "akak", "kak" };\r
-\r
-  int runktdep = 1;\r
- // double ktrng[8] = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0};\r
+
+/*********************************************************************
+ *                                                                   *
+ * ConfigFemtoAnalysis.C - configuration macro for the femtoscopic   *
+ * analysis, meant as a QA process for two-particle effects          *
+ *                                                                   *
+ * Author: Adam Kisiel (Adam.Kisiel@cern.ch)                         *
+ *                                                                   *
+ *********************************************************************/
+
+#if !defined(__CINT__) || defined(__MAKECINT_)
+#include "AliFemtoManager.h"
+#include "AliFemtoEventReaderESDChain.h"
+#include "AliFemtoEventReaderESDChainKine.h"
+#include "AliFemtoEventReaderAODChain.h"
+#include "AliFemtoSimpleAnalysis.h"
+#include "AliFemtoBasicEventCut.h"
+#include "AliFemtoESDTrackCut.h"
+#include "AliFemtoCorrFctn.h"
+#include "AliFemtoCutMonitorParticleYPt.h"
+#include "AliFemtoCutMonitorParticleVertPos.h"
+#include "AliFemtoCutMonitorParticleMomRes.h"
+#include "AliFemtoCutMonitorParticlePID.h"
+#include "AliFemtoCutMonitorEventMult.h"
+#include "AliFemtoCutMonitorEventVertex.h"
+#include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
+#include "AliFemtoPairCutAntiGamma.h"
+#include "AliFemtoPairCutRadialDistance.h"
+#include "AliFemtoQinvCorrFctn.h"
+#include "AliFemtoShareQualityCorrFctn.h"
+#include "AliFemtoTPCInnerCorrFctn.h"
+#include "AliFemtoVertexMultAnalysis.h"
+#include "AliFemtoCorrFctn3DSpherical.h"
+#include "AliFemtoChi2CorrFctn.h"
+#include "AliFemtoCorrFctnTPCNcls.h"
+#include "AliFemtoBPLCMS3DCorrFctn.h"
+#include "AliFemtoCorrFctn3DLCMSSym.h"
+#include "AliFemtoModelBPLCMSCorrFctn.h"
+#include "AliFemtoModelCorrFctn3DSpherical.h"
+#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
+#include "AliFemtoModelGausRinvFreezeOutGenerator.h"
+#include "AliFemtoModelManager.h"
+#include "AliFemtoModelWeightGeneratorBasic.h"
+#include "AliFemtoModelWeightGeneratorLednicky.h"
+#include "AliFemtoCorrFctnDirectYlm.h"
+#include "AliFemtoModelCorrFctnDirectYlm.h"
+#include "AliFemtoModelCorrFctnSource.h"
+#include "AliFemtoModelCorrFctn.h"
+#include "AliFemtoCutMonitorParticlePtPDG.h"
+#include "AliFemtoKTPairCut.h"
+#include "AliFemtoCutMonitorCollections.h"
+#endif
+//_  6 fm  pure kaons _(0.2-0.4)______________________________________________________________________
+AliFemtoManager* ConfigFemtoAnalysis() {
+
+  double PionMass = 0.13956995;
+  double KaonMass = 0.493677;
+  double ProtonMass = 0.938272013;
+
+  //multiplicity bins
+  int runmults[10] = {1, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+  int multbins[11] = {0.01, 200000, 400, 600, 900, 950, 500, 600, 700, 800, 900};
+
+  int runch[3] = {1, 1, 0};
+  const char *chrgs[3] = { "kk", "akak", "kak" };
+
+  int runktdep = 1;
+ // double ktrng[8] = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0};
   double ktrng[4] = {0.2, 0.4, 0.6, 1.0};
-\r
-  int run3d = 0; // Do 3D cartesian analysis?\r
-  int runshlcms = 1;\r
-\r
-  int nbinssh = 100;\r
-  double shqmax = 0.5;\r
-\r
-\r
-  AliFemtoEventReaderESDChainKine *Reader = new AliFemtoEventReaderESDChainKine();\r
+
+  int run3d = 0; // Do 3D cartesian analysis?
+  int runshlcms = 1;
+
+  int nbinssh = 100;
+  double shqmax = 0.5;
+
+
+  AliFemtoEventReaderESDChainKine *Reader = new AliFemtoEventReaderESDChainKine();
   Reader->SetKaonAnalysis(kTRUE);
-  //Reader->SetUseMultiplicity(AliFemtoEventReaderESDChainKine::kCentrality);\r
-  //Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);\r
-\r
-  // AliFemtoEventReaderKinematicsChain *Reader = new AliFemtoEventReaderKinematicsChain();\r
-  // Reader->SetUseMultiplicity(AliFemtoEventReaderKinematicsChain::kGlobalCount);\r
-  // // // Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);\r
-\r
-  AliFemtoModelGausRinvFreezeOutGenerator *tFreeze = new AliFemtoModelGausRinvFreezeOutGenerator();\r
-  tFreeze->SetSizeInv(5.0);\r
-\r
- // AliFemtoModelWeightGeneratorLednicky *tWeight = new AliFemtoModelWeightGeneratorLednicky();\r
+  //Reader->SetUseMultiplicity(AliFemtoEventReaderESDChainKine::kCentrality);
+  //Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);
+
+  // AliFemtoEventReaderKinematicsChain *Reader = new AliFemtoEventReaderKinematicsChain();
+  // Reader->SetUseMultiplicity(AliFemtoEventReaderKinematicsChain::kGlobalCount);
+  // // // Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);
+
+  AliFemtoModelGausRinvFreezeOutGenerator *tFreeze = new AliFemtoModelGausRinvFreezeOutGenerator();
+  tFreeze->SetSizeInv(5.0);
+
+ // AliFemtoModelWeightGeneratorLednicky *tWeight = new AliFemtoModelWeightGeneratorLednicky();
   // tWeight->SetPairType(AliFemtoModelWeightGenerator::KaonPlusKaonPlus());
    
     AliFemtoModelWeightGeneratorBasic *tWeight = new AliFemtoModelWeightGeneratorBasic();
         tWeight->SetPairType(AliFemtoModelWeightGenerator::KaonPlusKaonPlus());
    
- //////// tWeight->SetPairType(5);\r
-  //tWeight->SetPairType(1);\r
-\r
-  AliFemtoModelManager *tModelManager = new AliFemtoModelManager();\r
-  tModelManager->AcceptFreezeOutGenerator(tFreeze);\r
-  tModelManager->AcceptWeightGenerator(tWeight);\r
-  tModelManager->CreateCopyHiddenInfo(kFALSE);\r
-  //tModelManager->CreateCopyHiddenInfo(kTRUE);\r
-\r
-  AliFemtoManager* Manager=new AliFemtoManager();\r
-  Manager->SetEventReader(Reader);\r
-\r
-  AliFemtoVertexMultAnalysis    *anetaphitpc[30];\r
-  AliFemtoBasicEventCut         *mecetaphitpc[30];\r
-  AliFemtoCutMonitorEventMult   *cutPassEvMetaphitpc[30];\r
-  AliFemtoCutMonitorEventMult   *cutFailEvMetaphitpc[30];\r
-  AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[30];\r
-  AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[30];\r
-  AliFemtoCutMonitorCollections   *cutPassColletaphitpc[30];\r
-  AliFemtoCutMonitorCollections   *cutFailColletaphitpc[30];\r
-  AliFemtoESDTrackCut           *dtc1etaphitpc[30];\r
-  AliFemtoESDTrackCut           *dtc2etaphitpc[30];\r
-  AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[30];\r
-  AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[30];\r
-  AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[30];\r
-  AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[30];\r
-  AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[30];\r
-  AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[30];\r
-  AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[30];\r
-  AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[30];\r
-  //  AliFemtoPairCutAntiGamma      *sqpcetaphitpc[30];\r
-  AliFemtoShareQualityPairCut      *sqpcetaphitpc[30];\r
-  //AliFemtoPairCutRadialDistance      *sqpcetaphitpc[30];\r
-  AliFemtoModelCorrFctnDirectYlm     *cylmetaphitpc[30];\r
-  AliFemtoCorrFctnDEtaDPhi      *cdedpetaphi[30];\r
-  AliFemtoChi2CorrFctn          *cchiqinvetaphitpc[30];\r
-  AliFemtoKTPairCut             *ktpcuts[30*7];\r
-  AliFemtoModelCorrFctnDirectYlm     *cylmkttpc[30*7];\r
-  //AliFemtoModelQinvCorrFctn          *cqinvkttpc[30*7];\r
-  AliFemtoModelCorrFctnSource          *cqinvkttpc[30*7];\r
-//  AliFemtoModelCorrFctn          *cqinvkttpc[30*7];\r
-  //AliFemtoModelCorrFctn3DLCMSSym     *cq3dlcmskttpc[30*7];\r
-  AliFemtoCorrFctnTPCNcls       *cqinvnclstpc[30];\r
-  AliFemtoShareQualityCorrFctn  *cqinvsqtpc[30*10];\r
-  AliFemtoChi2CorrFctn          *cqinvchi2tpc[30];\r
-  AliFemtoTPCInnerCorrFctn      *cqinvinnertpc[30*10];\r
-\r
-\r
-  // *** Begin pion-pion analysis ***\r
-  int aniter = 0;\r
-\r
-  for (int imult=0; imult<10; imult++) {\r
-    if (runmults[imult]) {\r
-      for (int ichg=0; ichg<3; ichg++) {\r
-        if (runch[ichg]) {\r
-          aniter = ichg*10+imult;\r
-\r
-          anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(10, -8.0, 8.0, 4, multbins[imult], multbins[imult+1]);\r
-          anetaphitpc[aniter]->SetNumEventsToMix(5);\r
-          anetaphitpc[aniter]->SetMinSizePartCollection(1);\r
-\r
-          mecetaphitpc[aniter] = new AliFemtoBasicEventCut();\r
-          mecetaphitpc[aniter]->SetEventMult(0,200000);\r
-          mecetaphitpc[aniter]->SetVertZPos(-8,8);\r
-\r
-\r
-          dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();\r
-          dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();\r
+ //////// tWeight->SetPairType(5);
+  //tWeight->SetPairType(1);
+
+  AliFemtoModelManager *tModelManager = new AliFemtoModelManager();
+  tModelManager->AcceptFreezeOutGenerator(tFreeze);
+  tModelManager->AcceptWeightGenerator(tWeight);
+  tModelManager->CreateCopyHiddenInfo(kFALSE);
+  //tModelManager->CreateCopyHiddenInfo(kTRUE);
+
+  AliFemtoManager* Manager=new AliFemtoManager();
+  Manager->SetEventReader(Reader);
+
+  AliFemtoVertexMultAnalysis    *anetaphitpc[30];
+  AliFemtoBasicEventCut         *mecetaphitpc[30];
+  AliFemtoCutMonitorEventMult   *cutPassEvMetaphitpc[30];
+  AliFemtoCutMonitorEventMult   *cutFailEvMetaphitpc[30];
+  AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[30];
+  AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[30];
+  AliFemtoCutMonitorCollections   *cutPassColletaphitpc[30];
+  AliFemtoCutMonitorCollections   *cutFailColletaphitpc[30];
+  AliFemtoESDTrackCut           *dtc1etaphitpc[30];
+  AliFemtoESDTrackCut           *dtc2etaphitpc[30];
+  AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[30];
+  AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[30];
+  AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[30];
+  AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[30];
+  AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[30];
+  AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[30];
+  AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[30];
+  AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[30];
+  //  AliFemtoPairCutAntiGamma      *sqpcetaphitpc[30];
+  AliFemtoShareQualityPairCut      *sqpcetaphitpc[30];
+  //AliFemtoPairCutRadialDistance      *sqpcetaphitpc[30];
+  AliFemtoModelCorrFctnDirectYlm     *cylmetaphitpc[30];
+  AliFemtoCorrFctnDEtaDPhi      *cdedpetaphi[30];
+  AliFemtoChi2CorrFctn          *cchiqinvetaphitpc[30];
+  AliFemtoKTPairCut             *ktpcuts[30*7];
+  AliFemtoModelCorrFctnDirectYlm     *cylmkttpc[30*7];
+  //AliFemtoModelQinvCorrFctn          *cqinvkttpc[30*7];
+  AliFemtoModelCorrFctnSource          *cqinvkttpc[30*7];
+//  AliFemtoModelCorrFctn          *cqinvkttpc[30*7];
+  //AliFemtoModelCorrFctn3DLCMSSym     *cq3dlcmskttpc[30*7];
+  AliFemtoCorrFctnTPCNcls       *cqinvnclstpc[30];
+  AliFemtoShareQualityCorrFctn  *cqinvsqtpc[30*10];
+  AliFemtoChi2CorrFctn          *cqinvchi2tpc[30];
+  AliFemtoTPCInnerCorrFctn      *cqinvinnertpc[30*10];
+
+
+  // *** Begin pion-pion analysis ***
+  int aniter = 0;
+
+  for (int imult=0; imult<10; imult++) {
+    if (runmults[imult]) {
+      for (int ichg=0; ichg<3; ichg++) {
+        if (runch[ichg]) {
+          aniter = ichg*10+imult;
+
+          anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(10, -8.0, 8.0, 4, multbins[imult], multbins[imult+1]);
+          anetaphitpc[aniter]->SetNumEventsToMix(5);
+          anetaphitpc[aniter]->SetMinSizePartCollection(1);
+
+          mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
+          mecetaphitpc[aniter]->SetEventMult(0,200000);
+          mecetaphitpc[aniter]->SetVertZPos(-8,8);
+
+
+          dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
+          dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
 
 
     
-          if (ichg == 0) {\r
-            dtc1etaphitpc[aniter]->SetCharge(1.0);\r
-            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);\r
-          }\r
-          else if (ichg == 1) {\r
-            dtc1etaphitpc[aniter]->SetCharge(-1.0);\r
-            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);\r
-          }\r
-          else if (ichg == 2) {\r
-            dtc1etaphitpc[aniter]->SetCharge(-1.0);\r
-            dtc2etaphitpc[aniter]->SetCharge(1.0);\r
-            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);\r
-            dtc2etaphitpc[aniter]->SetPt(0.2,1.5);\r
-          }\r
-\r
-          dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);\r
-          dtc1etaphitpc[aniter]->SetMass(KaonMass);\r
-//test          dtc1etaphitpc[aniter]->SetMostProbableKaon();\r
-//test          dtc1etaphitpc[aniter]->SetNsigma(3.0);\r
-          //dtc1etaphitpc[aniter]->SetNsigma(2.0);\r
-//test          dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);\r
-          //dtc1etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);\r
-\r
-          if (ichg == 2) {\r
-            dtc2etaphitpc[aniter]->SetEta(-0.8,0.8);\r
-            dtc2etaphitpc[aniter]->SetMass(KaonMass);\r
-//test            dtc2etaphitpc[aniter]->SetMostProbableKaon();\r
-//test            dtc2etaphitpc[aniter]->SetNsigma(3.0);\r
-            //dtc2etaphitpc[aniter]->SetNsigma(2.0);\r
-//test            dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);\r
-            //dtc2etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);\r
-\r
-          }\r
-\r
-\r
-           dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);\r
-           dtc1etaphitpc[aniter]->SetminTPCncls(80);\r
-           dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);\r
-           dtc1etaphitpc[aniter]->SetLabel(kFALSE);\r
-           dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);\r
-           dtc1etaphitpc[aniter]->SetMaxImpactXY(2.4); // 2.4 0.1\r
-           dtc1etaphitpc[aniter]->SetMaxImpactZ(3.2); // 2.0 0.1\r
+          if (ichg == 0) {
+            dtc1etaphitpc[aniter]->SetCharge(1.0);
+            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
+          }
+          else if (ichg == 1) {
+            dtc1etaphitpc[aniter]->SetCharge(-1.0);
+            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
+          }
+          else if (ichg == 2) {
+            dtc1etaphitpc[aniter]->SetCharge(-1.0);
+            dtc2etaphitpc[aniter]->SetCharge(1.0);
+            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
+            dtc2etaphitpc[aniter]->SetPt(0.2,1.5);
+          }
+
+          dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
+          dtc1etaphitpc[aniter]->SetMass(KaonMass);
+//test          dtc1etaphitpc[aniter]->SetMostProbableKaon();
+//test          dtc1etaphitpc[aniter]->SetNsigma(3.0);
+          //dtc1etaphitpc[aniter]->SetNsigma(2.0);
+//test          dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
+          //dtc1etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);
+
+          if (ichg == 2) {
+            dtc2etaphitpc[aniter]->SetEta(-0.8,0.8);
+            dtc2etaphitpc[aniter]->SetMass(KaonMass);
+//test            dtc2etaphitpc[aniter]->SetMostProbableKaon();
+//test            dtc2etaphitpc[aniter]->SetNsigma(3.0);
+            //dtc2etaphitpc[aniter]->SetNsigma(2.0);
+//test            dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
+            //dtc2etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);
+
+          }
+
+
+           dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
+           dtc1etaphitpc[aniter]->SetminTPCncls(80);
+           dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
+           dtc1etaphitpc[aniter]->SetLabel(kFALSE);
+           dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
+           dtc1etaphitpc[aniter]->SetMaxImpactXY(2.4); // 2.4 0.1
+           dtc1etaphitpc[aniter]->SetMaxImpactZ(3.2); // 2.0 0.1
 
            AliFemtoCutMonitorParticleYPt *cutPassYPtkp = new AliFemtoCutMonitorParticleYPt("cutPasskp", 0.493677);
            AliFemtoCutMonitorParticleYPt *cutFailYPtkp = new AliFemtoCutMonitorParticleYPt("cutFailkp", 0.493677);
@@ -221,64 +221,64 @@ AliFemtoManager* ConfigFemtoAnalysis() {
            dtc1etaphitpc[aniter]->AddCutMonitor(cutPassMRkp, cutFailMRkp);
 
 
-          sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();\r
-          //sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();\r
-          sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);\r
-          sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);\r
-          sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);\r
-
-\r
-          anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);\r
-          if (ichg == 2) {\r
-            anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);\r
-            anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);\r
-          }\r
-          else {\r
-            anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);\r
-            anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);\r
-          }\r
-\r
-          anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);\r
-\r
-\r
-          //Qinv (without kT bins)\r
-          cqinvkttpc[aniter] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);\r
-          //cqinvkttpc[aniter] = new AliFemtoModelCorrFctn(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);\r
-        //  ktpcuts[0] = new AliFemtoKTPairCut(0.2,0.4);\r
-         // cqinvkttpc[aniter]->SetPairSelectionCut(ktpcuts[0]);\r
-          cqinvkttpc[aniter]->ConnectToManager(tModelManager);\r
-          anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[aniter]);\r
-\r
-\r
-          if (runktdep) {\r
-            int ktm;\r
-            for (int ikt=0; ikt<3; ikt++) {\r
-              ktm = aniter*3 + ikt;\r
-              ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);\r
-
-              cqinvkttpc[ktm] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0, shqmax);\r
-              cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);\r
-              cqinvkttpc[ktm]->ConnectToManager(tModelManager);\r
-              anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);\r
-
-              //if (run3d) {\r
-              //               cq3dlcmskttpc[ktm] = new AliFemtoCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),60,(imult>3)?((imult>6)?((imult>7)?0.6:0.4):0.25):0.15);\r
-              //cq3dlcmskttpc[ktm] = new AliFemtoModelCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),60,0.5);\r
-              //cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);\r
-              //anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);\r
-              //}\r
-            }\r
-          }\r
-\r
-          //cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),39, 39);\r
-          //anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);\r
-\r
-          Manager->AddAnalysis(anetaphitpc[aniter]);\r
-        }\r
-      }\r
-    }\r
-  }\r
-  // *** End pion-pion analysis\r
-\r
-  return Manager;\r
-}\r
+          sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();
+          //sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
+          sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
+          sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
+          sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
+
+
+          anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
+          if (ichg == 2) {
+            anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
+            anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
+          }
+          else {
+            anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
+            anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
+          }
+
+          anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
+
+
+          //Qinv (without kT bins)
+          cqinvkttpc[aniter] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
+          //cqinvkttpc[aniter] = new AliFemtoModelCorrFctn(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
+        //  ktpcuts[0] = new AliFemtoKTPairCut(0.2,0.4);
+         // cqinvkttpc[aniter]->SetPairSelectionCut(ktpcuts[0]);
+          cqinvkttpc[aniter]->ConnectToManager(tModelManager);
+          anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[aniter]);
+
+
+          if (runktdep) {
+            int ktm;
+            for (int ikt=0; ikt<3; ikt++) {
+              ktm = aniter*3 + ikt;
+              ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
+
+              cqinvkttpc[ktm] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0, shqmax);
+              cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
+              cqinvkttpc[ktm]->ConnectToManager(tModelManager);
+              anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
+
+              //if (run3d) {
+              //               cq3dlcmskttpc[ktm] = new AliFemtoCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),60,(imult>3)?((imult>6)?((imult>7)?0.6:0.4):0.25):0.15);
+              //cq3dlcmskttpc[ktm] = new AliFemtoModelCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),60,0.5);
+              //cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
+              //anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
+              //}
+            }
+          }
+
+          //cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),39, 39);
+          //anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
+
+          Manager->AddAnalysis(anetaphitpc[aniter]);
+        }
+      }
+    }
+  }
+  // *** End pion-pion analysis
+
+  return Manager;
+}
index 8c795eb..b3fb58d 100644 (file)
-/**************************************************************************\r
- * Author : Benjamin Dönigus (benjamin.doenigus@cern.ch)                  *\r
- *                                                                        *\r
- * Contributors are mentioned in the code where appropriate.              *\r
- *                                                                        *\r
- * Permission to use, copy, modify and distribute this software and its   *\r
- * documentation strictly for non-commercial purposes is hereby granted   *\r
- * without fee, provided that the above copyright notice appears in all   *\r
- * copies and that both the copyright notice and this permission notice   *\r
- * appear in the supporting documentation. The authors make no claims     *\r
- * about the suitability of this software for any purpose. It is          *\r
- * provided "as is" without express or implied warranty.                  *\r
- **************************************************************************/\r
-\r
-//-----------------------------------------------------------------\r
-//                 AliAnalysisTaskHdibaryonLPpi class\r
-//          used to search for the H-Dibaryon in weak \r
-//          (Lambda Proton Pion) and strong (Lambda Lambda) decays\r
-//-----------------------------------------------------------------\r
-\r
-#include "Riostream.h"\r
-#include "TROOT.h"\r
-#include "TChain.h"\r
-#include "TStyle.h"\r
-#include "TSystem.h"\r
-#include "TTree.h"\r
-#include "TH1F.h"\r
-#include "TH1.h"\r
-#include "TH2D.h"\r
-#include "TH3F.h"\r
-#include "TCanvas.h"\r
-#include "TParticle.h"\r
-#include "TNtuple.h"\r
-#include "TObjString.h"\r
-#include "TLorentzVector.h"\r
-#include "TDatabasePDG.h"\r
-\r
-#include "AliAnalysisTask.h"\r
-#include "AliAnalysisManager.h"\r
-#include "AliMCEventHandler.h"\r
-#include "AliMCEvent.h"\r
-\r
-#include "AliESD.h"\r
-#include "AliESDpid.h"\r
-#include "AliESDEvent.h"\r
-#include "AliGenEventHeader.h"\r
-#include "AliESDInputHandler.h"\r
-#include "AliLog.h"\r
-#include "AliStack.h"\r
-#include "AliHeader.h"\r
-\r
-#include "AliKFParticle.h"\r
-#include "AliKFVertex.h"\r
-#include "AliESDtrackCuts.h"\r
-#include "AliESDv0Cuts.h"\r
-#include "AliAnalysisTaskHdibaryonLPpi.h"\r
-\r
-#include "TFile.h"\r
-#include "TH2F.h"\r
-#include "TF1.h"\r
-#include "TCanvas.h"\r
-#include "TList.h"\r
-#include "TParallelCoord.h"\r
-\r
-#include "AliMCParticle.h"\r
-#include "AliGenPythiaEventHeader.h"\r
-\r
-#include "AliPID.h"\r
-#include "AliESDtrack.h"\r
-#include "AliCentrality.h"\r
-\r
-#include "THnSparse.h"\r
-\r
-#include "AliVertexerTracks.h"\r
-\r
-using namespace std;\r
-ClassImp(AliAnalysisTaskHdibaryonLPpi)\r
-//________________________________________________________________________\r
-AliAnalysisTaskHdibaryonLPpi::AliAnalysisTaskHdibaryonLPpi() : AliAnalysisTaskSE()/*AliAnalysisTask(name, ""), fMCEvent(0)*/, fESD(0),   fESDtrackCutsV0(0),\r
-  fESDCutsV0(0),\r
-  fEsdTrackCuts(0),\r
-  fBin(0),\r
-  fEvent(0x0),\r
-  fHistList(0),  \r
-  fHistMassDPi(0), \r
-  fHistMassLPi(0),\r
-  fHistMassLambdaPi(0),\r
-  fHistMassLambda(0),\r
-  fHistMassLambdaPPi(0),\r
-  fHistMassLambdaP(0),\r
-  fHistMassLambdaK(0),\r
-  fHistMassK0onl(0),\r
-  fHistMassK0offl(0),\r
-  fHistMassK0onlC(0),\r
-  fHistMassK0offlC(0),\r
-  fHistMassPQonl(0), \r
-  fHistMassPQoffl(0),\r
-  fHistDC(0),\r
-  fHistArmenterosPodolanski(0),\r
-  fHistArmenterosPodolanskiCut(0), \r
-  fHistHDibaryonInvaMassGen(0),\r
-  fHistHDibaryonInvaMassGenRes(0),\r
-  fHistAntiHDibaryonInvaMassGen(0),\r
-  fHistHDibaryonInvaMassAso(0),\r
-  fHistHDibaryonInvaMassAsoReso(0),\r
-  fHistAntiHDibaryonInvaMassAso(0),\r
-  fHistCheck(0),\r
-  fHistHPointingAngle(0),\r
-  fHistMassH(0),\r
-  fHistMassLambdaFromH(0),\r
-  fHistMassLambdaFromHtLorentz(0),\r
-  fHistMassPpi(0),\r
-  fHistMassPpiReso(0),\r
-  fHistMassLpi(0),\r
-  fHistMassLP(0),\r
-  fHistProtonPIDBb(0),\r
-  fHistPionPIDBb(0),\r
-  fHistProtonPIDLambda(0),\r
-  fHistPionPIDLambda(0),\r
-  fHistMCdcaPvtxDvtx(0),\r
-  fHistMCdcaPvtxLvtx(0),\r
-  fHistMCdcaDvtxLvtx(0),\r
-  fHistMCangleLH(0),\r
-  fHistMCdecayAngle(0),\r
-  fHistMCpointingAngle(0),\r
-  fHistMCap(0),\r
-  fHistMCdcaPvtxDvtxReso(0),\r
-  fHistMCdcaPvtxLvtxReso(0),\r
-  fHistMCdcaDvtxLvtxReso(0),\r
-  fHistMCangleLHReso(0),\r
-  fHistMCdecayAngleReso(0),\r
-  fHistMCpointingAngleReso(0),\r
-  fHistMCapReso(0),\r
-  fHistCentrality(0),\r
-  fHistCentralityAC(0), \r
-  fHistMultiplicity(0),\r
-  fHistHilf1(0),\r
-  fHistHilf2(0), \r
-  fHistHilf3(0),\r
-  fHistHilf4(0),\r
-  fHistHilf5(0), \r
-  fHistHilf6(0),\r
-  fHistPtvsEtaGen(0),\r
-  fHistPtvsEtaAso(0),\r
-  fHistPtvsYGen(0),\r
-  fHistPtvsYAso(0), \r
-  fHistRap(0),\r
-  fHistCount(0),\r
-  fPIDtpcESD(0),\r
-  fHistTriggerStat(0),\r
-  fHistTriggerStatAfterEventSelection(0), \r
-  fHistMassHcentMult(0),\r
-  fHistNdim(0)\r
-{\r
-  // DefaultConstructor\r
-\r
-}\r
-\r
-//________________________________________________________________________\r
-AliAnalysisTaskHdibaryonLPpi::AliAnalysisTaskHdibaryonLPpi(const char *name) : AliAnalysisTaskSE(name)/*AliAnalysisTask(name, ""), fMCEvent(0)*/, fESD(0),   fESDtrackCutsV0(0),\r
-  fESDCutsV0(0),\r
-  fEsdTrackCuts(0),\r
-  fBin(0),\r
-  fEvent(0x0),\r
-  fHistList(0),  \r
-  fHistMassDPi(0), \r
-  fHistMassLPi(0),\r
-  fHistMassLambdaPi(0),\r
-  fHistMassLambda(0),\r
-  fHistMassLambdaPPi(0),\r
-  fHistMassLambdaP(0),\r
-  fHistMassLambdaK(0),\r
-  fHistMassK0onl(0),\r
-  fHistMassK0offl(0),\r
-  fHistMassK0onlC(0),\r
-  fHistMassK0offlC(0),\r
-  fHistMassPQonl(0), \r
-  fHistMassPQoffl(0),\r
-  fHistDC(0),\r
-  fHistArmenterosPodolanski(0),\r
-  fHistArmenterosPodolanskiCut(0), \r
-  fHistHDibaryonInvaMassGen(0),\r
-  fHistHDibaryonInvaMassGenRes(0),\r
-  fHistAntiHDibaryonInvaMassGen(0),\r
-  fHistHDibaryonInvaMassAso(0),\r
-  fHistHDibaryonInvaMassAsoReso(0),\r
-  fHistAntiHDibaryonInvaMassAso(0),\r
-  fHistCheck(0),\r
-  fHistHPointingAngle(0),\r
-  fHistMassH(0),\r
-  fHistMassLambdaFromH(0),\r
-  fHistMassLambdaFromHtLorentz(0),\r
-  fHistMassPpi(0),\r
-  fHistMassPpiReso(0),\r
-  fHistMassLpi(0),\r
-  fHistMassLP(0),\r
-  fHistProtonPIDBb(0),\r
-  fHistPionPIDBb(0),\r
-  fHistProtonPIDLambda(0),\r
-  fHistPionPIDLambda(0),\r
-  fHistMCdcaPvtxDvtx(0),\r
-  fHistMCdcaPvtxLvtx(0),\r
-  fHistMCdcaDvtxLvtx(0),\r
-  fHistMCangleLH(0),\r
-  fHistMCdecayAngle(0),\r
-  fHistMCpointingAngle(0),\r
-  fHistMCap(0),\r
-  fHistMCdcaPvtxDvtxReso(0),\r
-  fHistMCdcaPvtxLvtxReso(0),\r
-  fHistMCdcaDvtxLvtxReso(0),\r
-  fHistMCangleLHReso(0),\r
-  fHistMCdecayAngleReso(0),\r
-  fHistMCpointingAngleReso(0),\r
-  fHistMCapReso(0),\r
-  fHistCentrality(0),\r
-  fHistCentralityAC(0), \r
-  fHistMultiplicity(0), \r
-  fHistHilf1(0),\r
-  fHistHilf2(0), \r
-  fHistHilf3(0),\r
-  fHistHilf4(0),\r
-  fHistHilf5(0), \r
-  fHistHilf6(0),\r
-  fHistPtvsEtaGen(0),\r
-  fHistPtvsEtaAso(0),\r
-  fHistPtvsYGen(0),\r
-  fHistPtvsYAso(0), \r
-  fHistRap(0),\r
-  fHistCount(0),\r
-  fPIDtpcESD(0),\r
-  fHistTriggerStat(0),\r
-  fHistTriggerStatAfterEventSelection(0),\r
-  fHistMassHcentMult(0),\r
-  fHistNdim(0)\r
-\r
-{\r
-  // Constructor\r
-\r
-  // Define input and output slots here\r
-  // Input from a TChain\r
-  DefineInput(0, TChain::Class());\r
-  // Output to TList container\r
-  DefineOutput(1, TList::Class()); //full\r
-\r
-  //MC info contol\r
-  if (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())SetHasMC();\r
-\r
-  //V0 cuts\r
-\r
-  fESDtrackCutsV0 = new AliESDtrackCuts("AliESDtrackCutsV0","AliESDtrackCutsV0");\r
-  fESDtrackCutsV0->SetAcceptKinkDaughters(kFALSE);\r
-  fESDtrackCutsV0->SetMinNClustersTPC(80);\r
-  fESDtrackCutsV0->SetMaxChi2PerClusterTPC(5);\r
-  fESDtrackCutsV0->SetRequireTPCRefit(kTRUE);\r
+/**************************************************************************
+ * Author : Benjamin Dönigus (benjamin.doenigus@cern.ch)                  *
+ *                                                                        *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+//-----------------------------------------------------------------
+//                 AliAnalysisTaskHdibaryonLPpi class
+//          used to search for the H-Dibaryon in weak 
+//          (Lambda Proton Pion) and strong (Lambda Lambda) decays
+//-----------------------------------------------------------------
+
+#include "Riostream.h"
+#include "TROOT.h"
+#include "TChain.h"
+#include "TStyle.h"
+#include "TSystem.h"
+#include "TTree.h"
+#include "TH1F.h"
+#include "TH1.h"
+#include "TH2D.h"
+#include "TH3F.h"
+#include "TCanvas.h"
+#include "TParticle.h"
+#include "TNtuple.h"
+#include "TObjString.h"
+#include "TLorentzVector.h"
+#include "TDatabasePDG.h"
+
+#include "AliAnalysisTask.h"
+#include "AliAnalysisManager.h"
+#include "AliMCEventHandler.h"
+#include "AliMCEvent.h"
+
+#include "AliESD.h"
+#include "AliESDpid.h"
+#include "AliESDEvent.h"
+#include "AliGenEventHeader.h"
+#include "AliESDInputHandler.h"
+#include "AliLog.h"
+#include "AliStack.h"
+#include "AliHeader.h"
+
+#include "AliKFParticle.h"
+#include "AliKFVertex.h"
+#include "AliESDtrackCuts.h"
+#include "AliESDv0Cuts.h"
+#include "AliAnalysisTaskHdibaryonLPpi.h"
+
+#include "TFile.h"
+#include "TH2F.h"
+#include "TF1.h"
+#include "TCanvas.h"
+#include "TList.h"
+#include "TParallelCoord.h"
+
+#include "AliMCParticle.h"
+#include "AliGenPythiaEventHeader.h"
+
+#include "AliPID.h"
+#include "AliESDtrack.h"
+#include "AliCentrality.h"
+
+#include "THnSparse.h"
+
+#include "AliVertexerTracks.h"
+
+using namespace std;
+ClassImp(AliAnalysisTaskHdibaryonLPpi)
+//________________________________________________________________________
+AliAnalysisTaskHdibaryonLPpi::AliAnalysisTaskHdibaryonLPpi() : AliAnalysisTaskSE()/*AliAnalysisTask(name, ""), fMCEvent(0)*/, fESD(0),   fESDtrackCutsV0(0),
+  fESDCutsV0(0),
+  fEsdTrackCuts(0),
+  fBin(0),
+  fEvent(0x0),
+  fHistList(0),  
+  fHistMassDPi(0), 
+  fHistMassLPi(0),
+  fHistMassLambdaPi(0),
+  fHistMassLambda(0),
+  fHistMassLambdaPPi(0),
+  fHistMassLambdaP(0),
+  fHistMassLambdaK(0),
+  fHistMassK0onl(0),
+  fHistMassK0offl(0),
+  fHistMassK0onlC(0),
+  fHistMassK0offlC(0),
+  fHistMassPQonl(0), 
+  fHistMassPQoffl(0),
+  fHistDC(0),
+  fHistArmenterosPodolanski(0),
+  fHistArmenterosPodolanskiCut(0), 
+  fHistHDibaryonInvaMassGen(0),
+  fHistHDibaryonInvaMassGenRes(0),
+  fHistAntiHDibaryonInvaMassGen(0),
+  fHistHDibaryonInvaMassAso(0),
+  fHistHDibaryonInvaMassAsoReso(0),
+  fHistAntiHDibaryonInvaMassAso(0),
+  fHistCheck(0),
+  fHistHPointingAngle(0),
+  fHistMassH(0),
+  fHistMassLambdaFromH(0),
+  fHistMassLambdaFromHtLorentz(0),
+  fHistMassPpi(0),
+  fHistMassPpiReso(0),
+  fHistMassLpi(0),
+  fHistMassLP(0),
+  fHistProtonPIDBb(0),
+  fHistPionPIDBb(0),
+  fHistProtonPIDLambda(0),
+  fHistPionPIDLambda(0),
+  fHistMCdcaPvtxDvtx(0),
+  fHistMCdcaPvtxLvtx(0),
+  fHistMCdcaDvtxLvtx(0),
+  fHistMCangleLH(0),
+  fHistMCdecayAngle(0),
+  fHistMCpointingAngle(0),
+  fHistMCap(0),
+  fHistMCdcaPvtxDvtxReso(0),
+  fHistMCdcaPvtxLvtxReso(0),
+  fHistMCdcaDvtxLvtxReso(0),
+  fHistMCangleLHReso(0),
+  fHistMCdecayAngleReso(0),
+  fHistMCpointingAngleReso(0),
+  fHistMCapReso(0),
+  fHistCentrality(0),
+  fHistCentralityAC(0), 
+  fHistMultiplicity(0),
+  fHistHilf1(0),
+  fHistHilf2(0), 
+  fHistHilf3(0),
+  fHistHilf4(0),
+  fHistHilf5(0), 
+  fHistHilf6(0),
+  fHistPtvsEtaGen(0),
+  fHistPtvsEtaAso(0),
+  fHistPtvsYGen(0),
+  fHistPtvsYAso(0), 
+  fHistRap(0),
+  fHistCount(0),
+  fPIDtpcESD(0),
+  fHistTriggerStat(0),
+  fHistTriggerStatAfterEventSelection(0), 
+  fHistMassHcentMult(0),
+  fHistNdim(0)
+{
+  // DefaultConstructor
+
+}
+
+//________________________________________________________________________
+AliAnalysisTaskHdibaryonLPpi::AliAnalysisTaskHdibaryonLPpi(const char *name) : AliAnalysisTaskSE(name)/*AliAnalysisTask(name, ""), fMCEvent(0)*/, fESD(0),   fESDtrackCutsV0(0),
+  fESDCutsV0(0),
+  fEsdTrackCuts(0),
+  fBin(0),
+  fEvent(0x0),
+  fHistList(0),  
+  fHistMassDPi(0), 
+  fHistMassLPi(0),
+  fHistMassLambdaPi(0),
+  fHistMassLambda(0),
+  fHistMassLambdaPPi(0),
+  fHistMassLambdaP(0),
+  fHistMassLambdaK(0),
+  fHistMassK0onl(0),
+  fHistMassK0offl(0),
+  fHistMassK0onlC(0),
+  fHistMassK0offlC(0),
+  fHistMassPQonl(0), 
+  fHistMassPQoffl(0),
+  fHistDC(0),
+  fHistArmenterosPodolanski(0),
+  fHistArmenterosPodolanskiCut(0), 
+  fHistHDibaryonInvaMassGen(0),
+  fHistHDibaryonInvaMassGenRes(0),
+  fHistAntiHDibaryonInvaMassGen(0),
+  fHistHDibaryonInvaMassAso(0),
+  fHistHDibaryonInvaMassAsoReso(0),
+  fHistAntiHDibaryonInvaMassAso(0),
+  fHistCheck(0),
+  fHistHPointingAngle(0),
+  fHistMassH(0),
+  fHistMassLambdaFromH(0),
+  fHistMassLambdaFromHtLorentz(0),
+  fHistMassPpi(0),
+  fHistMassPpiReso(0),
+  fHistMassLpi(0),
+  fHistMassLP(0),
+  fHistProtonPIDBb(0),
+  fHistPionPIDBb(0),
+  fHistProtonPIDLambda(0),
+  fHistPionPIDLambda(0),
+  fHistMCdcaPvtxDvtx(0),
+  fHistMCdcaPvtxLvtx(0),
+  fHistMCdcaDvtxLvtx(0),
+  fHistMCangleLH(0),
+  fHistMCdecayAngle(0),
+  fHistMCpointingAngle(0),
+  fHistMCap(0),
+  fHistMCdcaPvtxDvtxReso(0),
+  fHistMCdcaPvtxLvtxReso(0),
+  fHistMCdcaDvtxLvtxReso(0),
+  fHistMCangleLHReso(0),
+  fHistMCdecayAngleReso(0),
+  fHistMCpointingAngleReso(0),
+  fHistMCapReso(0),
+  fHistCentrality(0),
+  fHistCentralityAC(0), 
+  fHistMultiplicity(0), 
+  fHistHilf1(0),
+  fHistHilf2(0), 
+  fHistHilf3(0),
+  fHistHilf4(0),
+  fHistHilf5(0), 
+  fHistHilf6(0),
+  fHistPtvsEtaGen(0),
+  fHistPtvsEtaAso(0),
+  fHistPtvsYGen(0),
+  fHistPtvsYAso(0), 
+  fHistRap(0),
+  fHistCount(0),
+  fPIDtpcESD(0),
+  fHistTriggerStat(0),
+  fHistTriggerStatAfterEventSelection(0),
+  fHistMassHcentMult(0),
+  fHistNdim(0)
+
+{
+  // Constructor
+
+  // Define input and output slots here
+  // Input from a TChain
+  DefineInput(0, TChain::Class());
+  // Output to TList container
+  DefineOutput(1, TList::Class()); //full
+
+  //MC info contol
+  if (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())SetHasMC();
+
+  //V0 cuts
+
+  fESDtrackCutsV0 = new AliESDtrackCuts("AliESDtrackCutsV0","AliESDtrackCutsV0");
+  fESDtrackCutsV0->SetAcceptKinkDaughters(kFALSE);
+  fESDtrackCutsV0->SetMinNClustersTPC(80);
+  fESDtrackCutsV0->SetMaxChi2PerClusterTPC(5);
+  fESDtrackCutsV0->SetRequireTPCRefit(kTRUE);
   fESDtrackCutsV0->SetEtaRange(-0.9,0.9);
-  fESDtrackCutsV0->SetPtRange(0.2,1.5);\r
-  fESDtrackCutsV0->SetMinDCAToVertexXY(2);\r //war inzwischen 1 & 3
-  fESDtrackCutsV0->SetMinDCAToVertexZ(2);\r //war inzwischen 1 & 3
-\r
-  fESDCutsV0 = new AliESDv0Cuts("AliESDCutsV0","AliESDCutsV0");\r
+  fESDtrackCutsV0->SetPtRange(0.2,1.5);
+  fESDtrackCutsV0->SetMinDCAToVertexXY(2); //war inzwischen 1 & 3
+  fESDtrackCutsV0->SetMinDCAToVertexZ(2); //war inzwischen 1 & 3
+
+  fESDCutsV0 = new AliESDv0Cuts("AliESDCutsV0","AliESDCutsV0");
   fESDCutsV0->SetMaxDcaV0Daughters(1.0);
-  fESDCutsV0->SetMinDcaNegToVertex(2);\r //1.5
+  fESDCutsV0->SetMinDcaNegToVertex(2); //1.5
     fESDCutsV0->SetMinDcaPosToVertex(2); //1.5
-\r
-  //ESD Track cuts\r
-  fEsdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts");    \r
-  fEsdTrackCuts->SetMinNClustersTPC(80);\r
-  fEsdTrackCuts->SetAcceptKinkDaughters(kFALSE);\r
-  fEsdTrackCuts->SetMaxChi2PerClusterTPC(5);\r
-  fEsdTrackCuts->SetRequireTPCRefit(kTRUE);\r
-  fEsdTrackCuts->SetEtaRange(-0.9,0.9);\r
-}\r
-\r
-//____________________________________________________________\r
-AliAnalysisTaskHdibaryonLPpi::~AliAnalysisTaskHdibaryonLPpi(){\r
-  //\r
-  // Destructor\r
-  //\r
-  if(fHistList){ \r
-    fHistList->Clear();\r
-    delete fHistList;\r
-  }\r
-  if(fEsdTrackCuts) delete fEsdTrackCuts;\r
-  if(fESDtrackCutsV0) delete fESDtrackCutsV0;\r
-  if(fESDCutsV0) delete fESDCutsV0;\r
-\r
-}\r
-\r
-//________________________________________________________________________\r
-void AliAnalysisTaskHdibaryonLPpi::UserCreateOutputObjects()\r
-{\r
-  // Create histograms\r
-  // Called once\r
-\r
- fHistList = new TList();\r
- fHistList->SetOwner();\r
-\r
- fHistMassDPi = new TH1F("fHistMassDPi", "Invariant mass distribution p+#pi^{-} ", 500, 1.0, 1.25);\r
- fHistMassDPi->GetXaxis()->SetTitle("Invariant mass p+#pi^{-} (GeV/c^{2})");\r
- fHistMassDPi->GetYaxis()->SetTitle("Entries");\r
- fHistMassDPi->SetMarkerStyle(kFullCircle); \r
-\r
- fHistMassLPi = new TH1F("fHistMassLPi", "Offline Invariant mass distribution p+#pi^{-} ", 500, 1.0, 1.25);\r
- fHistMassLPi->GetXaxis()->SetTitle("Invariant mass p+#pi^{-} (GeV/c^{2})");\r
- fHistMassLPi->GetYaxis()->SetTitle("Entries");\r
- fHistMassLPi->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassLambdaPi = new TH1F("fHistMassLambdaPi", "Invariant mass distribution #Lambda+#pi^{-} ", 500, 1.2, 1.5);\r
- fHistMassLambdaPi->GetXaxis()->SetTitle("Invariant mass #Lambda+#pi^{-} (GeV/c^{2})");\r
- fHistMassLambdaPi->GetYaxis()->SetTitle("Entries");\r
- fHistMassLambdaPi->SetMarkerStyle(kFullCircle);\r
\r
- fHistMassLambda = new TH1F("fHistMassLambda", "Invariant mass distribution of #Lambda for further analyis", 500, 1.0, 1.2);\r
- fHistMassLambda->GetXaxis()->SetTitle("Invariant mass p+#pi^{+} (GeV/c^{2})");\r
- fHistMassLambda->GetYaxis()->SetTitle("Entries");\r
- fHistMassLambda->SetMarkerStyle(kFullCircle);\r
\r
- fHistMassLambdaPPi = new TH1F("fHistMassLambdaPPi", "Invariant mass distribution #Lambdap#pi^{-} ", 300, 2.2, 2.5);\r
- fHistMassLambdaPPi->GetXaxis()->SetTitle("Invariant mass #Lambdap#pi^{-} (GeV/c^{2})");\r
- fHistMassLambdaPPi->GetYaxis()->SetTitle("Entries");\r
- fHistMassLambdaPPi->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassLambdaP = new TH1F("fHistMassLambdaP", "Invariant mass distribution #Lambdap ", 300, 2.2, 2.5);\r
- fHistMassLambdaP->GetXaxis()->SetTitle("Invariant mass #Lambdap (GeV/c^{2})");\r
- fHistMassLambdaP->GetYaxis()->SetTitle("Entries");\r
- fHistMassLambdaP->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassLambdaK = new TH1F("fHistMassLambdaK", "Invariant mass distribution #LambdaK^{-} ", 300, 1.6, 1.9);\r
- fHistMassLambdaK->GetXaxis()->SetTitle("Invariant mass #LambdaK^{-} (GeV/c^{2})");\r
- fHistMassLambdaK->GetYaxis()->SetTitle("Entries");\r
- fHistMassLambdaK->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassK0onl = new TH1F("fHistMassK0onl", "Invariant mass distribution K_{s}^{0} online V0 finder ", 400, 0.2, 1.0);\r
- fHistMassK0onl->GetXaxis()->SetTitle("Invariant mass #pi^{+}#pi^{-} (GeV/c^{2})");\r
- fHistMassK0onl->GetYaxis()->SetTitle("Entries");\r
- fHistMassK0onl->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassK0offl = new TH1F("fHistMassK0offl", "Invariant mass distribution  K_{s}^{0} offline V0 finder ", 400, 0.2, 1.0);\r
- fHistMassK0offl->GetXaxis()->SetTitle("Invariant mass  #pi^{+}#pi^{-} (GeV/c^{2})");\r
- fHistMassK0offl->GetYaxis()->SetTitle("Entries");\r
- fHistMassK0offl->SetMarkerStyle(kFullCircle);\r
\r
- fHistMassK0onlC = new TH1F("fHistMassK0onlC", "Invariant mass distribution K_{s}^{0} online V0 finder ", 400, 0.2, 1.0);\r
- fHistMassK0onlC->GetXaxis()->SetTitle("Invariant mass #pi^{+}#pi^{-} (GeV/c^{2})");\r
- fHistMassK0onlC->GetYaxis()->SetTitle("Entries");\r
- fHistMassK0onlC->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassK0offlC = new TH1F("fHistMassK0offlC", "Invariant mass distribution  K_{s}^{0} offline V0 finder ", 400, 0.2, 1.0);\r
- fHistMassK0offlC->GetXaxis()->SetTitle("Invariant mass  #pi^{+}#pi^{-} (GeV/c^{2})");\r
- fHistMassK0offlC->GetYaxis()->SetTitle("Entries");\r
- fHistMassK0offlC->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassPQonl = new TH1F("fHistMassPQonl", "Invariant mass distribution K_{s}^{0}p using online V0 finder ", 500, 1.3, 2.3);\r
- fHistMassPQonl->GetXaxis()->SetTitle("Invariant mass K_{s}^{0}p (GeV/c^{2})");\r
- fHistMassPQonl->GetYaxis()->SetTitle("Entries");\r
- fHistMassPQonl->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassPQoffl = new TH1F("fHistMassPQoffl", "Invariant mass distribution  K_{s}^{0}p using offline V0 finder ", 500, 1.3, 2.3);\r
- fHistMassPQoffl->GetXaxis()->SetTitle("Invariant mass K_{s}^{0}p (GeV/c^{2})");\r
- fHistMassPQoffl->GetYaxis()->SetTitle("Entries");\r
- fHistMassPQoffl->SetMarkerStyle(kFullCircle);\r
-\r
- fHistDC = new TH1F("fHistDC", "Proper decay length", 500, 0.0, 25);\r
- fHistDC->GetXaxis()->SetTitle("c#tau (cm)");\r
- fHistDC->GetYaxis()->SetTitle("Entries");\r
- fHistDC->SetMarkerStyle(kFullCircle);\r
-  \r
- fHistArmenterosPodolanski = new TH2F("fHistArmenterosPodolanski", "Armenteros-Podolanski plot", 200,-1.0,1.0, 500,0,1);\r
- fHistArmenterosPodolanski->GetXaxis()->SetTitle("#alpha");\r
- fHistArmenterosPodolanski->GetYaxis()->SetTitle("q_{t}");\r
- fHistArmenterosPodolanski->SetMarkerStyle(kFullCircle);\r
-\r
- fHistArmenterosPodolanskiCut = new TH2F("fHistArmenterosPodolanskiCut", "Armenteros-Podolanski plot after cut", 200,-1.0,1.0, 500,0,1);\r
- fHistArmenterosPodolanskiCut->GetXaxis()->SetTitle("#alpha");\r
- fHistArmenterosPodolanskiCut->GetYaxis()->SetTitle("q_{t}");\r
- fHistArmenterosPodolanskiCut->SetMarkerStyle(kFullCircle);\r
-\r
- fHistHDibaryonInvaMassGen = new TH1F("fHistHDibaryonInvaMassGen", "Generated  #Lambda p #pi^{-}", 200, 2.1, 2.3);\r
- fHistHDibaryonInvaMassGen->GetYaxis()->SetTitle("Counts");\r
- fHistHDibaryonInvaMassGen->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");\r
-\r
- fHistHDibaryonInvaMassGenRes = new TH1F("fHistHDibaryonInvaMassGenRes", "Generated  #Lambda p #pi^{-} with particles in rapidity!!!", 200, 2.1, 2.3);\r
- fHistHDibaryonInvaMassGenRes->GetYaxis()->SetTitle("Counts");\r
- fHistHDibaryonInvaMassGenRes->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");\r
-\r
- fHistAntiHDibaryonInvaMassGen = new TH1F("fHistAntiHDibaryonInvaMassGen", "Generated  #bar{#Lambda} #bar{p} #pi^{+}", 200, 2.1, 2.3);\r
- fHistAntiHDibaryonInvaMassGen->GetYaxis()->SetTitle("Counts");\r
- fHistAntiHDibaryonInvaMassGen->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");\r
-\r
- fHistHDibaryonInvaMassAso = new TH1F("fHistHDibaryonInvaMassAso", "Associated  #Lambda p #pi^{-}", 200, 2.1, 2.3);\r
- fHistHDibaryonInvaMassAso->GetYaxis()->SetTitle("Counts");\r
- fHistHDibaryonInvaMassAso->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");\r
-\r
- fHistHDibaryonInvaMassAsoReso = new TH1F("fHistHDibaryonInvaMassAsoReso", "Associated  #Lambda p #pi^{-}", 200, 2.1, 2.3);\r
- fHistHDibaryonInvaMassAsoReso->GetYaxis()->SetTitle("Counts");\r
- fHistHDibaryonInvaMassAsoReso->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");\r
-\r
- fHistAntiHDibaryonInvaMassAso = new TH1F("fHistAntiHDibaryonInvaMassAso", "Associated  #bar{#Lambda} #bar{p} #pi^{+}", 200, 2.1, 2.3);\r
- fHistAntiHDibaryonInvaMassAso->GetYaxis()->SetTitle("Counts");\r
- fHistAntiHDibaryonInvaMassAso->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");\r
-\r
- fHistCheck = new TH2F("fHistCheck", "Check online/offline", 200, -0.5, 1.5, 200, -0.5, 1.5);\r
- fHistCheck->GetXaxis()->SetTitle("offline");\r
- fHistCheck->GetYaxis()->SetTitle("online");\r
- fHistCheck->SetMarkerStyle(kFullCircle);\r
-\r
- fHistHPointingAngle= new TH1F("fHistHPointingAngle", "Pointing angle distribution for #Lambdap#pi^{-}", 200, 0., 2*TMath::Pi());\r
- fHistHPointingAngle->GetXaxis()->SetTitle("Pointing angle distribution for #Lambdap#pi^{-}");\r
- fHistHPointingAngle->GetYaxis()->SetTitle("Entries");\r
- fHistHPointingAngle->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassH= new TH1F("fHistMassH", "Invariant mass distribution #Lambdap#pi^{-} (GeV/c^{2}) after pointing angle cut", 3000, 2.2, 2.5);\r
- fHistMassH->GetXaxis()->SetTitle("Invariant mass #Lambdap#pi^{-} (GeV/c^{2})");\r
- fHistMassH->GetYaxis()->SetTitle("Entries");\r
- fHistMassH->SetMarkerStyle(kFullCircle);\r
\r
- fHistMassLambdaFromH= new TH1F("fHistMassLambdaFromH", "Invariant mass distribution #Lambda (GeV/c^{2}) asking for Mother to be H", 300, 1.0, 1.3);\r
- fHistMassLambdaFromH->GetXaxis()->SetTitle("Invariant mass #Lambda (GeV/c^{2})");\r
- fHistMassLambdaFromH->GetYaxis()->SetTitle("Entries");\r
- fHistMassLambdaFromH->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassLambdaFromHtLorentz= new TH1F("fHistMassLambdaFromHtLorentz", "Invariant mass distribution #Lambda (GeV/c^{2}) asking for Mother to be H", 300, 1.0, 1.3);\r
- fHistMassLambdaFromHtLorentz->GetXaxis()->SetTitle("Invariant mass #Lambda (GeV/c^{2})");\r
- fHistMassLambdaFromHtLorentz->GetYaxis()->SetTitle("Entries");\r
- fHistMassLambdaFromHtLorentz->SetMarkerStyle(kFullCircle);\r
\r
- fHistMassPpi= new TH1F("fHistMassPpi", "Invariant mass distribution of the p#pi^{-} used for combing with Lambda", 300, 1.0, 1.6);\r
- fHistMassPpi->GetXaxis()->SetTitle("Invariant mass p#pi^{-} (GeV/c^{2})");\r
- fHistMassPpi->GetYaxis()->SetTitle("Entries");\r
- fHistMassPpi->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassPpiReso= new TH1F("fHistMassPpiReso", "Invariant mass distribution of the p#pi^{-} used for combing with Lambda", 300, 1.0, 1.6);\r
- fHistMassPpiReso->GetXaxis()->SetTitle("Invariant mass p#pi^{-} (GeV/c^{2})");\r
- fHistMassPpiReso->GetYaxis()->SetTitle("Entries");\r
- fHistMassPpiReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassLpi= new TH1F("fHistMassLpi", "Invariant mass distribution of the #Lambda#pi^{-} used for combing with p", 300, 1.1, 1.7);\r
- fHistMassLpi->GetXaxis()->SetTitle("Invariant mass #Lambda#pi^{-} (GeV/c^{2})");\r
- fHistMassLpi->GetYaxis()->SetTitle("Entries");\r
- fHistMassLpi->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMassLP= new TH1F("fHistMassLP", "Invariant mass distribution of the #Lambda p used for combing with #pi^{-}", 300, 2.0, 2.3);\r
- fHistMassLP->GetXaxis()->SetTitle("Invariant mass #Lambda p (GeV/c^{2})");\r
- fHistMassLP->GetYaxis()->SetTitle("Entries");\r
- fHistMassLP->SetMarkerStyle(kFullCircle);\r
-\r
- fHistProtonPIDBb = new TH2F("fHistProtonPIDBb", "dE/dx after p PID", 100, 0., 10, 100, 0, 100);\r
- fHistProtonPIDBb->GetYaxis()->SetTitle("TPC Signal");\r
- fHistProtonPIDBb->GetXaxis()->SetTitle("P (GeV/c)");\r
- fHistProtonPIDBb->SetOption("scat");\r
- fHistProtonPIDBb->SetMarkerStyle(kFullCircle);\r
-\r
- fHistPionPIDBb = new TH2F("fHistPionPIDBb", "dE/dx after K PID", 100, 0., 10, 100, 0, 100);\r
- fHistPionPIDBb->GetYaxis()->SetTitle("TPC Signal");\r
- fHistPionPIDBb->GetXaxis()->SetTitle("P (GeV/c)");\r
- fHistPionPIDBb->SetOption("scat");\r
- fHistPionPIDBb->SetMarkerStyle(kFullCircle);\r
-\r
- fHistProtonPIDLambda = new TH2F("fHistProtonPIDLambda", "dE/dx after p PID from V0", 100, 0., 10, 100, 0, 100);\r
- fHistProtonPIDLambda->GetYaxis()->SetTitle("TPC Signal");\r
- fHistProtonPIDLambda->GetXaxis()->SetTitle("P (GeV/c)");\r
- fHistProtonPIDLambda->SetOption("scat");\r
- fHistProtonPIDLambda->SetMarkerStyle(kFullCircle);\r
-\r
- fHistPionPIDLambda = new TH2F("fHistPionPIDLambda", "dE/dx after #pi PID from V0", 100, 0, 10, 100, 0, 100);\r
- fHistPionPIDLambda->GetYaxis()->SetTitle("TPC Signal");\r
- fHistPionPIDLambda->GetXaxis()->SetTitle("P (GeV/c)");\r
- fHistPionPIDLambda->SetOption("scat");\r
- fHistPionPIDLambda->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdcaPvtxDvtx= new TH1F("fHistMCdcaPvtxDvtx", "MC True DCA Primary Vertex - H Decay Vertex", 300, -0.1, 11.9);\r
- fHistMCdcaPvtxDvtx->GetXaxis()->SetTitle("dca prim. vtx- decay vtx (cm)");\r
- fHistMCdcaPvtxDvtx->GetYaxis()->SetTitle("Entries");\r
- fHistMCdcaPvtxDvtx->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdcaPvtxLvtx= new TH1F("fHistMCdcaPvtxLvtx", "MC True DCA Primary Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);\r
- fHistMCdcaPvtxLvtx->GetXaxis()->SetTitle("dca prim. vtx-#Lambda decay vtx (cm)");\r
- fHistMCdcaPvtxLvtx->GetYaxis()->SetTitle("Entries");\r
- fHistMCdcaPvtxLvtx->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdcaDvtxLvtx= new TH1F("fHistMCdcaDvtxLvtx", "MC True DCA H Decay Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);\r
- fHistMCdcaDvtxLvtx->GetXaxis()->SetTitle("dca H deacy vtx-#Lambda decay vtx (cm)");\r
- fHistMCdcaDvtxLvtx->GetYaxis()->SetTitle("Entries");\r
- fHistMCdcaDvtxLvtx->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCangleLH= new TH1F("fHistMCangleLH", "MC True Angle between #Lambda and H", 300, 0., 2*TMath::Pi());\r
- fHistMCangleLH->GetXaxis()->SetTitle("Angle (#Lambda-H)");\r
- fHistMCangleLH->GetYaxis()->SetTitle("Entries");\r
- fHistMCangleLH->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdecayAngle= new TH1F("fHistMCdecayAngle", "MC True Angle between decay products", 300, 0., 2*TMath::Pi());\r
- fHistMCdecayAngle->GetXaxis()->SetTitle("Angle (#Lambda-p#pi)");\r
- fHistMCdecayAngle->GetYaxis()->SetTitle("Entries");\r
- fHistMCdecayAngle->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCpointingAngle= new TH1F("fHistMCpointingAngle", "MC True Pointing Angle", 3000, 0., 2*TMath::Pi());\r
- fHistMCpointingAngle->GetXaxis()->SetTitle("Pointing Angle");\r
- fHistMCpointingAngle->GetYaxis()->SetTitle("Entries");\r
- fHistMCpointingAngle->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCap = new TH2F("fHistMCap", "True MC Armenteros-Podolanski", 200,-1.0,1.0, 500,0,1);\r
- fHistMCap->GetYaxis()->SetTitle("#alpha");\r
- fHistMCap->GetXaxis()->SetTitle("q_{t}");\r
- fHistMCap->SetOption("scat");\r
- fHistMCap->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdcaPvtxDvtxReso= new TH1F("fHistMCdcaPvtxDvtxReso", "MC True DCA Primary Vertex - H Decay Vertex", 300, -0.1, 11.9);\r
- fHistMCdcaPvtxDvtxReso->GetXaxis()->SetTitle("dca prim. vtx- decay vtx (cm)");\r
- fHistMCdcaPvtxDvtxReso->GetYaxis()->SetTitle("Entries");\r
- fHistMCdcaPvtxDvtxReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdcaPvtxLvtxReso= new TH1F("fHistMCdcaPvtxLvtxReso", "MC True DCA Primary Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);\r
- fHistMCdcaPvtxLvtxReso->GetXaxis()->SetTitle("dca prim. vtx-#Lambda decay vtx (cm)");\r
- fHistMCdcaPvtxLvtxReso->GetYaxis()->SetTitle("Entries");\r
- fHistMCdcaPvtxLvtxReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdcaDvtxLvtxReso= new TH1F("fHistMCdcaDvtxLvtxReso", "MC True DCA H Decay Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);\r
- fHistMCdcaDvtxLvtxReso->GetXaxis()->SetTitle("dca H deacy vtx-#Lambda decay vtx (cm)");\r
- fHistMCdcaDvtxLvtxReso->GetYaxis()->SetTitle("Entries");\r
- fHistMCdcaDvtxLvtxReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCangleLHReso= new TH1F("fHistMCangleLHReso", "MC True Angle between #Lambda and H", 300, 0., 2*TMath::Pi());\r
- fHistMCangleLHReso->GetXaxis()->SetTitle("Angle (#Lambda-H)");\r
- fHistMCangleLHReso->GetYaxis()->SetTitle("Entries");\r
- fHistMCangleLHReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCdecayAngleReso= new TH1F("fHistMCdecayAngleReso", "MC True Angle between decay products", 300, 0., 2*TMath::Pi());\r
- fHistMCdecayAngleReso->GetXaxis()->SetTitle("Angle (#Lambda-p#pi)");\r
- fHistMCdecayAngleReso->GetYaxis()->SetTitle("Entries");\r
- fHistMCdecayAngleReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCpointingAngleReso= new TH1F("fHistMCpointingAngleReso", "MC True Pointing Angle", 300, 0., 2*TMath::Pi());\r
- fHistMCpointingAngleReso->GetXaxis()->SetTitle("Pointing Angle");\r
- fHistMCpointingAngleReso->GetYaxis()->SetTitle("Entries");\r
- fHistMCpointingAngleReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistMCapReso = new TH2F("fHistMCapReso", "True MC Armenteros-Podolanski", 200,-1.0,1.0, 500,0,1);\r
- fHistMCapReso->GetYaxis()->SetTitle("#alpha");\r
- fHistMCapReso->GetXaxis()->SetTitle("q_{t}");\r
- fHistMCapReso->SetOption("scat");\r
- fHistMCapReso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistCentrality = new TH1F("Centrality ", "Centrality", 11, -0.5, 10.5);\r
- fHistCentrality ->GetXaxis()->SetTitle("Centrality");\r
- fHistCentrality ->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistCentralityAC = new TH1F("CentralityAC ", "CentralityAC", 11, -0.5, 10.5);\r
- fHistCentralityAC ->GetXaxis()->SetTitle("Centrality");\r
- fHistCentralityAC ->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistMultiplicity = new TH1F("Multiplicity ", "Multiplicity", 100, 0, 20000);\r
- fHistMultiplicity ->GetXaxis()->SetTitle("Centrality");\r
- fHistMultiplicity ->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistList->Add(fHistMassDPi);\r
- fHistList->Add(fHistMassLPi);\r
- fHistList->Add(fHistMassLambdaPi);\r
- fHistList->Add(fHistMassLambda);\r
- fHistList->Add(fHistMassLambdaPPi);\r
- fHistList->Add(fHistMassLambdaP);\r
- fHistList->Add(fHistMassLambdaK);\r
- fHistList->Add(fHistMassK0onl);\r
- fHistList->Add(fHistMassK0offl);\r
- fHistList->Add(fHistMassK0onlC);\r
- fHistList->Add(fHistMassK0offlC);\r
- fHistList->Add(fHistMassPQonl); \r
- fHistList->Add(fHistMassPQoffl);\r
- fHistList->Add(fHistDC); \r
- fHistList->Add(fHistArmenterosPodolanski);\r
- fHistList->Add(fHistArmenterosPodolanskiCut);\r
- fHistList->Add(fHistHDibaryonInvaMassGen);\r
- fHistList->Add(fHistHDibaryonInvaMassGenRes);\r
- fHistList->Add(fHistAntiHDibaryonInvaMassGen);\r
- fHistList->Add(fHistHDibaryonInvaMassAso);\r
- fHistList->Add(fHistHDibaryonInvaMassAsoReso);\r
- fHistList->Add(fHistAntiHDibaryonInvaMassAso);\r
- fHistList->Add(fHistCheck);\r
- fHistList->Add(fHistHPointingAngle);\r
- fHistList->Add(fHistMassH);\r
- fHistList->Add(fHistMassPpi);\r
- fHistList->Add(fHistMassPpiReso);\r
- fHistList->Add(fHistMassLpi);\r
- fHistList->Add(fHistMassLP);\r
- fHistList->Add(fHistMassLambdaFromH);\r
- fHistList->Add(fHistMassLambdaFromHtLorentz);\r
- fHistList->Add(fHistProtonPIDBb);\r
- fHistList->Add(fHistPionPIDBb);\r
- fHistList->Add(fHistProtonPIDLambda);\r
- fHistList->Add(fHistPionPIDLambda);\r
- fHistList->Add(fHistMCdcaPvtxDvtx);\r
- fHistList->Add(fHistMCdcaPvtxLvtx);\r
- fHistList->Add(fHistMCdcaDvtxLvtx);\r
- fHistList->Add(fHistMCangleLH);\r
- fHistList->Add(fHistMCdecayAngle);\r
- fHistList->Add(fHistMCpointingAngle);\r
- fHistList->Add(fHistMCap);\r
- fHistList->Add(fHistMCdcaPvtxDvtxReso);\r
- fHistList->Add(fHistMCdcaPvtxLvtxReso);\r
- fHistList->Add(fHistMCdcaDvtxLvtxReso);\r
- fHistList->Add(fHistMCangleLHReso);\r
- fHistList->Add(fHistMCdecayAngleReso);\r
- fHistList->Add(fHistMCpointingAngleReso);\r
- fHistList->Add(fHistMCapReso);\r
- fHistList->Add(fHistCentrality);\r
- fHistList->Add(fHistCentralityAC);\r
- fHistList->Add(fHistMultiplicity);\r
\r
- fHistHilf1= new TH1F("fHistHilf1", "HD", 300, 0., 10);\r
- fHistHilf1->GetXaxis()->SetTitle("");\r
- fHistHilf1->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistHilf2= new TH1F("fHistHilf2", "HD", 300, 0., 10);\r
- fHistHilf2->GetXaxis()->SetTitle("");\r
- fHistHilf2->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistHilf3= new TH1F("fHistHilf3", "HD", 300, 0., 10);\r
- fHistHilf3->GetXaxis()->SetTitle("");\r
- fHistHilf3->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistHilf4= new TH1F("fHistHilf4", "HD", 300, 0., 10);\r
- fHistHilf4->GetXaxis()->SetTitle("");\r
- fHistHilf4->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistHilf5= new TH1F("fHistHilf5", "HD", 300, 0., 10);\r
- fHistHilf5->GetXaxis()->SetTitle("");\r
- fHistHilf5->GetYaxis()->SetTitle("Entries");\r
-\r
- fHistHilf6= new TH1F("fHistHilf6", "HD", 300, 0., 10);\r
- fHistHilf6->GetXaxis()->SetTitle("");\r
- fHistHilf6->GetYaxis()->SetTitle("Entries");\r
\r
- fHistPtvsEtaGen = new TH2F("fHistPtvsEtaGen", "p_{t} vs #eta from generated H", 200,0.0,10.0, 200,-1,1);\r
- fHistPtvsEtaGen->GetXaxis()->SetTitle("p_{t} (GeV/c)");\r
- fHistPtvsEtaGen->GetYaxis()->SetTitle("#eta");\r
- fHistPtvsEtaGen->SetOption("scat");\r
- fHistPtvsEtaGen->SetMarkerStyle(kFullCircle);\r
-\r
- fHistPtvsEtaAso = new TH2F("fHistPtvsEtaAso", "p_{t} vs #eta from associated H", 200,0.0,10.0, 200,-1,1);\r
- fHistPtvsEtaAso->GetYaxis()->SetTitle("p_{t} (GeV/c)");\r
- fHistPtvsEtaAso->GetXaxis()->SetTitle("#eta");\r
- fHistPtvsEtaAso->SetOption("scat");\r
- fHistPtvsEtaAso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistPtvsYGen = new TH2F("fHistPtvsYGen", "p_{t} vs rapidity from generated H", 200,0.0,10.0, 200,-1,1);\r
- fHistPtvsYGen->GetXaxis()->SetTitle("p_{t} (GeV/c)");\r
- fHistPtvsYGen->GetYaxis()->SetTitle("y");\r
- fHistPtvsYGen->SetOption("scat");\r
- fHistPtvsYGen->SetMarkerStyle(kFullCircle);\r
-\r
- fHistPtvsYAso = new TH2F("fHistPtvsYAso", "p_{t} vs rapidity from associated H", 200,0.0,10.0, 200,-1,1);\r
- fHistPtvsYAso->GetXaxis()->SetTitle("p_{t} (GeV/c)");\r
- fHistPtvsYAso->GetYaxis()->SetTitle("y");\r
- fHistPtvsYAso->SetOption("scat");\r
- fHistPtvsYAso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistRap= new TH1F("fHistRap", "Rapidity", 400, -2., 2);\r
- fHistRap->GetXaxis()->SetTitle("Y");\r
- fHistRap->GetYaxis()->SetTitle("Entries");\r
- fHistPtvsEtaAso->SetMarkerStyle(kFullCircle);\r
-\r
- fHistList->Add(fHistHilf1);\r
- fHistList->Add(fHistHilf2); \r
- fHistList->Add(fHistHilf3);\r
- fHistList->Add(fHistHilf4);\r
- fHistList->Add(fHistHilf5); \r
- fHistList->Add(fHistHilf6);\r
- fHistList->Add(fHistPtvsEtaGen);\r
- fHistList->Add(fHistPtvsEtaAso);\r
- fHistList->Add(fHistPtvsYGen);\r
- fHistList->Add(fHistPtvsYAso);\r
- fHistList->Add(fHistRap);\r
-\r
- fHistCount = new TH1F("fHistCount","test",17,0,17);\r
- fHistCount->GetXaxis()->SetBinLabel(1,"Events");\r
- fHistCount->GetXaxis()->SetBinLabel(2,"MC All");\r
- fHistCount->GetXaxis()->SetBinLabel(3,"MC from Primary Vtx");\r
- fHistCount->GetXaxis()->SetBinLabel(4,"Horst");\r
- fHistCount->GetXaxis()->SetBinLabel(5,"Lambda Candidates");\r
- fHistCount->GetXaxis()->SetBinLabel(6,"Sigma Candidates");\r
- fHistCount->GetXaxis()->SetBinLabel(7,"Horst");\r
- fHistCount->GetXaxis()->SetBinLabel(8,"Horst");\r
- fHistCount->GetXaxis()->SetBinLabel(9,"Horst");\r
- fHistCount->GetXaxis()->SetBinLabel(10,"MC All #bar{Lambda}(1520)s");\r
- fHistCount->GetXaxis()->SetBinLabel(11,"");\r
- fHistCount->GetXaxis()->SetBinLabel(12,"H-Dibaryon");\r
- fHistCount->GetXaxis()->SetBinLabel(13,"Hypertriton 2-Body");\r
- fHistCount->GetXaxis()->SetBinLabel(14,"Hypertriton 3-Body");\r
- fHistCount->GetXaxis()->SetBinLabel(15,"");\r
- fHistCount->GetXaxis()->SetBinLabel(16,"");\r
- fHistCount->GetXaxis()->SetBinLabel(17,"Lambdas!!!");\r
- fHistCount->SetStats(0);\r
- fHistCount->SetFillColor(38);\r
- fHistList->Add(fHistCount);\r
-\r
- //trigger statistics histogram\r
-  fHistTriggerStat = new TH1F("fHistTriggerStat","Trigger statistics", 4,-0.5, 3.5);\r
-  const Char_t* aTriggerNames[] = { "kMB", "kCentral", "kSemiCentral" };\r
-  for ( Int_t ii=0; ii < 3; ii++ )\r
-    fHistTriggerStat->GetXaxis()->SetBinLabel(ii+1, aTriggerNames[ii]);\r
-\r
-  fHistTriggerStatAfterEventSelection = new TH1F("fHistTriggerStatAfterEventSelection","Trigger statistics after event selection", 4,-0.5, 3.5);\r
-  for ( Int_t ii=0; ii < 3; ii++ )\r
-    fHistTriggerStatAfterEventSelection->GetXaxis()->SetBinLabel(ii+1, aTriggerNames[ii]);\r
-  fHistList->Add(fHistTriggerStat);\r
-  fHistList->Add(fHistTriggerStatAfterEventSelection);\r
-\r
-  fHistMassHcentMult =  new TH3F("fHistMassHcentMult", "Inv. Mass vs. centrality vs. multiplicity", 100, 2.2, 2.3, 5, 0, 4, 300, 0, 6000);\r
-  fHistMassHcentMult->GetXaxis()->SetTitle("Invariant mass #Lambdap#pi^{-} (GeV/c^{2})"); // inv. mass\r
-  fHistMassHcentMult->GetYaxis()->SetTitle("Centrality"); // triggertype\r
-  fHistMassHcentMult->GetZaxis()->SetTitle("Multiplicity"); // refTPC\r
-  fHistList->Add(fHistMassHcentMult);\r
-\r
-  \r
-  const Double_t kz = 2*TMath::Pi();\r
-  Int_t binsD01[16]={  300, 200, 100, 100, 100, 100, 100, 100, 200, 200, 200, 200, 400, 200, 200, 3};\r
-  Double_t xminD01[16]={2.0, 1.0, 0., -1, 0., 0., 0., 0., 0., 0., 0., 0., -1, 0.,  0., 0};\r
-  Double_t xmaxD01[16]={2.3, 1.2, kz, 1, 1, 10, 10, 5, 5, 5, 5, 100, 1, 100,  4000, 1};\r
-  \r
-\r
-  fHistNdim = new THnSparseF("fHistNdim","THnS;InvMass, InvMassLambda, pointingAngle, armPoAlpha, armPoQt, pTL, pTH, d0p, d0n, dcaHd, dca, decayL, cosPA, centr, multi, mcinf;InvMassH", 16,binsD01,xminD01,xmaxD01);\r
-  fHistList->Add(fHistNdim);\r
-\r
- AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();\r
- AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());\r
- fPIDtpcESD = inputHandler->GetPIDResponse();\r
-\r
-// Post output data (if histograms are not used later, PostData is at least called here)\r
-  PostData(1, fHistList);\r
-\r
-}\r
-\r
- //________________________________________________________________________\r
-void AliAnalysisTaskHdibaryonLPpi::UserExec(Option_t *) \r
-{ \r
-  // Main loop\r
-  // Called for each event\r
-\r
-  //define improtant masses\r
-  Double_t pionMass     = 0.13957;\r
-  Double_t protonMass   = 0.93827;\r
-\r
-  //define PDGCodes\r
-  Long_t pdgPionPlus          =         211;\r
-  Long_t pdgPionMinus         =        -211;\r
-  Long_t pdgProton            =        2212;\r
-  Long_t pdgAntiProton        =       -2212;\r
-  Long_t pdgLambda            =        3122;\r
-  Long_t pdgAntiLambda        =       -3122;\r
-  Long_t pdgHDibaryon         =  1020000020;\r
-  Long_t pdgAntiHDibaryon     = -1020000020;\r
-\r
-  AliStack* stack(NULL);\r
-  if(HasMC()){\r
-\r
-    if(!fMCEvent)return;\r
-\r
-    AliHeader *head = fMCEvent->Header();\r
-    if(!head)return;\r
-    AliGenPythiaEventHeader *pyheader = (AliGenPythiaEventHeader*)head->GenEventHeader();\r
-    if(!pyheader)return;\r
-    \r
-    if(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()){\r
-      if(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->MCEvent()){\r
-       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->TreeK()))return;\r
-       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->TreeTR()))return;\r
-       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->Init("local")))return;\r
-       stack = static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->MCEvent()->Stack();\r
-       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->MCEvent()->Stack()->TreeK()))return;\r
-      }\r
-    }\r
-       \r
-    if(!stack)return;\r
-  }\r
-\r
-  // -------------------------------------------------------\r
-  // Loop for Inv. Mass via ESD tracks\r
-  // -------------------------------------------------------\r
-\r
-  fHistCount->Fill(0);\r
-  \r
-  fESD = dynamic_cast<AliESDEvent *>(fInputEvent);\r
-\r
+
+  //ESD Track cuts
+  fEsdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts");    
+  fEsdTrackCuts->SetMinNClustersTPC(80);
+  fEsdTrackCuts->SetAcceptKinkDaughters(kFALSE);
+  fEsdTrackCuts->SetMaxChi2PerClusterTPC(5);
+  fEsdTrackCuts->SetRequireTPCRefit(kTRUE);
+  fEsdTrackCuts->SetEtaRange(-0.9,0.9);
+}
+
+//____________________________________________________________
+AliAnalysisTaskHdibaryonLPpi::~AliAnalysisTaskHdibaryonLPpi(){
+  //
+  // Destructor
+  //
+  if(fHistList){ 
+    fHistList->Clear();
+    delete fHistList;
+  }
+  if(fEsdTrackCuts) delete fEsdTrackCuts;
+  if(fESDtrackCutsV0) delete fESDtrackCutsV0;
+  if(fESDCutsV0) delete fESDCutsV0;
+
+}
+
+//________________________________________________________________________
+void AliAnalysisTaskHdibaryonLPpi::UserCreateOutputObjects()
+{
+  // Create histograms
+  // Called once
+
+ fHistList = new TList();
+ fHistList->SetOwner();
+
+ fHistMassDPi = new TH1F("fHistMassDPi", "Invariant mass distribution p+#pi^{-} ", 500, 1.0, 1.25);
+ fHistMassDPi->GetXaxis()->SetTitle("Invariant mass p+#pi^{-} (GeV/c^{2})");
+ fHistMassDPi->GetYaxis()->SetTitle("Entries");
+ fHistMassDPi->SetMarkerStyle(kFullCircle); 
+
+ fHistMassLPi = new TH1F("fHistMassLPi", "Offline Invariant mass distribution p+#pi^{-} ", 500, 1.0, 1.25);
+ fHistMassLPi->GetXaxis()->SetTitle("Invariant mass p+#pi^{-} (GeV/c^{2})");
+ fHistMassLPi->GetYaxis()->SetTitle("Entries");
+ fHistMassLPi->SetMarkerStyle(kFullCircle);
+
+ fHistMassLambdaPi = new TH1F("fHistMassLambdaPi", "Invariant mass distribution #Lambda+#pi^{-} ", 500, 1.2, 1.5);
+ fHistMassLambdaPi->GetXaxis()->SetTitle("Invariant mass #Lambda+#pi^{-} (GeV/c^{2})");
+ fHistMassLambdaPi->GetYaxis()->SetTitle("Entries");
+ fHistMassLambdaPi->SetMarkerStyle(kFullCircle);
+ fHistMassLambda = new TH1F("fHistMassLambda", "Invariant mass distribution of #Lambda for further analyis", 500, 1.0, 1.2);
+ fHistMassLambda->GetXaxis()->SetTitle("Invariant mass p+#pi^{+} (GeV/c^{2})");
+ fHistMassLambda->GetYaxis()->SetTitle("Entries");
+ fHistMassLambda->SetMarkerStyle(kFullCircle);
+ fHistMassLambdaPPi = new TH1F("fHistMassLambdaPPi", "Invariant mass distribution #Lambdap#pi^{-} ", 300, 2.2, 2.5);
+ fHistMassLambdaPPi->GetXaxis()->SetTitle("Invariant mass #Lambdap#pi^{-} (GeV/c^{2})");
+ fHistMassLambdaPPi->GetYaxis()->SetTitle("Entries");
+ fHistMassLambdaPPi->SetMarkerStyle(kFullCircle);
+
+ fHistMassLambdaP = new TH1F("fHistMassLambdaP", "Invariant mass distribution #Lambdap ", 300, 2.2, 2.5);
+ fHistMassLambdaP->GetXaxis()->SetTitle("Invariant mass #Lambdap (GeV/c^{2})");
+ fHistMassLambdaP->GetYaxis()->SetTitle("Entries");
+ fHistMassLambdaP->SetMarkerStyle(kFullCircle);
+
+ fHistMassLambdaK = new TH1F("fHistMassLambdaK", "Invariant mass distribution #LambdaK^{-} ", 300, 1.6, 1.9);
+ fHistMassLambdaK->GetXaxis()->SetTitle("Invariant mass #LambdaK^{-} (GeV/c^{2})");
+ fHistMassLambdaK->GetYaxis()->SetTitle("Entries");
+ fHistMassLambdaK->SetMarkerStyle(kFullCircle);
+
+ fHistMassK0onl = new TH1F("fHistMassK0onl", "Invariant mass distribution K_{s}^{0} online V0 finder ", 400, 0.2, 1.0);
+ fHistMassK0onl->GetXaxis()->SetTitle("Invariant mass #pi^{+}#pi^{-} (GeV/c^{2})");
+ fHistMassK0onl->GetYaxis()->SetTitle("Entries");
+ fHistMassK0onl->SetMarkerStyle(kFullCircle);
+
+ fHistMassK0offl = new TH1F("fHistMassK0offl", "Invariant mass distribution  K_{s}^{0} offline V0 finder ", 400, 0.2, 1.0);
+ fHistMassK0offl->GetXaxis()->SetTitle("Invariant mass  #pi^{+}#pi^{-} (GeV/c^{2})");
+ fHistMassK0offl->GetYaxis()->SetTitle("Entries");
+ fHistMassK0offl->SetMarkerStyle(kFullCircle);
+ fHistMassK0onlC = new TH1F("fHistMassK0onlC", "Invariant mass distribution K_{s}^{0} online V0 finder ", 400, 0.2, 1.0);
+ fHistMassK0onlC->GetXaxis()->SetTitle("Invariant mass #pi^{+}#pi^{-} (GeV/c^{2})");
+ fHistMassK0onlC->GetYaxis()->SetTitle("Entries");
+ fHistMassK0onlC->SetMarkerStyle(kFullCircle);
+
+ fHistMassK0offlC = new TH1F("fHistMassK0offlC", "Invariant mass distribution  K_{s}^{0} offline V0 finder ", 400, 0.2, 1.0);
+ fHistMassK0offlC->GetXaxis()->SetTitle("Invariant mass  #pi^{+}#pi^{-} (GeV/c^{2})");
+ fHistMassK0offlC->GetYaxis()->SetTitle("Entries");
+ fHistMassK0offlC->SetMarkerStyle(kFullCircle);
+
+ fHistMassPQonl = new TH1F("fHistMassPQonl", "Invariant mass distribution K_{s}^{0}p using online V0 finder ", 500, 1.3, 2.3);
+ fHistMassPQonl->GetXaxis()->SetTitle("Invariant mass K_{s}^{0}p (GeV/c^{2})");
+ fHistMassPQonl->GetYaxis()->SetTitle("Entries");
+ fHistMassPQonl->SetMarkerStyle(kFullCircle);
+
+ fHistMassPQoffl = new TH1F("fHistMassPQoffl", "Invariant mass distribution  K_{s}^{0}p using offline V0 finder ", 500, 1.3, 2.3);
+ fHistMassPQoffl->GetXaxis()->SetTitle("Invariant mass K_{s}^{0}p (GeV/c^{2})");
+ fHistMassPQoffl->GetYaxis()->SetTitle("Entries");
+ fHistMassPQoffl->SetMarkerStyle(kFullCircle);
+
+ fHistDC = new TH1F("fHistDC", "Proper decay length", 500, 0.0, 25);
+ fHistDC->GetXaxis()->SetTitle("c#tau (cm)");
+ fHistDC->GetYaxis()->SetTitle("Entries");
+ fHistDC->SetMarkerStyle(kFullCircle);
+  
+ fHistArmenterosPodolanski = new TH2F("fHistArmenterosPodolanski", "Armenteros-Podolanski plot", 200,-1.0,1.0, 500,0,1);
+ fHistArmenterosPodolanski->GetXaxis()->SetTitle("#alpha");
+ fHistArmenterosPodolanski->GetYaxis()->SetTitle("q_{t}");
+ fHistArmenterosPodolanski->SetMarkerStyle(kFullCircle);
+
+ fHistArmenterosPodolanskiCut = new TH2F("fHistArmenterosPodolanskiCut", "Armenteros-Podolanski plot after cut", 200,-1.0,1.0, 500,0,1);
+ fHistArmenterosPodolanskiCut->GetXaxis()->SetTitle("#alpha");
+ fHistArmenterosPodolanskiCut->GetYaxis()->SetTitle("q_{t}");
+ fHistArmenterosPodolanskiCut->SetMarkerStyle(kFullCircle);
+
+ fHistHDibaryonInvaMassGen = new TH1F("fHistHDibaryonInvaMassGen", "Generated  #Lambda p #pi^{-}", 200, 2.1, 2.3);
+ fHistHDibaryonInvaMassGen->GetYaxis()->SetTitle("Counts");
+ fHistHDibaryonInvaMassGen->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");
+
+ fHistHDibaryonInvaMassGenRes = new TH1F("fHistHDibaryonInvaMassGenRes", "Generated  #Lambda p #pi^{-} with particles in rapidity!!!", 200, 2.1, 2.3);
+ fHistHDibaryonInvaMassGenRes->GetYaxis()->SetTitle("Counts");
+ fHistHDibaryonInvaMassGenRes->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");
+
+ fHistAntiHDibaryonInvaMassGen = new TH1F("fHistAntiHDibaryonInvaMassGen", "Generated  #bar{#Lambda} #bar{p} #pi^{+}", 200, 2.1, 2.3);
+ fHistAntiHDibaryonInvaMassGen->GetYaxis()->SetTitle("Counts");
+ fHistAntiHDibaryonInvaMassGen->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");
+
+ fHistHDibaryonInvaMassAso = new TH1F("fHistHDibaryonInvaMassAso", "Associated  #Lambda p #pi^{-}", 200, 2.1, 2.3);
+ fHistHDibaryonInvaMassAso->GetYaxis()->SetTitle("Counts");
+ fHistHDibaryonInvaMassAso->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");
+
+ fHistHDibaryonInvaMassAsoReso = new TH1F("fHistHDibaryonInvaMassAsoReso", "Associated  #Lambda p #pi^{-}", 200, 2.1, 2.3);
+ fHistHDibaryonInvaMassAsoReso->GetYaxis()->SetTitle("Counts");
+ fHistHDibaryonInvaMassAsoReso->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");
+
+ fHistAntiHDibaryonInvaMassAso = new TH1F("fHistAntiHDibaryonInvaMassAso", "Associated  #bar{#Lambda} #bar{p} #pi^{+}", 200, 2.1, 2.3);
+ fHistAntiHDibaryonInvaMassAso->GetYaxis()->SetTitle("Counts");
+ fHistAntiHDibaryonInvaMassAso->GetXaxis()->SetTitle("Invariant mass  (GeV/c^{2})");
+
+ fHistCheck = new TH2F("fHistCheck", "Check online/offline", 200, -0.5, 1.5, 200, -0.5, 1.5);
+ fHistCheck->GetXaxis()->SetTitle("offline");
+ fHistCheck->GetYaxis()->SetTitle("online");
+ fHistCheck->SetMarkerStyle(kFullCircle);
+
+ fHistHPointingAngle= new TH1F("fHistHPointingAngle", "Pointing angle distribution for #Lambdap#pi^{-}", 200, 0., 2*TMath::Pi());
+ fHistHPointingAngle->GetXaxis()->SetTitle("Pointing angle distribution for #Lambdap#pi^{-}");
+ fHistHPointingAngle->GetYaxis()->SetTitle("Entries");
+ fHistHPointingAngle->SetMarkerStyle(kFullCircle);
+
+ fHistMassH= new TH1F("fHistMassH", "Invariant mass distribution #Lambdap#pi^{-} (GeV/c^{2}) after pointing angle cut", 3000, 2.2, 2.5);
+ fHistMassH->GetXaxis()->SetTitle("Invariant mass #Lambdap#pi^{-} (GeV/c^{2})");
+ fHistMassH->GetYaxis()->SetTitle("Entries");
+ fHistMassH->SetMarkerStyle(kFullCircle);
+ fHistMassLambdaFromH= new TH1F("fHistMassLambdaFromH", "Invariant mass distribution #Lambda (GeV/c^{2}) asking for Mother to be H", 300, 1.0, 1.3);
+ fHistMassLambdaFromH->GetXaxis()->SetTitle("Invariant mass #Lambda (GeV/c^{2})");
+ fHistMassLambdaFromH->GetYaxis()->SetTitle("Entries");
+ fHistMassLambdaFromH->SetMarkerStyle(kFullCircle);
+
+ fHistMassLambdaFromHtLorentz= new TH1F("fHistMassLambdaFromHtLorentz", "Invariant mass distribution #Lambda (GeV/c^{2}) asking for Mother to be H", 300, 1.0, 1.3);
+ fHistMassLambdaFromHtLorentz->GetXaxis()->SetTitle("Invariant mass #Lambda (GeV/c^{2})");
+ fHistMassLambdaFromHtLorentz->GetYaxis()->SetTitle("Entries");
+ fHistMassLambdaFromHtLorentz->SetMarkerStyle(kFullCircle);
+ fHistMassPpi= new TH1F("fHistMassPpi", "Invariant mass distribution of the p#pi^{-} used for combing with Lambda", 300, 1.0, 1.6);
+ fHistMassPpi->GetXaxis()->SetTitle("Invariant mass p#pi^{-} (GeV/c^{2})");
+ fHistMassPpi->GetYaxis()->SetTitle("Entries");
+ fHistMassPpi->SetMarkerStyle(kFullCircle);
+
+ fHistMassPpiReso= new TH1F("fHistMassPpiReso", "Invariant mass distribution of the p#pi^{-} used for combing with Lambda", 300, 1.0, 1.6);
+ fHistMassPpiReso->GetXaxis()->SetTitle("Invariant mass p#pi^{-} (GeV/c^{2})");
+ fHistMassPpiReso->GetYaxis()->SetTitle("Entries");
+ fHistMassPpiReso->SetMarkerStyle(kFullCircle);
+
+ fHistMassLpi= new TH1F("fHistMassLpi", "Invariant mass distribution of the #Lambda#pi^{-} used for combing with p", 300, 1.1, 1.7);
+ fHistMassLpi->GetXaxis()->SetTitle("Invariant mass #Lambda#pi^{-} (GeV/c^{2})");
+ fHistMassLpi->GetYaxis()->SetTitle("Entries");
+ fHistMassLpi->SetMarkerStyle(kFullCircle);
+
+ fHistMassLP= new TH1F("fHistMassLP", "Invariant mass distribution of the #Lambda p used for combing with #pi^{-}", 300, 2.0, 2.3);
+ fHistMassLP->GetXaxis()->SetTitle("Invariant mass #Lambda p (GeV/c^{2})");
+ fHistMassLP->GetYaxis()->SetTitle("Entries");
+ fHistMassLP->SetMarkerStyle(kFullCircle);
+
+ fHistProtonPIDBb = new TH2F("fHistProtonPIDBb", "dE/dx after p PID", 100, 0., 10, 100, 0, 100);
+ fHistProtonPIDBb->GetYaxis()->SetTitle("TPC Signal");
+ fHistProtonPIDBb->GetXaxis()->SetTitle("P (GeV/c)");
+ fHistProtonPIDBb->SetOption("scat");
+ fHistProtonPIDBb->SetMarkerStyle(kFullCircle);
+
+ fHistPionPIDBb = new TH2F("fHistPionPIDBb", "dE/dx after K PID", 100, 0., 10, 100, 0, 100);
+ fHistPionPIDBb->GetYaxis()->SetTitle("TPC Signal");
+ fHistPionPIDBb->GetXaxis()->SetTitle("P (GeV/c)");
+ fHistPionPIDBb->SetOption("scat");
+ fHistPionPIDBb->SetMarkerStyle(kFullCircle);
+
+ fHistProtonPIDLambda = new TH2F("fHistProtonPIDLambda", "dE/dx after p PID from V0", 100, 0., 10, 100, 0, 100);
+ fHistProtonPIDLambda->GetYaxis()->SetTitle("TPC Signal");
+ fHistProtonPIDLambda->GetXaxis()->SetTitle("P (GeV/c)");
+ fHistProtonPIDLambda->SetOption("scat");
+ fHistProtonPIDLambda->SetMarkerStyle(kFullCircle);
+
+ fHistPionPIDLambda = new TH2F("fHistPionPIDLambda", "dE/dx after #pi PID from V0", 100, 0, 10, 100, 0, 100);
+ fHistPionPIDLambda->GetYaxis()->SetTitle("TPC Signal");
+ fHistPionPIDLambda->GetXaxis()->SetTitle("P (GeV/c)");
+ fHistPionPIDLambda->SetOption("scat");
+ fHistPionPIDLambda->SetMarkerStyle(kFullCircle);
+
+ fHistMCdcaPvtxDvtx= new TH1F("fHistMCdcaPvtxDvtx", "MC True DCA Primary Vertex - H Decay Vertex", 300, -0.1, 11.9);
+ fHistMCdcaPvtxDvtx->GetXaxis()->SetTitle("dca prim. vtx- decay vtx (cm)");
+ fHistMCdcaPvtxDvtx->GetYaxis()->SetTitle("Entries");
+ fHistMCdcaPvtxDvtx->SetMarkerStyle(kFullCircle);
+
+ fHistMCdcaPvtxLvtx= new TH1F("fHistMCdcaPvtxLvtx", "MC True DCA Primary Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);
+ fHistMCdcaPvtxLvtx->GetXaxis()->SetTitle("dca prim. vtx-#Lambda decay vtx (cm)");
+ fHistMCdcaPvtxLvtx->GetYaxis()->SetTitle("Entries");
+ fHistMCdcaPvtxLvtx->SetMarkerStyle(kFullCircle);
+
+ fHistMCdcaDvtxLvtx= new TH1F("fHistMCdcaDvtxLvtx", "MC True DCA H Decay Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);
+ fHistMCdcaDvtxLvtx->GetXaxis()->SetTitle("dca H deacy vtx-#Lambda decay vtx (cm)");
+ fHistMCdcaDvtxLvtx->GetYaxis()->SetTitle("Entries");
+ fHistMCdcaDvtxLvtx->SetMarkerStyle(kFullCircle);
+
+ fHistMCangleLH= new TH1F("fHistMCangleLH", "MC True Angle between #Lambda and H", 300, 0., 2*TMath::Pi());
+ fHistMCangleLH->GetXaxis()->SetTitle("Angle (#Lambda-H)");
+ fHistMCangleLH->GetYaxis()->SetTitle("Entries");
+ fHistMCangleLH->SetMarkerStyle(kFullCircle);
+
+ fHistMCdecayAngle= new TH1F("fHistMCdecayAngle", "MC True Angle between decay products", 300, 0., 2*TMath::Pi());
+ fHistMCdecayAngle->GetXaxis()->SetTitle("Angle (#Lambda-p#pi)");
+ fHistMCdecayAngle->GetYaxis()->SetTitle("Entries");
+ fHistMCdecayAngle->SetMarkerStyle(kFullCircle);
+
+ fHistMCpointingAngle= new TH1F("fHistMCpointingAngle", "MC True Pointing Angle", 3000, 0., 2*TMath::Pi());
+ fHistMCpointingAngle->GetXaxis()->SetTitle("Pointing Angle");
+ fHistMCpointingAngle->GetYaxis()->SetTitle("Entries");
+ fHistMCpointingAngle->SetMarkerStyle(kFullCircle);
+
+ fHistMCap = new TH2F("fHistMCap", "True MC Armenteros-Podolanski", 200,-1.0,1.0, 500,0,1);
+ fHistMCap->GetYaxis()->SetTitle("#alpha");
+ fHistMCap->GetXaxis()->SetTitle("q_{t}");
+ fHistMCap->SetOption("scat");
+ fHistMCap->SetMarkerStyle(kFullCircle);
+
+ fHistMCdcaPvtxDvtxReso= new TH1F("fHistMCdcaPvtxDvtxReso", "MC True DCA Primary Vertex - H Decay Vertex", 300, -0.1, 11.9);
+ fHistMCdcaPvtxDvtxReso->GetXaxis()->SetTitle("dca prim. vtx- decay vtx (cm)");
+ fHistMCdcaPvtxDvtxReso->GetYaxis()->SetTitle("Entries");
+ fHistMCdcaPvtxDvtxReso->SetMarkerStyle(kFullCircle);
+
+ fHistMCdcaPvtxLvtxReso= new TH1F("fHistMCdcaPvtxLvtxReso", "MC True DCA Primary Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);
+ fHistMCdcaPvtxLvtxReso->GetXaxis()->SetTitle("dca prim. vtx-#Lambda decay vtx (cm)");
+ fHistMCdcaPvtxLvtxReso->GetYaxis()->SetTitle("Entries");
+ fHistMCdcaPvtxLvtxReso->SetMarkerStyle(kFullCircle);
+
+ fHistMCdcaDvtxLvtxReso= new TH1F("fHistMCdcaDvtxLvtxReso", "MC True DCA H Decay Vertex - Lambda Decay Vertex", 300, -0.1, 11.9);
+ fHistMCdcaDvtxLvtxReso->GetXaxis()->SetTitle("dca H deacy vtx-#Lambda decay vtx (cm)");
+ fHistMCdcaDvtxLvtxReso->GetYaxis()->SetTitle("Entries");
+ fHistMCdcaDvtxLvtxReso->SetMarkerStyle(kFullCircle);
+
+ fHistMCangleLHReso= new TH1F("fHistMCangleLHReso", "MC True Angle between #Lambda and H", 300, 0., 2*TMath::Pi());
+ fHistMCangleLHReso->GetXaxis()->SetTitle("Angle (#Lambda-H)");
+ fHistMCangleLHReso->GetYaxis()->SetTitle("Entries");
+ fHistMCangleLHReso->SetMarkerStyle(kFullCircle);
+
+ fHistMCdecayAngleReso= new TH1F("fHistMCdecayAngleReso", "MC True Angle between decay products", 300, 0., 2*TMath::Pi());
+ fHistMCdecayAngleReso->GetXaxis()->SetTitle("Angle (#Lambda-p#pi)");
+ fHistMCdecayAngleReso->GetYaxis()->SetTitle("Entries");
+ fHistMCdecayAngleReso->SetMarkerStyle(kFullCircle);
+
+ fHistMCpointingAngleReso= new TH1F("fHistMCpointingAngleReso", "MC True Pointing Angle", 300, 0., 2*TMath::Pi());
+ fHistMCpointingAngleReso->GetXaxis()->SetTitle("Pointing Angle");
+ fHistMCpointingAngleReso->GetYaxis()->SetTitle("Entries");
+ fHistMCpointingAngleReso->SetMarkerStyle(kFullCircle);
+
+ fHistMCapReso = new TH2F("fHistMCapReso", "True MC Armenteros-Podolanski", 200,-1.0,1.0, 500,0,1);
+ fHistMCapReso->GetYaxis()->SetTitle("#alpha");
+ fHistMCapReso->GetXaxis()->SetTitle("q_{t}");
+ fHistMCapReso->SetOption("scat");
+ fHistMCapReso->SetMarkerStyle(kFullCircle);
+
+ fHistCentrality = new TH1F("Centrality ", "Centrality", 11, -0.5, 10.5);
+ fHistCentrality ->GetXaxis()->SetTitle("Centrality");
+ fHistCentrality ->GetYaxis()->SetTitle("Entries");
+
+ fHistCentralityAC = new TH1F("CentralityAC ", "CentralityAC", 11, -0.5, 10.5);
+ fHistCentralityAC ->GetXaxis()->SetTitle("Centrality");
+ fHistCentralityAC ->GetYaxis()->SetTitle("Entries");
+
+ fHistMultiplicity = new TH1F("Multiplicity ", "Multiplicity", 100, 0, 20000);
+ fHistMultiplicity ->GetXaxis()->SetTitle("Centrality");
+ fHistMultiplicity ->GetYaxis()->SetTitle("Entries");
+
+ fHistList->Add(fHistMassDPi);
+ fHistList->Add(fHistMassLPi);
+ fHistList->Add(fHistMassLambdaPi);
+ fHistList->Add(fHistMassLambda);
+ fHistList->Add(fHistMassLambdaPPi);
+ fHistList->Add(fHistMassLambdaP);
+ fHistList->Add(fHistMassLambdaK);
+ fHistList->Add(fHistMassK0onl);
+ fHistList->Add(fHistMassK0offl);
+ fHistList->Add(fHistMassK0onlC);
+ fHistList->Add(fHistMassK0offlC);
+ fHistList->Add(fHistMassPQonl); 
+ fHistList->Add(fHistMassPQoffl);
+ fHistList->Add(fHistDC); 
+ fHistList->Add(fHistArmenterosPodolanski);
+ fHistList->Add(fHistArmenterosPodolanskiCut);
+ fHistList->Add(fHistHDibaryonInvaMassGen);
+ fHistList->Add(fHistHDibaryonInvaMassGenRes);
+ fHistList->Add(fHistAntiHDibaryonInvaMassGen);
+ fHistList->Add(fHistHDibaryonInvaMassAso);
+ fHistList->Add(fHistHDibaryonInvaMassAsoReso);
+ fHistList->Add(fHistAntiHDibaryonInvaMassAso);
+ fHistList->Add(fHistCheck);
+ fHistList->Add(fHistHPointingAngle);
+ fHistList->Add(fHistMassH);
+ fHistList->Add(fHistMassPpi);
+ fHistList->Add(fHistMassPpiReso);
+ fHistList->Add(fHistMassLpi);
+ fHistList->Add(fHistMassLP);
+ fHistList->Add(fHistMassLambdaFromH);
+ fHistList->Add(fHistMassLambdaFromHtLorentz);
+ fHistList->Add(fHistProtonPIDBb);
+ fHistList->Add(fHistPionPIDBb);
+ fHistList->Add(fHistProtonPIDLambda);
+ fHistList->Add(fHistPionPIDLambda);
+ fHistList->Add(fHistMCdcaPvtxDvtx);
+ fHistList->Add(fHistMCdcaPvtxLvtx);
+ fHistList->Add(fHistMCdcaDvtxLvtx);
+ fHistList->Add(fHistMCangleLH);
+ fHistList->Add(fHistMCdecayAngle);
+ fHistList->Add(fHistMCpointingAngle);
+ fHistList->Add(fHistMCap);
+ fHistList->Add(fHistMCdcaPvtxDvtxReso);
+ fHistList->Add(fHistMCdcaPvtxLvtxReso);
+ fHistList->Add(fHistMCdcaDvtxLvtxReso);
+ fHistList->Add(fHistMCangleLHReso);
+ fHistList->Add(fHistMCdecayAngleReso);
+ fHistList->Add(fHistMCpointingAngleReso);
+ fHistList->Add(fHistMCapReso);
+ fHistList->Add(fHistCentrality);
+ fHistList->Add(fHistCentralityAC);
+ fHistList->Add(fHistMultiplicity);
+ fHistHilf1= new TH1F("fHistHilf1", "HD", 300, 0., 10);
+ fHistHilf1->GetXaxis()->SetTitle("");
+ fHistHilf1->GetYaxis()->SetTitle("Entries");
+
+ fHistHilf2= new TH1F("fHistHilf2", "HD", 300, 0., 10);
+ fHistHilf2->GetXaxis()->SetTitle("");
+ fHistHilf2->GetYaxis()->SetTitle("Entries");
+
+ fHistHilf3= new TH1F("fHistHilf3", "HD", 300, 0., 10);
+ fHistHilf3->GetXaxis()->SetTitle("");
+ fHistHilf3->GetYaxis()->SetTitle("Entries");
+
+ fHistHilf4= new TH1F("fHistHilf4", "HD", 300, 0., 10);
+ fHistHilf4->GetXaxis()->SetTitle("");
+ fHistHilf4->GetYaxis()->SetTitle("Entries");
+
+ fHistHilf5= new TH1F("fHistHilf5", "HD", 300, 0., 10);
+ fHistHilf5->GetXaxis()->SetTitle("");
+ fHistHilf5->GetYaxis()->SetTitle("Entries");
+
+ fHistHilf6= new TH1F("fHistHilf6", "HD", 300, 0., 10);
+ fHistHilf6->GetXaxis()->SetTitle("");
+ fHistHilf6->GetYaxis()->SetTitle("Entries");
+ fHistPtvsEtaGen = new TH2F("fHistPtvsEtaGen", "p_{t} vs #eta from generated H", 200,0.0,10.0, 200,-1,1);
+ fHistPtvsEtaGen->GetXaxis()->SetTitle("p_{t} (GeV/c)");
+ fHistPtvsEtaGen->GetYaxis()->SetTitle("#eta");
+ fHistPtvsEtaGen->SetOption("scat");
+ fHistPtvsEtaGen->SetMarkerStyle(kFullCircle);
+
+ fHistPtvsEtaAso = new TH2F("fHistPtvsEtaAso", "p_{t} vs #eta from associated H", 200,0.0,10.0, 200,-1,1);
+ fHistPtvsEtaAso->GetYaxis()->SetTitle("p_{t} (GeV/c)");
+ fHistPtvsEtaAso->GetXaxis()->SetTitle("#eta");
+ fHistPtvsEtaAso->SetOption("scat");
+ fHistPtvsEtaAso->SetMarkerStyle(kFullCircle);
+
+ fHistPtvsYGen = new TH2F("fHistPtvsYGen", "p_{t} vs rapidity from generated H", 200,0.0,10.0, 200,-1,1);
+ fHistPtvsYGen->GetXaxis()->SetTitle("p_{t} (GeV/c)");
+ fHistPtvsYGen->GetYaxis()->SetTitle("y");
+ fHistPtvsYGen->SetOption("scat");
+ fHistPtvsYGen->SetMarkerStyle(kFullCircle);
+
+ fHistPtvsYAso = new TH2F("fHistPtvsYAso", "p_{t} vs rapidity from associated H", 200,0.0,10.0, 200,-1,1);
+ fHistPtvsYAso->GetXaxis()->SetTitle("p_{t} (GeV/c)");
+ fHistPtvsYAso->GetYaxis()->SetTitle("y");
+ fHistPtvsYAso->SetOption("scat");
+ fHistPtvsYAso->SetMarkerStyle(kFullCircle);
+
+ fHistRap= new TH1F("fHistRap", "Rapidity", 400, -2., 2);
+ fHistRap->GetXaxis()->SetTitle("Y");
+ fHistRap->GetYaxis()->SetTitle("Entries");
+ fHistPtvsEtaAso->SetMarkerStyle(kFullCircle);
+
+ fHistList->Add(fHistHilf1);
+ fHistList->Add(fHistHilf2); 
+ fHistList->Add(fHistHilf3);
+ fHistList->Add(fHistHilf4);
+ fHistList->Add(fHistHilf5); 
+ fHistList->Add(fHistHilf6);
+ fHistList->Add(fHistPtvsEtaGen);
+ fHistList->Add(fHistPtvsEtaAso);
+ fHistList->Add(fHistPtvsYGen);
+ fHistList->Add(fHistPtvsYAso);
+ fHistList->Add(fHistRap);
+
+ fHistCount = new TH1F("fHistCount","test",17,0,17);
+ fHistCount->GetXaxis()->SetBinLabel(1,"Events");
+ fHistCount->GetXaxis()->SetBinLabel(2,"MC All");
+ fHistCount->GetXaxis()->SetBinLabel(3,"MC from Primary Vtx");
+ fHistCount->GetXaxis()->SetBinLabel(4,"Horst");
+ fHistCount->GetXaxis()->SetBinLabel(5,"Lambda Candidates");
+ fHistCount->GetXaxis()->SetBinLabel(6,"Sigma Candidates");
+ fHistCount->GetXaxis()->SetBinLabel(7,"Horst");
+ fHistCount->GetXaxis()->SetBinLabel(8,"Horst");
+ fHistCount->GetXaxis()->SetBinLabel(9,"Horst");
+ fHistCount->GetXaxis()->SetBinLabel(10,"MC All #bar{Lambda}(1520)s");
+ fHistCount->GetXaxis()->SetBinLabel(11,"");
+ fHistCount->GetXaxis()->SetBinLabel(12,"H-Dibaryon");
+ fHistCount->GetXaxis()->SetBinLabel(13,"Hypertriton 2-Body");
+ fHistCount->GetXaxis()->SetBinLabel(14,"Hypertriton 3-Body");
+ fHistCount->GetXaxis()->SetBinLabel(15,"");
+ fHistCount->GetXaxis()->SetBinLabel(16,"");
+ fHistCount->GetXaxis()->SetBinLabel(17,"Lambdas!!!");
+ fHistCount->SetStats(0);
+ fHistCount->SetFillColor(38);
+ fHistList->Add(fHistCount);
+
+ //trigger statistics histogram
+  fHistTriggerStat = new TH1F("fHistTriggerStat","Trigger statistics", 4,-0.5, 3.5);
+  const Char_t* aTriggerNames[] = { "kMB", "kCentral", "kSemiCentral" };
+  for ( Int_t ii=0; ii < 3; ii++ )
+    fHistTriggerStat->GetXaxis()->SetBinLabel(ii+1, aTriggerNames[ii]);
+
+  fHistTriggerStatAfterEventSelection = new TH1F("fHistTriggerStatAfterEventSelection","Trigger statistics after event selection", 4,-0.5, 3.5);
+  for ( Int_t ii=0; ii < 3; ii++ )
+    fHistTriggerStatAfterEventSelection->GetXaxis()->SetBinLabel(ii+1, aTriggerNames[ii]);
+  fHistList->Add(fHistTriggerStat);
+  fHistList->Add(fHistTriggerStatAfterEventSelection);
+
+  fHistMassHcentMult =  new TH3F("fHistMassHcentMult", "Inv. Mass vs. centrality vs. multiplicity", 100, 2.2, 2.3, 5, 0, 4, 300, 0, 6000);
+  fHistMassHcentMult->GetXaxis()->SetTitle("Invariant mass #Lambdap#pi^{-} (GeV/c^{2})"); // inv. mass
+  fHistMassHcentMult->GetYaxis()->SetTitle("Centrality"); // triggertype
+  fHistMassHcentMult->GetZaxis()->SetTitle("Multiplicity"); // refTPC
+  fHistList->Add(fHistMassHcentMult);
+
+  
+  const Double_t kz = 2*TMath::Pi();
+  Int_t binsD01[16]={  300, 200, 100, 100, 100, 100, 100, 100, 200, 200, 200, 200, 400, 200, 200, 3};
+  Double_t xminD01[16]={2.0, 1.0, 0., -1, 0., 0., 0., 0., 0., 0., 0., 0., -1, 0.,  0., 0};
+  Double_t xmaxD01[16]={2.3, 1.2, kz, 1, 1, 10, 10, 5, 5, 5, 5, 100, 1, 100,  4000, 1};
+  
+
+  fHistNdim = new THnSparseF("fHistNdim","THnS;InvMass, InvMassLambda, pointingAngle, armPoAlpha, armPoQt, pTL, pTH, d0p, d0n, dcaHd, dca, decayL, cosPA, centr, multi, mcinf;InvMassH", 16,binsD01,xminD01,xmaxD01);
+  fHistList->Add(fHistNdim);
+
+ AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
+ AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
+ fPIDtpcESD = inputHandler->GetPIDResponse();
+
+// Post output data (if histograms are not used later, PostData is at least called here)
+  PostData(1, fHistList);
+
+}
+
+ //________________________________________________________________________
+void AliAnalysisTaskHdibaryonLPpi::UserExec(Option_t *) 
+{ 
+  // Main loop
+  // Called for each event
+
+  //define improtant masses
+  Double_t pionMass     = 0.13957;
+  Double_t protonMass   = 0.93827;
+
+  //define PDGCodes
+  Long_t pdgPionPlus          =         211;
+  Long_t pdgPionMinus         =        -211;
+  Long_t pdgProton            =        2212;
+  Long_t pdgAntiProton        =       -2212;
+  Long_t pdgLambda            =        3122;
+  Long_t pdgAntiLambda        =       -3122;
+  Long_t pdgHDibaryon         =  1020000020;
+  Long_t pdgAntiHDibaryon     = -1020000020;
+
+  AliStack* stack(NULL);
+  if(HasMC()){
+
+    if(!fMCEvent)return;
+
+    AliHeader *head = fMCEvent->Header();
+    if(!head)return;
+    AliGenPythiaEventHeader *pyheader = (AliGenPythiaEventHeader*)head->GenEventHeader();
+    if(!pyheader)return;
+    
+    if(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()){
+      if(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->MCEvent()){
+       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->TreeK()))return;
+       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->TreeTR()))return;
+       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->Init("local")))return;
+       stack = static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->MCEvent()->Stack();
+       if(!(static_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler())->MCEvent()->Stack()->TreeK()))return;
+      }
+    }
+       
+    if(!stack)return;
+  }
+
+  // -------------------------------------------------------
+  // Loop for Inv. Mass via ESD tracks
+  // -------------------------------------------------------
+
+  fHistCount->Fill(0);
+  
+  fESD = dynamic_cast<AliESDEvent *>(fInputEvent);
+
     if (!fESD) {
       //Printf("ERROR: fESD not available");
       return;
     }
 
 
-  const AliESDVertex *vertex = fESD->GetPrimaryVertexTracks();\r
-  if (vertex->GetNContributors()<1) \r
-    {\r
-      // SPD vertex\r
-      vertex = fESD->GetPrimaryVertexSPD();\r
-       if(vertex->GetNContributors()<1) {\r
-        return;\r
-       }  \r
-    }\r
-  if (TMath::Abs(vertex->GetZv()) > 10) return;\r
-  \r
-  Int_t centrality = -5;\r
-  Double_t centrPerc = -5;\r
-\r
-  if (fESD->GetEventSpecie() == 4) \r
-    { // PbPb\r
-      AliCentrality *esdCentrality = fESD->GetCentrality();\r
-      centrality = esdCentrality->GetCentralityClass10("V0M"); // centrality percentile determined with V0\r
-      centrPerc = esdCentrality->GetCentralityPercentile("V0M");\r
-//      if (centrality < 0. || centrality > 8.) return; //0 bis 80 %\r
-      if (centrality > 8) return; //0 bis 80 %\r
-      //  cout<<"Centrality: "<< centrality << endl;\r
-    }\r
-\r
-\r
-  fHistCentrality->Fill(centrality);   \r
-\r
-  //*****************//  \r
-  //*   Centrality  *//\r
-  //*****************//\r
\r
-  //  Float_t percentile=centrality->GetCentralityPercentile("V0M");\r
-\r
-  Bool_t isSelectedCentral = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kCentral);\r
-  Bool_t isSelectedSemiCentral = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kSemiCentral);\r
-  Bool_t isSelectedMB = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kMB);\r
\r
-  Int_t triggertype = 17;\r
-\r
-  if(isSelectedCentral){\r
-    fHistTriggerStat->Fill(1);\r
-    triggertype=1;\r
-  }\r
-\r
-  if(isSelectedSemiCentral){\r
-    fHistTriggerStat->Fill(2);\r
-    triggertype=2;\r
-  }\r
-\r
-  if(isSelectedMB){\r
-    fHistTriggerStat->Fill(0);\r
-    triggertype=3;\r
-  }\r
-\r
-  //  if(isSelectedCentral || isSelectedSemiCentral || isSelectedMB){\r
-\r
-  //*******************************\r
-\r
-  Int_t runNumber = 0;\r
-  //  itrk = 0;\r
-  runNumber = fESD->GetRunNumber();\r
-/*  \r
-  if (!fPIDtpcESD) fPIDtpcESD = ((AliESDInputHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->GetESDpid();\r
-  if (!fPIDtpcESD) {\r
-    fPIDtpcESD = new AliESDpid(); // HACK FOR MC PBPB --> PLEASE REMOVE AS SOON AS POSSIBLE\r
-    fPIDtpcESD->GetTPCResponse().SetBetheBlochParameters(1.28778e+00/50., 3.13539e+01, TMath::Exp(-3.16327e+01), 1.87901e+00, 6.41583e+00);\r
-  }\r
-*/\r
-  Double_t pionK=1;\r
-  Double_t pK=1;\r
\r
-  TObjArray* listCrossV0   = fESDCutsV0->GetAcceptedV0s(fESD);  \r
-  Int_t nGoodV0s      = listCrossV0->GetEntries();\r
-\r
-  const AliESDVertex *esdVer = fESD->GetPrimaryVertex();\r
-  AliESDVertex *esdVer1 = new AliESDVertex(*esdVer);\r
-   \r
-  AliVertexerTracks *vertexer = new AliVertexerTracks(fESD->GetMagneticField());\r
-  TObjArray *trkArray = new TObjArray(2);\r
-  AliVertexerTracks *vertexer1 = new AliVertexerTracks(fESD->GetMagneticField());\r
-  TObjArray *trkArray1 = new TObjArray(2);\r
-\r
-  AliKFParticle::SetField(fESD->GetMagneticField());\r
-       \r
-  AliKFVertex primVtx(*(fESD->GetPrimaryVertex()));\r
-\r
-  Int_t refMultTpc = AliESDtrackCuts::GetReferenceMultiplicity(fESD, kTRUE);\r
-  //cout<<"Multiplicity: "<< refMultTpc << endl;\r
-  fHistMultiplicity->Fill(refMultTpc);\r
-  \r
-  Double_t mn[3] = {0,0,0};\r
-  Double_t mp[3] = {0,0,0}; \r
-  Double_t mm[3] = {0,0,0};\r
-  Double_t dd[3] = {0,0,0};\r
-  Double_t dd1[3] = {0,0,0};\r
-  const Double_t cProtonMass=TDatabasePDG::Instance()->GetParticle(2212)->Mass();\r
-  const Double_t cPionMass=TDatabasePDG::Instance()->GetParticle(211)->Mass();\r
-  const Double_t cElectronMass=TDatabasePDG::Instance()->GetParticle(11)->Mass();\r
-  const Double_t cLambdaMass=TDatabasePDG::Instance()->GetParticle(3122)->Mass();\r
-  Double_t decayLength=0;\r
-  Double_t decayLengthH=0;\r
-\r
-  //V0 Loop for Lambda and Anti-Lambda\r
-  for(Int_t iV0MI = 0; iV0MI < nGoodV0s ; iV0MI++) {\r
-    AliESDv0 * fV0MIs = fESD->GetV0(iV0MI);\r
-    Int_t    lOnFlyStatus = 0;\r
-\r
-    lOnFlyStatus = fV0MIs->GetOnFlyStatus();\r
-    Double_t lInvMassLambda=0;\r
-    Double_t lInvMassLambdaPi=0;\r
-    Double_t lPtLambda=0;\r
-    Double_t lPzLambda=0;\r
-    Double_t lPLambda=0;\r
-    Int_t onl=0;\r
-    Int_t offl=0;\r
-\r
-    TLorentzVector posE;\r
-    TLorentzVector negE;\r
-    TLorentzVector photon;\r
-\r
-    if (lOnFlyStatus){\r
-      onl=1; \r
-      //      return;\r
-    }\r
-    if (!lOnFlyStatus){\r
-      offl=1;\r
-      //return;\r
-    }\r
-\r
-    //    fHistMultiplicity->Fill(refMultTpc); \r
-    fHistCentralityAC->Fill(centrality);\r
-    fHistCheck->Fill(offl,onl);\r
-\r
-    AliESDtrack* trackPosTest = fESD->GetTrack(fV0MIs->GetPindex());\r
-    AliESDtrack* trackNegTest = fESD->GetTrack(fV0MIs->GetNindex());\r
-\r
-    //    if (!\r
-    if (!fEsdTrackCuts->AcceptTrack(trackPosTest)) continue;\r
-    if (!fESDtrackCutsV0->AcceptTrack(trackNegTest)) continue;\r
-\r
-\r
-      //PID via specific energy loss in the TPC\r
-      //define the arrays for the Bethe-Bloch-Parameters\r
-      Double_t parProton[5]   = {0,0,0,0,0};\r
-\r
-      if(runNumber < 166500) //LHC10h\r
-       {\r
-         parProton[0]   = 1.45802; // ALEPH parameters for protons (pass2)\r
-         parProton[1]   = 27.4992;\r
-         parProton[2]   = 4.00313e-15;\r
-         parProton[3]   = 2.48485;\r
-         parProton[4]   = 8.31768; \r
-       }\r
-      \r
-      if(runNumber > 166500) //LHC11h\r
-       { \r
-         parProton[0]   = 1.11243; // ALEPH parameters for protons (pass2)\r
-         parProton[1]   = 26.1144;\r
-         parProton[2]   = 4.00313e-15;\r
-         parProton[3]   = 2.72969;\r
-         parProton[4]   = 9.15038; \r
-       }\r
\r
-      //Get the total momentum for each track at the inner readout of the TPC\r
-      Double_t ptotN = trackNegTest->GetInnerParam()->GetP();\r
-      Double_t ptotP = trackPosTest->GetInnerParam()->GetP();\r
-\r
-      //define expected signals for the various species\r
-      Double_t expSignalPionP = 0;\r
-      Double_t expSignalPionN = 0;\r
-      Double_t expSignalProtonN = 0;\r
-      Double_t expSignalProtonP = 0;\r
-\r
-      //for data\r
-      if(!HasMC())\r
-       {\r
-         expSignalProtonN = AliExternalTrackParam::BetheBlochAleph(ptotN/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);\r
-         expSignalProtonP = AliExternalTrackParam::BetheBlochAleph(ptotP/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);\r
-       }\r
-    \r
-      //for MC\r
-      if(HasMC())\r
-       {\r
-         expSignalPionP = 0.7*AliExternalTrackParam::BetheBlochAleph(ptotP/(pionMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]); \r
-         expSignalPionN = 0.7*AliExternalTrackParam::BetheBlochAleph(ptotN/(pionMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]); \r
-\r
-         expSignalProtonN = 0.65*AliExternalTrackParam::BetheBlochAleph(ptotN/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);\r
-         expSignalProtonP = 0.65*AliExternalTrackParam::BetheBlochAleph(ptotP/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);   \r
-       }\r
-\r
-      // PID cut on the nuclei (proton, deuteron, triton, helium3)\r
-      Bool_t corrParticle = kFALSE;\r
-\r
-      Bool_t posProton = kFALSE;\r
-      //data\r
-      if(!HasMC())\r
-       {\r
-         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackPosTest, AliPID::kProton)) < 3)\r
-           {\r
-             posProton = kTRUE;\r
-             corrParticle = kTRUE;\r
-           }\r
-       }\r
-      //MC\r
-      if(HasMC())\r
-       {\r
-         if(//trackPosTest->GetTPCsignal() < 1200 && \r
-            TMath::Abs(trackPosTest->GetTPCsignal() - expSignalProtonP)/expSignalProtonP < 0.4)\r
-           {\r
-             posProton = kTRUE;\r
-             corrParticle = kTRUE;\r
-           }\r
-       }\r
-\r
-      Bool_t negProton = kFALSE;\r
-      //data\r
-      if(!HasMC())\r
-       {\r
-         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackNegTest, AliPID::kProton)) < 3)\r
-           {\r
-             negProton = kTRUE;\r
-             corrParticle = kTRUE;\r
-           }\r
-       }\r
-     //MC\r
-      if(HasMC())\r
-       {\r
-         if(//trackNegTest->GetTPCsignal() < 1200 && \r
-            TMath::Abs(trackNegTest->GetTPCsignal() - expSignalProtonN)/expSignalProtonN < 0.4)\r
-           {\r
-             negProton = kTRUE;\r
-             corrParticle = kTRUE;\r
-           }\r
-       }\r
-\r
-     //PID cut for pions\r
-      //data: 3sigma cut on the pions\r
-\r
-      Bool_t negPion = kFALSE;     \r
-      Bool_t posPion = kFALSE;\r
-\r
-      if (!HasMC())\r
-       {\r
-         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackPosTest, AliPID::kPion)) < 4) posPion=kTRUE; //pos daughter has to be a pion\r
-         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackNegTest, AliPID::kPion)) < 4) negPion=kTRUE; // negative daughter has to be a pion\r
-       }\r
-      \r
-      //MC: like the nuclei via the specific energyloss in the TPC\r
-      if (HasMC())\r
-       {\r
-         if (TMath::Abs(trackPosTest->GetTPCsignal() - expSignalPionP)/expSignalPionP < 0.4) posPion=kTRUE;\r
-         if (TMath::Abs(trackNegTest->GetTPCsignal() - expSignalPionN)/expSignalPionN < 0.4) negPion=kTRUE;\r
-       }\r
-\r
-      if (!(posProton==kTRUE)) continue;\r
-      if (!(negPion==kTRUE)) continue;\r
-\r
-    //To avoid ghosts\r
-\r
-    if( !(trackPosTest->GetStatus() & AliESDtrack::kTPCrefit)){\r
-      continue;\r
-    }\r
-\r
-    if( !(trackNegTest->GetStatus() & AliESDtrack::kTPCrefit)){\r
-      continue;\r
-    }\r
-\r
-    if( trackPosTest->GetSign() >0 && trackNegTest->GetSign() <0){\r
-      fV0MIs->GetNPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter\r
-      fV0MIs->GetPPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter \r
-    }\r
-\r
-    if( trackPosTest->GetSign() <0 && trackNegTest->GetSign() >0){\r
-      fV0MIs->GetPPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter\r
-      fV0MIs->GetNPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter\r
-    }\r
-\r
-    fV0MIs->GetPxPyPz(mm[0],mm[1],mm[2]); //reconstructed cartesian momentum components of mother\r
-\r
-    TVector3 vecN(mn[0],mn[1],mn[2]);\r
-    TVector3 vecP(mp[0],mp[1],mp[2]);\r
-    TVector3 vecM(mm[0],mm[1],mm[2]);\r
-    \r
-    Double_t thetaP = acos((vecP * vecM)/(vecP.Mag() * vecM.Mag()));\r
-    Double_t thetaN = acos((vecN * vecM)/(vecN.Mag() * vecM.Mag()));\r
-    \r
-    Double_t alfa = ((vecP.Mag())*cos(thetaP)-(vecN.Mag())*cos(thetaN))/\r
-      ((vecP.Mag())*cos(thetaP)+(vecN.Mag())*cos(thetaN)) ;\r
-    Double_t qt = vecP.Mag()*sin(thetaP);\r
-\r
-    fHistArmenterosPodolanski->Fill(alfa,qt);    //Armenteros-Podolanski calculation\r
-\r
-    TLorentzVector k0;\r
-    TLorentzVector k0daugh1;\r
-    TLorentzVector k0daugh2;\r
-    TLorentzVector proton;\r
-    TLorentzVector pq;\r
-\r
-    k0daugh1.SetXYZM(mn[0],mn[1],mn[2],cPionMass);\r
-    k0daugh2.SetXYZM(mp[0],mp[1],mp[2],cPionMass);\r
-    k0=k0daugh1+k0daugh2;\r
-\r
-    fV0MIs->ChangeMassHypothesis(3122);\r
-    lInvMassLambda = fV0MIs->GetEffMass();\r
-    lPtLambda = fV0MIs->Pt();\r
-    lPzLambda = fV0MIs->Pz();\r
-    lPLambda = fV0MIs->P();\r
-\r
-    trkArray->AddAt(trackPosTest,0);\r
-    trkArray->AddAt(trackNegTest,1);\r
-\r
-    vertexer->SetVtxStart(esdVer1);\r
-    AliESDVertex *decayVertex = (AliESDVertex*)vertexer->VertexForSelectedESDTracks(trkArray);\r
-\r
-    dd[0]=fESD->GetPrimaryVertexSPD()->GetX()-decayVertex->GetX();\r
-    dd[1]=fESD->GetPrimaryVertexSPD()->GetY()-decayVertex->GetY();\r
-    dd[2]=fESD->GetPrimaryVertexSPD()->GetZ()-decayVertex->GetZ();\r
-    \r
-    decayLength=sqrt(dd[0]*dd[0]+dd[1]*dd[1]+dd[2]*dd[2]);\r
-\r
+  const AliESDVertex *vertex = fESD->GetPrimaryVertexTracks();
+  if (vertex->GetNContributors()<1) 
+    {
+      // SPD vertex
+      vertex = fESD->GetPrimaryVertexSPD();
+       if(vertex->GetNContributors()<1) {
+        return;
+       }  
+    }
+  if (TMath::Abs(vertex->GetZv()) > 10) return;
+  
+  Int_t centrality = -5;
+  Double_t centrPerc = -5;
+
+  if (fESD->GetEventSpecie() == 4) 
+    { // PbPb
+      AliCentrality *esdCentrality = fESD->GetCentrality();
+      centrality = esdCentrality->GetCentralityClass10("V0M"); // centrality percentile determined with V0
+      centrPerc = esdCentrality->GetCentralityPercentile("V0M");
+//      if (centrality < 0. || centrality > 8.) return; //0 bis 80 %
+      if (centrality > 8) return; //0 bis 80 %
+      //  cout<<"Centrality: "<< centrality << endl;
+    }
+
+
+  fHistCentrality->Fill(centrality);   
+
+  //*****************//  
+  //*   Centrality  *//
+  //*****************//
+  //  Float_t percentile=centrality->GetCentralityPercentile("V0M");
+
+  Bool_t isSelectedCentral = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kCentral);
+  Bool_t isSelectedSemiCentral = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kSemiCentral);
+  Bool_t isSelectedMB = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & AliVEvent::kMB);
+  Int_t triggertype = 17;
+
+  if(isSelectedCentral){
+    fHistTriggerStat->Fill(1);
+    triggertype=1;
+  }
+
+  if(isSelectedSemiCentral){
+    fHistTriggerStat->Fill(2);
+    triggertype=2;
+  }
+
+  if(isSelectedMB){
+    fHistTriggerStat->Fill(0);
+    triggertype=3;
+  }
+
+  //  if(isSelectedCentral || isSelectedSemiCentral || isSelectedMB){
+
+  //*******************************
+
+  Int_t runNumber = 0;
+  //  itrk = 0;
+  runNumber = fESD->GetRunNumber();
+/*  
+  if (!fPIDtpcESD) fPIDtpcESD = ((AliESDInputHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->GetESDpid();
+  if (!fPIDtpcESD) {
+    fPIDtpcESD = new AliESDpid(); // HACK FOR MC PBPB --> PLEASE REMOVE AS SOON AS POSSIBLE
+    fPIDtpcESD->GetTPCResponse().SetBetheBlochParameters(1.28778e+00/50., 3.13539e+01, TMath::Exp(-3.16327e+01), 1.87901e+00, 6.41583e+00);
+  }
+*/
+  Double_t pionK=1;
+  Double_t pK=1;
+  TObjArray* listCrossV0   = fESDCutsV0->GetAcceptedV0s(fESD);  
+  Int_t nGoodV0s      = listCrossV0->GetEntries();
+
+  const AliESDVertex *esdVer = fESD->GetPrimaryVertex();
+  AliESDVertex *esdVer1 = new AliESDVertex(*esdVer);
+   
+  AliVertexerTracks *vertexer = new AliVertexerTracks(fESD->GetMagneticField());
+  TObjArray *trkArray = new TObjArray(2);
+  AliVertexerTracks *vertexer1 = new AliVertexerTracks(fESD->GetMagneticField());
+  TObjArray *trkArray1 = new TObjArray(2);
+
+  AliKFParticle::SetField(fESD->GetMagneticField());
+       
+  AliKFVertex primVtx(*(fESD->GetPrimaryVertex()));
+
+  Int_t refMultTpc = AliESDtrackCuts::GetReferenceMultiplicity(fESD, kTRUE);
+  //cout<<"Multiplicity: "<< refMultTpc << endl;
+  fHistMultiplicity->Fill(refMultTpc);
+  
+  Double_t mn[3] = {0,0,0};
+  Double_t mp[3] = {0,0,0}; 
+  Double_t mm[3] = {0,0,0};
+  Double_t dd[3] = {0,0,0};
+  Double_t dd1[3] = {0,0,0};
+  const Double_t cProtonMass=TDatabasePDG::Instance()->GetParticle(2212)->Mass();
+  const Double_t cPionMass=TDatabasePDG::Instance()->GetParticle(211)->Mass();
+  const Double_t cElectronMass=TDatabasePDG::Instance()->GetParticle(11)->Mass();
+  const Double_t cLambdaMass=TDatabasePDG::Instance()->GetParticle(3122)->Mass();
+  Double_t decayLength=0;
+  Double_t decayLengthH=0;
+
+  //V0 Loop for Lambda and Anti-Lambda
+  for(Int_t iV0MI = 0; iV0MI < nGoodV0s ; iV0MI++) {
+    AliESDv0 * fV0MIs = fESD->GetV0(iV0MI);
+    Int_t    lOnFlyStatus = 0;
+
+    lOnFlyStatus = fV0MIs->GetOnFlyStatus();
+    Double_t lInvMassLambda=0;
+    Double_t lInvMassLambdaPi=0;
+    Double_t lPtLambda=0;
+    Double_t lPzLambda=0;
+    Double_t lPLambda=0;
+    Int_t onl=0;
+    Int_t offl=0;
+
+    TLorentzVector posE;
+    TLorentzVector negE;
+    TLorentzVector photon;
+
+    if (lOnFlyStatus){
+      onl=1; 
+      //      return;
+    }
+    if (!lOnFlyStatus){
+      offl=1;
+      //return;
+    }
+
+    //    fHistMultiplicity->Fill(refMultTpc); 
+    fHistCentralityAC->Fill(centrality);
+    fHistCheck->Fill(offl,onl);
+
+    AliESDtrack* trackPosTest = fESD->GetTrack(fV0MIs->GetPindex());
+    AliESDtrack* trackNegTest = fESD->GetTrack(fV0MIs->GetNindex());
+
+    //    if (!
+    if (!fEsdTrackCuts->AcceptTrack(trackPosTest)) continue;
+    if (!fESDtrackCutsV0->AcceptTrack(trackNegTest)) continue;
+
+
+      //PID via specific energy loss in the TPC
+      //define the arrays for the Bethe-Bloch-Parameters
+      Double_t parProton[5]   = {0,0,0,0,0};
+
+      if(runNumber < 166500) //LHC10h
+       {
+         parProton[0]   = 1.45802; // ALEPH parameters for protons (pass2)
+         parProton[1]   = 27.4992;
+         parProton[2]   = 4.00313e-15;
+         parProton[3]   = 2.48485;
+         parProton[4]   = 8.31768; 
+       }
+      
+      if(runNumber > 166500) //LHC11h
+       { 
+         parProton[0]   = 1.11243; // ALEPH parameters for protons (pass2)
+         parProton[1]   = 26.1144;
+         parProton[2]   = 4.00313e-15;
+         parProton[3]   = 2.72969;
+         parProton[4]   = 9.15038; 
+       }
+      //Get the total momentum for each track at the inner readout of the TPC
+      Double_t ptotN = trackNegTest->GetInnerParam()->GetP();
+      Double_t ptotP = trackPosTest->GetInnerParam()->GetP();
+
+      //define expected signals for the various species
+      Double_t expSignalPionP = 0;
+      Double_t expSignalPionN = 0;
+      Double_t expSignalProtonN = 0;
+      Double_t expSignalProtonP = 0;
+
+      //for data
+      if(!HasMC())
+       {
+         expSignalProtonN = AliExternalTrackParam::BetheBlochAleph(ptotN/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);
+         expSignalProtonP = AliExternalTrackParam::BetheBlochAleph(ptotP/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);
+       }
+    
+      //for MC
+      if(HasMC())
+       {
+         expSignalPionP = 0.7*AliExternalTrackParam::BetheBlochAleph(ptotP/(pionMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]); 
+         expSignalPionN = 0.7*AliExternalTrackParam::BetheBlochAleph(ptotN/(pionMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]); 
+
+         expSignalProtonN = 0.65*AliExternalTrackParam::BetheBlochAleph(ptotN/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);
+         expSignalProtonP = 0.65*AliExternalTrackParam::BetheBlochAleph(ptotP/(protonMass),parProton[0],parProton[1],parProton[2],parProton[3],parProton[4]);   
+       }
+
+      // PID cut on the nuclei (proton, deuteron, triton, helium3)
+      Bool_t corrParticle = kFALSE;
+
+      Bool_t posProton = kFALSE;
+      //data
+      if(!HasMC())
+       {
+         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackPosTest, AliPID::kProton)) < 3)
+           {
+             posProton = kTRUE;
+             corrParticle = kTRUE;
+           }
+       }
+      //MC
+      if(HasMC())
+       {
+         if(//trackPosTest->GetTPCsignal() < 1200 && 
+            TMath::Abs(trackPosTest->GetTPCsignal() - expSignalProtonP)/expSignalProtonP < 0.4)
+           {
+             posProton = kTRUE;
+             corrParticle = kTRUE;
+           }
+       }
+
+      Bool_t negProton = kFALSE;
+      //data
+      if(!HasMC())
+       {
+         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackNegTest, AliPID::kProton)) < 3)
+           {
+             negProton = kTRUE;
+             corrParticle = kTRUE;
+           }
+       }
+     //MC
+      if(HasMC())
+       {
+         if(//trackNegTest->GetTPCsignal() < 1200 && 
+            TMath::Abs(trackNegTest->GetTPCsignal() - expSignalProtonN)/expSignalProtonN < 0.4)
+           {
+             negProton = kTRUE;
+             corrParticle = kTRUE;
+           }
+       }
+
+     //PID cut for pions
+      //data: 3sigma cut on the pions
+
+      Bool_t negPion = kFALSE;     
+      Bool_t posPion = kFALSE;
+
+      if (!HasMC())
+       {
+         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackPosTest, AliPID::kPion)) < 4) posPion=kTRUE; //pos daughter has to be a pion
+         if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackNegTest, AliPID::kPion)) < 4) negPion=kTRUE; // negative daughter has to be a pion
+       }
+      
+      //MC: like the nuclei via the specific energyloss in the TPC
+      if (HasMC())
+       {
+         if (TMath::Abs(trackPosTest->GetTPCsignal() - expSignalPionP)/expSignalPionP < 0.4) posPion=kTRUE;
+         if (TMath::Abs(trackNegTest->GetTPCsignal() - expSignalPionN)/expSignalPionN < 0.4) negPion=kTRUE;
+       }
+
+      if (!(posProton==kTRUE)) continue;
+      if (!(negPion==kTRUE)) continue;
+
+    //To avoid ghosts
+
+    if( !(trackPosTest->GetStatus() & AliESDtrack::kTPCrefit)){
+      continue;
+    }
+
+    if( !(trackNegTest->GetStatus() & AliESDtrack::kTPCrefit)){
+      continue;
+    }
+
+    if( trackPosTest->GetSign() >0 && trackNegTest->GetSign() <0){
+      fV0MIs->GetNPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter
+      fV0MIs->GetPPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter 
+    }
+
+    if( trackPosTest->GetSign() <0 && trackNegTest->GetSign() >0){
+      fV0MIs->GetPPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter
+      fV0MIs->GetNPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter
+    }
+
+    fV0MIs->GetPxPyPz(mm[0],mm[1],mm[2]); //reconstructed cartesian momentum components of mother
+
+    TVector3 vecN(mn[0],mn[1],mn[2]);
+    TVector3 vecP(mp[0],mp[1],mp[2]);
+    TVector3 vecM(mm[0],mm[1],mm[2]);
+    
+    Double_t thetaP = acos((vecP * vecM)/(vecP.Mag() * vecM.Mag()));
+    Double_t thetaN = acos((vecN * vecM)/(vecN.Mag() * vecM.Mag()));
+    
+    Double_t alfa = ((vecP.Mag())*cos(thetaP)-(vecN.Mag())*cos(thetaN))/
+      ((vecP.Mag())*cos(thetaP)+(vecN.Mag())*cos(thetaN)) ;
+    Double_t qt = vecP.Mag()*sin(thetaP);
+
+    fHistArmenterosPodolanski->Fill(alfa,qt);    //Armenteros-Podolanski calculation
+
+    TLorentzVector k0;
+    TLorentzVector k0daugh1;
+    TLorentzVector k0daugh2;
+    TLorentzVector proton;
+    TLorentzVector pq;
+
+    k0daugh1.SetXYZM(mn[0],mn[1],mn[2],cPionMass);
+    k0daugh2.SetXYZM(mp[0],mp[1],mp[2],cPionMass);
+    k0=k0daugh1+k0daugh2;
+
+    fV0MIs->ChangeMassHypothesis(3122);
+    lInvMassLambda = fV0MIs->GetEffMass();
+    lPtLambda = fV0MIs->Pt();
+    lPzLambda = fV0MIs->Pz();
+    lPLambda = fV0MIs->P();
+
+    trkArray->AddAt(trackPosTest,0);
+    trkArray->AddAt(trackNegTest,1);
+
+    vertexer->SetVtxStart(esdVer1);
+    AliESDVertex *decayVertex = (AliESDVertex*)vertexer->VertexForSelectedESDTracks(trkArray);
+
+    dd[0]=fESD->GetPrimaryVertexSPD()->GetX()-decayVertex->GetX();
+    dd[1]=fESD->GetPrimaryVertexSPD()->GetY()-decayVertex->GetY();
+    dd[2]=fESD->GetPrimaryVertexSPD()->GetZ()-decayVertex->GetZ();
+    
+    decayLength=sqrt(dd[0]*dd[0]+dd[1]*dd[1]+dd[2]*dd[2]);
+
       if (decayVertex == NULL) cout << "Lambda decay vtx pointer NULL" << endl;
-    if (decayVertex) delete decayVertex;\r
-\r
-    TLorentzVector negPio1;\r
-    TLorentzVector posProt1;        \r
-    TLorentzVector posP;\r
-    TLorentzVector posProt;\r
-    TLorentzVector negK;\r
-    TLorentzVector negPio;\r
-    TLorentzVector negPi;\r
-    TLorentzVector omega;\r
-    TLorentzVector threeSum;\r
-    TLorentzVector fourSum;\r
-    TLorentzVector ppK;\r
-    TLorentzVector posPiK;\r
-    TLorentzVector negPiK;\r
-    TLorentzVector kaon;\r
-    TLorentzVector lambda;\r
-    TLorentzVector lambdaH;\r
-    TLorentzVector hDibaryon;\r
-    TVector3 h;\r
-    TVector3 h1;\r
-\r
-    Int_t mcStatus=0;\r
-\r
-    h.SetXYZ(-dd[0],-dd[1],-dd[2]);\r
-\r
-\r
-    if (offl==1)fHistMassDPi->Fill(lInvMassLambda);\r
-\r
-    if (onl==1){\r
-      fHistMassLPi->Fill(lInvMassLambda);\r
-      \r
-      negE.SetXYZM(mn[0],mn[1],mn[2],cElectronMass);\r
-      posE.SetXYZM(mp[0],mp[1],mp[2],cElectronMass);\r
-      photon=posE+negE;\r
-      \r
-      negPiK.SetXYZM(mn[0],mn[1],mn[2],cPionMass);\r
-      posPiK.SetXYZM(mp[0],mp[1],mp[2],cPionMass);\r
-      kaon=posPiK+negPiK;\r
-      \r
-      negPi.SetXYZM(mn[0],mn[1],mn[2],cPionMass);\r
-      posP.SetXYZM(mp[0],mp[1],mp[2],cProtonMass);\r
-      lambda=negPi+posP;\r
-      lambdaH.SetXYZM(mm[0],mm[1],mm[2],cLambdaMass);\r
-\r
-      if (lInvMassLambda>1.1113&&lInvMassLambda<1.1202){\r
-       \r
-       if (!HasMC()){\r
-         if (qt<-2.21*alfa*alfa+2.945*alfa-0.887) continue;\r
-         if (qt>-2.21*alfa*alfa+2.945*alfa-0.873) continue;\r
-         if (photon.M()<0.005) continue;\r
-         if (kaon.M()>0.495 && kaon.M()<0.500 ) continue;\r
-       }\r
-       \r
-\r
-       fHistMassLambda->Fill(lInvMassLambda);\r
-       //\r
-       Bool_t isCorrectlyAssociatedLambda = kFALSE;\r
-       Bool_t isPartialCorrectlyAssociatedLambda = kFALSE;\r
-       Int_t labelAssociatedH=-1; \r
-       Int_t labelLambda=-1;\r
-       //\r
-       if (HasMC()) {\r
-         Int_t labelPosTest = trackPosTest->GetLabel();\r
-         TParticle *tparticleDaughter = stack->Particle(TMath::Abs(labelPosTest));\r
-         Int_t labelMother = tparticleDaughter->GetFirstMother();\r
-         TParticle *tparticleMother = stack->Particle(TMath::Abs(labelMother));\r
-         \r
-         Int_t labelOma = tparticleMother->GetFirstMother();\r
-         TParticle *tparticleOma = stack->Particle(TMath::Abs(labelOma));\r
-\r
-         if ((tparticleOma->GetPdgCode() < 0) && TMath::Abs(tparticleMother->GetPdgCode())==pdgLambda){// check mother to be Lambda \r
-           Int_t labelFirstDaughter  = tparticleMother->GetDaughter(1);// Proton\r
-           Int_t labelThirdDaughter  = tparticleMother->GetDaughter(0);// Pion\r
-           \r
-           TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));\r
-           TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));\r
-           \r
-           if((tparticleFirstDaughter->GetPdgCode() == pdgProton && tparticleThirdDaughter->GetPdgCode()== pdgPionMinus) || \r
-              (tparticleFirstDaughter->GetPdgCode() == pdgPionMinus && tparticleThirdDaughter->GetPdgCode()== pdgProton)){ //daughter PDGs\r
-             isPartialCorrectlyAssociatedLambda = kTRUE;\r
-             labelLambda=labelMother;\r
-           }\r
-         }\r
-         \r
-         //H-Dibaryon\r
-         if(tparticleOma->GetPdgCode() == pdgHDibaryon){ //check grandmother to be H PDG\r
-           if (TMath::Abs(tparticleMother->GetPdgCode())==pdgLambda){// check mother to be Lambda \r
-             Int_t labelFirstDaughter  = tparticleMother->GetDaughter(1);// Proton\r
-             Int_t labelThirdDaughter  = tparticleMother->GetDaughter(0);// Pion\r
-                       \r
-             TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));\r
-             TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));\r
-             \r
-             if((tparticleFirstDaughter->GetPdgCode() == pdgProton && tparticleThirdDaughter->GetPdgCode()== pdgPionMinus) || \r
-                (tparticleFirstDaughter->GetPdgCode() == pdgPionMinus && tparticleThirdDaughter->GetPdgCode()== pdgProton)){ //daughter PDGs\r
-               isCorrectlyAssociatedLambda = kTRUE;\r
-               labelAssociatedH=labelOma;\r
-               fHistMassLambdaFromH->Fill(lInvMassLambda);\r
-               fHistMassLambdaFromHtLorentz->Fill(lambda.M());\r
-             }\r
-           }\r
-         }\r
-       }\r
-\r
-       fHistProtonPIDLambda->Fill(trackPosTest->GetInnerParam()->GetP(), trackPosTest->GetTPCsignal());\r
-       fHistPionPIDLambda->Fill(trackNegTest->GetInnerParam()->GetP(), trackNegTest->GetTPCsignal());\r
-\r
-       //---------------------------------------------------------\r
-       // Proton track loop\r
-       //---------------------------------------------------------\r
-       fHistArmenterosPodolanskiCut->Fill(alfa,qt);\r
-\r
-       for (Int_t iTracksP = 0; iTracksP < fESD->GetNumberOfTracks(); iTracksP++) {\r
-         AliESDtrack* trackP = dynamic_cast<AliESDtrack*> (fESD->GetTrack(iTracksP));\r
-         if (trackP->GetSign()<0) continue;\r
-         \r
-         if (iTracksP==fV0MIs->GetPindex())continue;\r
-         if (iTracksP==fV0MIs->GetNindex())continue;\r
-               \r
-         if (!fEsdTrackCuts->AcceptTrack(trackP)) continue;\r
-               \r
-         AliKFParticle protonKF( *(trackP), 2212);\r
-\r
-         if (!trackP->GetInnerParam()) continue;\r
-         \r
-         if (HasMC()) {\r
-           pK=0.65;\r
+    if (decayVertex) delete decayVertex;
+
+    TLorentzVector negPio1;
+    TLorentzVector posProt1;        
+    TLorentzVector posP;
+    TLorentzVector posProt;
+    TLorentzVector negK;
+    TLorentzVector negPio;
+    TLorentzVector negPi;
+    TLorentzVector omega;
+    TLorentzVector threeSum;
+    TLorentzVector fourSum;
+    TLorentzVector ppK;
+    TLorentzVector posPiK;
+    TLorentzVector negPiK;
+    TLorentzVector kaon;
+    TLorentzVector lambda;
+    TLorentzVector lambdaH;
+    TLorentzVector hDibaryon;
+    TVector3 h;
+    TVector3 h1;
+
+    Int_t mcStatus=0;
+
+    h.SetXYZ(-dd[0],-dd[1],-dd[2]);
+
+
+    if (offl==1)fHistMassDPi->Fill(lInvMassLambda);
+
+    if (onl==1){
+      fHistMassLPi->Fill(lInvMassLambda);
+      
+      negE.SetXYZM(mn[0],mn[1],mn[2],cElectronMass);
+      posE.SetXYZM(mp[0],mp[1],mp[2],cElectronMass);
+      photon=posE+negE;
+      
+      negPiK.SetXYZM(mn[0],mn[1],mn[2],cPionMass);
+      posPiK.SetXYZM(mp[0],mp[1],mp[2],cPionMass);
+      kaon=posPiK+negPiK;
+      
+      negPi.SetXYZM(mn[0],mn[1],mn[2],cPionMass);
+      posP.SetXYZM(mp[0],mp[1],mp[2],cProtonMass);
+      lambda=negPi+posP;
+      lambdaH.SetXYZM(mm[0],mm[1],mm[2],cLambdaMass);
+
+      if (lInvMassLambda>1.1113&&lInvMassLambda<1.1202){
+       
+       if (!HasMC()){
+         if (qt<-2.21*alfa*alfa+2.945*alfa-0.887) continue;
+         if (qt>-2.21*alfa*alfa+2.945*alfa-0.873) continue;
+         if (photon.M()<0.005) continue;
+         if (kaon.M()>0.495 && kaon.M()<0.500 ) continue;
+       }
+       
+
+       fHistMassLambda->Fill(lInvMassLambda);
+       //
+       Bool_t isCorrectlyAssociatedLambda = kFALSE;
+       Bool_t isPartialCorrectlyAssociatedLambda = kFALSE;
+       Int_t labelAssociatedH=-1; 
+       Int_t labelLambda=-1;
+       //
+       if (HasMC()) {
+         Int_t labelPosTest = trackPosTest->GetLabel();
+         TParticle *tparticleDaughter = stack->Particle(TMath::Abs(labelPosTest));
+         Int_t labelMother = tparticleDaughter->GetFirstMother();
+         TParticle *tparticleMother = stack->Particle(TMath::Abs(labelMother));
+         
+         Int_t labelOma = tparticleMother->GetFirstMother();
+         TParticle *tparticleOma = stack->Particle(TMath::Abs(labelOma));
+
+         if ((tparticleOma->GetPdgCode() < 0) && TMath::Abs(tparticleMother->GetPdgCode())==pdgLambda){// check mother to be Lambda 
+           Int_t labelFirstDaughter  = tparticleMother->GetDaughter(1);// Proton
+           Int_t labelThirdDaughter  = tparticleMother->GetDaughter(0);// Pion
+           
+           TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));
+           TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));
+           
+           if((tparticleFirstDaughter->GetPdgCode() == pdgProton && tparticleThirdDaughter->GetPdgCode()== pdgPionMinus) || 
+              (tparticleFirstDaughter->GetPdgCode() == pdgPionMinus && tparticleThirdDaughter->GetPdgCode()== pdgProton)){ //daughter PDGs
+             isPartialCorrectlyAssociatedLambda = kTRUE;
+             labelLambda=labelMother;
+           }
+         }
+         
+         //H-Dibaryon
+         if(tparticleOma->GetPdgCode() == pdgHDibaryon){ //check grandmother to be H PDG
+           if (TMath::Abs(tparticleMother->GetPdgCode())==pdgLambda){// check mother to be Lambda 
+             Int_t labelFirstDaughter  = tparticleMother->GetDaughter(1);// Proton
+             Int_t labelThirdDaughter  = tparticleMother->GetDaughter(0);// Pion
+                       
+             TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));
+             TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));
+             
+             if((tparticleFirstDaughter->GetPdgCode() == pdgProton && tparticleThirdDaughter->GetPdgCode()== pdgPionMinus) || 
+                (tparticleFirstDaughter->GetPdgCode() == pdgPionMinus && tparticleThirdDaughter->GetPdgCode()== pdgProton)){ //daughter PDGs
+               isCorrectlyAssociatedLambda = kTRUE;
+               labelAssociatedH=labelOma;
+               fHistMassLambdaFromH->Fill(lInvMassLambda);
+               fHistMassLambdaFromHtLorentz->Fill(lambda.M());
+             }
+           }
+         }
+       }
+
+       fHistProtonPIDLambda->Fill(trackPosTest->GetInnerParam()->GetP(), trackPosTest->GetTPCsignal());
+       fHistPionPIDLambda->Fill(trackNegTest->GetInnerParam()->GetP(), trackNegTest->GetTPCsignal());
+
+       //---------------------------------------------------------
+       // Proton track loop
+       //---------------------------------------------------------
+       fHistArmenterosPodolanskiCut->Fill(alfa,qt);
+
+       for (Int_t iTracksP = 0; iTracksP < fESD->GetNumberOfTracks(); iTracksP++) {
+         AliESDtrack* trackP = dynamic_cast<AliESDtrack*> (fESD->GetTrack(iTracksP));
+         if (trackP->GetSign()<0) continue;
+         
+         if (iTracksP==fV0MIs->GetPindex())continue;
+         if (iTracksP==fV0MIs->GetNindex())continue;
+               
+         if (!fEsdTrackCuts->AcceptTrack(trackP)) continue;
+               
+         AliKFParticle protonKF( *(trackP), 2212);
+
+         if (!trackP->GetInnerParam()) continue;
+         
+         if (HasMC()) {
+           pK=0.65;
          }
 
-         if (!HasMC()){\r
-           if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackP, AliPID::kProton)) > 3) continue;\r
-         }             \r
-         \r
-         fHistProtonPIDBb->Fill(trackP->GetInnerParam()->GetP(), trackP->GetTPCsignal());\r
-         \r
-         posProt.SetXYZM(trackP->Px(),trackP->Py(),trackP->Pz(),cProtonMass);\r
-\r
-         //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\r
-         //Pion Track loop!!!!!!!!!!!!!!!!!!!!!\r
-         //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\r
-         for (Int_t iTracksN = iTracksP+1; iTracksN < fESD->GetNumberOfTracks(); iTracksN++) {\r
-           AliESDtrack* trackN = dynamic_cast<AliESDtrack*> (fESD->GetTrack(iTracksN));\r
-\r
-           if (iTracksN==fV0MIs->GetPindex())continue;\r
-           if (iTracksN==fV0MIs->GetNindex())continue;\r
-           if (trackN->GetSign()>0) continue;\r
-         \r
-           if (!fEsdTrackCuts->AcceptTrack(trackN)) continue;\r
-           if (!fESDtrackCutsV0->AcceptTrack(trackN)) continue;\r
+         if (!HasMC()){
+           if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackP, AliPID::kProton)) > 3) continue;
+         }             
+         
+         fHistProtonPIDBb->Fill(trackP->GetInnerParam()->GetP(), trackP->GetTPCsignal());
+         
+         posProt.SetXYZM(trackP->Px(),trackP->Py(),trackP->Pz(),cProtonMass);
+
+         //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+         //Pion Track loop!!!!!!!!!!!!!!!!!!!!!
+         //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+         for (Int_t iTracksN = iTracksP+1; iTracksN < fESD->GetNumberOfTracks(); iTracksN++) {
+           AliESDtrack* trackN = dynamic_cast<AliESDtrack*> (fESD->GetTrack(iTracksN));
+
+           if (iTracksN==fV0MIs->GetPindex())continue;
+           if (iTracksN==fV0MIs->GetNindex())continue;
+           if (trackN->GetSign()>0) continue;
+         
+           if (!fEsdTrackCuts->AcceptTrack(trackN)) continue;
+           if (!fESDtrackCutsV0->AcceptTrack(trackN)) continue;
 
            Double_t bz = fESD->GetMagneticField();
           
@@ -1265,323 +1265,323 @@ void AliAnalysisTaskHdibaryonLPpi::UserExec(Option_t *)
            Double_t dca = trackN->GetDCA(trackP,bz,xthiss,xpp);
            if (dca>0.5) continue;
 
-           negPi.SetXYZM(mn[0],mn[1],mn[2],cPionMass);\r
-           posP.SetXYZM(mp[0],mp[1],mp[2],cProtonMass);\r
-           negPio.SetXYZM(trackN->Px(),trackN->Py(),trackN->Pz(),cPionMass);\r
-             \r
-           threeSum=negPi+posP+negPio;\r
-           lInvMassLambdaPi=threeSum.M();         \r
-           \r
-           fHistDC->Fill(decayLength*lInvMassLambda/lPLambda);\r
-           \r
-           AliKFParticle posPionKF( *(trackN) ,-211);\r
-           \r
-           if (!trackN->GetInnerParam()) continue;\r
-           if (HasMC()) {\r
-             pionK=0.7;\r
-           }\r
-           \r
-         //  if (!HasMC()){\r
-             if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackN, AliPID::kPion)) > 3) continue;\r
-         //  }\r
-           fHistPionPIDBb->Fill(trackN->GetInnerParam()->GetP(), trackN->GetTPCsignal());\r
-\r
-           trkArray1->AddAt(trackP,0);\r
-           trkArray1->AddAt(trackN,1);\r
-           \r
-           vertexer1->SetVtxStart(esdVer1);\r
-           AliESDVertex *decayVertex1 = (AliESDVertex*)vertexer1->VertexForSelectedESDTracks(trkArray1);\r
-           \r
-           dd1[0]=fESD->GetPrimaryVertexSPD()->GetX()-decayVertex1->GetX();\r
-           dd1[1]=fESD->GetPrimaryVertexSPD()->GetY()-decayVertex1->GetY();\r
-           dd1[2]=fESD->GetPrimaryVertexSPD()->GetZ()-decayVertex1->GetZ();\r
-           \r
-           decayLengthH=sqrt(dd1[0]*dd1[0]+dd1[1]*dd1[1]+dd1[2]*dd1[2]);\r
-\r
-             //            Double_t bz = fESD->GetMagneticField();\r
-           \r
-           trackP->PropagateToDCA(decayVertex1, bz, 10);\r
-           trackN->PropagateToDCA(decayVertex1, bz, 10);\r
-\r
-             //            Double_t xthiss(0.0);\r
-             //            Double_t xpp(0.0);\r
-             //            Double_t dca = trackN->GetDCA(trackP,bz,xthiss,xpp);\r
-\r
+           negPi.SetXYZM(mn[0],mn[1],mn[2],cPionMass);
+           posP.SetXYZM(mp[0],mp[1],mp[2],cProtonMass);
+           negPio.SetXYZM(trackN->Px(),trackN->Py(),trackN->Pz(),cPionMass);
+             
+           threeSum=negPi+posP+negPio;
+           lInvMassLambdaPi=threeSum.M();         
+           
+           fHistDC->Fill(decayLength*lInvMassLambda/lPLambda);
+           
+           AliKFParticle posPionKF( *(trackN) ,-211);
+           
+           if (!trackN->GetInnerParam()) continue;
+           if (HasMC()) {
+             pionK=0.7;
+           }
+           
+         //  if (!HasMC()){
+             if (TMath::Abs(fPIDtpcESD->NumberOfSigmasTPC(trackN, AliPID::kPion)) > 3) continue;
+         //  }
+           fHistPionPIDBb->Fill(trackN->GetInnerParam()->GetP(), trackN->GetTPCsignal());
+
+           trkArray1->AddAt(trackP,0);
+           trkArray1->AddAt(trackN,1);
+           
+           vertexer1->SetVtxStart(esdVer1);
+           AliESDVertex *decayVertex1 = (AliESDVertex*)vertexer1->VertexForSelectedESDTracks(trkArray1);
+           
+           dd1[0]=fESD->GetPrimaryVertexSPD()->GetX()-decayVertex1->GetX();
+           dd1[1]=fESD->GetPrimaryVertexSPD()->GetY()-decayVertex1->GetY();
+           dd1[2]=fESD->GetPrimaryVertexSPD()->GetZ()-decayVertex1->GetZ();
+           
+           decayLengthH=sqrt(dd1[0]*dd1[0]+dd1[1]*dd1[1]+dd1[2]*dd1[2]);
+
+             //            Double_t bz = fESD->GetMagneticField();
+           
+           trackP->PropagateToDCA(decayVertex1, bz, 10);
+           trackN->PropagateToDCA(decayVertex1, bz, 10);
+
+             //            Double_t xthiss(0.0);
+             //            Double_t xpp(0.0);
+             //            Double_t dca = trackN->GetDCA(trackP,bz,xthiss,xpp);
+
              if (decayVertex1 == NULL) cout << "Secondary decay vtx pointer NULL" << endl;
-           if (decayVertex1) delete decayVertex1;\r
-           h1.SetXYZ(-dd1[0],-dd1[1],-dd1[2]);\r
-\r
-           //      if (dca>1) continue;\r
-             //            if (dca>0.1) continue;\r
-\r
-           fourSum=threeSum+posProt;\r
-\r
-           posProt1.SetXYZM(trackP->Px(),trackP->Py(),trackP->Pz(),cProtonMass);\r
-           negPio1.SetXYZM(trackN->Px(),trackN->Py(),trackN->Pz(),cPionMass);\r
-           hDibaryon=lambdaH+posProt1+negPio1;\r
-           Double_t hPointingAngle = hDibaryon.Angle(h);\r
-           Double_t pointingAngleH = hDibaryon.Angle(h1);\r
-           Double_t decayAngleH = h.Angle(h1);\r
-           TVector3 vecDist(dd[0]-dd1[0],dd[1]-dd1[1],dd[2]-dd1[2]);\r
-           fHistMassLambdaPPi->Fill(hDibaryon.M());\r
-           fHistHPointingAngle->Fill(pointingAngleH);\r
-\r
-            fHistMassHcentMult->Fill(hDibaryon.M(),triggertype,refMultTpc);\r
-\r
-           Double_t rapidity = hDibaryon.Rapidity();\r
-           if(rapidity > 1.0 || rapidity < -1.0) continue;\r
-\r
-             //Double_t vec[16]={hDibaryon.M(), lInvMassLambda, pointingAngleH, alfa, qt, lPtLambda, hDibaryon.Pt(), posPionKF.GetDistanceFromVertex(primVtx), protonKF.GetDistanceFromVertex(primVtx), dca, protonKF.GetDistanceFromVertex(posPionKF), TMath::Cos(pointingAngleH), centrPerc, refMultTpc, mcStatus};\r
-             //fHistNdim->Fill(vec);\r
-\r
-           fHistRap->Fill(rapidity);\r
-           //if (pointingAngleH > 0.1) continue;\r
-           if (pointingAngleH > 0.05) continue;\r
-\r
-           ///////////////////////////\r
-           //MC part for Associated H\r
-           ///////////////////////////\r
-\r
-           if (HasMC() && isCorrectlyAssociatedLambda) {\r
-             Int_t labelP = trackP->GetLabel();\r
-             TParticle *tparticleDaughter = stack->Particle(TMath::Abs(labelP));\r
-             Int_t labelMother = tparticleDaughter->GetFirstMother();\r
-             TParticle *tparticleMother = stack->Particle(TMath::Abs(labelMother));\r
-             \r
-             Int_t labelProton = trackP->GetLabel();\r
-             Int_t labelPion = trackN->GetLabel();\r
-             \r
-             //H-Dibaryon\r
-             if(tparticleMother->GetPdgCode() == pdgHDibaryon && labelAssociatedH==labelMother){ //check mother PDG\r
-               Int_t labelFirstDaughter  = tparticleMother->GetDaughter(0);\r
-               Int_t labelThirdDaughter  = tparticleMother->GetDaughter(1);\r
-               Int_t labelSecondDaughter = labelFirstDaughter +1;\r
-\r
-               TParticle *tparticleSecondDaughter = stack->Particle(TMath::Abs(labelSecondDaughter));\r
-               TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));\r
-               \r
-               TLorentzVector ppi;\r
-               TLorentzVector lpi;\r
-               TLorentzVector lP;\r
-               \r
-               ppi=posProt+negPio;\r
-               lpi=lambdaH+negPio;\r
-               lP=lambdaH+posProt;\r
-\r
-               if((tparticleThirdDaughter->GetPdgCode() == pdgPionMinus && labelPion==labelThirdDaughter)&&(tparticleSecondDaughter->GetPdgCode() == pdgProton||tparticleSecondDaughter->GetPdgCode() == pdgPionMinus) && labelProton==labelSecondDaughter) fHistMassPpi->Fill(ppi.M());\r
-               \r
-               if(tparticleThirdDaughter->GetPdgCode() == pdgPionMinus && labelPion==labelThirdDaughter) fHistMassLpi->Fill(lpi.M());\r
-               \r
-               if(tparticleSecondDaughter->GetPdgCode() == pdgProton && labelProton==labelSecondDaughter) fHistMassLP->Fill(lP.M());\r
-               \r
-               if(tparticleSecondDaughter->GetPdgCode() == pdgProton && labelProton==labelSecondDaughter){//check second daughter PDG\r
-                 if(tparticleThirdDaughter->GetPdgCode() == pdgPionMinus && labelPion==labelThirdDaughter){//check second daughter PDG\r
-\r
-                   fHistHDibaryonInvaMassAso->Fill(hDibaryon.M()); \r
-                   \r
-                   Double_t distance01=vecDist.Mag();\r
-                   fHistMCdcaPvtxDvtx->Fill(decayLengthH);\r
-                   fHistMCdcaPvtxLvtx->Fill(decayLength);\r
-                   fHistMCdcaDvtxLvtx->Fill(distance01);\r
-                   fHistMCangleLH->Fill(hPointingAngle);\r
-                   fHistMCdecayAngle->Fill(decayAngleH);\r
-                   fHistMCpointingAngle->Fill(pointingAngleH);\r
-                   fHistMCap->Fill(alfa,qt);\r
-\r
-                   fHistHilf1->Fill(posPionKF.GetDistanceFromVertex(primVtx));\r
-                   fHistHilf2->Fill(protonKF.GetDistanceFromVertex(primVtx));\r
-                   fHistHilf3->Fill(protonKF.GetDistanceFromVertex(posPionKF));\r
-                   fHistHilf6->Fill(dca);\r
-                   fHistPtvsYAso->Fill(hDibaryon.Pt(),hDibaryon.Rapidity());\r
-                   fHistPtvsEtaAso->Fill(hDibaryon.Pt(),hDibaryon.Eta());\r
-                   mcStatus=1;\r
-                 }//end check for third daughter PDG\r
-               }//end check second daughter PDG\r
-             }//end H-Dibaryon\r
-           }//end MC\r
-           \r
-           //      cout<<"Trigger: "<<triggertype<<endl;\r
-           fHistMassH->Fill(hDibaryon.M());\r
-           fHistMassHcentMult->Fill(hDibaryon.M(),triggertype,refMultTpc);\r
-           ppK=lambdaH+posProt;\r
-           fHistMassLambdaP->Fill(ppK.M());\r
-\r
-             //fHistNdim = new THnSparseF("fHistNdim","THnS;InvMass, InvMassLambda, pointingAngle, armPoAlpha, armPoQt, pTL, pTH, d0p, d0n, dcaHd, dca, decayL, cosPA, centr, multi, mcinf;InvMassH", 16,binsD01,xminD01,xmaxD01);\r
-\r
-           Double_t vec[16]={hDibaryon.M(), lInvMassLambda, pointingAngleH, alfa, qt, lPtLambda, hDibaryon.Pt(), posPionKF.GetDistanceFromVertex(primVtx), protonKF.GetDistanceFromVertex(primVtx), dca, protonKF.GetDistanceFromVertex(posPionKF), TMath::Cos(pointingAngleH), centrPerc, static_cast<Double_t>(refMultTpc), static_cast<Double_t>(mcStatus)};\r
-                     fHistNdim->Fill(vec);\r
-\r
-         }\r
-       }\r
-      }\r
-    }    \r
-  }\r
-  \r
-  //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\r
-  //Pure MC Part!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\r
-  //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\r
-\r
-  // Monte Carlo for genenerated particles\r
-  if (HasMC()) //MC loop  \r
-    {\r
-\r
-      Int_t stackN = 0;\r
-\r
-      Double_t momentumPionGen[3]={0,0,0};\r
-      Double_t momentumNucleonGen[3]={0,0,0};\r
-      Double_t momentumLambdaGen[3]={0,0,0};\r
-\r
-      Double_t energyPionGen = 0;\r
-      Double_t energyNucleonGen = 0;\r
-      Double_t energyLambdaGen = 0;\r
-\r
-      Double_t transversMomentumMotherGen = 0;\r
-      Double_t longitudinalMomentumMotherGen = 0;\r
-      Double_t totalEnergyMotherGen = 0;\r
-      \r
-      Double_t rapidityGen = 2;\r
-\r
-      for(stackN = 0; stackN < stack->GetNtrack(); stackN++) //loop over stack\r
-       {\r
-\r
-         TParticle *tparticleMother = stack->Particle(stackN);\r
-\r
-         if(tparticleMother->GetPdgCode() == pdgLambda) fHistCount->Fill(16); \r
-\r
-         //H-Dibaryon\r
-         if(tparticleMother->GetPdgCode() == pdgHDibaryon) //check mother PDG\r
-           {\r
-             Int_t labelFirstDaughter  = tparticleMother->GetDaughter(0);\r
-             Int_t labelThirdDaughter  = tparticleMother->GetDaughter(1);\r
-             Int_t labelSecondDaughter = labelFirstDaughter +1;\r
-\r
-             TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));\r
-             TParticle *tparticleSecondDaughter = stack->Particle(TMath::Abs(labelSecondDaughter));\r
-             TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));\r
-            \r
-             if(tparticleFirstDaughter->GetPdgCode() == pdgLambda) //check first daughter PDG\r
-               {\r
-                 if(tparticleSecondDaughter->GetPdgCode() == pdgProton)//check second daughter PDG\r
-                   {\r
-                     if(tparticleThirdDaughter->GetPdgCode() == pdgPionMinus)//check second daughter PDG\r
-                       {                                     \r
-                         momentumLambdaGen[0] = tparticleFirstDaughter->Px();\r
-                         momentumLambdaGen[1] = tparticleFirstDaughter->Py();\r
-                         momentumLambdaGen[2] = tparticleFirstDaughter->Pz();\r
-               \r
-                         momentumNucleonGen[0] = tparticleSecondDaughter->Px();\r
-                         momentumNucleonGen[1] = tparticleSecondDaughter->Py();\r
-                         momentumNucleonGen[2] = tparticleSecondDaughter->Pz();\r
-\r
-                         momentumPionGen[0] = tparticleThirdDaughter->Px();\r
-                         momentumPionGen[1] = tparticleThirdDaughter->Py();\r
-                         momentumPionGen[2] = tparticleThirdDaughter->Pz();\r
-\r
-                         TLorentzVector lorentzVectorLambda;\r
-                         TLorentzVector lorentzVectorProton;\r
-                         TLorentzVector lorentzVectorPion;\r
-                         TLorentzVector lorentzVectorHDibaryon;\r
-                         \r
-                         lorentzVectorLambda.SetXYZM(momentumLambdaGen[0],momentumLambdaGen[1],momentumLambdaGen[2],1.115);\r
-                         lorentzVectorProton.SetXYZM(momentumNucleonGen[0],momentumNucleonGen[1],momentumNucleonGen[2],protonMass);\r
-                         lorentzVectorPion.SetXYZM(momentumPionGen[0],momentumPionGen[1],momentumPionGen[2],pionMass);\r
-                         \r
-                         lorentzVectorHDibaryon = lorentzVectorLambda + lorentzVectorProton + lorentzVectorPion;\r
-                         rapidityGen=lorentzVectorHDibaryon.Rapidity();\r
-                         transversMomentumMotherGen = lorentzVectorHDibaryon.Pt();\r
-                         longitudinalMomentumMotherGen = lorentzVectorHDibaryon.Pz();\r
-                         totalEnergyMotherGen = lorentzVectorHDibaryon.Energy();\r
-\r
-                         if(rapidityGen > 1.0 || rapidityGen < -1 ) continue;\r
-                         //lorentzVectorLambda\r
-                         fHistHDibaryonInvaMassGen->Fill(lorentzVectorHDibaryon.M()); \r
-                         if (lorentzVectorLambda.Rapidity()  > 1.0 || lorentzVectorLambda.Rapidity() < -1) continue;\r
-                         if (lorentzVectorProton.Rapidity()  > 1.0 ||      lorentzVectorProton.Rapidity() < -1) continue;\r
-\r
-                         if (lorentzVectorPion.Rapidity()  > 1.0 ||  lorentzVectorPion.Rapidity() < -1) continue;\r
-                         fHistHDibaryonInvaMassGenRes->Fill(lorentzVectorHDibaryon.M());\r
-                         fHistPtvsEtaGen->Fill(lorentzVectorHDibaryon.Pt(),lorentzVectorHDibaryon.Eta());\r
-                         fHistPtvsYGen->Fill(lorentzVectorHDibaryon.Pt(),lorentzVectorHDibaryon.Rapidity());\r
-                         fHistPtvsEtaGen->Fill(lorentzVectorHDibaryon.Pt(),lorentzVectorHDibaryon.Eta());\r
-                         fHistCount->Fill(11);\r
-                       }//end of check third daughter PDG\r
-                   }//end of check second daughter PDG\r
-               }//end of check first daughter PDG\r
-           }//end of H-Dibaryon\r
-\r
-         //Anti-H-Dibaryon\r
-         if(tparticleMother->GetPdgCode() == pdgAntiHDibaryon) //check mother PDG\r
-           {\r
-             Int_t labelFirstDaughter  = tparticleMother->GetDaughter(0);\r
-             Int_t labelThirdDaughter  = tparticleMother->GetDaughter(1);\r
-             Int_t labelSecondDaughter = labelFirstDaughter +1;\r
-\r
-             TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));\r
-             TParticle *tparticleSecondDaughter = stack->Particle(TMath::Abs(labelSecondDaughter));\r
-             TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));\r
-\r
-             if(tparticleFirstDaughter->GetPdgCode() == pdgAntiLambda) //check first daughter PDG\r
-               {\r
-                 if(tparticleSecondDaughter->GetPdgCode() == pdgAntiProton)//check second daughter PDG\r
-                   {\r
-                     if(tparticleThirdDaughter->GetPdgCode() == pdgPionPlus)//check second daughter PDG\r
-                       {                \r
-                         momentumLambdaGen[0] = tparticleFirstDaughter->Px();\r
-                         momentumLambdaGen[1] = tparticleFirstDaughter->Py();\r
-                         momentumLambdaGen[2] = tparticleFirstDaughter->Pz();\r
-               \r
-                         momentumNucleonGen[0] = tparticleSecondDaughter->Px();\r
-                         momentumNucleonGen[1] = tparticleSecondDaughter->Py();\r
-                         momentumNucleonGen[2] = tparticleSecondDaughter->Pz();\r
-\r
-                         momentumPionGen[0] = tparticleThirdDaughter->Px();\r
-                         momentumPionGen[1] = tparticleThirdDaughter->Py();\r
-                         momentumPionGen[2] = tparticleThirdDaughter->Pz();\r
-                         \r
-                         energyLambdaGen  = tparticleFirstDaughter->Energy();\r
-                         energyNucleonGen = tparticleSecondDaughter->Energy();\r
-                         energyPionGen    = tparticleThirdDaughter->Energy();\r
-                          \r
-                         TLorentzVector lorentzVectorLambda;\r
-                         TLorentzVector lorentzVectorProton;\r
-                         TLorentzVector lorentzVectorPion;\r
-                         TLorentzVector lorentzVectorHDibaryon;\r
-                         \r
-                         lorentzVectorLambda.SetXYZM(momentumLambdaGen[0],momentumLambdaGen[1],momentumLambdaGen[2],1.115);\r
-                         lorentzVectorProton.SetXYZM(momentumNucleonGen[0],momentumNucleonGen[1],momentumNucleonGen[2],protonMass);\r
-                         lorentzVectorPion.SetXYZM(momentumPionGen[0],momentumPionGen[1],momentumPionGen[2],pionMass);\r
-                         \r
-                         lorentzVectorHDibaryon = lorentzVectorLambda + lorentzVectorProton + lorentzVectorPion;\r
-\r
-                         rapidityGen=lorentzVectorHDibaryon.Rapidity();\r
-                         if(rapidityGen > 1.0 || rapidityGen < -1 ) continue;\r
-                         fHistAntiHDibaryonInvaMassGen->Fill(lorentzVectorHDibaryon.M()); \r
-                       }//end of check third daughter PDG\r
-                   }//end of check second daughter PDG\r
-               }//end of check first daughter PDG\r
-           }//end of Anti-H-Dibaryon\r
-       }      \r
-    }//end MC\r
-\r
-  // Post output data.\r
-  PostData(1,fHistList);\r
-  //PostData(0,fHistList);\r
-\r
+           if (decayVertex1) delete decayVertex1;
+           h1.SetXYZ(-dd1[0],-dd1[1],-dd1[2]);
+
+           //      if (dca>1) continue;
+             //            if (dca>0.1) continue;
+
+           fourSum=threeSum+posProt;
+
+           posProt1.SetXYZM(trackP->Px(),trackP->Py(),trackP->Pz(),cProtonMass);
+           negPio1.SetXYZM(trackN->Px(),trackN->Py(),trackN->Pz(),cPionMass);
+           hDibaryon=lambdaH+posProt1+negPio1;
+           Double_t hPointingAngle = hDibaryon.Angle(h);
+           Double_t pointingAngleH = hDibaryon.Angle(h1);
+           Double_t decayAngleH = h.Angle(h1);
+           TVector3 vecDist(dd[0]-dd1[0],dd[1]-dd1[1],dd[2]-dd1[2]);
+           fHistMassLambdaPPi->Fill(hDibaryon.M());
+           fHistHPointingAngle->Fill(pointingAngleH);
+
+            fHistMassHcentMult->Fill(hDibaryon.M(),triggertype,refMultTpc);
+
+           Double_t rapidity = hDibaryon.Rapidity();
+           if(rapidity > 1.0 || rapidity < -1.0) continue;
+
+             //Double_t vec[16]={hDibaryon.M(), lInvMassLambda, pointingAngleH, alfa, qt, lPtLambda, hDibaryon.Pt(), posPionKF.GetDistanceFromVertex(primVtx), protonKF.GetDistanceFromVertex(primVtx), dca, protonKF.GetDistanceFromVertex(posPionKF), TMath::Cos(pointingAngleH), centrPerc, refMultTpc, mcStatus};
+             //fHistNdim->Fill(vec);
+
+           fHistRap->Fill(rapidity);
+           //if (pointingAngleH > 0.1) continue;
+           if (pointingAngleH > 0.05) continue;
+
+           ///////////////////////////
+           //MC part for Associated H
+           ///////////////////////////
+
+           if (HasMC() && isCorrectlyAssociatedLambda) {
+             Int_t labelP = trackP->GetLabel();
+             TParticle *tparticleDaughter = stack->Particle(TMath::Abs(labelP));
+             Int_t labelMother = tparticleDaughter->GetFirstMother();
+             TParticle *tparticleMother = stack->Particle(TMath::Abs(labelMother));
+             
+             Int_t labelProton = trackP->GetLabel();
+             Int_t labelPion = trackN->GetLabel();
+             
+             //H-Dibaryon
+             if(tparticleMother->GetPdgCode() == pdgHDibaryon && labelAssociatedH==labelMother){ //check mother PDG
+               Int_t labelFirstDaughter  = tparticleMother->GetDaughter(0);
+               Int_t labelThirdDaughter  = tparticleMother->GetDaughter(1);
+               Int_t labelSecondDaughter = labelFirstDaughter +1;
+
+               TParticle *tparticleSecondDaughter = stack->Particle(TMath::Abs(labelSecondDaughter));
+               TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));
+               
+               TLorentzVector ppi;
+               TLorentzVector lpi;
+               TLorentzVector lP;
+               
+               ppi=posProt+negPio;
+               lpi=lambdaH+negPio;
+               lP=lambdaH+posProt;
+
+               if((tparticleThirdDaughter->GetPdgCode() == pdgPionMinus && labelPion==labelThirdDaughter)&&(tparticleSecondDaughter->GetPdgCode() == pdgProton||tparticleSecondDaughter->GetPdgCode() == pdgPionMinus) && labelProton==labelSecondDaughter) fHistMassPpi->Fill(ppi.M());
+               
+               if(tparticleThirdDaughter->GetPdgCode() == pdgPionMinus && labelPion==labelThirdDaughter) fHistMassLpi->Fill(lpi.M());
+               
+               if(tparticleSecondDaughter->GetPdgCode() == pdgProton && labelProton==labelSecondDaughter) fHistMassLP->Fill(lP.M());
+               
+               if(tparticleSecondDaughter->GetPdgCode() == pdgProton && labelProton==labelSecondDaughter){//check second daughter PDG
+                 if(tparticleThirdDaughter->GetPdgCode() == pdgPionMinus && labelPion==labelThirdDaughter){//check second daughter PDG
+
+                   fHistHDibaryonInvaMassAso->Fill(hDibaryon.M()); 
+                   
+                   Double_t distance01=vecDist.Mag();
+                   fHistMCdcaPvtxDvtx->Fill(decayLengthH);
+                   fHistMCdcaPvtxLvtx->Fill(decayLength);
+                   fHistMCdcaDvtxLvtx->Fill(distance01);
+                   fHistMCangleLH->Fill(hPointingAngle);
+                   fHistMCdecayAngle->Fill(decayAngleH);
+                   fHistMCpointingAngle->Fill(pointingAngleH);
+                   fHistMCap->Fill(alfa,qt);
+
+                   fHistHilf1->Fill(posPionKF.GetDistanceFromVertex(primVtx));
+                   fHistHilf2->Fill(protonKF.GetDistanceFromVertex(primVtx));
+                   fHistHilf3->Fill(protonKF.GetDistanceFromVertex(posPionKF));
+                   fHistHilf6->Fill(dca);
+                   fHistPtvsYAso->Fill(hDibaryon.Pt(),hDibaryon.Rapidity());
+                   fHistPtvsEtaAso->Fill(hDibaryon.Pt(),hDibaryon.Eta());
+                   mcStatus=1;
+                 }//end check for third daughter PDG
+               }//end check second daughter PDG
+             }//end H-Dibaryon
+           }//end MC
+           
+           //      cout<<"Trigger: "<<triggertype<<endl;
+           fHistMassH->Fill(hDibaryon.M());
+           fHistMassHcentMult->Fill(hDibaryon.M(),triggertype,refMultTpc);
+           ppK=lambdaH+posProt;
+           fHistMassLambdaP->Fill(ppK.M());
+
+             //fHistNdim = new THnSparseF("fHistNdim","THnS;InvMass, InvMassLambda, pointingAngle, armPoAlpha, armPoQt, pTL, pTH, d0p, d0n, dcaHd, dca, decayL, cosPA, centr, multi, mcinf;InvMassH", 16,binsD01,xminD01,xmaxD01);
+
+           Double_t vec[16]={hDibaryon.M(), lInvMassLambda, pointingAngleH, alfa, qt, lPtLambda, hDibaryon.Pt(), posPionKF.GetDistanceFromVertex(primVtx), protonKF.GetDistanceFromVertex(primVtx), dca, protonKF.GetDistanceFromVertex(posPionKF), TMath::Cos(pointingAngleH), centrPerc, static_cast<Double_t>(refMultTpc), static_cast<Double_t>(mcStatus)};
+                     fHistNdim->Fill(vec);
+
+         }
+       }
+      }
+    }    
+  }
+  
+  //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+  //Pure MC Part!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+  //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+  // Monte Carlo for genenerated particles
+  if (HasMC()) //MC loop  
+    {
+
+      Int_t stackN = 0;
+
+      Double_t momentumPionGen[3]={0,0,0};
+      Double_t momentumNucleonGen[3]={0,0,0};
+      Double_t momentumLambdaGen[3]={0,0,0};
+
+      Double_t energyPionGen = 0;
+      Double_t energyNucleonGen = 0;
+      Double_t energyLambdaGen = 0;
+
+      Double_t transversMomentumMotherGen = 0;
+      Double_t longitudinalMomentumMotherGen = 0;
+      Double_t totalEnergyMotherGen = 0;
+      
+      Double_t rapidityGen = 2;
+
+      for(stackN = 0; stackN < stack->GetNtrack(); stackN++) //loop over stack
+       {
+
+         TParticle *tparticleMother = stack->Particle(stackN);
+
+         if(tparticleMother->GetPdgCode() == pdgLambda) fHistCount->Fill(16); 
+
+         //H-Dibaryon
+         if(tparticleMother->GetPdgCode() == pdgHDibaryon) //check mother PDG
+           {
+             Int_t labelFirstDaughter  = tparticleMother->GetDaughter(0);
+             Int_t labelThirdDaughter  = tparticleMother->GetDaughter(1);
+             Int_t labelSecondDaughter = labelFirstDaughter +1;
+
+             TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));
+             TParticle *tparticleSecondDaughter = stack->Particle(TMath::Abs(labelSecondDaughter));
+             TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));
+            
+             if(tparticleFirstDaughter->GetPdgCode() == pdgLambda) //check first daughter PDG
+               {
+                 if(tparticleSecondDaughter->GetPdgCode() == pdgProton)//check second daughter PDG
+                   {
+                     if(tparticleThirdDaughter->GetPdgCode() == pdgPionMinus)//check second daughter PDG
+                       {                                     
+                         momentumLambdaGen[0] = tparticleFirstDaughter->Px();
+                         momentumLambdaGen[1] = tparticleFirstDaughter->Py();
+                         momentumLambdaGen[2] = tparticleFirstDaughter->Pz();
+               
+                         momentumNucleonGen[0] = tparticleSecondDaughter->Px();
+                         momentumNucleonGen[1] = tparticleSecondDaughter->Py();
+                         momentumNucleonGen[2] = tparticleSecondDaughter->Pz();
+
+                         momentumPionGen[0] = tparticleThirdDaughter->Px();
+                         momentumPionGen[1] = tparticleThirdDaughter->Py();
+                         momentumPionGen[2] = tparticleThirdDaughter->Pz();
+
+                         TLorentzVector lorentzVectorLambda;
+                         TLorentzVector lorentzVectorProton;
+                         TLorentzVector lorentzVectorPion;
+                         TLorentzVector lorentzVectorHDibaryon;
+                         
+                         lorentzVectorLambda.SetXYZM(momentumLambdaGen[0],momentumLambdaGen[1],momentumLambdaGen[2],1.115);
+                         lorentzVectorProton.SetXYZM(momentumNucleonGen[0],momentumNucleonGen[1],momentumNucleonGen[2],protonMass);
+                         lorentzVectorPion.SetXYZM(momentumPionGen[0],momentumPionGen[1],momentumPionGen[2],pionMass);
+                         
+                         lorentzVectorHDibaryon = lorentzVectorLambda + lorentzVectorProton + lorentzVectorPion;
+                         rapidityGen=lorentzVectorHDibaryon.Rapidity();
+                         transversMomentumMotherGen = lorentzVectorHDibaryon.Pt();
+                         longitudinalMomentumMotherGen = lorentzVectorHDibaryon.Pz();
+                         totalEnergyMotherGen = lorentzVectorHDibaryon.Energy();
+
+                         if(rapidityGen > 1.0 || rapidityGen < -1 ) continue;
+                         //lorentzVectorLambda
+                         fHistHDibaryonInvaMassGen->Fill(lorentzVectorHDibaryon.M()); 
+                         if (lorentzVectorLambda.Rapidity()  > 1.0 || lorentzVectorLambda.Rapidity() < -1) continue;
+                         if (lorentzVectorProton.Rapidity()  > 1.0 ||      lorentzVectorProton.Rapidity() < -1) continue;
+
+                         if (lorentzVectorPion.Rapidity()  > 1.0 ||  lorentzVectorPion.Rapidity() < -1) continue;
+                         fHistHDibaryonInvaMassGenRes->Fill(lorentzVectorHDibaryon.M());
+                         fHistPtvsEtaGen->Fill(lorentzVectorHDibaryon.Pt(),lorentzVectorHDibaryon.Eta());
+                         fHistPtvsYGen->Fill(lorentzVectorHDibaryon.Pt(),lorentzVectorHDibaryon.Rapidity());
+                         fHistPtvsEtaGen->Fill(lorentzVectorHDibaryon.Pt(),lorentzVectorHDibaryon.Eta());
+                         fHistCount->Fill(11);
+                       }//end of check third daughter PDG
+                   }//end of check second daughter PDG
+               }//end of check first daughter PDG
+           }//end of H-Dibaryon
+
+         //Anti-H-Dibaryon
+         if(tparticleMother->GetPdgCode() == pdgAntiHDibaryon) //check mother PDG
+           {
+             Int_t labelFirstDaughter  = tparticleMother->GetDaughter(0);
+             Int_t labelThirdDaughter  = tparticleMother->GetDaughter(1);
+             Int_t labelSecondDaughter = labelFirstDaughter +1;
+
+             TParticle *tparticleFirstDaughter  = stack->Particle(TMath::Abs(labelFirstDaughter));
+             TParticle *tparticleSecondDaughter = stack->Particle(TMath::Abs(labelSecondDaughter));
+             TParticle *tparticleThirdDaughter  = stack->Particle(TMath::Abs(labelThirdDaughter));
+
+             if(tparticleFirstDaughter->GetPdgCode() == pdgAntiLambda) //check first daughter PDG
+               {
+                 if(tparticleSecondDaughter->GetPdgCode() == pdgAntiProton)//check second daughter PDG
+                   {
+                     if(tparticleThirdDaughter->GetPdgCode() == pdgPionPlus)//check second daughter PDG
+                       {                
+                         momentumLambdaGen[0] = tparticleFirstDaughter->Px();
+                         momentumLambdaGen[1] = tparticleFirstDaughter->Py();
+                         momentumLambdaGen[2] = tparticleFirstDaughter->Pz();
+               
+                         momentumNucleonGen[0] = tparticleSecondDaughter->Px();
+                         momentumNucleonGen[1] = tparticleSecondDaughter->Py();
+                         momentumNucleonGen[2] = tparticleSecondDaughter->Pz();
+
+                         momentumPionGen[0] = tparticleThirdDaughter->Px();
+                         momentumPionGen[1] = tparticleThirdDaughter->Py();
+                         momentumPionGen[2] = tparticleThirdDaughter->Pz();
+                         
+                         energyLambdaGen  = tparticleFirstDaughter->Energy();
+                         energyNucleonGen = tparticleSecondDaughter->Energy();
+                         energyPionGen    = tparticleThirdDaughter->Energy();
+                          
+                         TLorentzVector lorentzVectorLambda;
+                         TLorentzVector lorentzVectorProton;
+                         TLorentzVector lorentzVectorPion;
+                         TLorentzVector lorentzVectorHDibaryon;
+                         
+                         lorentzVectorLambda.SetXYZM(momentumLambdaGen[0],momentumLambdaGen[1],momentumLambdaGen[2],1.115);
+                         lorentzVectorProton.SetXYZM(momentumNucleonGen[0],momentumNucleonGen[1],momentumNucleonGen[2],protonMass);
+                         lorentzVectorPion.SetXYZM(momentumPionGen[0],momentumPionGen[1],momentumPionGen[2],pionMass);
+                         
+                         lorentzVectorHDibaryon = lorentzVectorLambda + lorentzVectorProton + lorentzVectorPion;
+
+                         rapidityGen=lorentzVectorHDibaryon.Rapidity();
+                         if(rapidityGen > 1.0 || rapidityGen < -1 ) continue;
+                         fHistAntiHDibaryonInvaMassGen->Fill(lorentzVectorHDibaryon.M()); 
+                       }//end of check third daughter PDG
+                   }//end of check second daughter PDG
+               }//end of check first daughter PDG
+           }//end of Anti-H-Dibaryon
+       }      
+    }//end MC
+
+  // Post output data.
+  PostData(1,fHistList);
+  //PostData(0,fHistList);
+
     if (listCrossV0 == NULL) return;
 
-  if (listCrossV0) delete listCrossV0;\r
-  if (esdVer1) delete esdVer1;\r
-  if (vertexer) delete vertexer;\r
-  if (vertexer1) delete vertexer1;\r
-  if (trkArray) delete trkArray;\r
-  if (trkArray1) delete trkArray1;\r
-}\r
-\r
-//________________________________________________________________________\r
-void AliAnalysisTaskHdibaryonLPpi::Terminate(Option_t *) \r
-{\r
-  // Draw result to the screen\r
-  // Called once at the end of the query\r
-\r
-}\r
-\r
+  if (listCrossV0) delete listCrossV0;
+  if (esdVer1) delete esdVer1;
+  if (vertexer) delete vertexer;
+  if (vertexer1) delete vertexer1;
+  if (trkArray) delete trkArray;
+  if (trkArray1) delete trkArray1;
+}
+
+//________________________________________________________________________
+void AliAnalysisTaskHdibaryonLPpi::Terminate(Option_t *) 
+{
+  // Draw result to the screen
+  // Called once at the end of the query
+
+}
+