1 / 34

BESIII Simulation Status

BESIII Simulation Status. Z.Y. Deng M. He H.M. Liu Y.J. Mao Z.Y. You Y.Yuan BES Annual Meeting June 1, 2005. Outline. Introduction Current status Performance Future work Summary. BOOST Project ( B ESIII O bject O riented S imulation T ool). Event GENBES Generator.

nathaniel-ewing
Download Presentation

BESIII Simulation Status

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. BESIII Simulation Status Z.Y. Deng M. He H.M. Liu Y.J. Mao Z.Y. You Y.Yuan BES Annual Meeting June 1, 2005

  2. Outline • Introduction • Current status • Performance • Future work • Summary

  3. BOOST Project (BESIII Object Oriented Simulation Tool) Event GENBES Generator Geometry Geant4 Tracking Detector Digitization Response Raw data MC truth HepEvt format Hit objects

  4. BOOST Project (BESIII Object Oriented Simulation Tool) • Components • Generator • Particles and physics processes • Magnetic field • Material • MC truth • Data I/O • User interface • Geometry • Hits recording • Digitization common parts sub-detector parts

  5. BOOST code structure boost/ bes/ common part src/ source c++ codes include/header files gen/ … event generators phy/ … physics processes tru/ … MC truth gmk/ … common makefiles mac/ … common card files dat/ … geometry data files exe/ … main program for execution doc/ … documentations mdc/ …main drift chamber tof/ …time of fight emc/ …electromagnetic calorimeter muc/ …muon chamber

  6. BOOST Project (BESIII Object Oriented Simulation Tool) • Developers • Liu Huaimin: common parts • Yuan Ye: MDC • Deng Ziyan: TOF & common parts • Fu Chengdong/He Miao: EMC • Peking University: MUC & xml geometry

  7. BOOST Project (BESIII Object Oriented Simulation Tool) • Phase 1: Framework prototype, simple geometry with hit • Finished, BOOST prototype realized in May,2003 • Phase 2: Detailed detector with simple digitization • Finished, now digits and MC truth output can be used for reconstruction • Phase 3: Detailed digitization, data/MC comparisons • In progress, need information from beam test and data. A long-term work!

  8. Geant4 and BOOST releases • Geant4 • Release 6.1 March, 2004 • Release 6.2 June, 2004 • Release 7.0 December, 2004 • Patch of 7.0 February, 2005 • BOOST (mantained by cvs) • ‘boost-1-0’ (March, 2004) -----Geant4.6.1 • ‘boost-1-1’ (June, 2004) -----Geant4.6.2 • …… • ‘boost-2-0’ (February, 2005) ----Geant4.7.0 • ‘boost-2-1’ (May, 2005) ----- Geant4.7.0+p01

  9. BOOST Current Statusgeometry • Detailed BESIII geometry is constructed in BOOST • Geometry added since last annual meeting • Accurate stereo cells in MDC • New End-cap EMC • SCM (Super-Conducting Magnet ) • BESIII Geometry in XML format • Geometry in BOOST now can also be constructed with XML • Simulation, reconstruction, event display will share the same BESIII geometry • See detailed information in Mao Yajun’s report

  10. BOOST Current StatusMDC stereo cell • Accurate description • New class (J4TwistedTubs) introduced • all stereo wire cells are constructed • Display for checking • Twisted tube can’t be displayed directly • hit positions in stereo cell used • 3D scatter plot with ROOT package

  11. BOOST Current Status end-capEMC geometry • Quadrangular crystal • Irregular boxoid • Can’t be described by existing classes in Geant4 releases • New class (G4IrregBox) introduced • Pentangular crystal • Divided into 2 boxoids • 2 boxoids share the same copy number quadrangular crystal pentangular crystal

  12. BOOST Current Statusend-capEMC geometry • construct all crystals in one sector • Place other sectors through rotation and reflection • Rotate east end-cap to place west end-cap end-cap EMC rotation reflection sector

  13. BOOST Current Statusgeometry From boost-2-1 end-cap:TOF, EMC, MUC barrel: MDC, TOF, EMC, SCM, MUC

  14. BOOST Current Statusgenerator • genbes: (BESII generator system) • Generate events in an ASCII file which supports HEPEVT interface with Geant4 • The same user interface as BESII • genbes.cards , genbes.user • All BESII generators can be used in BOOST • tester: (single-particle generator for debugging) • Shoot particles of given type • to a given direction • with given kinetic energy • Users can randomize above quantities

  15. BOOST Current Statusgenerator: howl Hits in MDC Hits in barrel TOF (no threshold cut) Hits in end-cap TOF Hits in end-cap MUC J/ψKs Ksπ+π-π+π-

  16. BOOST Current Statusgenerator: rhopi gamma Hits in barrel MUC (π- punch through) J/ψρπ π+π-  Shower in EMC

  17. BOOST Current Statusprocesses & magnetic field & material • Physics processes • BesPhysicsList constructed with Geant4 classes • Many hadronic models installed, LHEP, LHEP-GN, QGSP-GN, …… • More research on them needed to choose suitable ones for BESIII physics • Magnetic field • 1 Tesla uniform magnetic field defined inside SCM • Material • All materials needed in BOOST defined in XML

  18. BOOST Current Statushits and digitization • Hits recording • Finished last year • digitization • MDC • Drift distance-> drift time • dE/dx • Wire resolution added • Simple background added • TOF • Detailed digitization • Light emission & light propagation • PMT response & threshold discrimination • EMC • waveform added • Simulate main amplifier output • Provide risetime to eliminate noise • No noise, so not used now • MUC • Digit information added

  19. BOOST Current StatusMC truth • Purpose • To provide truth information for reconstruction debugging • Similar to MCMADE in BESII • What to save • True state of the particles and their association with detector response • Particle • particles from generator • uninterrupted decayed daughters • Detector response • hits in sub-detectors associated with particles

  20. Track and vertex (associated with particle) Track PDGcode charge original vertex terminal vertex four-momentum decayed daughters Vertex parent track position time truth hits (associated with track) Considering demand of each reconstruction developer MDC Hit history TOF First-hit EMC Total energy loss MUC Hit history BOOST Current StatusMC truth

  21. BOOST Current StatusData I/O • Ascii data I/O • It’s temporary, but useful in software (SIM + REC) development • Ascii data format fixed for the moment • MC truth (TRUTH) & Raw data (DIGI) included in one file • EVHEAD • TRACKTRUTH, VERTEXTRUTH • MDCTRUTH, MDCDIGI • TOFTRUTH, TOFDIGI • EMCTRUTH, EMCDIGI • MUCTRUTH, MUCDIGI • Users can turn on/off TRUTH and DIGI part of each sub-detector in run cards

  22. BOOST Current StatusUser interface • Friendly to both users and developers • source ~dengzy/.boostenv • cvs co boost • edit boost.cards & boost.user • more controls in run cards • Run background/terminal/interactive/PBS job freely • boost -b/-t/-i/-q • use ‘boost –h’ for detail • Reconstruction developers can generate needed data files themselves

  23. From LHC experience G4 direct comparison with G3 G4 is slower (1.5~2) than G3 G4 tracking particle to 0 energy Compared with SIMBES SIMBES (G3/BESII), BOOST (G4/BESIII) How to compare Similar cuts CUTE CUTM (10KeV), range cuts (1mm) Same platform (Besfarm7) Physics events used (dimuon, rhopi, bhabha) BOOST PerformanceSpeed

  24. BOOST PerformanceSpeed • BOOST 1.5 times slower than SIMBES • More to do about speed improvement • Threshold cut in digitization • Speed of hadronic models to be compared • Stable (no crash for 100,000 events)

  25. BOOST PerformanceMDC Reconstruction result using BOOST data sample P distribution (1GeV/c mu-) σ p=5.8MeV P distribution (1GeV/c e-) σp =6.1MeV

  26. BOOST PerformanceTOF 1GeV/c e- vertical incidence total energy loss in one scintillator (5cm width)

  27. BOOST PerformanceTOF T : 1/sqrt(Q) time-walk forward backward e- 1GeV/c Z=0 : forward TDC resolution = 121ps backward TDC resolution = 122ps

  28. BOOST PerformanceTOF e μ π k p • Randomly generate particles: e-, mu-, pi-, kaon-, proton • Momentum: uniform between 0-1.6 GeV/c • Only intrinsic time resolution

  29. BOOST PerformanceEMC total energy loss of pi-e-

  30. BOOST PerformanceMUC Probability of single particle(μπ) entering MUC

  31. Different developing environment BOOST On besfarm Configured with gmake Based on GEANT4 BOSS On koala Configured with CMT Based on GAUDI Integration Status BOOSTAlg algorithm/packages created in BOSS Key functions of BOOST realized in BOSS Interface of saving data to TDS is ok Codes imported into cvs BOOST Integration to BOSSframework

  32. Future Work • Move developing environment from besfarm to lxplus • Magnetic field mapping • Physics validation, optimization • Data I/O • Objectivity? • Raw data in online format, MC-Truth format? • Background, noise? • Detailed detector response

  33. Summary • Detailed BESIII detector is built in BOOST • All BESII generators can be used in BOOST • Raw data/MC-truth data are successfully used in reconstruction • Friendly user interface • BOOST integrated into BOSS • Prototype of TOF detailed digitization is realized • More to do about detector response !

  34. Thank you ! 

More Related