Made more robust

diff --git a/EMCAL/AliEMCALJetFinderAlgoOmni.cxx b/EMCAL/AliEMCALJetFinderAlgoOmni.cxx
index 6c473c66988650891586f72d3f9b7f10d5db57cf..3f25e4d6796f7b1706f9ba8c4eb43ab50866b986 100644 (file)
--- a/EMCAL/AliEMCALJetFinderAlgoOmni.cxx
+++ b/EMCAL/AliEMCALJetFinderAlgoOmni.cxx
@@ -20,6 +20,9 @@
 /*
  
 $Log$
 /*
  
 $Log$

+Revision 1.10  2004/03/06 01:56:09  mhorner
+Pythai comp code
+

 Revision 1.9  2004/02/23 20:37:32  mhorner
 changed geometry call
 
 Revision 1.9  2004/02/23 20:37:32  mhorner
 changed geometry call
 


@@ -327,12 +330,12 @@ if (fDebug>0) Info("AliEMCALJetFinderAlgoOmni","Beginning Default Constructor");
             for(Int_t k=numTracks; k < fNumUnits; k++)
             {
                     fUnit[k].SetUnitID(k);
             for(Int_t k=numTracks; k < fNumUnits; k++)
             {
                     fUnit[k].SetUnitID(k);

-                    fUnit[k].SetUnitFlag(kOutJet);
+                    fUnit[k].SetUnitFlag(kBelowMinEt);

                     fUnit[k].SetUnitEta(0.0);
                     fUnit[k].SetUnitPhi(0.0);
                     fUnit[k].SetUnitEnergy(0.0);                    
                     fUnitNoCuts[k].SetUnitID(k);
                     fUnit[k].SetUnitEta(0.0);
                     fUnit[k].SetUnitPhi(0.0);
                     fUnit[k].SetUnitEnergy(0.0);                    
                     fUnitNoCuts[k].SetUnitID(k);

-                    fUnitNoCuts[k].SetUnitFlag(kOutJet);
+                    fUnitNoCuts[k].SetUnitFlag(kBelowMinEt);

                     fUnitNoCuts[k].SetUnitEta(0.0);
                     fUnitNoCuts[k].SetUnitPhi(0.0);
                     fUnitNoCuts[k].SetUnitEnergy(0.0);                   
                     fUnitNoCuts[k].SetUnitEta(0.0);
                     fUnitNoCuts[k].SetUnitPhi(0.0);
                     fUnitNoCuts[k].SetUnitEnergy(0.0);                   


@@ -495,8 +498,11 @@ if (fDebug>0) Info("AliEMCALJetFinderAlgoOmni","Beginning Default Constructor");
      fJetEtaSum += fEnergy * fDEta;
      fJetPhiSum += fEnergy * fDPhi;
      fJetESum += fEnergy;
      fJetEtaSum += fEnergy * fDEta;
      fJetPhiSum += fEnergy * fDPhi;
      fJetESum += fEnergy;

-     fJetEta = fEtaInit + (fJetEtaSum / fJetESum);
-     fJetPhi = fPhiInit + (fJetPhiSum / fJetESum);
+     if (fJetESum >1.0e-7)
+     {
+            fJetEta = fEtaInit + (fJetEtaSum / fJetESum);
+            fJetPhi = fPhiInit + (fJetPhiSum / fJetESum);
+     }

    }
 
 
    }
 
 


@@ -783,7 +789,7 @@ if (fDebug>0) Info("AliEMCALJetFinderAlgoOmni","Beginning Default Constructor");
 
                      //Find the distance cone centre is from initiator cell
                      if (fDebug>10) Info("FindJets","Checking if cone move too large");
 
                      //Find the distance cone centre is from initiator cell
                      if (fDebug>10) Info("FindJets","Checking if cone move too large");

-                     if ( ((fJetEtaSum)*(fJetEtaSum))+((fJetPhiSum)*(fJetPhiSum)) >1.e-4)
+                     if ( ((fJetEtaSum)*(fJetEtaSum))+((fJetPhiSum)*(fJetPhiSum)) >1.e-7)

                      {
                              fDistI = TMath::Sqrt( ((fJetEtaSum/fJetESum)*(fJetEtaSum/fJetESum)) + ((fJetPhiSum/fJetESum)*
                                                                                                    (fJetPhiSum/fJetESum)));
                      {
                              fDistI = TMath::Sqrt( ((fJetEtaSum/fJetESum)*(fJetEtaSum/fJetESum)) + ((fJetPhiSum/fJetESum)*
                                                                                                    (fJetPhiSum/fJetESum)));