1 / 17

Speed up ALICE simulation with TFluka

Speed up ALICE simulation with TFluka. Ernesto Lopez Torres. Problem. HIJING Central event takes ~ 60 h with TFluka ~ 4 h with Geant3. Why TFluka takes so long time?. First ideas: Energy transport cuts on region basics is not allow for hadrons, so we have 1MeV cut for all volumes

zarita
Download Presentation

Speed up ALICE simulation with TFluka

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Speed up ALICE simulation with TFluka Ernesto Lopez Torres

  2. Problem • HIJING Central event takes ~ 60 h with TFluka ~ 4 h with Geant3

  3. Why TFluka takes so long time? First ideas: • Energy transport cuts on region basics is not allow for hadrons, so we have 1MeV cut for all volumes • There is not cut on number of track steps (for GEANT we have set max. 30 000 steps) • More detail physics?

  4. Fix on the CVS Where is Fluka spending time? Some primary particle create a lot secondaries (105 - 106) • Small bug on HMPID: by error someone comment the initialization of the array with the quantum efficiency properties. QEff = 1031 instead of = 1 T0 HMPID

  5. Where is Fluka spending time? TFluka HIJING Pb-Pb 5.5TeV central event 56140 particles(Laptop Centrino 2GHz 1Gb RAM) Total transport time: 236389s (65.6 h) # Primaries with (transp. time > 30s): 657 (1.1%) with a total time: 151067 s = 41.9632 h (63.8%)

  6. Where is Fluka spending time? from AliZDCv2::CreateMaterials() // Thresholds for showering in the ZDCs i = 1; //tantalum gMC->Gstpar(idtmed[i], "CUTGAM", .001); gMC->Gstpar(idtmed[i], "CUTELE", .001); gMC->Gstpar(idtmed[i], "CUTNEU", .01); gMC->Gstpar(idtmed[i], "CUTHAD", .01); …… i = 5; //lead gMC->Gstpar(idtmed[i], "CUTGAM", .001); gMC->Gstpar(idtmed[i], "CUTELE", .001); gMC->Gstpar(idtmed[i], "CUTNEU", .01); gMC->Gstpar(idtmed[i], "CUTHAD", .01); …. // Avoid too detailed showering along the beam line i = 7; //iron with energy loss (ZIRON) gMC->Gstpar(idtmed[i], "CUTGAM", .1); gMC->Gstpar(idtmed[i], "CUTELE", .1); gMC->Gstpar(idtmed[i], "CUTNEU", 1.); gMC->Gstpar(idtmed[i], "CUTHAD", 1.); • Small hadronic cut (1MeV) on ZDC calorimeter (PIPE, etc) allow the develop of detail hadronic showers. • Those hadrons escape the volume and develop the EM part also in ALIC

  7. Where is Fluka spending time? AliMC::Stepping() timer.Stop() TH3F->Fill(x,y,z,timer.CpuTime()); ... timer.Start() } ZDC (ZEM) , PIPE, Compensator DIPO =>ALIC ZEM T0

  8. ALIC Where is Fluka spending time? (ZDC (ZEM) , PIPE, Compensator DIPO) ->ALIC 1100 cm pi- P=92Gev #secondaries=17626 transport time=134s (2.3 min)

  9. Where is Fluka spending time? ZDC and PIPE => ALIC ZDC ZNST ZNSL ZPST TH3F->Fill(x,y,z,timer.CpuTime());

  10. Where is Fluka spending time? ZDC and PIPE -> ALIC

  11. To reduce the transport time in ALIC we introduced new volumes with high Energy cuts for gammas an electrons (0.1GeV) ALIZ DCAR PIAR ALIC Shielding around dipole compensator

  12. We change from to

  13. Final Results # Primaries with (transp. time > 30s): 33 (0.05%) with a total time: 8497 s = 2.3 h (25.5%)

  14. TODO • Check how the changes could affect the ZDC performance.

  15. END

  16. Where is Fluka spending time? ALIC ZDC

  17. Final FLUKA GEANT

More Related