else man->SetDefaultStorage("alien://Folder=/alice/data/2010/OCDB");
man->SetRun(137161);
AliCDBEntry* obj = man->Get(AliCDBPath("GRP", "Geometry", "Data"));
+ if(!obj) AliFatal("Unable to load geometry from CDB!");
AliGeomManager::SetGeometry((TGeoManager*) obj->GetObject());
AliGeomManager::GetNalignable("ITS");
AliGeomManager::ApplyAlignObjsFromCDB("ITS");
delete[] partners;
delete[] minDists;
+ delete[] associatedLay1;
for (Int_t i=0; i<fNClustersLay1; i++)
if (blacklist[i])
if (idt==80)
Printf("First Cl2 LoadClArr %f %f %f ",cluGlo[0],cluGlo[1],cluGlo[2]);*/
new (clArr[nclLayer++]) AliITSRecPoint(*cluster);
- nClustersInChip[ seg.GetChipFromLocal(0,cluster->GetDetLocalZ()) ]++;
+ if (seg.GetChipFromLocal(0,cluster->GetDetLocalZ())>=0)
+ nClustersInChip[ seg.GetChipFromLocal(0,cluster->GetDetLocalZ()) ]++;
}
for(Int_t ifChip=5;ifChip--;) if (nClustersInChip[ifChip]) fNFiredChips[il]++;
}
AliITSsegmentationSPD seg;
AliITSRecPointContainer* rpcont=AliITSRecPointContainer::Instance();
- TClonesArray* itsClusters=rpcont->FetchClusters(0,itsClusterTree);
+ TClonesArray* itsClusters=NULL;
+ rpcont->FetchClusters(0,itsClusterTree);
if(!rpcont->IsSPDActive()){
AliWarning("No SPD rec points found, multiplicity not calculated");
return;
// find the chip for the current cluster
Float_t locz = cluster->GetDetLocalZ();
Int_t iChip = seg.GetChipFromLocal(0,locz);
- nClustersInChip[iChip]++;
+ if (iChip>=0)
+ nClustersInChip[iChip]++;
}// end of cluster loop