Changes in modules/Calorimeter.cc [4b9a2dc:82575a3] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Calorimeter.cc
r4b9a2dc r82575a3 135 135 } 136 136 */ 137 138 // read min E value for towers to be saved139 fEcalEnergyMin = GetDouble("EcalTowerMinEnergy", 0.0);140 fHcalEnergyMin = GetDouble("HcalTowerMinEnergy", 0.0);141 142 fEcalSigmaMin = GetDouble("EcalTowerMinSignificance", 0.0);143 fHcalSigmaMin = GetDouble("HcalTowerMinSignificance", 0.0);144 145 146 137 // read resolution formulas 147 138 fECalResolutionFormula->Compile(GetString("ECalResolutionFormula", "0")); … … 423 414 424 415 if(!fTower) return; 425 // cout<<"----------------------"<<endl;426 // cout<<"Finalize Tower"<<endl;427 // cout<<""<<endl;428 429 416 430 417 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy); … … 443 430 hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma); 444 431 hcalTime = (fTowerHCalWeightTime < 1.0E-09 ) ? 0 : fTowerHCalTime/fTowerHCalWeightTime; 445 446 447 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, ecalEnergy);448 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy);449 450 ecalEnergy = (ecalEnergy < fEcalEnergyMin || ecalEnergy < fEcalSigmaMin*ecalSigma) ? 0 : ecalEnergy;451 hcalEnergy = (hcalEnergy < fHcalEnergyMin || hcalEnergy < fHcalSigmaMin*hcalSigma) ? 0 : hcalEnergy;452 432 453 433 energy = ecalEnergy + hcalEnergy; … … 473 453 fTower->Edges[3] = fTowerEdges[3]; 474 454 475 if( energy > 0.0 ) 455 456 // fill calorimeter towers 457 if(energy > 0.0) 476 458 { 477 459 if(fTowerPhotonHits > 0 && fTowerTrackHits == 0) … … 493 475 494 476 ecalEnergy -= fTrackECalEnergy; 495 if(ecalEnergy < fEcalEnergyMin || ecalEnergy < fEcalSigmaMin*fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, ecalEnergy)) ecalEnergy = 0.0;477 if(ecalEnergy < 0.0) ecalEnergy = 0.0; 496 478 497 479 hcalEnergy -= fTrackHCalEnergy; 498 if(hcalEnergy < fHcalEnergyMin || hcalEnergy < fHcalSigmaMin*fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy)) hcalEnergy = 0.0;480 if(hcalEnergy < 0.0) hcalEnergy = 0.0; 499 481 500 482 energy = ecalEnergy + hcalEnergy; 501 483 484 485 // save ECAL and/or HCAL energy excess as an energy flow tower 502 486 if(ecalEnergy > 0.0) 503 487 { … … 513 497 fEFlowPhotonOutputArray->Add(tower); 514 498 } 499 515 500 if(hcalEnergy > 0.0) 516 501 {
Note:
See TracChangeset
for help on using the changeset viewer.