Changeset 8624f58 in git
- Timestamp:
- Dec 19, 2014, 2:11:47 PM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- ca29ad7
- Parents:
- c04eb8a
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/delphes_card_ATLAS.tcl
rc04eb8a r8624f58 207 207 set PhotonOutputArray photons 208 208 209 set E calTowerMinEnergy 0.50210 set H calTowerMinEnergy 1.00211 212 set E calTowerMinSignificance 1.0213 set H calTowerMinSignificance 1.0209 set ECalMinEnergy 0.50 210 set HCalMinEnergy 1.00 211 212 set ECalMinSignificance 1.0 213 set HCalMinSignificance 1.0 214 214 215 215 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_ATLAS_PileUp.tcl
rc04eb8a r8624f58 249 249 set PhotonOutputArray photons 250 250 251 set E calTowerMinEnergy 0.50252 set H calTowerMinEnergy 1.00253 254 set E calTowerMinSignificance 1.0255 set H calTowerMinSignificance 1.0251 set ECalMinEnergy 0.50 252 set HCalMinEnergy 1.00 253 254 set ECalMinSignificance 1.0 255 set HCalMinSignificance 1.0 256 256 257 257 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_CMS.tcl
rc04eb8a r8624f58 212 212 set PhotonOutputArray photons 213 213 214 set E calTowerMinEnergy 0.5215 set H calTowerMinEnergy 1.0216 217 set E calTowerMinSignificance 1.0218 set H calTowerMinSignificance 1.0214 set ECalMinEnergy 0.5 215 set HCalMinEnergy 1.0 216 217 set ECalMinSignificance 1.0 218 set HCalMinSignificance 1.0 219 219 220 220 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_CMS_NoFastJet.tcl
rc04eb8a r8624f58 186 186 set PhotonOutputArray photons 187 187 188 set E calTowerMinEnergy 0.5189 set H calTowerMinEnergy 1.0190 191 set E calTowerMinSignificance 1.0192 set H calTowerMinSignificance 1.0188 set ECalMinEnergy 0.5 189 set HCalMinEnergy 1.0 190 191 set ECalMinSignificance 1.0 192 set HCalMinSignificance 1.0 193 193 194 194 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_CMS_PileUp.tcl
rc04eb8a r8624f58 255 255 set PhotonOutputArray photons 256 256 257 set E calTowerMinEnergy 0.5258 set H calTowerMinEnergy 1.0259 260 set E calTowerMinSignificance 1.0261 set H calTowerMinSignificance 1.0257 set ECalMinEnergy 0.5 258 set HCalMinEnergy 1.0 259 260 set ECalMinSignificance 1.0 261 set HCalMinSignificance 1.0 262 262 263 263 set EFlowTrackOutputArray eflowTracks -
cards/delphes_card_FCC_basic.tcl
rc04eb8a r8624f58 18 18 ImpactParameterSmearing 19 19 20 E cal21 H cal20 ECal 21 HCal 22 22 23 23 TowerMerger … … 219 219 ############# 220 220 221 module SimpleCalorimeter E cal {221 module SimpleCalorimeter ECal { 222 222 set ParticleInputArray ParticlePropagator/stableParticles 223 223 set TrackInputArray ImpactParameterSmearing/tracks … … 281 281 ############# 282 282 283 module SimpleCalorimeter H cal {283 module SimpleCalorimeter HCal { 284 284 set ParticleInputArray ParticlePropagator/stableParticles 285 285 set TrackInputArray ImpactParameterSmearing/tracks … … 346 346 module Merger TowerMerger { 347 347 # add InputArray InputArray 348 add InputArray E cal/ecalTowers349 add InputArray H cal/hcalTowers348 add InputArray ECal/ecalTowers 349 add InputArray HCal/hcalTowers 350 350 set OutputArray towers 351 351 } … … 358 358 # add InputArray InputArray 359 359 add InputArray ImpactParameterSmearing/tracks 360 add InputArray E cal/eflowPhotons361 add InputArray H cal/eflowNeutralHadrons360 add InputArray ECal/eflowPhotons 361 add InputArray HCal/eflowNeutralHadrons 362 362 set OutputArray eflow 363 363 } … … 512 512 513 513 add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track 514 add Branch H cal/eflowNeutralHadrons NeutralHadron Tower515 add Branch E cal/eflowPhotons Photon Photon514 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower 515 add Branch ECal/eflowPhotons Photon Photon 516 516 517 517 add Branch ElectronEnergySmearing/electrons Electron Electron -
cards/delphes_card_LHCb.tcl
rc04eb8a r8624f58 20 20 IdentificationMap 21 21 22 E cal23 H cal22 ECal 23 HCal 24 24 25 25 TreeWriter … … 217 217 ############# 218 218 219 module SimpleCalorimeter E cal {219 module SimpleCalorimeter ECal { 220 220 set ParticleInputArray ParticlePropagator/stableParticles 221 221 set TrackInputArray IdentificationMap/tracks … … 292 292 ############# 293 293 294 module SimpleCalorimeter H cal {294 module SimpleCalorimeter HCal { 295 295 set ParticleInputArray ParticlePropagator/stableParticles 296 296 set TrackInputArray IdentificationMap/tracks … … 369 369 370 370 add Branch IdentificationMap/tracks Track Track 371 add Branch H cal/eflowNeutralHadrons NeutralHadron Tower372 add Branch E cal/eflowPhotons Photon Photon373 374 } 375 371 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower 372 add Branch ECal/eflowPhotons Photon Photon 373 374 } 375 -
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.