for(Int_t j=0; j<nOfPlanesInRpPerCoordinate; ++j){ Int_t nClusters = rps->numberOfClusters(i,Planes[coordinate][j]); nOfClusters[i][Planes[coordinate][j]]->Fill(nClusters); if(nClusters < maxNumberOfClusterPerPlane) for(Int_t k=0; k < nClusters; ++k){ Int_t lenCluster = rps->lengthCluster(i,Planes[coordinate][j],k); clusterLength[i][Planes[coordinate][j]]->Fill(lenCluster); if(lenCluster <= maxClusterLength && lenCluster>0){ Int_t enCluster = rps->energyCluster(i,Planes[coordinate][j],k); clusterEnergy[i][Planes[coordinate][j]][lenCluster-1]->Fill(enCluster); if(lenCluster==1) clusterEnergy_vs_strip[i][Planes[coordinate][j]]->Fill(1e3*rps->positionCluster(i,Planes[coordinate][j],k)/Pitch[coordinate], enCluster); } } }