- Timestamp:
- Nov 20, 2008, 11:28:51 AM (16 years ago)
- Location:
- trunk
- Files:
-
- 2 edited
-
Delphes.cpp (modified) (3 diffs)
-
src/SmearUtil.cc (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/Delphes.cpp
r47 r49 252 252 case 4: { 253 253 jet_def = fastjet::JetDefinition(fastjet::kt_algorithm, DET->CONERADIUS); 254 //jet_defs[4] = fastjet::JetDefinition(fastjet::cambridge_algorithm,jet_radius); 255 //jet_defs[5] = fastjet::JetDefinition(fastjet::antikt_algorithm,jet_radius); 254 } 255 break; 256 257 case 5: { 258 jet_def = fastjet::JetDefinition(fastjet::cambridge_algorithm,DET->CONERADIUS); 259 } 260 break; 261 262 case 6: { 263 jet_def = fastjet::JetDefinition(fastjet::antikt_algorithm,DET->CONERADIUS); 256 264 } 257 265 break; … … 375 383 TrackCentral.push_back(genMomentum); 376 384 } 377 378 } // switch385 386 } // switch 379 387 380 388 // Forward particles in CASTOR ? … … 504 512 // Tau jet identification : 1! track and electromagnetic collimation 505 513 if(fabs(JET.Eta()) < (DET->MAX_TRACKER - DET->TAU_CONE_TRACKS)) { 506 507 508 509 510 511 512 513 514 515 516 514 double Energie_tau_central = DET->EnergySmallCone(towers,JET.Eta(),JET.Phi()); 515 if( 516 ( Energie_tau_central/JET.E() > DET->TAU_EM_COLLIMATION ) && 517 ( DET->NumTracks(TrackCentral,DET->PT_TRACK_TAU,JET.Eta(),JET.Phi()) == 1 ) && 518 ( JET.Pt() > DET->TAUJET_pt) 519 ) { 520 elementTauJet = (TRootTauJet*) branchTauJet->NewEntry(); 521 elementTauJet->Set(JET); 522 } // if tau jet 523 } // if JET.eta < tracker - tau_cone : Tau jet identification 524 517 525 if(JET.Pt() > DET->JET_pt) 518 526 { -
trunk/src/SmearUtil.cc
r48 r49 353 353 f_out<<"#*************** *"<<"\n"; 354 354 f_out<<"* *"<<"\n"; 355 f_out<<"* Six algorithms are currently available: *"<<"\n"; 356 f_out<<"* - 1) CDF cone algorithm, *"<<"\n"; 357 f_out<<"* - 2) CDF MidPoint algorithm, *"<<"\n"; 358 f_out<<"* - 3) SIScone algorithm, *"<<"\n"; 359 f_out<<"* - 4) kt algorithm, *"<<"\n"; 360 f_out<<"* - 5) Cambrigde/Aachen algorithm, *"<<"\n"; 361 f_out<<"* - 6) Anti-kt algorithm. *"<<"\n"; 362 f_out<<"* *"<<"\n"; 363 f_out<<"* You have chosen *"<<"\n"; 355 364 switch(JETALGO) { 356 365 default: 357 366 case 1: { 358 f_out<<"* CDF JetClu jet algorithm with:*"<<"\n";359 360 361 362 363 364 365 366 367 368 369 370 367 f_out<<"* CDF JetClu jet algorithm with parameters: *"<<"\n"; 368 f_out << left << setw(40) <<"* - Seed threshold: "<<"" 369 << left << setw(10) <<SEEDTHRESHOLD <<""<< right << setw(20)<<"! not in datacard *"<<"\n"; 370 f_out << left << setw(40) <<"* - Cone radius: "<<"" 371 << left << setw(10) <<CONERADIUS <<""<< right << setw(20)<<"*"<<"\n"; 372 f_out << left << setw(40) <<"* - Adjacency cut: "<<"" 373 << left << setw(10) <<C_ADJACENCYCUT <<""<< right << setw(20)<<"! not in datacard *"<<"\n"; 374 f_out << left << setw(40) <<"* - Max iterations: "<<"" 375 << left << setw(10) <<C_MAXITERATIONS <<""<< right << setw(20)<<"! not in datacard *"<<"\n"; 376 f_out << left << setw(40) <<"* - Iratch: "<<"" 377 << left << setw(10) <<C_IRATCH <<""<< right << setw(20)<<"! not in datacard *"<<"\n"; 378 f_out << left << setw(40) <<"* - Overlap threshold: "<<"" 379 << left << setw(10) <<OVERLAPTHRESHOLD <<""<< right << setw(20)<<"! not in datacard *"<<"\n"; 371 380 } 372 381 break; 373 382 case 2: { 374 f_out<<"* CDF midpoint jet algorithm with:*"<<"\n";375 376 377 378 379 380 381 382 383 384 385 386 383 f_out<<"* CDF midpoint jet algorithm with parameters: *"<<"\n"; 384 f_out << left << setw(40) <<"* - Seed threshold: "<<"" 385 << left << setw(20) <<SEEDTHRESHOLD <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 386 f_out << left << setw(40) <<"* - Cone radius: "<<"" 387 << left << setw(20) <<CONERADIUS <<""<< right << setw(10)<<"*"<<"\n"; 388 f_out << left << setw(40) <<"* - Cone area fraction:"<<"" 389 << left << setw(20) <<M_CONEAREAFRACTION <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 390 f_out << left << setw(40) <<"* - Maximum pair size: "<<"" 391 << left << setw(20) <<M_MAXPAIRSIZE <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 392 f_out << left << setw(40) <<"* - Max iterations: "<<"" 393 << left << setw(20) <<M_MAXITERATIONS <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 394 f_out << left << setw(40) <<"* - Overlap threshold: "<<"" 395 << left << setw(20) <<OVERLAPTHRESHOLD <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 387 396 } 388 397 break; 389 398 case 3: { 390 f_out<<"* SISCone jet algorithm with:*"<<"\n";391 392 393 394 395 396 397 398 399 f_out <<"* SISCone jet algorithm with parameters: *"<<"\n"; 400 f_out << left << setw(40) <<"* - Cone radius: "<<"" 401 << left << setw(20) <<CONERADIUS <<""<< right << setw(10)<<"*"<<"\n"; 402 f_out << left << setw(40) <<"* - Overlap threshold: "<<"" 403 << left << setw(20) <<OVERLAPTHRESHOLD <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 404 f_out << left << setw(40) <<"* - Number pass max: "<<"" 405 << left << setw(20) <<NPASS <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 406 f_out << left << setw(40) <<"* - Minimum pT for protojet: "<<"" 407 << left << setw(20) <<PROTOJET_PTMIN <<""<< right << setw(10)<<"! not in datacard *"<<"\n"; 399 408 } 400 409 break; 401 410 case 4: { 402 f_out<<"* KT jet algorithm with:*"<<"\n";403 404 411 f_out <<"* KT jet algorithm with parameters: *"<<"\n"; 412 f_out << left << setw(40) <<"* - Cone radius: "<<"" 413 << left << setw(20) <<CONERADIUS <<""<< right << setw(10)<<"*"<<"\n"; 405 414 } 406 415 break; 416 case 5: { 417 f_out <<"* Cambridge/Aachen jet algorithm with parameters: *"<<"\n"; 418 f_out << left << setw(40) <<"* - Cone radius: "<<"" 419 << left << setw(20) <<CONERADIUS <<""<< right << setw(10)<<"*"<<"\n"; 420 } 421 break; 422 case 6: { 423 f_out <<"* Anti-kt jet algorithm with parameters: *"<<"\n"; 424 f_out << left << setw(40) <<"* - Cone radius: "<<"" 425 << left << setw(20) <<CONERADIUS <<""<< right << setw(10)<<"*"<<"\n"; 426 } 427 break; 428 429 407 430 } 408 431 f_out<<"* *"<<"\n";
Note:
See TracChangeset
for help on using the changeset viewer.
