- adding check that digits contain data from HG channel, if not

we might get high energy clusters from noisy LG channels

diff --git a/HLT/CALO/AliHLTCaloClusterizerComponent.cxx b/HLT/CALO/AliHLTCaloClusterizerComponent.cxx
index 0e8b7640d8cf1b7cbe970ab6909e5708501ccf85..a7fb1416d229f5d441da5a5380a01c91d457ccd1 100644 (file)
--- a/HLT/CALO/AliHLTCaloClusterizerComponent.cxx
+++ b/HLT/CALO/AliHLTCaloClusterizerComponent.cxx
@@ -114,9 +114,14 @@ AliHLTCaloClusterizerComponent::DoEvent(const AliHLTComponentEventData& evtData,
             // Could be changed if this is a bottle neck.
             for (Int_t i = 0; i < nDigits; i++)
             {
+              // If we have a digit based on a low gain channel, but there has been no high gain channel, 
+              // we shouldn't use it since we are then very sensitive to noise (e.g. for PHOS 1 LG ADC count = 40 MeV)
+              if(digitDataPtr->fHgPresent)
+              {
                 fDigitsPointerArray[digCount] = digitDataPtr;
                 digCount++;
                 digitDataPtr++;
+              }
             }
         }
     }
@@ -282,6 +287,7 @@ AliHLTCaloClusterizerComponent::DoInit(int argc, const char** argv )
             fAnalyserPtr->SetRecoParamHandler(fRecoParamsPtr);
             fClusterizerPtr->SetEmcClusteringThreshold(fRecoParamsPtr->GetRecPointThreshold());
             fClusterizerPtr->SetEmcMinEnergyThreshold(fRecoParamsPtr->GetRecPointMemberThreshold());
+           HLTInfo("Setting thresholds for clusterizer: %f, %f", fRecoParamsPtr->GetRecPointThreshold(), fRecoParamsPtr->GetRecPointMemberThreshold());
         }
     }
     //