- Timestamp:
- Oct 18, 2021, 12:05:26 PM (3 years ago)
- Branches:
- master
- Children:
- 0b8551f, ebf40fd, fc2ef70
- Parents:
- 78e1fb3
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/DualReadoutCalorimeter.cc
r78e1fb3 r9a7ea36 425 425 } 426 426 427 428 427 // in Dual Readout we do not care if tracks are ECAL of HCAL 429 428 if(fECalTrackFractions[number] > 1.0E-9 || fHCalTrackFractions[number] > 1.0E-9) … … 444 443 fTrackSigma += (track->TrackResolution)*energyGuess*(track->TrackResolution)*energyGuess; 445 444 fTowerTrackArray->Add(track); 446 447 445 } 448 446 else … … 472 470 fTowerTime += (ecalEnergy + hcalEnergy) * position.T(); //sigma_t ~ 1/sqrt(E) 473 471 fTowerTimeWeight += ecalEnergy + hcalEnergy; 474 //fTowerTime += (hcalEnergy) * position.T(); //sigma_t ~ 1/sqrt(E)475 //fTowerTimeWeight += hcalEnergy;476 //fTowerTime += position.T(); //sigma_t ~ 1/sqrt(E)477 //fTowerTimeWeight += 1;478 479 //cout<<" tower particle PID, pt, eta, phi, l, tof: "<<particle->PID<<", "<<momentum.E()<<", "<<momentum.Eta()<<", "<<momentum.Phi()<<", "<<position.Vect().Mag()<<", "<<position.T()/2.99792458E2<<endl;480 //cout<<" tower particle time, weight: "<<fTowerTime/2.99792458E2<<", "<<fTowerTimeWeight<<endl;481 472 482 473 fTower->AddCandidate(particle); … … 515 506 energy = fECalTowerEnergy; 516 507 sigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, energy); 517 //cout<<"em energy"<<energy<<", sigma: "<<sigma<<endl;518 508 } 519 509 … … 523 513 energy = fECalTowerEnergy + fHCalTowerEnergy; 524 514 sigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, energy); 525 //cout<<"had energy: "<<energy<<", sigma: "<<sigma<<endl;526 515 } 527 516
Note:
See TracChangeset
for help on using the changeset viewer.