// p-p collisions
////////////////////////////////////////////////////////////////
-const Int_t AliITSVertexer3D::fgkMaxNumOfClDefault = 500;
+const Int_t AliITSVertexer3D::fgkMaxNumOfClDefault = 1000;
ClassImp(AliITSVertexer3D)
if(nolines>=2){
rc=Prepare3DVertex(1);
if(fPileupAlgo == 2 && rc == 0) FindVertex3DIterative();
- else if(fPileupAlgo<2 && rc == 0) FindVertex3D(itsClusterTree);
+ else if(fPileupAlgo!=2 && rc == 0) FindVertex3D(itsClusterTree);
if(rc!=0) fVert3D.SetNContributors(0); // exclude this vertex
}
}
switch(fPileupAlgo){
case 0: PileupFromZ(); break;
case 1: FindOther3DVertices(itsClusterTree); break;
+ case 3: break; // no pileup algo
default: AliError("Wrong pileup algorithm"); break;
}
if(fNoVertices==1){
Int_t nolines = 0;
// Loop on modules of layer 1
- Int_t firstL1 = AliITSgeomTGeo::GetModuleIndex(1,1,1);
+ Int_t firstL1 = TMath::Max(0,AliITSgeomTGeo::GetModuleIndex(1,1,1));
Int_t lastL1 = AliITSgeomTGeo::GetModuleIndex(2,1,1)-1;
for(Int_t modul1= firstL1; modul1<=lastL1;modul1++){ // Loop on modules of layer 1
if(!fUseModule[modul1]) continue;
Int_t ladmod=fLadders[ladder-1]+ladl2;
if(ladmod>AliITSgeomTGeo::GetNLadders(2)) ladmod=ladmod-AliITSgeomTGeo::GetNLadders(2);
Int_t modul2=AliITSgeomTGeo::GetModuleIndex(2,ladmod,k+1);
+ if(modul2<0)continue;
if(!fUseModule[modul2]) continue;
itsRec=rpcont->UncheckedGetClusters(modul2);
Int_t nrecp2 = itsRec->GetEntries();