Changeset 5eda6767 in git
- Timestamp:
- Mar 18, 2021, 6:28:10 PM (4 years ago)
- Branches:
- master
- Children:
- cb80e6f
- Parents:
- d1678fd
- Location:
- modules
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Calorimeter.cc
rd1678fd r5eda6767 231 231 fItParticleInputArray->Reset(); 232 232 number = -1; 233 fTowerRmax=0.; 233 234 while((particle = static_cast<Candidate *>(fItParticleInputArray->Next()))) 234 235 { 235 236 const TLorentzVector &particlePosition = particle->Position; 236 237 ++number; 238 239 // compute maximum radius (needed in FinalizeTower to assess whether barrel or endcap tower) 240 if (particlePosition.Perp() > fTowerRmax) 241 fTowerRmax=particlePosition.Perp(); 237 242 238 243 pdgCode = TMath::Abs(particle->PID); … … 518 523 } 519 524 520 r = TMath::Sqrt(fTower->Position.X()*fTower->Position.X()+fTower->Position.Y()*fTower->Position.Y()); 525 // check whether barrel or endcap tower 526 if (fTower->Position.Perp() < fTowerRmax && TMath::Abs(eta) > 0.) 527 r = fTower->Position.Z()/TMath::SinH(eta); 528 else 529 r = fTower->Position.Pt(); 521 530 522 531 if(sumWeight > 0.0) -
modules/Calorimeter.h
rd1678fd r5eda6767 60 60 Double_t fTimingEnergyMin; 61 61 Bool_t fElectronsFromTrack; 62 Double_t fTowerRmax; 62 63 63 64 Int_t fTowerTrackHits, fTowerPhotonHits; -
modules/DualReadoutCalorimeter.cc
rd1678fd r5eda6767 247 247 fItParticleInputArray->Reset(); 248 248 number = -1; 249 fTowerRmax=0.; 249 250 while((particle = static_cast<Candidate*>(fItParticleInputArray->Next()))) 250 251 { 251 252 const TLorentzVector &particlePosition = particle->Position; 252 253 ++number; 254 255 // compute maximum radius (needed in FinalizeTower to assess whether barrel or endcap tower) 256 if (particlePosition.Perp() > fTowerRmax) 257 fTowerRmax=particlePosition.Perp(); 253 258 254 259 pdgCode = TMath::Abs(particle->PID); … … 558 563 } 559 564 560 r = TMath::Sqrt(fTower->Position.X()*fTower->Position.X()+fTower->Position.Y()*fTower->Position.Y()); 565 // check whether barrel or endcap tower 566 if (fTower->Position.Perp() < fTowerRmax && TMath::Abs(eta) > 0.) 567 r = fTower->Position.Z()/TMath::SinH(eta); 568 else 569 r = fTower->Position.Pt(); 561 570 562 571 if(sumWeight > 0.0) -
modules/DualReadoutCalorimeter.h
rd1678fd r5eda6767 60 60 Double_t fECalTrackEnergy, fHCalTrackEnergy; 61 61 Double_t fTrackEnergy; 62 Double_t fTowerRmax; 62 63 63 64 Double_t fTimingEnergyMin; -
modules/SimpleCalorimeter.cc
rd1678fd r5eda6767 208 208 fItParticleInputArray->Reset(); 209 209 number = -1; 210 fTowerRmax=0.; 210 211 while((particle = static_cast<Candidate *>(fItParticleInputArray->Next()))) 211 212 { 212 213 const TLorentzVector &particlePosition = particle->Position; 213 214 ++number; 215 216 // compute maximum radius (needed in FinalizeTower to assess whether barrel or endcap tower) 217 if (particlePosition.Perp() > fTowerRmax) 218 fTowerRmax=particlePosition.Perp(); 214 219 215 220 pdgCode = TMath::Abs(particle->PID); … … 444 449 445 450 pt = energy / TMath::CosH(eta); 446 r = TMath::Sqrt(fTower->Position.X()*fTower->Position.X()+fTower->Position.Y()*fTower->Position.Y()); 447 451 452 // check whether barrel or endcap tower 453 if (fTower->Position.Perp() < fTowerRmax && TMath::Abs(eta) > 0.) 454 r = fTower->Position.Z()/TMath::SinH(eta); 455 else 456 r = fTower->Position.Pt(); 457 448 458 fTower->Position.SetPtEtaPhiE(r, eta, phi, time); 449 459 fTower->Momentum.SetPtEtaPhiE(pt, eta, phi, energy); -
modules/SimpleCalorimeter.h
rd1678fd r5eda6767 62 62 Double_t fTrackTime; 63 63 64 Double_t fTowerRmax; 65 64 66 Double_t fTowerTimeWeight; 65 67 Double_t fTrackTimeWeight;
Note:
See TracChangeset
for help on using the changeset viewer.