Changeset db316cd in git
- Timestamp:
- May 24, 2021, 3:31:22 PM (3 years ago)
- Branches:
- master
- Children:
- 7e6c201, b03c0d0, f76741a
- Parents:
- f8e61b2 (diff), fa0b54e (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the(diff)
links above to see all the changes relative to each parent. - Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/delphes_card_IDEA.tcl
rf8e61b2 rdb316cd 68 68 TauTagging 69 69 70 UniqueObjectFinder71 72 ScalarHT73 70 TreeWriter 74 71 } … … 80 77 module TruthVertexFinder TruthVertexFinder { 81 78 82 ## below this distance two vertices are assumed to be the same79 ## below this distance two vertices are assumed to be merged 83 80 set Resolution 1E-06 84 81 … … 116 113 set InputArray ParticlePropagator/chargedHadrons 117 114 set OutputArray chargedHadrons 115 set UseMomentumVector true 116 118 117 # We use only one efficiency, we set only 0 effincency out of eta bounds: 119 118 120 119 set EfficiencyFormula { 121 120 (abs(eta) > 3.0) * (0.000) + 122 ( energy>= 0.5) * (abs(eta) <= 3.0) * (0.997) +123 ( energy < 0.5 && energy>= 0.3) * (abs(eta) <= 3.0) * (0.65) +124 ( energy< 0.3) * (abs(eta) <= 3.0) * (0.06)121 (pt >= 0.5) * (abs(eta) <= 3.0) * (0.997) + 122 (pt < 0.5 && pt >= 0.3) * (abs(eta) <= 3.0) * (0.65) + 123 (pt < 0.3) * (abs(eta) <= 3.0) * (0.06) 125 124 } 126 125 } … … 135 134 set InputArray ParticlePropagator/electrons 136 135 set OutputArray electrons 137 136 set UseMomentumVector true 138 137 139 138 # Current full simulation with CLICdet provides for electrons: 140 139 set EfficiencyFormula { 141 140 (abs(eta) > 3.0) * (0.000) + 142 ( energy>= 0.5) * (abs(eta) <= 3.0) * (0.997) +143 ( energy < 0.5 && energy>= 0.3) * (abs(eta) <= 3.0) * (0.65) +144 ( energy< 0.3) * (abs(eta) <= 3.0) * (0.06)141 (pt >= 0.5) * (abs(eta) <= 3.0) * (0.997) + 142 (pt < 0.5 && pt >= 0.3) * (abs(eta) <= 3.0) * (0.65) + 143 (pt < 0.3) * (abs(eta) <= 3.0) * (0.06) 145 144 } 146 145 } … … 154 153 set InputArray ParticlePropagator/muons 155 154 set OutputArray muons 155 set UseMomentumVector true 156 156 157 157 # Current full simulation with CLICdet provides for muons: 158 158 set EfficiencyFormula { 159 159 (abs(eta) > 3.0) * (0.000) + 160 ( energy>= 0.5) * (abs(eta) <= 3.0) * (0.997) +161 ( energy < 0.5 && energy>= 0.3) * (abs(eta) <= 3.0) * (0.65) +162 ( energy< 0.3) * (abs(eta) <= 3.0) * (0.06)160 (pt >= 0.5) * (abs(eta) <= 3.0) * (0.997) + 161 (pt < 0.5 && pt >= 0.3) * (abs(eta) <= 3.0) * (0.65) + 162 (pt < 0.3) * (abs(eta) <= 3.0) * (0.06) 163 163 } 164 164 } … … 417 417 418 418 # 0: assume vertex time tV from MC Truth (ideal case) 419 # 1: assume vertex time tV =0419 # 1: assume vertex time tV = 0 420 420 # 2: calculate vertex time as vertex TOF, assuming tPV=0 421 421 422 set VertexTimeMode 2 422 423 } 424 423 } 425 424 426 425 ############## … … 721 720 module FastJetFinder GenJetFinder { 722 721 set InputArray NeutrinoFilter/filteredParticles 723 724 722 set OutputArray jets 725 723 726 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt727 set JetAlgorithm 6728 set Parameter R 0.4724 set JetAlgorithm 10 725 set ParameterR 1.5 726 set ParameterP -1.0 729 727 set JetPTMin 1.0 728 730 729 } 731 730 … … 752 751 753 752 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 754 set JetAlgorithm 6 755 set ParameterR 0.4 753 set JetAlgorithm 10 754 set ParameterR 1.5 755 set ParameterP -1.0 756 756 set JetPTMin 1.0 757 758 757 759 } 758 760 … … 766 768 767 769 # scale formula for jets 768 set ScaleFormula {1.0 8}770 set ScaleFormula {1.00} 769 771 } 770 772 … … 826 828 827 829 828 #####################################################829 # Find uniquely identified photons/electrons/tau/jets830 #####################################################831 832 module UniqueObjectFinder UniqueObjectFinder {833 834 set UseUniqueID true835 # earlier arrays take precedence over later ones836 # add InputArray InputArray OutputArray837 add InputArray PhotonIsolation/photons photons838 add InputArray ElectronIsolation/electrons electrons839 add InputArray MuonIsolation/muons muons840 add InputArray JetEnergyScale/jets jets841 }842 843 830 844 831 ################## … … 854 841 855 842 add Branch Delphes/allParticles Particle GenParticle 843 add Branch TruthVertexFinder/vertices GenVertex Vertex 856 844 857 845 add Branch TrackMerger/tracks Track Track … … 862 850 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower 863 851 852 add Branch EFlowMerger/eflow ParticleFlowCandidate ParticleFlowCandidate 853 864 854 add Branch Calorimeter/photons CaloPhoton Photon 865 855 add Branch PhotonEfficiency/photons PhotonEff Photon … … 869 859 add Branch GenMissingET/momentum GenMissingET MissingET 870 860 871 add Branch UniqueObjectFinder/jets Jet Jet 872 add Branch UniqueObjectFinder/electrons Electron Electron 873 add Branch UniqueObjectFinder/photons Photon Photon 874 add Branch UniqueObjectFinder/muons Muon Muon 875 876 add Branch JetEnergyScale/jets AntiKtJet Jet 861 add Branch JetEnergyScale/jets Jet Jet 862 add Branch ElectronIsolation/electrons Electron Electron 863 add Branch PhotonIsolation/photons Photon Photon 864 add Branch MuonIsolation/muons Muon Muon 877 865 878 866 add Branch MissingET/momentum MissingET MissingET 879 add Branch ScalarHT/energy ScalarHT ScalarHT880 867 881 868 # add Info InfoName InfoValue -
modules/TimeOfFlight.cc
rf8e61b2 rdb316cd 164 164 } 165 165 166 167 //ti = ti - t_truth;168 169 166 p = candidateMomentum.P(); 170 167 … … 190 187 candidate->InitialPosition.SetT(ti * 1.0E3 * c_light); 191 188 192 // update particle mass based on TOF-based PID 193 candidate->Momentum.SetVectM(candidateMomentum.Vect(), mass);189 // update particle mass based on TOF-based PID (commented for now, assume this calculation is done offline) 190 //candidate->Momentum.SetVectM(candidateMomentum.Vect(), mass); 194 191 195 192 candidate->AddCandidate(mother); -
modules/TimeSmearing.cc
rf8e61b2 rdb316cd 103 103 while((candidate = static_cast<Candidate *>(fItTrackInputArray->Next()))) 104 104 { 105 // converting to meters 105 106 106 const TLorentzVector &candidateFinalPosition = candidate->Position; 107 107 const TLorentzVector &candidateMomentum = candidate->Momentum;
Note:
See TracChangeset
for help on using the changeset viewer.