name = "DCAz"; _dcaz = createHisto1F(name,name, 500, -5.0, 5.0, "dcaZ","counts");
name = "DCAxy"; _dcaxy = createHisto1F(name,name, 500, -5.0, 5.0, "dcaXY","counts");
- //name = "Nclus1"; _Ncluster1 = createHisto1F(name,name, 200, 0, 200, "Ncluster1","counts");
+ name = "Nclus1"; _Ncluster1 = createHisto1F(name,name, 200, 0, 200, "Ncluster1","counts");
//name = "Nclus2"; _Ncluster2 = createHisto1F(name,name, 200, 0, 200, "Ncluster2","counts");
if (_singlesOnly)
//double b[2];
//double bCov[3];
const AliAODVertex* vertex;
+ int nClus;
bool bitOK;
AliAnalysisManager* manager = AliAnalysisManager::GetAnalysisManager();
eta = t->Eta();
//dcaXY = t->DCA();
//dcaZ = t->ZAtDCA();
-
+ nClus = t->GetTPCNcls();
+
+ if ( nClus<_nClusterMin ) continue;
+
+ _Ncluster1->Fill(nClus);
+
//for Global tracks
Double_t nsigmaelectron = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
Double_t nsigmapion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kPion));