New code for full visualization of the vzero raw data. ctrl-alt-letf-button shows...

[u/mrichter/AliRoot.git] / PHOS / AliPHOSClusterizerv1.cxx

diff --git a/PHOS/AliPHOSClusterizerv1.cxx b/PHOS/AliPHOSClusterizerv1.cxx
index 33098168647765d97783799e9ade37d28da9eb91..5f5d221c391dbda69a0743a75f729512dd7d4622 100644 (file)
--- a/PHOS/AliPHOSClusterizerv1.cxx
+++ b/PHOS/AliPHOSClusterizerv1.cxx
@@ -384,8 +384,10 @@ void AliPHOSClusterizerv1::InitParameters()
   const AliPHOSRecoParam* recoParam = AliPHOSReconstructor::GetRecoParam();
   if(!recoParam) AliFatal("Reconstruction parameters are not set!");
 
-  fCpvClusteringThreshold  = recoParam->GetEMCClusteringThreshold();
-  fEmcClusteringThreshold  = recoParam->GetCPVClusteringThreshold();
+  recoParam->Print();
+
+  fEmcClusteringThreshold  = recoParam->GetEMCClusteringThreshold();
+  fCpvClusteringThreshold  = recoParam->GetCPVClusteringThreshold();
   
   fEmcLocMaxCut            = recoParam->GetEMCLocalMaxCut();
   fCpvLocMaxCut            = recoParam->GetCPVLocalMaxCut();
@@ -541,14 +543,20 @@ void AliPHOSClusterizerv1::MakeClusters()
   fNumberOfEmcClusters     = 0 ;
 
   //Mark all digits as unused yet
+  const Int_t maxNDigits = 1500;
   Int_t nDigits=fDigitsArr->GetEntriesFast() ;
+  if (nDigits > maxNDigits) {
+    AliWarning(Form("Skip event with too high digit occupancy: nDigits=%d",nDigits));
+    return;
+  }
+
   for(Int_t i=0; i<nDigits; i++){
     fDigitsUsed[i]=0 ;
   }
   Int_t iFirst = 0 ; //first index of digit which potentially can be a part of cluster
                      //e.g. first digit in this module, first CPV digit etc.
   AliPHOSDigit * digit ; 
-  TArrayI clusterdigitslist(1500) ;   
+  TArrayI clusterdigitslist(maxNDigits) ;   
   AliPHOSRecPoint * clu = 0 ; 
   for(Int_t i=0; i<nDigits; i++){
     if(fDigitsUsed[i])