Changeset 38bf1ae in git
- Timestamp:
- Dec 21, 2014, 12:12:45 PM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- f10999a
- Parents:
- 6cdc544
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/delphes_card_ATLAS.tcl
r6cdc544 r38bf1ae 207 207 set PhotonOutputArray photons 208 208 209 set ECal MinEnergy0.50210 set HCal MinEnergy1.00211 212 set ECal MinSignificance1.0213 set HCal MinSignificance1.0209 set ECalEnergyMin 0.50 210 set HCalEnergyMin 1.00 211 212 set ECalEnergySignificanceMin 1.0 213 set HCalEnergySignificanceMin 1.0 214 214 215 215 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_ATLAS_PileUp.tcl
r6cdc544 r38bf1ae 249 249 set PhotonOutputArray photons 250 250 251 set ECal MinEnergy0.50252 set HCal MinEnergy1.00253 254 set ECal MinSignificance1.0255 set HCal MinSignificance1.0251 set ECalEnergyMin 0.50 252 set HCalEnergyMin 1.00 253 254 set ECalEnergySignificanceMin 1.0 255 set HCalEnergySignificanceMin 1.0 256 256 257 257 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_CMS.tcl
r6cdc544 r38bf1ae 212 212 set PhotonOutputArray photons 213 213 214 set ECal MinEnergy0.5215 set HCal MinEnergy1.0216 217 set ECal MinSignificance1.0218 set HCal MinSignificance1.0214 set ECalEnergyMin 0.5 215 set HCalEnergyMin 1.0 216 217 set ECalEnergySignificanceMin 1.0 218 set HCalEnergySignificanceMin 1.0 219 219 220 220 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_CMS_NoFastJet.tcl
r6cdc544 r38bf1ae 186 186 set PhotonOutputArray photons 187 187 188 set ECal MinEnergy0.5189 set HCal MinEnergy1.0190 191 set ECal MinSignificance1.0192 set HCal MinSignificance1.0188 set ECalEnergyMin 0.5 189 set HCalEnergyMin 1.0 190 191 set ECalEnergySignificanceMin 1.0 192 set HCalEnergySignificanceMin 1.0 193 193 194 194 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_CMS_PileUp.tcl
r6cdc544 r38bf1ae 255 255 set PhotonOutputArray photons 256 256 257 set ECal MinEnergy0.5258 set HCal MinEnergy1.0259 260 set ECal MinSignificance1.0261 set HCal MinSignificance1.0257 set ECalEnergyMin 0.5 258 set HCalEnergyMin 1.0 259 260 set ECalEnergySignificanceMin 1.0 261 set HCalEnergySignificanceMin 1.0 262 262 263 263 set EFlowTrackOutputArray eflowTracks -
modules/Calorimeter.cc
r6cdc544 r38bf1ae 151 151 152 152 // read min E value for towers to be saved 153 fECalEnergyMin = GetDouble("ECal MinEnergy", 0.0);154 fHCalEnergyMin = GetDouble("HCal MinEnergy", 0.0);155 156 fECal SigmaMin = GetDouble("ECalMinSignificance", 0.0);157 fHCal SigmaMin = GetDouble("HCalMinSignificance", 0.0);153 fECalEnergyMin = GetDouble("ECalEnergyMin", 0.0); 154 fHCalEnergyMin = GetDouble("HCalEnergyMin", 0.0); 155 156 fECalEnergySignificanceMin = GetDouble("ECalEnergySignificanceMin", 0.0); 157 fHCalEnergySignificanceMin = GetDouble("HCalEnergySignificanceMin", 0.0); 158 158 159 159 // switch on or off the dithering of the center of calorimeter towers … … 178 178 fEFlowPhotonOutputArray = ExportArray(GetString("EFlowPhotonOutputArray", "eflowPhotons")); 179 179 fEFlowNeutralHadronOutputArray = ExportArray(GetString("EFlowNeutralHadronOutputArray", "eflowNeutralHadrons")); 180 181 180 } 182 181 … … 440 439 441 440 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerECalEnergy); 441 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy); 442 442 443 443 ecalEnergy = LogNormal(fTowerECalEnergy, ecalSigma); 444 ecalTime = (fTowerECalTimeWeight < 1.0E-09 ) ? 0 : fTowerECalTime/fTowerECalTimeWeight;445 446 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, fTowerHCalEnergy);447 448 444 hcalEnergy = LogNormal(fTowerHCalEnergy, hcalSigma); 449 hcalTime = (fTowerHCalTimeWeight < 1.0E-09 ) ? 0 : fTowerHCalTime/fTowerHCalTimeWeight; 450 445 446 ecalTime = (fTowerECalTimeWeight < 1.0E-09 ) ? 0.0 : fTowerECalTime/fTowerECalTimeWeight; 447 hcalTime = (fTowerHCalTimeWeight < 1.0E-09 ) ? 0.0 : fTowerHCalTime/fTowerHCalTimeWeight; 451 448 452 449 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, ecalEnergy); 453 450 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy); 454 451 455 ecalEnergy = (ecalEnergy < fECalEnergyMin || ecalEnergy < fECalSigmaMin*ecalSigma) ? 0 : ecalEnergy;456 hcalEnergy = (hcalEnergy < fHCalEnergyMin || hcalEnergy < fHCalSigmaMin*hcalSigma) ? 0 : hcalEnergy;452 if(ecalEnergy < fECalEnergyMin || ecalEnergy < fECalEnergySignificanceMin*ecalSigma) ecalEnergy = 0.0; 453 if(hcalEnergy < fHCalEnergyMin || hcalEnergy < fHCalEnergySignificanceMin*hcalSigma) hcalEnergy = 0.0; 457 454 458 455 energy = ecalEnergy + hcalEnergy; … … 461 458 if(fDitherTowerCenter) 462 459 { 463 464 460 eta = gRandom->Uniform(fTowerEdges[0], fTowerEdges[1]); 461 phi = gRandom->Uniform(fTowerEdges[2], fTowerEdges[3]); 465 462 } 466 463 else 467 464 { 468 469 465 eta = fTowerEta; 466 phi = fTowerPhi; 470 467 } 471 468 … … 482 479 fTower->Edges[3] = fTowerEdges[3]; 483 480 484 if( energy > 0.0)481 if(energy > 0.0) 485 482 { 486 483 if(fTowerPhotonHits > 0 && fTowerTrackHits == 0) … … 502 499 503 500 ecalEnergy -= fTrackECalEnergy; 504 if(ecalEnergy < fECalEnergyMin || ecalEnergy < fECalSigmaMin*fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, ecalEnergy)) ecalEnergy = 0.0;505 506 501 hcalEnergy -= fTrackHCalEnergy; 507 if(hcalEnergy < fHCalEnergyMin || hcalEnergy < fHCalSigmaMin*fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy)) hcalEnergy = 0.0; 502 503 ecalSigma = fECalResolutionFormula->Eval(0.0, fTowerEta, 0.0, ecalEnergy); 504 hcalSigma = fHCalResolutionFormula->Eval(0.0, fTowerEta, 0.0, hcalEnergy); 505 506 if(ecalEnergy < fECalEnergyMin || ecalEnergy < fECalEnergySignificanceMin*ecalSigma) ecalEnergy = 0.0; 507 if(hcalEnergy < fHCalEnergyMin || hcalEnergy < fHCalEnergySignificanceMin*hcalSigma) hcalEnergy = 0.0; 508 508 509 509 energy = ecalEnergy + hcalEnergy; … … 518 518 tower->Momentum.SetPtEtaPhiE(pt, eta, phi, ecalEnergy); 519 519 tower->Eem = ecalEnergy; 520 tower->Ehad = 0 ;520 tower->Ehad = 0.0; 521 521 522 522 fEFlowPhotonOutputArray->Add(tower); … … 530 530 531 531 tower->Momentum.SetPtEtaPhiE(pt, eta, phi, hcalEnergy); 532 tower->Eem = 0 ;532 tower->Eem = 0.0; 533 533 tower->Ehad = hcalEnergy; 534 534 … … 548 548 a = TMath::Log(mean) - 0.5*b*b; 549 549 550 return TMath::Exp(a + b*gRandom->Gaus(0 , 1));550 return TMath::Exp(a + b*gRandom->Gaus(0.0, 1.0)); 551 551 } 552 552 else -
modules/Calorimeter.h
r6cdc544 r38bf1ae 71 71 Double_t fHCalEnergyMin; 72 72 73 Double_t fECal SigmaMin;74 Double_t fHCal SigmaMin;73 Double_t fECalEnergySignificanceMin; 74 Double_t fHCalEnergySignificanceMin; 75 75 76 76 Bool_t fDitherTowerCenter;
Note:
See TracChangeset
for help on using the changeset viewer.