Changeset 1230 for trunk/examples/extended/hadronic/Hadr01/src
- Timestamp:
- Jan 8, 2010, 3:02:48 PM (15 years ago)
- Location:
- trunk/examples/extended/hadronic/Hadr01/src
- Files:
-
- 15 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/hadronic/Hadr01/src/CheckVolumeSD.cc
r807 r1230 25 25 // 26 26 // $Id: CheckVolumeSD.cc,v 1.3 2006/06/29 17:24:02 gunter Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/DetectorConstruction.cc
r807 r1230 25 25 // 26 26 // $Id: DetectorConstruction.cc,v 1.6 2006/10/04 09:56:03 vnivanch Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/DetectorMessenger.cc
r807 r1230 25 25 // 26 26 // $Id: DetectorMessenger.cc,v 1.7 2007/05/16 11:43:30 vnivanch Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/EventAction.cc
r807 r1230 25 25 // 26 26 // $Id: EventAction.cc,v 1.4 2007/05/16 11:43:30 vnivanch Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/EventActionMessenger.cc
r807 r1230 25 25 // 26 26 // $Id: EventActionMessenger.cc,v 1.3 2006/06/29 17:24:11 gunter Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/Histo.cc
r807 r1230 51 51 verbose = ver; 52 52 histName = "histo"; 53 histType = " hbook";53 histType = "root"; 54 54 option = "--noErrors uncompress"; 55 55 nHisto = 0; 56 56 defaultAct = 1; 57 tupleName = "t uple.hbook";57 tupleName = "tree.root"; 58 58 tupleId = "100"; 59 59 tupleList = ""; -
trunk/examples/extended/hadronic/Hadr01/src/HistoManager.cc
r807 r1230 24 24 // ******************************************************************** 25 25 // 26 // $Id: HistoManager.cc,v 1.1 4 2007/05/24 13:52:31vnivanch Exp $27 // GEANT4 tag $Name: $26 // $Id: HistoManager.cc,v 1.16 2009/09/02 10:21:32 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 //--------------------------------------------------------------------------- … … 90 90 nSlices = 300; 91 91 nBinsE = 100; 92 nHisto = 2 2;92 nHisto = 25; 93 93 length = 300.*mm; 94 94 edepMax = 1.0*GeV; … … 133 133 histo->add1D("20","Log10 Energy (MeV) of pi+",nBinsE,-4.,6.,1.0); 134 134 histo->add1D("21","Log10 Energy (MeV) of pi-",nBinsE,-4.,6.,1.0); 135 histo->add1D("22","Energy deposition (GeV) in the target", 136 nBinsE,0.0,edepMax,GeV); 135 histo->add1D("22","Energy deposition in the target normalized to beam energy", 136 110,0.0,1.1,1.0); 137 histo->add1D("23","EM energy deposition in the target normalized to beam energy", 138 110,0.0,1.1,1.0); 139 histo->add1D("24","Pion energy deposition in the target normalized to beam energy", 140 110,0.0,1.1,1.0); 141 histo->add1D("25","Proton energy deposition in the target normalized to beam energy", 142 110,0.0,1.1,1.0); 137 143 } 138 144 … … 256 262 { 257 263 edepEvt = 0.0; 264 edepEM = 0.0; 265 edepPI = 0.0; 266 edepP = 0.0; 258 267 } 259 268 … … 264 273 edepSum += edepEvt; 265 274 edepSum2 += edepEvt*edepEvt; 266 histo->fill(21,edepEvt,1.0); 275 histo->fill(21,edepEvt/primaryKineticEnergy,1.0); 276 histo->fill(22,edepEM/primaryKineticEnergy,1.0); 277 histo->fill(23,edepPI/primaryKineticEnergy,1.0); 278 histo->fill(24,edepP/primaryKineticEnergy,1.0); 267 279 } 268 280 … … 372 384 edepEvt += edep; 373 385 histo->fill(0,z,edep); 386 const G4ParticleDefinition* pd = currentDef; 387 388 if(pd == G4Gamma::Gamma() || pd == G4Electron::Electron() 389 || pd == G4Positron::Positron()) { 390 edepEM += edep; 391 } else if ( pd == G4PionPlus::PionPlus() || pd == G4PionMinus::PionMinus()) { 392 edepPI += edep; 393 } else if ( pd == G4Proton::Proton() || pd == G4AntiProton::AntiProton()) { 394 edepP += edep; 395 } 374 396 375 397 if(1 < verbose) -
trunk/examples/extended/hadronic/Hadr01/src/HistoMessenger.cc
r807 r1230 26 26 // 27 27 // $Id: HistoMessenger.cc,v 1.5 2006/06/29 17:24:18 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/PhysicsList.cc
r807 r1230 25 25 // 26 26 // 27 // $Id: PhysicsList.cc,v 1. 24 2007/12/12 12:00:24 vnivanch Exp $28 // GEANT4 tag $Name: $27 // $Id: PhysicsList.cc,v 1.31 2009/08/28 09:21:34 vnivanch Exp $ 28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 ///////////////////////////////////////////////////////////////////////// … … 45 45 46 46 #include "G4DecayPhysics.hh" 47 #include "G4EmStandardPhysics_option3.hh" 47 48 #include "G4EmStandardPhysics_option2.hh" 48 49 #include "G4EmStandardPhysics_option1.hh" … … 52 53 #include "G4HadronQElasticPhysics.hh" 53 54 #include "G4HadronHElasticPhysics.hh" 55 #include "G4ChargeExchangePhysics.hh" 54 56 #include "G4NeutronTrackingCut.hh" 55 57 #include "G4QStoppingPhysics.hh" … … 60 62 #include "G4EmProcessOptions.hh" 61 63 62 #include "HadronPhysicsFTF C.hh"64 #include "HadronPhysicsFTFP_BERT_TRV.hh" 63 65 #include "HadronPhysicsFTFP.hh" 64 66 #include "HadronPhysicsFTFP_BERT.hh" 67 #include "HadronPhysicsFTF_BIC.hh" 65 68 #include "HadronPhysicsLHEP.hh" 66 #include "HadronPhysicsLHEP_BERT.hh"67 69 #include "HadronPhysicsLHEP_EMV.hh" 68 #include "HadronPhysicsLHEP_PRECO_HP.hh"69 70 #include "G4HadronInelasticQBBC.hh" 70 71 #include "HadronPhysicsQGSC.hh" 71 72 #include "HadronPhysicsQGSC_BERT.hh" 72 #include "HadronPhysicsQGSC_EFLOW.hh" 73 #include "HadronPhysicsQGSC_CHIPS.hh" 74 #include "HadronPhysicsQGSC_QGSC.hh" 73 75 #include "HadronPhysicsQGSP.hh" 74 76 #include "HadronPhysicsQGSP_BERT.hh" … … 76 78 #include "HadronPhysicsQGSP_BIC.hh" 77 79 #include "HadronPhysicsQGSP_BIC_HP.hh" 80 #include "HadronPhysicsQGSP_FTFP_BERT.hh" 81 #include "HadronPhysicsQGS_BIC.hh" 78 82 79 83 #include "G4HadronInelasticQLHEP.hh" 80 84 #include "G4IonPhysics.hh" 81 #include "G4HadronProcessStore.hh"82 85 83 86 #include "G4LossTableManager.hh" … … 100 103 cutForPositron = defaultCutValue; 101 104 dump = false; 105 verboseLevel = 1; 102 106 103 107 pMessenger = new PhysicsListMessenger(this); … … 139 143 hadronPhys[i]->ConstructProcess(); 140 144 } 141 142 if(dump) G4HadronProcessStore::Instance()->Dump(1);143 145 } 144 146 … … 155 157 emPhysicsList = new G4EmStandardPhysics_option2(); 156 158 159 } else if (name == "emstandard_opt3") { 160 161 delete emPhysicsList; 162 emPhysicsList = new G4EmStandardPhysics_option3(); 163 157 164 } else if (name == "emstandard_opt1") { 158 165 … … 160 167 emPhysicsList = new G4EmStandardPhysics_option1(); 161 168 162 } else if (name == "FTFC") { 163 164 SetBuilderList1(); 165 hadronPhys.push_back( new HadronPhysicsFTFC("hadron",true)); 166 dump = true; 169 } else if (name == "FTFP_BERT_EMV") { 170 171 AddPhysicsList("emstandard_opt1"); 172 AddPhysicsList("FTFP_BERT"); 173 174 } else if (name == "FTFP_BERT_EMX") { 175 176 AddPhysicsList("emstandard_opt2"); 177 AddPhysicsList("FTFP_BERT"); 167 178 168 179 } else if (name == "FTFP") { … … 178 189 dump = true; 179 190 180 } else if (name == "FTFP_EMV") { 181 182 AddPhysicsList("emstandard_opt1"); 183 AddPhysicsList("FTFP"); 191 } else if (name == "FTFP_BERT_TRV") { 192 193 SetBuilderList1(); 194 hadronPhys.push_back( new HadronPhysicsFTFP_BERT_TRV("hadron",true)); 195 dump = true; 196 197 } else if (name == "FTF_BIC") { 198 199 SetBuilderList0(); 200 hadronPhys.push_back( new HadronPhysicsFTF_BIC("hadron",true)); 201 dump = true; 184 202 185 203 } else if (name == "LHEP") { … … 187 205 SetBuilderList2(); 188 206 hadronPhys.push_back( new HadronPhysicsLHEP("hadron")); 189 dump = true;190 191 } else if (name == "LHEP_BERT") {192 193 SetBuilderList3();194 hadronPhys.push_back( new HadronPhysicsLHEP_BERT("hadron"));195 207 dump = true; 196 208 … … 202 214 dump = true; 203 215 204 } else if (name == "LHEP_PRECO_HP") {205 206 SetBuilderList3(true);207 hadronPhys.push_back( new HadronPhysicsLHEP_PRECO_HP("hadron"));208 dump = true;209 210 216 } else if (name == "QBBC") { 211 217 212 218 SetBuilderList0(); 213 219 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel, 214 false,false,false,false,true)); 220 false,true,false,false,false)); 221 222 } else if (name == "QBBCG") { 223 224 SetBuilderList0(); 225 hadronPhys.push_back( new G4ChargeExchangePhysics(verboseLevel,false)); 226 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel, 227 false,true,false,false,true)); 228 } else if (name == "QBBCF") { 229 230 SetBuilderList0(); 231 hadronPhys.push_back( new G4ChargeExchangePhysics(verboseLevel,false)); 232 hadronPhys.push_back( new G4HadronInelasticQBBC("QBBC",verboseLevel, 233 false,true,false,false,false)); 215 234 216 235 } else if (name == "QBBC_DEL") { … … 243 262 dump = true; 244 263 245 } else if (name == "QGSC_ EFLOW") {264 } else if (name == "QGSC_CHIPS") { 246 265 247 266 SetBuilderList4(); 248 hadronPhys.push_back( new HadronPhysicsQGSC_EFLOW("hadron",true)); 267 hadronPhys.push_back( new HadronPhysicsQGSC_CHIPS("hadron",true)); 268 dump = true; 269 270 } else if (name == "QGSC_QGSC") { 271 272 SetBuilderList4(); 273 hadronPhys.push_back( new HadronPhysicsQGSC_QGSC("hadron",true)); 249 274 dump = true; 250 275 … … 266 291 dump = true; 267 292 293 } else if (name == "QGSP_FTFP_BERT") { 294 295 SetBuilderList1(); 296 hadronPhys.push_back( new HadronPhysicsQGSP_FTFP_BERT("hadron",true)); 297 dump = true; 298 268 299 } else if (name == "QGSP_BERT_EMV") { 269 300 … … 271 302 AddPhysicsList("QGSP_BERT"); 272 303 304 } else if (name == "QGSP_BERT_EMX") { 305 306 AddPhysicsList("emstandard_opt2"); 307 AddPhysicsList("QGSP_BERT"); 308 273 309 } else if (name == "QGSP_BERT_HP") { 274 310 … … 282 318 dump = true; 283 319 320 } else if (name == "QGS_BIC") { 321 322 SetBuilderList0(); 323 hadronPhys.push_back( new HadronPhysicsQGS_BIC("hadron",true)); 324 dump = true; 325 284 326 } else if (name == "QGSP_BIC_HP") { 285 327 … … 288 330 dump = true; 289 331 290 } else if (name == "QGSP_DIF") { 291 292 SetBuilderList1(); 293 HadronPhysicsQGSP * hp = new HadronPhysicsQGSP("hadron"); 294 hp->SetQuasiElastic(true); 295 hp->SetProjectileDiffraction(true); 296 hadronPhys.push_back(hp); 297 dump = true; 298 299 } else if (name == "QGSP_EMV") { 300 301 AddPhysicsList("emstandard_opt1"); 302 AddPhysicsList("QGSP"); 303 dump = true; 304 305 } else if (name == "QGSP_EMX") { 306 307 AddPhysicsList("emstandard_opt2"); 308 AddPhysicsList("QGSP"); 309 dump = true; 310 311 } else if (name == "QGSP_NQE") { 312 313 SetBuilderList1(); 314 hadronPhys.push_back( new HadronPhysicsQGSP("hadron",false)); 315 dump = true; 316 317 } else if (name == "QGSP_QEL") { 318 319 SetBuilderList4(); 320 hadronPhys.push_back( new HadronPhysicsQGSP("hadron",true)); 321 dump = true; 332 } else if (name == "QGSP_BIC_EMY") { 333 334 AddPhysicsList("emstandard_opt3"); 335 AddPhysicsList("QGSP_BIC"); 322 336 323 337 } else { … … 338 352 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel)); 339 353 hadronPhys.push_back( new G4IonBinaryCascadePhysics("ionBIC")); 340 hadronPhys.push_back( new G4NeutronTrackingCut(" Neutron tracking cut"));354 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel)); 341 355 } 342 356 … … 350 364 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel)); 351 365 hadronPhys.push_back( new G4IonPhysics("ion")); 352 hadronPhys.push_back( new G4NeutronTrackingCut(" Neutron tracking cut"));366 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel)); 353 367 } 354 368 … … 382 396 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel)); 383 397 hadronPhys.push_back( new G4IonPhysics("ion")); 384 hadronPhys.push_back( new G4NeutronTrackingCut(" Neutron tracking cut"));398 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel)); 385 399 } 386 400 … … 393 407 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel)); 394 408 hadronPhys.push_back( new G4IonBinaryCascadePhysics("ionBIC")); 395 hadronPhys.push_back( new G4NeutronTrackingCut(" Neutron tracking cut"));409 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel)); 396 410 } 397 411 … … 404 418 hadronPhys.push_back( new G4QStoppingPhysics("stopping",verboseLevel)); 405 419 hadronPhys.push_back( new G4IonBinaryCascadePhysics("ionBIC")); 406 hadronPhys.push_back( new G4NeutronTrackingCut(" Neutron tracking cut"));420 hadronPhys.push_back( new G4NeutronTrackingCut("nTackingCut",verboseLevel)); 407 421 } 408 422 … … 452 466 void PhysicsList::List() 453 467 { 454 G4cout << "### PhysicsLists available: FTF C FTFP FTFP_BERT FTFP_EMV LHEP LHEP_BERT LHEP_EMV"468 G4cout << "### PhysicsLists available: FTFP FTFP_BERT FTFP_BERT_EMV FTFP_BERT_EMX FTFP_BERT_TRV FTF_BIC LHEP LHEP_EMV" 455 469 << G4endl; 456 G4cout << " LHEP_PRECO_HPQBBC QBBC_DEL QBBC_HEL QBBC_HP QGSC "470 G4cout << " QBBC QBBC_DEL QBBC_HEL QBBC_HP QGSC " 457 471 << G4endl; 458 G4cout << " QGSC_BERT QGSC_ EFLOWQGSC_EMV QGSP QGSP_BERT QGSP_BER_EMV "472 G4cout << " QGSC_BERT QGSC_CHIPS QGSC_QGSC QGSC_EMV QGSP QGSP_BERT QGSP_BER_EMV " 459 473 << G4endl; 460 G4cout << " QGSP_BERT_ HP QGSP_BIC QGSP_BIC_HP QGSP_DIF"474 G4cout << " QGSP_BERT_EMX QGSP_BERT_HP QGSP_BIC QGSP_BIC_HP QGSP_FTFP_BERT " 461 475 << G4endl; 462 G4cout << " QGSP_EMV QGSP_EMX QGSP_NQE QGSP_QEL "463 << G4endl;464 476 } 465 477 -
trunk/examples/extended/hadronic/Hadr01/src/PhysicsListMessenger.cc
r807 r1230 26 26 // 27 27 // $Id: PhysicsListMessenger.cc,v 1.4 2006/08/11 14:38:11 vnivanch Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // -
trunk/examples/extended/hadronic/Hadr01/src/PrimaryGeneratorAction.cc
r807 r1230 25 25 // 26 26 // $Id: PrimaryGeneratorAction.cc,v 1.4 2006/11/15 14:58:10 vnivanch Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/RunAction.cc
r807 r1230 40 40 #include "G4UImanager.hh" 41 41 #include "G4VVisManager.hh" 42 #include "G4HadronProcessStore.hh"43 42 #include "G4NistManager.hh" 44 43 #include "G4Element.hh" -
trunk/examples/extended/hadronic/Hadr01/src/StackingAction.cc
r807 r1230 25 25 // 26 26 // $Id: StackingAction.cc,v 1.5 2006/10/04 09:56:03 vnivanch Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/StackingMessenger.cc
r807 r1230 25 25 // 26 26 // $Id: StackingMessenger.cc,v 1.3 2006/06/29 17:24:32 gunter Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 ///////////////////////////////////////////////////////////////////////// -
trunk/examples/extended/hadronic/Hadr01/src/TargetSD.cc
r807 r1230 25 25 // 26 26 // $Id: TargetSD.cc,v 1.4 2006/06/29 17:24:34 gunter Exp $ 27 // GEANT4 tag $Name: $27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 // 29 29 /////////////////////////////////////////////////////////////////////////
Note: See TracChangeset
for help on using the changeset viewer.