--- PamCut/CaloCuts/CaloNotIntCut/CaloNotIntCut.cpp 2009/10/27 18:20:18 1.2 +++ PamCut/CaloCuts/CaloNotIntCut/CaloNotIntCut.cpp 2009/11/24 14:42:31 1.4 @@ -18,20 +18,29 @@ // energy released along the event axis in the calo / total energy release in the calo - float qratio = ( _xCaloAxis->GetQaxis() + _yCaloAxis->GetQaxis() ) / event->GetCaloLevel2()->qtot; + Double_t qtot = event->GetCaloLevel2()->qtot; + if ( qtot>0 ) { + + Double_t qratio = ( _xCaloAxis->GetQaxis() + _yCaloAxis->GetQaxis() ) / qtot; - if (! (qratio > _qRatioMin) ) { - return 1; - } + if (! (qratio > _qRatioMin) ) { + return 1; + } + } } if (_pamTrack) { // PAMELA track - float qratio = ( _pamTrack->GetCaloTrack()->qtrack ) / event->GetCaloLevel2()->qtot; + Double_t qtot = event->GetCaloLevel2()->qtot; + if ( qtot>0 ) { + + Double_t qratio = ( (*_pamTrack)->GetCaloTrack()->qtrack ) / qtot; - if (! (qratio > _qRatioMin) ) { - return 2; + if (! (qratio > _qRatioMin) ) { + return 2; + } + } }