Changeset 31 in svn

Timestamp:

Nov 13, 2008, 5:15:34 PM (16 years ago)

Author:

severine ovyn

Message:

add isolation

Location:

trunk

Files:

• interface/SmearUtil.h (modified) (1 diff)
• src/SmearUtil.cc (modified) (3 diffs)

trunk/interface/SmearUtil.h

@@ -128 +128 @@
   bool Btaggedjet(const TLorentzVector &JET, const TSimpleArray<TRootGenParticle> &subarray);
 
+  //******************************isolation criteria**************************************
+  bool Isolation(Float_t phi,Float_t eta,const vector<TLorentzVector> &tracks,float PT_TRACK2);
+
 };
 

trunk/src/SmearUtil.cc

@@ -180 +180 @@
                                          pow(ELG_Sfwd*sqrt(energy),2) ) ); 
   }
+float theta=electron.Theta();
+float phi=electron.Phi();
+  float px = energyS * sin(theta) * cos(phi);
+  float py = energyS * sin(theta) * sin(phi);
+  float pz = energyS * cos(theta);
+  //cout<<" px "<<px<<" py "<<py<<" pz "<<pz<<endl;
   electron.SetPtEtaPhiE(energyS/cosh(electron.Eta()), electron.Eta(), electron.Phi(), energyS);
-  
+  //cout<<"nous "<<electron.Px()<<"  "<<electron.Py()<<"  "<<electron.Pz()<<endl;
   if(electron.E() < 0)electron.SetPxPyPzE(0,0,0,0); // no negative values after smearing !
 }
@@ -212 +218 @@
   // some of the particles decay mostly in the ecal, some mostly in the hcal
   
-float energyS1,energyS2;
+  float energyS1,energyS2;
   if(fabs(hadron.Eta()) < MAX_CALO_CEN) {
    energyS1 = gRandom->Gaus(energy_hcal, sqrt(
@@ -302 +308 @@
 }
 
+//***********************Isolation criteria***********************
+//****************************************************************
+bool RESOLution::Isolation(Float_t phi,Float_t eta,const vector<TLorentzVector> &tracks,float PT_TRACK2)
+{
+   bool isolated = false;
+   Float_t deltar=5000.; // Initial value; should be high; no further repercussion
+   // loop on all final charged particles, with p_t >2, close enough from the electron
+   for(unsigned int i=0; i < tracks.size(); i++) 
+      {
+       if(tracks[i].Pt() < PT_TRACK2)continue;
+       Float_t genDeltaR = DeltaR(phi,eta,tracks[i].Phi(),tracks[i].Eta());    // slower to evaluate
+         if(
+             (genDeltaR > deltar) ||
+             (genDeltaR==0)
+           ) continue ;
+           deltar=genDeltaR;
+      }
+   if(deltar > 0.5)isolated = true; // returns the closest distance
+   return isolated;
+}
+
+
 //**************************** Returns the delta Phi ****************************
 float DeltaPhi(const float phi1, const float phi2) {