Changeset 29b722a in git
- Timestamp:
- Dec 9, 2021, 8:42:12 PM (3 years ago)
- Children:
- a09b75f
- Parents:
- d612dec
- Files:
-
- 2 added
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
Makefile
rd612dec r29b722a 522 522 classes/DelphesFactory.h \ 523 523 classes/SortableObject.h 524 tmp/classes/DelphesCscClusterFormula.$(ObjSuf): \ 525 classes/DelphesCscClusterFormula.$(SrcSuf) \ 526 classes/DelphesCscClusterFormula.h 524 527 tmp/classes/DelphesCylindricalFormula.$(ObjSuf): \ 525 528 classes/DelphesCylindricalFormula.$(SrcSuf) \ … … 751 754 classes/DelphesClasses.h \ 752 755 classes/DelphesFactory.h \ 753 classes/Delphes Formula.h \756 classes/DelphesCscClusterFormula.h \ 754 757 external/ExRootAnalysis/ExRootClassifier.h \ 755 758 external/ExRootAnalysis/ExRootFilter.h \ … … 759 762 modules/CscClusterId.h \ 760 763 classes/DelphesClasses.h \ 761 762 classes/DelphesFormula.h \764 classes/DelphesFactory.h \ 765 classes/DelphesCscClusterFormula.h \ 763 766 external/ExRootAnalysis/ExRootClassifier.h \ 764 767 external/ExRootAnalysis/ExRootFilter.h \ … … 1194 1197 DELPHES_OBJ += \ 1195 1198 tmp/classes/DelphesClasses.$(ObjSuf) \ 1199 tmp/classes/DelphesCscClusterFormula.$(ObjSuf) \ 1196 1200 tmp/classes/DelphesCylindricalFormula.$(ObjSuf) \ 1197 1201 tmp/classes/DelphesFactory.$(ObjSuf) \ … … 1801 1805 ifeq ($(HAS_PYTHIA8),true) 1802 1806 FASTJET_OBJ += \ 1803 1807 1804 1808 endif 1805 1809 … … 1848 1852 ifeq ($(HAS_PYTHIA8),true) 1849 1853 DISPLAY_OBJ += \ 1850 1854 1851 1855 endif 1852 1856 … … 2594 2598 2595 2599 ### 2596 2597 -
cards/delphes_card_CMS_CSCCluster.tcl
rd612dec r29b722a 65 65 TreeWriter 66 66 } 67 67 68 68 69 ################################# … … 408 409 ################# 409 410 410 module PdgCodeFilter HiggsFilter {411 set InputArray Delphes/allParticles412 set OutputArray higgs413 set Invert true414 add PdgCode {25}415 }416 417 418 419 411 420 412 # filter out LLPs that decay in CSC … … 426 418 # DecayRegion = 1: select LLP that decays in CSC volume 427 419 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation 428 set RequireStatus true429 set Status 2430 add PdgCode {9000005} 420 set RequireStatus false 421 add PdgCode {1500001} 422 431 423 432 424 } … … 436 428 set OutputArray LLP 437 429 set DecayRegion 0 438 set RequireStatus true 439 set Status 2 430 set RequireStatus false 440 431 441 432 # DecayRegion = 0: no cuts on decay region 442 433 # DecayRegion = 1: select LLP that decays in CSC volume 443 434 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation 444 add PdgCode { 9000005}435 add PdgCode {1500001} 445 436 446 437 } … … 455 446 set EfficiencyFormula { 456 447 457 (pt > 3900 && eta < 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0049) + 458 (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0130) + 459 (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.0346) + 460 (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.0623) + 461 (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.0919) + 462 (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.1086) + 463 (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.1292) + 464 (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.1106) + 465 (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0249) + 466 (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.0285) + 467 (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.0501) + 468 (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.0841) + 469 (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.1021) + 470 (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.1129) + 471 (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.1141) + 472 (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.1370) + 473 (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.0282) + 474 (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.0445) + 475 (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.0643) + 476 (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.0903) + 477 (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.0998) + 478 (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.1420) + 479 (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.1429) + 480 (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.0882) + 481 (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.0594) + 482 (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.0521) + 483 (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.0605) + 484 (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.0791) + 485 (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.1117) + 486 (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.0862) + 487 (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.0698) + 488 (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.0500) + 489 (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.0758) + 490 (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.0414) + 491 (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.0755) + 492 (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.1027) + 493 (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.0440) + 494 (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.0811) + 495 (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.1538) + 496 (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.0833) + 497 (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.0300) + 498 (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.0609) + 499 (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.0745) + 500 (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.0610) + 501 (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.1224) + 502 (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.1667) + 503 (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.0000) + 504 (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.0000) + 505 (energy >= 150.0 && 0.0 == phi)*(0.0282) + 506 (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.0809) + 507 (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.0352) + 508 (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.0984) + 509 (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.0968) + 510 (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.1282) + 511 (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.2105) + 512 (energy >= 150.0 && phi >= 200.0)*(0.0769)) + 513 (eta > 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0184) + 514 (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0772) + 515 (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.2086) + 516 (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.3091) + 517 (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.3867) + 518 (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.4500) + 519 (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.4746) + 520 (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.4906) + 521 (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0955) + 522 (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.1461) + 523 (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.2594) + 524 (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.3556) + 525 (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.4165) + 526 (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.4693) + 527 (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.5054) + 528 (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.5219) + 529 (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.1472) + 530 (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.1970) + 531 (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.2974) + 532 (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.3783) + 533 (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.4335) + 534 (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.4736) + 535 (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.4937) + 536 (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.5077) + 537 (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.2053) + 538 (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.2314) + 539 (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.3114) + 540 (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.3799) + 541 (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.4420) + 542 (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.4502) + 543 (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.5348) + 544 (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.5115) + 545 (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.2198) + 546 (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.2404) + 547 (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.3295) + 548 (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.3932) + 549 (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.4327) + 550 (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.4377) + 551 (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.5175) + 552 (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.6087) + 553 (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.2147) + 554 (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.2605) + 555 (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.3442) + 556 (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.3622) + 557 (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.4407) + 558 (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.5168) + 559 (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.5056) + 560 (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.4559) + 561 (energy >= 150.0 && 0.0 == phi)*(0.2824) + 562 (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.2447) + 563 (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.3519) + 564 (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.3772) + 565 (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.4447) + 566 (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.4703) + 567 (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.4460) + 568 (energy >= 150.0 && phi >= 200.0)*(0.4400)) + 569 (pt < 2700 && eta < 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0002) + 570 (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0001) + 571 (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.0006) + 572 (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.0014) + 573 (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.0025) + 574 (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.0046) + 575 (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.0060) + 576 (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.0136) + 577 (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0000) + 578 (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.0000) + 579 (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.0006) + 580 (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.0015) + 581 (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.0033) + 582 (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.0051) + 583 (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.0098) + 584 (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.0146) + 585 (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.0000) + 586 (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.0001) + 587 (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.0003) + 588 (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.0015) + 589 (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.0038) + 590 (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.0052) + 591 (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.0114) + 592 (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.0181) + 593 (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.0000) + 594 (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.0001) + 595 (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.0005) + 596 (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.0022) + 597 (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.0067) + 598 (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.0047) + 599 (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.0113) + 600 (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.0145) + 601 (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.0000) + 602 (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.0001) + 603 (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.0003) + 604 (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.0016) + 605 (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.0110) + 606 (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.0029) + 607 (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.0138) + 608 (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.0000) + 609 (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.0000) + 610 (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.0000) + 611 (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.0000) + 612 (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.0026) + 613 (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.0047) + 614 (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.0085) + 615 (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.0152) + 616 (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.0164) + 617 (energy >= 150.0 && 0.0 == phi)*(0.0000) + 618 (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.0000) + 619 (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.0000) + 620 (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.0000) + 621 (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.0000) + 622 (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.0000) + 623 (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.0080) + 624 (energy >= 150.0 && phi >= 200.0)*(0.0143)) 625 } 448 (decayR > 3900 && decayZ < 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0049) + 449 (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0130) + 450 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0346) + 451 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0623) + 452 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0919) + 453 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1086) + 454 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1292) + 455 (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.1106) + 456 (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0249) + 457 (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0285) + 458 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0501) + 459 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0841) + 460 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1021) + 461 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1129) + 462 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1141) + 463 (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.1370) + 464 (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.0282) + 465 (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0445) + 466 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0643) + 467 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0903) + 468 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0998) + 469 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1420) + 470 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1429) + 471 (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.0882) + 472 (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.0594) + 473 (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0521) + 474 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0605) + 475 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0791) + 476 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1117) + 477 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0862) + 478 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0698) + 479 (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.0500) + 480 (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.0758) + 481 (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0414) + 482 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0755) + 483 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.1027) + 484 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0440) + 485 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0811) + 486 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1538) + 487 (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.0833) + 488 (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.0300) + 489 (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0609) + 490 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0745) + 491 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0610) + 492 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1224) + 493 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1667) + 494 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0000) + 495 (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.0000) + 496 (Eem >= 150.0 && 0.0 == Ehad)*(0.0282) + 497 (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0809) + 498 (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0352) + 499 (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0984) + 500 (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0968) + 501 (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1282) + 502 (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.2105) + 503 (Eem >= 150.0 && Ehad >= 200.0)*(0.0769)) + 504 (decayZ > 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0184) + 505 (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0772) + 506 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2086) + 507 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3091) + 508 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.3867) + 509 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4500) + 510 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4746) + 511 (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.4906) + 512 (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0955) + 513 (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.1461) + 514 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2594) + 515 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3556) + 516 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4165) + 517 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4693) + 518 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5054) + 519 (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.5219) + 520 (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.1472) + 521 (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.1970) + 522 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2974) + 523 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3783) + 524 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4335) + 525 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4736) + 526 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4937) + 527 (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.5077) + 528 (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.2053) + 529 (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2314) + 530 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3114) + 531 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3799) + 532 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4420) + 533 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4502) + 534 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5348) + 535 (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.5115) + 536 (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.2198) + 537 (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2404) + 538 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3295) + 539 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3932) + 540 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4327) + 541 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4377) + 542 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5175) + 543 (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.6087) + 544 (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.2147) + 545 (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2605) + 546 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3442) + 547 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3622) + 548 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4407) + 549 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.5168) + 550 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5056) + 551 (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.4559) + 552 (Eem >= 150.0 && 0.0 == Ehad)*(0.2824) + 553 (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2447) + 554 (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3519) + 555 (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3772) + 556 (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4447) + 557 (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4703) + 558 (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4460) + 559 (Eem >= 150.0 && Ehad >= 200.0)*(0.4400)) + 560 (decayR < 2700 && decayZ < 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0002) + 561 (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) + 562 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0006) + 563 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0014) + 564 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0025) + 565 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0046) + 566 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0060) + 567 (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.0136) + 568 (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0000) + 569 (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) + 570 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0006) + 571 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0015) + 572 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0033) + 573 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0051) + 574 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0098) + 575 (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.0146) + 576 (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.0000) + 577 (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) + 578 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0003) + 579 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0015) + 580 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0038) + 581 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0052) + 582 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0114) + 583 (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.0181) + 584 (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.0000) + 585 (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) + 586 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0005) + 587 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0022) + 588 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0067) + 589 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0047) + 590 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0113) + 591 (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.0145) + 592 (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.0000) + 593 (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) + 594 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0003) + 595 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0016) + 596 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0110) + 597 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0029) + 598 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0138) + 599 (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.0000) + 600 (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.0000) + 601 (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) + 602 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0000) + 603 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0026) + 604 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0047) + 605 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0085) + 606 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0152) + 607 (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.0164) + 608 (Eem >= 150.0 && 0.0 == Ehad)*(0.0000) + 609 (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) + 610 (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0000) + 611 (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0000) + 612 (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0000) + 613 (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0000) + 614 (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0080) + 615 (Eem >= 150.0 && Ehad >= 200.0)*(0.0143)) } 626 616 } 627 617 … … 632 622 # efficiency formula for Csc Cluster, as a function of LLP decay vertex in R, Z and hadronic and EM energy 633 623 set EfficiencyFormula { 634 ( pt > 3900 && eta < 6710) * ((0.0 == phi)*(0.0656) +635 ( phi > 0.0 && phi< 25.0)*(0.0777) +636 ( phi >= 25.0 && phi< 50.0)*(0.1607) +637 ( phi >= 50.0 && phi< 75.0)*(0.2294) +638 ( phi >= 75.0 && phi< 100.0)*(0.3146) +639 ( phi >= 100.0 && phi< 150.0)*(0.3107) +640 ( phi >= 150.0 && phi< 200.0)*(0.3177) +641 ( phi>= 200.0)*(0.3229)) +642 ( eta > 6710) * ((0.0 == phi)*(0.2987) +643 ( phi > 0.0 && phi< 25.0)*(0.3100) +644 ( phi >= 25.0 && phi< 50.0)*(0.4476) +645 ( phi >= 50.0 && phi< 75.0)*(0.5335) +646 ( phi >= 75.0 && phi< 100.0)*(0.5961) +647 ( phi >= 100.0 && phi< 150.0)*(0.6368) +648 ( phi >= 150.0 && phi< 200.0)*(0.6814) +649 ( phi>= 200.0)*(0.6998)) +650 ( pt < 2700 && eta < 6710) * ((0.0 == phi)*(0.8604) +651 ( phi > 0.0 && phi< 25.0)*(0.3335) +652 ( phi >= 25.0 && phi< 50.0)*(0.2457) +653 ( phi >= 50.0 && phi< 75.0)*(0.1831) +654 ( phi >= 75.0 && phi< 100.0)*(0.2100) +655 ( phi >= 100.0 && phi< 150.0)*(0.2443) +656 ( phi >= 150.0 && phi< 200.0)*(0.2532) +657 ( phi>= 200.0)*(0.2404))624 (decayR > 3900 && decayZ < 6710) * ((0.0 == Ehad)*(0.0656) + 625 (Ehad > 0.0 && Ehad < 25.0)*(0.0777) + 626 (Ehad >= 25.0 && Ehad < 50.0)*(0.1607) + 627 (Ehad >= 50.0 && Ehad < 75.0)*(0.2294) + 628 (Ehad >= 75.0 && Ehad < 100.0)*(0.3146) + 629 (Ehad >= 100.0 && Ehad < 150.0)*(0.3107) + 630 (Ehad >= 150.0 && Ehad < 200.0)*(0.3177) + 631 (Ehad >= 200.0)*(0.3229)) + 632 (decayZ > 6710) * ((0.0 == Ehad)*(0.2987) + 633 (Ehad > 0.0 && Ehad < 25.0)*(0.3100) + 634 (Ehad >= 25.0 && Ehad < 50.0)*(0.4476) + 635 (Ehad >= 50.0 && Ehad < 75.0)*(0.5335) + 636 (Ehad >= 75.0 && Ehad < 100.0)*(0.5961) + 637 (Ehad >= 100.0 && Ehad < 150.0)*(0.6368) + 638 (Ehad >= 150.0 && Ehad < 200.0)*(0.6814) + 639 (Ehad >= 200.0)*(0.6998)) + 640 (decayR < 2700 && decayZ < 6710) * ((0.0 == Ehad)*(0.8604) + 641 (Ehad > 0.0 && Ehad < 25.0)*(0.3335) + 642 (Ehad >= 25.0 && Ehad < 50.0)*(0.2457) + 643 (Ehad >= 50.0 && Ehad < 75.0)*(0.1831) + 644 (Ehad >= 75.0 && Ehad < 100.0)*(0.2100) + 645 (Ehad >= 100.0 && Ehad < 150.0)*(0.2443) + 646 (Ehad >= 150.0 && Ehad < 200.0)*(0.2532) + 647 (Ehad >= 200.0)*(0.2404)) 658 648 } 659 649 } … … 1055 1045 #add Branch GenJetFinder/jets GenJet Jet 1056 1046 #add Branch GenMissingET/momentum GenMissingET MissingET 1057 add Branch HiggsFilter/higgs higgs GenParticle1058 1047 1059 1048 … … 1071 1060 add Branch CSCFilter/LLP Cscllp CscCluster 1072 1061 1073 add Branch ClusterEfficiency/cluster CscCluster 130CscCluster1074 } 1062 add Branch ClusterEfficiency/cluster CscCluster CscCluster 1063 } -
classes/DelphesHepMC2Reader.cc
rd612dec r29b722a 438 438 { 439 439 Candidate *candidate; 440 Candidate *candidateDaughter; 440 441 map<int, pair<int, int> >::iterator itMotherMap; 441 442 map<int, pair<int, int> >::iterator itDaughterMap; … … 446 447 candidate = static_cast<Candidate *>(allParticleOutputArray->At(i)); 447 448 449 448 450 if(candidate->M1 > 0) 449 451 { … … 477 479 candidate->D1 = -1; 478 480 candidate->D2 = -1; 481 const TLorentzVector &decayPosition = candidate->Position; 482 candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position 479 483 } 480 484 else … … 482 486 candidate->D1 = itDaughterMap->second.first; 483 487 candidate->D2 = itDaughterMap->second.second; 484 } 485 } 486 } 487 } 488 489 //--------------------------------------------------------------------------- 488 candidateDaughter = static_cast<Candidate *>(allParticleOutputArray->At(candidate->D1)); 489 const TLorentzVector &decayPosition = candidateDaughter->Position; 490 candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position 491 492 } 493 } 494 } 495 } 496 497 //--------------------------------------------------------------------------- -
classes/DelphesHepMC3Reader.cc
rd612dec r29b722a 458 458 TObjArray *array; 459 459 Candidate *candidate; 460 Candidate *candidateDaughter; 460 461 TParticlePDG *pdgParticle; 461 462 int pdgCode; … … 573 574 candidate->D1 = -1; 574 575 candidate->D2 = -1; 576 const TLorentzVector &decayPosition = candidate->Position; 577 candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position 575 578 } 576 579 else … … 578 581 candidate->D1 = itDaughterMap->second.first; 579 582 candidate->D2 = itDaughterMap->second.second; 580 } 581 } 582 } 583 } 584 585 //--------------------------------------------------------------------------- 583 candidateDaughter = static_cast<Candidate *>(allParticleOutputArray->At(candidate->D1)); 584 const TLorentzVector &decayPosition = candidateDaughter->Position; 585 candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position 586 } 587 } 588 } 589 } 590 591 //--------------------------------------------------------------------------- -
modules/CscClusterEfficiency.cc
rd612dec r29b722a 29 29 #include "classes/DelphesClasses.h" 30 30 #include "classes/DelphesFactory.h" 31 #include "classes/Delphes Formula.h"31 #include "classes/DelphesCscClusterFormula.h" 32 32 33 33 #include "ExRootAnalysis/ExRootClassifier.h" … … 55 55 fFormula(0), fItInputArray(0) 56 56 { 57 fFormula = new Delphes Formula;57 fFormula = new DelphesCscClusterFormula; 58 58 } 59 59 … … 70 70 { 71 71 // read CscClusterEfficiency formula 72 72 cout<<"initialize"<<endl; 73 73 fFormula->Compile(GetString("EfficiencyFormula", "1.0")); 74 74 … … 106 106 Eem = candidate->Eem; 107 107 // apply an efficency formula 108 //pt, eta, phi, energy109 108 if(gRandom->Uniform() > fFormula->Eval(decayR, decayZ, Ehad, Eem)) continue; 109 110 110 111 111 fOutputArray->Add(candidate); -
modules/CscClusterEfficiency.h
rd612dec r29b722a 32 32 class TIterator; 33 33 class TObjArray; 34 class Delphes Formula;34 class DelphesCscClusterFormula; 35 35 36 36 class CscClusterEfficiency: public DelphesModule … … 45 45 46 46 private: 47 Delphes Formula *fFormula; //!47 DelphesCscClusterFormula *fFormula; //! 48 48 49 49 TIterator *fItInputArray; //! -
modules/CscClusterId.cc
rd612dec r29b722a 29 29 #include "classes/DelphesClasses.h" 30 30 #include "classes/DelphesFactory.h" 31 #include "classes/Delphes Formula.h"31 #include "classes/DelphesCscClusterFormula.h" 32 32 33 33 #include "ExRootAnalysis/ExRootClassifier.h" … … 55 55 fFormula(0), fItInputArray(0) 56 56 { 57 fFormula = new Delphes Formula;57 fFormula = new DelphesCscClusterFormula; 58 58 } 59 59 … … 113 113 eta = (cosTheta == 1.0 ? signPz * 999.9 : momentum.Eta()); 114 114 115 NStationEff = fFormula->Eval(decayR, decayZ, Ehad); //pt is used as argument in Delphes Formula115 NStationEff = fFormula->Eval(decayR, decayZ, Ehad); //pt is used as argument in DelphesCscClusterFormula 116 116 117 117 // assign average station for the cluster -
modules/CscClusterId.h
rd612dec r29b722a 32 32 class TIterator; 33 33 class TObjArray; 34 class Delphes Formula;34 class DelphesCscClusterFormula; 35 35 36 36 class CscClusterId: public DelphesModule … … 45 45 46 46 private: 47 Delphes Formula *fFormula; //!47 DelphesCscClusterFormula *fFormula; //! 48 48 49 49 TIterator *fItInputArray; //!
Note:
See TracChangeset
for help on using the changeset viewer.