--- UserCode/MitPhysics/Utils/interface/IsolationTools.h	2011/05/12 23:16:11	1.14
+++ UserCode/MitPhysics/Utils/interface/IsolationTools.h	2011/11/02 20:12:03	1.21
@@ -1,5 +1,5 @@
 //--------------------------------------------------------------------------------------------------
-// $Id: IsolationTools.h,v 1.14 2011/05/12 23:16:11 ceballos Exp $
+// $Id: IsolationTools.h,v 1.21 2011/11/02 20:12:03 ceballos Exp $
 //
 // IsolationTools
 //
@@ -13,6 +13,7 @@
 
 #include <TMath.h>
 #include "MitAna/DataTree/interface/Track.h"
+#include "MitAna/DataTree/interface/Photon.h"
 #include "MitAna/DataTree/interface/BasicCluster.h"
 #include "MitAna/DataTree/interface/SuperCluster.h"
 #include "MitAna/DataTree/interface/CaloTower.h"
@@ -39,21 +40,21 @@ namespace mithep
       static Double_t CaloTowerEmIsolation(const ThreeVector *p, Double_t extRadius, 
                                            Double_t intRadius, Double_t etLow, 
                                            const mithep::Collection<mithep::CaloTower> *caloTowers);
-      static Double_t PFMuonIsolation(const Muon *p, const PFCandidateCol *PFCands, 
-                      	   	      const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
-			   	      Double_t extRadius, Double_t intRadius);
       static Double_t PFMuonIsolation(const Muon *p, const PFCandidateCol *PFCands, const Vertex *vertex, 
-				      const MuonCol *goodMuons, const ElectronCol *goodElectrons,
-				      Double_t  delta_z = 0.2, Double_t ptMin = 1.0, Double_t extRadius = 0.4,
-				      Double_t intRadius = 0.0, int isoType = 0, Double_t beta = 1.0);
+                      	   	      Double_t  delta_z = 0.1, Double_t ptMin = 1.0,
+			   	      Double_t extRadius = 0.4, Double_t intRadiusGamma = 0.07, Double_t intRadius = 0.0);
+      static Double_t PFMuonIsolation(const Muon *p, const PFCandidateCol *PFCands,
+                                      const MuonCol *goodMuons, const ElectronCol *goodElectrons, 
+                                      const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
+				      Double_t extRadius, Double_t intRadiusGamma, Double_t intRadius);
       static Double_t PFElectronIsolation(const Electron *p, const PFCandidateCol *PFCands, 
                       	     		  const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
 			     		  Double_t extRadius, Double_t intRadius);
-      static Double_t PFElectronIsolation(const Electron *p, const PFCandidateCol *PFCands, const Vertex *vertex, 
-                                          const MuonCol *goodMuons, const ElectronCol *goodElectrons, 
-                                          Double_t  delta_z = 0.2 , Double_t ptMin = 1.0,Double_t extRadius = 0.4,			     		  
-                                          Double_t intRadius = 0.0, int isoType = 0, Double_t beta = 1.0);
-      static Double_t BetaM(const TrackCol *tracks, const Muon *p, const Vertex *vertex, 
+      static Double_t PFElectronIsolation(const Electron *p, const PFCandidateCol *PFCands, 
+                      	     		  const MuonCol *goodMuons, const ElectronCol *goodElectrons,
+					  const Vertex *vertex, Double_t  delta_z, Double_t ptMin,
+			     		  Double_t extRadius, Double_t intRadius);
+       static Double_t BetaM(const TrackCol *tracks, const Muon *p, const Vertex *vertex, 
                             Double_t ptMin, Double_t  delta_z, Double_t extRadius,
 			    Double_t intRadius);
       static Double_t BetaE(const TrackCol *tracks, const Electron *p, const Vertex *vertex, 
@@ -72,6 +73,24 @@ namespace mithep
                                          UInt_t maxNExpectedHitsInner = 999,
                                          const mithep::DecayParticleCol *conversions = 0);
 
+      // methods for Hgg BaseLien Selection. These isoaltion are stupid, but what can we do.... ;(
+      static Double_t CiCTrackIsolation(const mithep::Photon*, 
+					const BaseVertex*, 
+					Double_t extRadius, 
+					Double_t intRadius, 
+					Double_t ptLow, 
+					Double_t etaStrip,
+					Double_t maxD0,
+					Double_t maxDZ,
+					const mithep::Collection<mithep::Track> *tracks,
+					unsigned int* worstVtxIdx = NULL,
+					const mithep::Collection<mithep::Vertex> *vtxs = NULL,
+					const mithep::Collection<mithep::Electron> *eles = NULL,
+					bool print=false,
+					double* ptmax=NULL, 
+					double* dRmax=NULL);
+
+
     ClassDef(IsolationTools, 0) // Isolation tools
   };
 }