Fork me on GitHub

Changeset 305 in svn


Ignore:
Timestamp:
Mar 9, 2009, 11:40:41 AM (16 years ago)
Author:
severine ovyn
Message:

add isolation to datacard

Location:
trunk
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/Delphes.cpp

    r304 r305  
    579579        elementElec->PhiCalo = electron[i].PhiCalo();
    580580        elementElec->Charge = sign(electron[i].PID());
    581         elementElec->IsolFlag = DET->Isolation(electron[i].Phi(),electron[i].Eta(),TrackCentral,2.0);//isolation based on tracks
     581        //elementElec->IsolFlag = DET->Isolation(electron[i].Phi(),electron[i].Eta(),TrackCentral,2.0,0.5);//isolation based on tracks
     582        elementElec->IsolFlag = DET->Isolation(electron[i].Phi(),electron[i].Eta(),TrackCentral,DET->ISOL_PT,DET->ISOL_Cone);//isolation based on tracks
    582583      }                                                                                 /////////////// HARDCODING
    583584      DET->SortedVector(muon);
    … …  
    588589        elementMu->EtaCalo = muon[i].EtaCalo();
    589590        elementMu->PhiCalo = muon[i].PhiCalo();
    590         elementMu->IsolFlag = DET->Isolation(muon[i].Phi(),muon[i].Eta(),TrackCentral,2.0); /////////////// HARDCODING
     591        //elementMu->IsolFlag = DET->Isolation(muon[i].Phi(),muon[i].Eta(),TrackCentral,2.0,0.5); /////////////// HARDCODING
     592        elementMu->IsolFlag = DET->Isolation(muon[i].Phi(),muon[i].Eta(),TrackCentral,DET->ISOL_PT,DET->ISOL_Cone); /////////////// HARDCODING
    591593      }
    592594      DET->SortedVector(gamma);

  • trunk/data/DetectorCard.dat

    r303 r305  
    5454PTCUT_taujet     10.0
    5555
     56ISOL_PT          2.0      //minimal pt of tracks for isolation criteria
     57ISOL_Cone        0.5      //Cone  for isolation criteria
     58
     59
    5660# General jet variable
    5761JET_coneradius   0.7            // generic jet radius ; not for tau's !!!

  • trunk/data/DetectorCard_ATLAS.dat

    r302 r305  
    5454PTCUT_taujet     10.0
    5555
     56ISOL_PT          2.0      //minimal pt of tracks for isolation criteria
     57ISOL_Cone        0.5      //Cone  for isolation criteria
     58
     59
    5660# General jet variable
    5761JET_coneradius   0.7            // generic jet radius ; not for tau's !!!

  • trunk/data/DetectorCard_CMS.dat

    r303 r305  
    5454PTCUT_taujet     10.0
    5555
     56ISOL_PT          2.0      //minimal pt of tracks for isolation criteria
     57ISOL_Cone        0.5      //Cone  for isolation criteria
     58
     59
    5660# General jet variable
    5761JET_coneradius   0.7            // generic jet radius ; not for tau's !!!

  • trunk/interface/SmearUtil.h

    r287 r305  
    152152  float PTCUT_gamma;
    153153  float PTCUT_taujet;
     154
     155  float ISOL_PT;      //minimal pt of tracks for isolation criteria
     156  float ISOL_Cone;    //Cone  for isolation criteria
    154157 
    155158  //General jet variable
    … …  
    228231  /// Lepton isolation
    229232  //bool Isolation(const float phi, const float eta,const vector<TLorentzVector> &tracks,float PT_TRACK2);
    230   bool Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float PT_TRACK2);
     233  bool Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float PT_TRACK2, const float isolCone);
    231234
    232235  //********************* returns a segmented value for eta and phi, for calo towers *****

  • trunk/src/SmearUtil.cc

    r287 r305  
    102102  PTCUT_taujet    = 10.0;
    103103
     104  ISOL_PT         = 2.0;      //minimal pt of tracks for isolation criteria
     105  ISOL_Cone       = 0.5;      //Cone  for isolation criteria
     106
     107
    104108  // General jet variable
    105109  JET_coneradius   = 0.7;            // generic jet radius ; not for tau's !!!
    … …  
    223227  PTCUT_taujet    = DET.PTCUT_taujet;
    224228
     229  ISOL_PT         = DET.ISOL_PT;      //minimal pt of tracks for isolation criteria
     230  ISOL_Cone       = DET.ISOL_Cone;      //Cone  for isolation criteria
     231
     232
    225233  // General jet variable
    226234  JET_coneradius   = DET.JET_coneradius;
    … …  
    336344  PTCUT_gamma     = DET.PTCUT_gamma;
    337345  PTCUT_taujet    = DET.PTCUT_taujet;
     346
     347  ISOL_PT         = DET.ISOL_PT;      //minimal pt of tracks for isolation criteria
     348  ISOL_Cone       = DET.ISOL_Cone;      //Cone  for isolation criteria
     349
    338350
    339351  // General jet variable
    … …  
    482494    else if(strstr(temp_string.c_str(),"PTCUT_taujet"))     {curstring >> varname >> value; PTCUT_taujet      = value;}
    483495
     496    else if(strstr(temp_string.c_str(),"ISOL_PT"))          {curstring >> varname >> value; ISOL_PT        = value;}
     497    else if(strstr(temp_string.c_str(),"ISOL_Cone"))        {curstring >> varname >> value; ISOL_Cone        = value;}
     498
     499
    484500    else if(strstr(temp_string.c_str(),"JET_coneradius"))   {curstring >> varname >> value; JET_coneradius    = value;}
    485501    else if(strstr(temp_string.c_str(),"JET_jetalgo"))      {curstring >> varname >> ivalue;JET_jetalgo       = ivalue;}
    … …  
    745761        << left << setw(20) <<PTCUT_gamma        <<""<< right << setw(10)<<"*"<<"\n";
    746762  f_out<<"*                                                                    *"<<"\n";
     763  f_out<<"#*******************                                                 *"<<"\n";
     764  f_out<<"# Isolation criteria                                                 *"<<"\n";
     765  f_out<<"#*******************                                                 *"<<"\n";
     766  f_out<<"*                                                                    *"<<"\n";
     767  f_out << left << setw(40) <<"* Minimum pT for tracks [GeV]: "<<""
     768        << left << setw(20) <<ISOL_PT               <<""<< right << setw(10)<<"*"<<"\n";
     769  f_out << left << setw(40) <<"* Cone for isolation criteria: "<<""
     770        << left << setw(20) <<ISOL_Cone             <<""<< right << setw(10)<<"*"<<"\n";
     771  f_out<<"*                                                                    *"<<"\n";
    747772  f_out<<"#***************                                                     *"<<"\n";
    748773  f_out<<"# Jet definition                                                     *"<<"\n";
    … …  
    10501075//****************************************************************
    10511076//bool RESOLution::Isolation(const float phi, const float eta,const vector<TLorentzVector> &tracks, const float pt_second_track)
    1052 bool RESOLution::Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float pt_second_track)
     1077bool RESOLution::Isolation(const float phi, const float eta,const vector<TRootTracks> &tracks, const float pt_second_track, const float isolCone)
    10531078{
    10541079   bool isolated = false;
    … …  
    10651090           deltar=genDeltaR;
    10661091      }
    1067    if(deltar > 0.5) isolated = true;
     1092   if(deltar > isolCone) isolated = true;
    10681093   return isolated;
    10691094}
Note: See TracChangeset for help on using the changeset viewer.