Changeset 77249aa in git for cards/FCC/FCChh_PileUpVtx.tcl
- Timestamp:
- Jan 26, 2020, 4:53:05 PM (5 years ago)
- Branches:
- Timing
- Children:
- 6777565
- Parents:
- f59b36de
- File:
-
- 1 edited
-
cards/FCC/FCChh_PileUpVtx.tcl (modified) (17 diffs)
Legend:
- Unmodified
- Added
- Removed
-
cards/FCC/FCChh_PileUpVtx.tcl
rf59b36de r77249aa 31 31 32 32 TrackMerger 33 34 33 TrackSmearing 35 34 36 35 ECal 37 36 HCal 37 38 TimeSmearingMIP 39 TimeSmearingPhotons 40 TimeSmearingNH 41 42 TimeSmearingEcal 43 TimeSmearingHcal 44 45 VertexFinderDA4D 46 HighMassVertexRecover 47 PileUpSubtractor4D 48 49 ElectronFilter 50 MuonFilter 51 ChargedHadronFilter 38 52 39 53 Calorimeter … … 41 55 EFlowFilter 42 56 43 TimeSmearingMIP44 TimeSmearingPhotons45 TimeSmearingNH46 47 VertexFinderDA4D48 PileUpSubtractor4D49 50 HighMassVertexRecover51 52 57 PhotonEfficiency 53 58 PhotonIsolation 54 59 55 ElectronFilter56 60 ElectronIsolation 57 58 ChargedHadronFilter59 61 60 62 MuonIsolation … … 256 258 set ActiveFraction 0.013 257 259 set ChargeCollectionEfficiency 0.75 258 260 259 261 # fixme: this number should probably be charge/energy dependent, or absolute number in MeV/cm? 260 262 set Resolution 0.15 … … 262 264 # active material properties (cf. http://pdg.lbl.gov/2014/AtomicNuclearProperties/properties8.dat) 263 265 set Z 14. 264 set A 28.0855 265 set rho 2.329 266 266 set A 28.0855 267 set rho 2.329 268 267 269 # material polarisation correction parameters 268 270 set a 0.1492 … … 270 272 set x0 0.2015 271 273 set x1 2.8716 272 set I 173.0 274 set I 173.0 273 275 set c0 4.4355 274 276 … … 298 300 set ApplyToPileUp true 299 301 300 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html301 source trackResolutionFCChh.tcl 302 # this is CMS resolution is d0 , dZ on top of the smearing in P applied before 303 source trackResolutionFCChh.tcl 302 304 # FIXME !!!! we need to add track resolution of FCC-hh baseline detector !!!!! 303 305 } … … 312 314 module SimpleCalorimeter ECal { 313 315 set ParticleInputArray ParticlePropagator/stableParticles 314 set TrackInputArray T imeSmearing/tracks316 set TrackInputArray TrackSmearing/tracks 315 317 316 318 set TowerOutputArray ecalTowers … … 469 471 } 470 472 471 #################472 # Electron filter473 #################474 475 module PdgCodeFilter ElectronFilter {476 set InputArray HCal/eflowTracks477 set OutputArray electrons478 set Invert true479 add PdgCode {11}480 add PdgCode {-11}481 }482 483 484 ######################485 # ChargedHadronFilter486 ######################487 488 module PdgCodeFilter ChargedHadronFilter {489 set InputArray HCal/eflowTracks490 set OutputArray chargedHadrons491 492 add PdgCode {11}493 add PdgCode {-11}494 add PdgCode {13}495 add PdgCode {-13}496 }497 498 499 ###################################################500 # Tower Merger (in case not using e-flow algorithm)501 ###################################################502 503 module Merger Calorimeter {504 # add InputArray InputArray505 add InputArray ECal/ecalTowers506 add InputArray HCal/hcalTowers507 add InputArray MuonMomentumSmearing/muons508 set OutputArray towers509 }510 511 ####################512 # Energy flow merger513 ####################514 515 module Merger EFlowMerger {516 # add InputArray InputArray517 add InputArray HCal/eflowTracks518 add InputArray ECal/eflowPhotons519 add InputArray HCal/eflowNeutralHadrons520 set OutputArray eflow521 }522 523 ######################524 # EFlowFilter525 ######################526 527 module PdgCodeFilter EFlowFilter {528 set InputArray EFlowMerger/eflow529 set OutputArray eflow530 531 add PdgCode {11}532 add PdgCode {-11}533 add PdgCode {13}534 add PdgCode {-13}535 }536 473 537 474 ######################################## … … 543 480 set OutputArray tracks 544 481 545 # assume 30 ps resolution for now 546 set TimeResolution {30E-12} 482 # assume constant 30 ps resolution for now 483 set TimeResolution { 484 (abs(eta) > 0.0 && abs(eta) <= 6.0)* 30E-12 485 } 547 486 } 548 487 … … 554 493 set InputArray ECal/eflowPhotons 555 494 set OutputArray photons 556 set TimeResolution {sqrt(20^2 + 150^2)/energy^2} 495 set TimeResolution { 496 (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2 497 } 557 498 } 558 499 … … 565 506 set OutputArray neutralhadrons 566 507 567 # assume 30 ps resolution for now 568 set TimeResolution {sqrt(20^2 + 150^2)/energy^2} 569 } 508 # assume resolution for now 509 set TimeResolution { 510 (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2 511 } 512 } 513 514 ######################################## 515 # Time Smearing Neutral Ecal 516 ######################################## 517 518 module TimeSmearing TimeSmearingEcal { 519 set InputArray ECal/ecalTowers 520 set OutputArray towers 521 set TimeResolution { 522 (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2 523 } 524 } 525 526 527 ######################################## 528 # Time Smearing Neutral Hcal 529 ######################################## 530 531 module TimeSmearing TimeSmearingHcal { 532 set InputArray HCal/hcalTowers 533 set OutputArray towers 534 set TimeResolution { 535 (abs(eta) > 0.0 && abs(eta) <= 6.0) * sqrt(20e-12^2 + 150e-12^2)/energy^2 536 } 537 } 538 539 570 540 571 541 … … 573 543 # Primary vertex reconstruction 574 544 ################################## 575 576 545 577 546 module VertexFinderDA4D VertexFinderDA4D { … … 602 571 } 603 572 604 ##########################605 # Track pile-up subtractor606 ##########################607 608 module PileUpSubtractor4D PileUpSubtractor4D {609 # add InputArray InputArray OutputArray610 611 add InputArray TimeSmearingMIP/tracks612 add InputArray TimeSmearingPhotons/photons613 add InputArray TimeSmearingNH/neutralhadrons614 615 set VertexInputArray VertexFinderDA4D/vertices616 617 set fChargedMinSignificance {3}618 set fNeutralMinSignificance {3}619 }620 621 573 ###################################### 622 574 # Heavy(slow) particles vertex recover … … 634 586 635 587 } 588 589 ########################## 590 # PileUpSubtractor4D 591 ########################## 592 593 module PileUpSubtractor4D PileUpSubtractor4D { 594 # add InputArray InputArray OutputArray 595 596 add InputArray TimeSmearingMIP/tracks tracks 597 add InputArray TimeSmearingPhotons/photons photons 598 add InputArray TimeSmearingNH/neutralhadrons neutralhadrons 599 add InputArray TimeSmearingEcal/towers towers 600 add InputArray TimeSmearingHcal/towers towers 601 602 set VertexInputArray VertexFinderDA4D/vertices 603 604 set fChargedMinSignificance 3.0 605 set fNeutralMinSignificance 3.0 606 } 607 608 609 ################# 610 # Electron filter 611 ################# 612 613 module PdgCodeFilter ElectronFilter { 614 set InputArray PileUpSubtractor4D/tracks 615 set OutputArray electrons 616 set Invert true 617 add PdgCode {11} 618 add PdgCode {-11} 619 } 620 621 ################# 622 # Muon filter 623 ################# 624 625 module PdgCodeFilter MuonFilter { 626 set InputArray PileUpSubtractor4D/tracks 627 set OutputArray muons 628 set Invert true 629 add PdgCode {13} 630 add PdgCode {-13} 631 } 632 633 634 ###################### 635 # ChargedHadronFilter 636 ###################### 637 638 module PdgCodeFilter ChargedHadronFilter { 639 set InputArray PileUpSubtractor4D/tracks 640 set OutputArray chargedHadrons 641 642 add PdgCode {11} 643 add PdgCode {-11} 644 add PdgCode {13} 645 add PdgCode {-13} 646 } 647 648 649 ################################################### 650 # Tower Merger (in case not using e-flow algorithm) 651 ################################################### 652 653 module Merger Calorimeter { 654 # add InputArray InputArray 655 add InputArray PileUpSubtractor4D/towers 656 add InputArray PileUpSubtractor4D/towers 657 add InputArray MuonFilter/muons 658 set OutputArray towers 659 } 660 661 #################### 662 # Energy flow merger 663 #################### 664 665 module Merger EFlowMerger { 666 # add InputArray InputArray 667 add InputArray PileUpSubtractor4D/tracks 668 add InputArray PileUpSubtractor4D/photons 669 add InputArray PileUpSubtractor4D/neutralhadrons 670 set OutputArray eflow 671 } 672 673 ############## 674 # EFlowFilter 675 ############## 676 677 module PdgCodeFilter EFlowFilter { 678 set InputArray EFlowMerger/eflow 679 set OutputArray eflow 680 681 add PdgCode {11} 682 add PdgCode {-11} 683 add PdgCode {13} 684 add PdgCode {-13} 685 } 686 687 636 688 637 689 ################### … … 644 696 set MomentumOutputArray momentum 645 697 } 646 647 648 698 649 699 ################## … … 1065 1115 add Branch GenMissingET/momentum GenMissingET MissingET 1066 1116 1067 add Branch T imeSmearing/tracks Track Track1117 add Branch TrackSmearing/tracks Track Track 1068 1118 add Branch Calorimeter/towers Tower Tower 1069 1119 1070 add Branch HCal/eflowTracks EFlowTrack Track1071 add Branch ECal/eflowPhotons EFlowPhoton Tower1072 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower1120 add Branch PileUpSubtractor4D/tracks EFlowTrack Track 1121 add Branch PileUpSubtractor4D/photons EFlowPhoton Tower 1122 add Branch PileUpSubtractor4D/neutralhadrons EFlowNeutralHadron Tower 1073 1123 1074 1124 add Branch UniqueObjectFinder/photons Photon Photon … … 1088 1138 add Branch HighMassVertexRecover/tracks Track Track 1089 1139 } 1090
Note:
See TracChangeset
for help on using the changeset viewer.
