- Timestamp:
- Dec 19, 2014, 2:11:47 PM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- ca29ad7
- Parents:
- c04eb8a
- Location:
- modules
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Calorimeter.cc
rc04eb8a r8624f58 145 145 fFractionMap[param[i*2].GetInt()] = make_pair(ecalFraction, hcalFraction); 146 146 } 147 147 148 /* 148 149 TFractionMap::iterator itFractionMap; … … 154 155 155 156 // read min E value for towers to be saved 156 fEcalEnergyMin = GetDouble("EcalTowerMinEnergy", 0.0); 157 fHcalEnergyMin = GetDouble("HcalTowerMinEnergy", 0.0); 158 159 fEcalSigmaMin = GetDouble("EcalTowerMinSignificance", 0.0); 160 fHcalSigmaMin = GetDouble("HcalTowerMinSignificance", 0.0); 161 157 fECalEnergyMin = GetDouble("ECalMinEnergy", 0.0); 158 fHCalEnergyMin = GetDouble("HCalMinEnergy", 0.0); 159 160 fECalSigmaMin = GetDouble("ECalMinSignificance", 0.0); 161 fHCalSigmaMin = GetDouble("HCalMinSignificance", 0.0); 162 163 // switch on or off the dithering of the center of calorimeter towers 164 fDitherTowerCenter = GetBool("DitherTowerCenter", true); 162 165 163 166 // read resolution formulas … … 180 183 fEFlowNeutralHadronOutputArray = ExportArray(GetString("EFlowNeutralHadronOutputArray", "eflowNeutralHadrons")); 181 184 182 fDitherTowerCenter = GetBool("DitherTowerCenter", true);183 185 } 184 186 … … 440 442 441 443 if(!fTower) return; 442 // cout<<"----------------------"<<endl;443 // cout<<"Finalize Tower"<<endl;444 // cout<<""<<endl;445 446 444 447 445 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy); 448 449 // ecalEnergy = gRandom->Gaus(fTowerECalEnergy, ecalSigma);450 // if(ecalEnergy < 0.0) ecalEnergy = 0.0;451 446 452 447 ecalEnergy = LogNormal(fTowerECalEnergy, ecalSigma); … … 455 450 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy); 456 451 457 // hcalEnergy = gRandom->Gaus(fTowerHCalEnergy, hcalSigma);458 // if(hcalEnergy < 0.0) hcalEnergy = 0.0;459 460 452 hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma); 461 453 hcalTime = (fTowerHCalWeightTime < 1.0E-09 ) ? 0 : fTowerHCalTime/fTowerHCalWeightTime; … … 465 457 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy); 466 458 467 ecalEnergy = (ecalEnergy < fE calEnergyMin || ecalEnergy < fEcalSigmaMin*ecalSigma) ? 0 : ecalEnergy;468 hcalEnergy = (hcalEnergy < fH calEnergyMin || hcalEnergy < fHcalSigmaMin*hcalSigma) ? 0 : hcalEnergy;459 ecalEnergy = (ecalEnergy < fECalEnergyMin || ecalEnergy < fECalSigmaMin*ecalSigma) ? 0 : ecalEnergy; 460 hcalEnergy = (hcalEnergy < fHCalEnergyMin || hcalEnergy < fHCalSigmaMin*hcalSigma) ? 0 : hcalEnergy; 469 461 470 462 energy = ecalEnergy + hcalEnergy; … … 484 476 pt = energy / TMath::CosH(eta); 485 477 486 // fTower->Position.SetXYZT(-time, 0.0, 0.0, time);487 478 fTower->Position.SetPtEtaPhiE(1.0, eta, phi, time); 488 479 fTower->Momentum.SetPtEtaPhiE(pt, eta, phi, energy); … … 515 506 516 507 ecalEnergy -= fTrackECalEnergy; 517 if(ecalEnergy < fE calEnergyMin || ecalEnergy < fEcalSigmaMin*fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, ecalEnergy)) ecalEnergy = 0.0;508 if(ecalEnergy < fECalEnergyMin || ecalEnergy < fECalSigmaMin*fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, ecalEnergy)) ecalEnergy = 0.0; 518 509 519 510 hcalEnergy -= fTrackHCalEnergy; 520 if(hcalEnergy < fH calEnergyMin || hcalEnergy < fHcalSigmaMin*fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy)) hcalEnergy = 0.0;511 if(hcalEnergy < fHCalEnergyMin || hcalEnergy < fHCalSigmaMin*fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy)) hcalEnergy = 0.0; 521 512 522 513 energy = ecalEnergy + hcalEnergy; … … 548 539 fEFlowNeutralHadronOutputArray->Add(tower); 549 540 } 550 551 552 553 554 541 } 555 542 -
modules/Calorimeter.h
rc04eb8a r8624f58 72 72 Int_t fTowerTrackHits, fTowerPhotonHits; 73 73 74 Double_t fECalEnergyMin; 75 Double_t fHCalEnergyMin; 76 77 Double_t fECalSigmaMin; 78 Double_t fHCalSigmaMin; 79 74 80 Bool_t fDitherTowerCenter; 75 76 Double_t fEcalEnergyMin;77 Double_t fHcalEnergyMin;78 79 Double_t fEcalSigmaMin;80 Double_t fHcalSigmaMin;81 81 82 82 TFractionMap fFractionMap; //!
Note:
See TracChangeset
for help on using the changeset viewer.