// Primary tracks from ESD/AOD
Float_t lPrimaryTrackMultiplicity= -1.;
+// Int_t lSPDTrackletsMultiplicity = -1;
+
AliCentrality* centrality = 0x0;
AliESDVZERO* esdV0 = 0x0;
AliAODVZERO* aodV0 = 0x0;
lPrimaryTrackMultiplicity = fESDtrackCuts->CountAcceptedTracks(lESDevent);
+// const AliMultiplicity *lAliMult = lESDevent->GetMultiplicity();
+// lSPDTrackletsMultiplicity = lAliMult->GetNumberOfTracklets();
+
// nNumberOfMCPrimaries = lMCstack->GetNprimary();
nNumberOfMCPrimaries = lMCstack->GetNtrack(); // MN: this stack->GetNtrack(); has to be used to check primaries because in HIJING decay products of D and B mesons are also primaries and produced in HIJING during transport, however this is not the number of primaries!
// Centrality determination
lPrimaryTrackMultiplicity = 0.;
for (Int_t itrack = 0; itrack<nTrackMultiplicity; itrack++) {
AliAODTrack* track = lAODevent->GetTrack(itrack);
- if (track->TestFilterBit(AliAODTrack::kTrkGlobal)) lPrimaryTrackMultiplicity++;
+ if (track->TestFilterBit(AliAODTrack::kTrkGlobalNoDCA)) lPrimaryTrackMultiplicity++;
}
+// lSPDTrackletsMultiplicity = lAODevent->GetTracklets()->GetNumberOfTracklets();
+
nNumberOfMCPrimaries = arrayMC->GetEntries();
Double_t cascadeMass = 0.;
-Int_t lSPDTrackletsMultiplicity = -1;
-
-
for (Int_t iXi = 0; iXi < ncascades; iXi++) {// This is the begining of the Cascade loop
// - Step 6 : Containers = Cascade cuts + PID
//-------------
- if (fAnalysisType == "ESD") {
- const AliMultiplicity *lAliMult = lESDevent->GetMultiplicity();
- lSPDTrackletsMultiplicity = lAliMult->GetNumberOfTracklets();
- }
- else if (fAnalysisType == "AOD") lSPDTrackletsMultiplicity = lAODevent->GetTracklets()->GetNumberOfTracklets();
-
-
// 6.3 - Filling the AliCFContainer (optimisation of topological selections + systematics)
Double_t lContainerCutVars[21] = {0.0};