- Timestamp:
- Nov 19, 2008, 2:06:15 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SmearUtil.cc
r33 r43 80 80 CONERADIUS = 0.7; // generic jet radius ; not for tau's !!! 81 81 JETALGO = 1; // 1 for Cone algorithm, 2 for MidPoint algorithm, 3 for SIScone algorithm, 4 for kt algorithm 82 83 //General jet parameters 84 SEEDTHRESHOLD = 1.0; 85 OVERLAPTHRESHOLD = 0.75; 86 82 87 // Define Cone algorithm. 83 C_SEEDTHRESHOLD = 1.0;84 88 C_ADJACENCYCUT = 2; 85 89 C_MAXITERATIONS = 100; 86 90 C_IRATCH = 1; 87 C_OVERLAPTHRESHOLD = 0.75;88 91 89 92 //Define MidPoint algorithm. 90 M_SEEDTHRESHOLD = 1.0;91 93 M_CONEAREAFRACTION = 0.25; 92 94 M_MAXPAIRSIZE = 2; 93 95 M_MAXITERATIONS = 100; 94 M_OVERLAPTHRESHOLD = 0.75;95 96 96 97 } … … 149 150 } 150 151 152 // General jet variables 153 SEEDTHRESHOLD = 1.0; 154 OVERLAPTHRESHOLD = 0.75; 155 151 156 // Define Cone algorithm. 152 C_SEEDTHRESHOLD = 1.0;153 157 C_ADJACENCYCUT = 2; 154 158 C_MAXITERATIONS = 100; 155 159 C_IRATCH = 1; 156 C_OVERLAPTHRESHOLD = 0.75;157 160 158 161 //Define MidPoint algorithm. 159 M_SEEDTHRESHOLD = 1.0;160 162 M_CONEAREAFRACTION = 0.25; 161 163 M_MAXPAIRSIZE = 2; 162 164 M_MAXITERATIONS = 100; 163 M_OVERLAPTHRESHOLD = 0.75;164 165 165 166 } … … 250 251 double Energie=0; 251 252 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; 253 254 if((DeltaR(phi,eta,towers[i].fourVector.phi(),towers[i].fourVector.eta()) < TAU_CONE_ENERGY)) { 254 255 Energie += towers[i].fourVector.E;
Note:
See TracChangeset
for help on using the changeset viewer.