Changeset 305 in svn for trunk/src/SmearUtil.cc

Timestamp:

Mar 9, 2009, 11:40:41 AM (16 years ago)

Author:

severine ovyn

Message:

add isolation to datacard

File:

• trunk/src/SmearUtil.cc (modified) (7 diffs)

trunk/src/SmearUtil.cc

@@ -102 +102 @@
   PTCUT_taujet    = 10.0;
 
+  ISOL_PT         = 2.0;      //minimal pt of tracks for isolation criteria
+  ISOL_Cone       = 0.5;      //Cone  for isolation criteria
+
+
   // General jet variable
   JET_coneradius   = 0.7;            // generic jet radius ; not for tau's !!!
@@ -223 +227 @@
   PTCUT_taujet    = DET.PTCUT_taujet;
 
+  ISOL_PT         = DET.ISOL_PT;      //minimal pt of tracks for isolation criteria
+  ISOL_Cone       = DET.ISOL_Cone;      //Cone  for isolation criteria
+
+
   // General jet variable
   JET_coneradius   = DET.JET_coneradius;
@@ -336 +344 @@
   PTCUT_gamma     = DET.PTCUT_gamma;
   PTCUT_taujet    = DET.PTCUT_taujet;
+
+  ISOL_PT         = DET.ISOL_PT;      //minimal pt of tracks for isolation criteria
+  ISOL_Cone       = DET.ISOL_Cone;      //Cone  for isolation criteria
+
 
   // General jet variable
@@ -482 +494 @@
     else if(strstr(temp_string.c_str(),"PTCUT_taujet"))     {curstring >> varname >> value; PTCUT_taujet      = value;}
 
+    else if(strstr(temp_string.c_str(),"ISOL_PT"))          {curstring >> varname >> value; ISOL_PT        = value;}
+    else if(strstr(temp_string.c_str(),"ISOL_Cone"))        {curstring >> varname >> value; ISOL_Cone        = value;}
+
+
     else if(strstr(temp_string.c_str(),"JET_coneradius"))   {curstring >> varname >> value; JET_coneradius    = value;}
     else if(strstr(temp_string.c_str(),"JET_jetalgo"))      {curstring >> varname >> ivalue;JET_jetalgo       = ivalue;}
@@ -745 +761 @@
         << left << setw(20) <<PTCUT_gamma        <<""<< right << setw(10)<<"*"<<"\n";
   f_out<<"*                                                                    *"<<"\n";
+  f_out<<"#*******************                                                 *"<<"\n";
+  f_out<<"# Isolation criteria                                                 *"<<"\n";
+  f_out<<"#*******************                                                 *"<<"\n";
+  f_out<<"*                                                                    *"<<"\n";
+  f_out << left << setw(40) <<"* Minimum pT for tracks [GeV]: "<<""
+        << left << setw(20) <<ISOL_PT               <<""<< right << setw(10)<<"*"<<"\n";
+  f_out << left << setw(40) <<"* Cone for isolation criteria: "<<""
+        << left << setw(20) <<ISOL_Cone             <<""<< right << setw(10)<<"*"<<"\n";
+  f_out<<"*                                                                    *"<<"\n";
   f_out<<"#***************                                                     *"<<"\n";
   f_out<<"# Jet definition                                                     *"<<"\n";
@@ -1050 +1075 @@
 //****************************************************************
 //bool RESOLution::Isolation(const float phi, const float eta,const vector<TLorentzVector> &tracks, const float pt_second_track)
-bool RESOLution::Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float pt_second_track)
+bool RESOLution::Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float pt_second_track, const float isolCone)
 {
    bool isolated = false;
@@ -1065 +1090 @@
            deltar=genDeltaR;
       }
-   if(deltar > 0.5) isolated = true;
+   if(deltar > isolCone) isolated = true;
    return isolated;
 }