Fork me on GitHub

Changeset 43 in svn for trunk/src


Ignore:
Timestamp:
Nov 19, 2008, 2:06:15 PM (16 years ago)
Author:
severine ovyn
Message:

general jet variable

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/SmearUtil.cc

    r33 r43  
    8080CONERADIUS   = 0.7;     // generic jet radius ; not for tau's !!!
    8181JETALGO      =   1;     // 1 for Cone algorithm, 2 for MidPoint algorithm, 3 for SIScone algorithm, 4 for kt algorithm
     82
     83//General jet parameters
     84SEEDTHRESHOLD     = 1.0;
     85OVERLAPTHRESHOLD = 0.75;
     86
    8287// Define Cone algorithm.
    83 C_SEEDTHRESHOLD     = 1.0;
    8488C_ADJACENCYCUT     = 2;
    8589C_MAXITERATIONS    = 100;
    8690C_IRATCH           = 1;
    87 C_OVERLAPTHRESHOLD = 0.75;
    8891
    8992//Define MidPoint algorithm.
    90 M_SEEDTHRESHOLD     =   1.0;
    9193M_CONEAREAFRACTION  =   0.25;
    9294M_MAXPAIRSIZE       =   2;
    9395M_MAXITERATIONS     =   100;
    94 M_OVERLAPTHRESHOLD  =  0.75;
    9596
    9697}
     
    149150  }
    150151
     152// General jet variables
     153 SEEDTHRESHOLD     = 1.0;
     154 OVERLAPTHRESHOLD = 0.75;
     155
    151156// Define Cone algorithm.
    152  C_SEEDTHRESHOLD     = 1.0;
    153157 C_ADJACENCYCUT     = 2;
    154158 C_MAXITERATIONS    = 100;
    155159 C_IRATCH           = 1;
    156  C_OVERLAPTHRESHOLD = 0.75;
    157160
    158161//Define MidPoint algorithm.
    159  M_SEEDTHRESHOLD     =   1.0;
    160162 M_CONEAREAFRACTION  =   0.25;
    161163 M_MAXPAIRSIZE       =   2;
    162164 M_MAXITERATIONS     =   100;
    163  M_OVERLAPTHRESHOLD  =  0.75;
    164165
    165166}
     
    250251  double Energie=0;
    251252  for(unsigned int i=0; i < towers.size(); i++) {
    252       if(towers[i].fourVector.pt() < M_SEEDTHRESHOLD) continue;
     253      if(towers[i].fourVector.pt() < SEEDTHRESHOLD) continue;
    253254      if((DeltaR(phi,eta,towers[i].fourVector.phi(),towers[i].fourVector.eta()) < TAU_CONE_ENERGY)) {
    254255          Energie += towers[i].fourVector.E;
Note: See TracChangeset for help on using the changeset viewer.