Changeset a09b75f in git
- Timestamp:
- Dec 9, 2021, 10:38:13 PM (3 years ago)
- Children:
- c61b5ce
- Parents:
- 29b722a
- Location:
- modules
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/CscClusterEfficiency.cc
r29b722a ra09b75f 70 70 { 71 71 // read CscClusterEfficiency formula 72 cout<<"initialize"<<endl;73 72 fFormula->Compile(GetString("EfficiencyFormula", "1.0")); 74 73 -
modules/CscClusterId.cc
r29b722a ra09b75f 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 DelphesCscClusterFormula115 NStationEff = fFormula->Eval(decayR, decayZ, Ehad); 116 116 117 117 // assign average station for the cluster -
modules/LLPFilter.cc
r29b722a ra09b75f 151 151 if (fDaughterNumber > 0) 152 152 { 153 if (candidate->D2-candidate->D1 != fDaughterNumber) continue;//require 3daughters153 if (candidate->D2-candidate->D1 != fDaughterNumber) continue;//require at least fDaughterNumber daughters 154 154 155 155 } … … 160 160 candidate->Eem = 0.0; 161 161 candidate->Ehad = 0.0; 162 163 162 fItParticleInputArray->Reset(); 164 163 … … 166 165 { 167 166 168 // check mother is LLP169 //check ID is 123456 or neutrinos or charged leptons170 167 daughterPdg = daughter->PID; 171 168 if (daughter->Status != 1)continue; 172 173 169 if (daughter->IsPU)continue; 174 175 170 if (abs(daughterPdg)==12 || abs(daughterPdg)==14 || abs(daughterPdg)==16 || abs(daughterPdg)==13)continue; // ignore neutrinos and muons 176 177 178 171 if (abs(daughterPdg) > 1000000) continue;//ignore BSM particles 172 179 173 const TLorentzVector &daughterProdPosition = daughter->Position; 180 174 const TLorentzVector &daughterMomentum = daughter->Momentum; 181 182 175 const TLorentzVector distance = daughterProdPosition - candidateDecayPosition; 183 176 184 if (sqrt(pow(distance.X(), 2) + pow(distance.Y(), 2) + pow(distance.Z(), 2))>100) continue; //if vertices are close, then matched, this 1 m only works for those decay in muon system, would be too large for tracker related177 // look for mother until find LLP or reach the top of the tree 185 178 tempCandidate = daughter; 186 179 while(tempCandidate->M1 != -1 && tempCandidate->M1 != index) 187 180 { 188 181 tempCandidate = static_cast<Candidate *>(fParticleInputArray->At(tempCandidate->M1)); 189 // if (tempCandidate->PID == pdgCode)break;190 191 182 } 192 183 if (tempCandidate->M1 == -1) continue; … … 194 185 if (abs(daughterPdg)==11 || abs(daughterPdg)==22 || abs(daughterPdg)==111)candidate->Eem += daughterMomentum.E(); 195 186 else candidate->Ehad += daughterMomentum.E(); 196 197 198 187 } 199 188
Note:
See TracChangeset
for help on using the changeset viewer.