Changeset 839deb7 in git for modules/PileUpJetID.cc
- Timestamp:
- Jun 29, 2015, 10:27:47 PM (9 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- c1ce3fe
- Parents:
- 5d2481f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/PileUpJetID.cc
r5d2481f r839deb7 154 154 float sumT40 = 0.; 155 155 float sumWeightsForT = 0.; 156 candidate->N times = 0;156 candidate->NTimeHits = 0; 157 157 158 158 float sumpt = 0.; … … 197 197 float tow_sumT = 0; 198 198 float tow_sumW = 0; 199 for (int i = 0 ; i < constituent->E cal_E_t.size() ; i++) {200 float w = TMath::Sqrt(constituent->E cal_E_t[i].first);199 for (int i = 0 ; i < constituent->ECalEnergyTimePairs.size() ; i++) { 200 float w = TMath::Sqrt(constituent->ECalEnergyTimePairs[i].first); 201 201 if (fAverageEachTower) { 202 tow_sumT += w*constituent->E cal_E_t[i].second;202 tow_sumT += w*constituent->ECalEnergyTimePairs[i].second; 203 203 tow_sumW += w; 204 204 } else { 205 sumT0 += w*constituent->E cal_E_t[i].second;206 sumT1 += w*gRandom->Gaus(constituent->E cal_E_t[i].second,0.001);207 sumT10 += w*gRandom->Gaus(constituent->E cal_E_t[i].second,0.010);208 sumT20 += w*gRandom->Gaus(constituent->E cal_E_t[i].second,0.020);209 sumT30 += w*gRandom->Gaus(constituent->E cal_E_t[i].second,0.030);210 sumT40 += w*gRandom->Gaus(constituent->E cal_E_t[i].second,0.040);205 sumT0 += w*constituent->ECalEnergyTimePairs[i].second; 206 sumT1 += w*gRandom->Gaus(constituent->ECalEnergyTimePairs[i].second,0.001); 207 sumT10 += w*gRandom->Gaus(constituent->ECalEnergyTimePairs[i].second,0.010); 208 sumT20 += w*gRandom->Gaus(constituent->ECalEnergyTimePairs[i].second,0.020); 209 sumT30 += w*gRandom->Gaus(constituent->ECalEnergyTimePairs[i].second,0.030); 210 sumT40 += w*gRandom->Gaus(constituent->ECalEnergyTimePairs[i].second,0.040); 211 211 sumWeightsForT += w; 212 candidate->N times++;212 candidate->NTimeHits++; 213 213 } 214 214 } … … 221 221 sumT40 += tow_sumW*gRandom->Gaus(tow_sumT/tow_sumW,0.0040); 222 222 sumWeightsForT += tow_sumW; 223 candidate->N times++;223 candidate->NTimeHits++; 224 224 } 225 225 }
Note:
See TracChangeset
for help on using the changeset viewer.