EIC software Alexander Kiselev NPPS Group Meeting August,23 2019 - PowerPoint PPT Presentation

EIC software Alexander Kiselev NPPS Group Meeting August,23 2019

Contents of this talk n Fast simulation tool: eic-smear n Software frameworks n GEMC n fun4all n EicRoot n Argonne EIC software initiative n PID consortium GEANT4 software (one slide) n Near-term future trend(s) Materials taken from presentations of T.Burton, M.Ungaro, E.Sichtermann, J.Repond, D.Lawrence, D.Romanov, Y.Furletova & others

eic-smear by Tom Burton (BNL TF group)

Overview Overview • C++ code, runs in ROOT • Build with configure/Make or CMake • libeicsmear.so to load in ROOT Tree code : Smearer : MC Build ROOT Perform fast generator tree containing detector ASCII output events smearing Large number PYTHIA DPMJet of EIC Monte PEPSI Milou gmc_trans Carlo generators with standard Rapgap LEPTO Djangoh ASCII format

Smearing ‣ Built-in standard smearers “ Smearer ” defines some provided with eic-smear element of performance ‣ Users can define own + acceptance smearers using inheritance “Detector” Smearer NOT a “physical Smearer ‣ Apply all smearers to an detector”: Smearer represents the MC event overall Smearer ‣ Yield smeared event performance in measuring Smearer ‣ Optionally recalculate a quantity. derived values e.g x, Q 2

How to use it Simple “Device” smearers define σ (X) • Write a ROOT script: via text string Smear::Detector createDetector() { // Resolution in momentum, sigma(P). // sigma(P) = 0.4%P + 0.3%P^2. Smear::Device tracking(“P”, “0.004 * P + 0.003 * pow(P, 2)”); // Add devices to a Detector. Smear::Detector detector; detector.AddDevice(tracking); Handles event return detector; loop, file I/O } • Smear your ROOT tree: root[0] SmearTree(createDetector(), “mc.root”, “smeared.root”); • “Standard” detector descriptions (like STAR or BeAST) exist See K.Kauder: talk at the EIC software meeting 07/10/2019

GEMC by Maurizio Ungaro (JLab)

GE ant4 ant4 M on onte C arlo rlo Arch Archit itect ecture re • Application independent geometry/digitization/fields: definitions stored in databases • Realistic hits treatment: electronic time window, voltage versus time signals. • Sensitive attributes assigned at run time: real calibration, survey tilts and displacements. • Plugins for generator formats (LUND, BEAGLE, easy expansion) • Plugins for output formats (TXT, CODA, JSon, easy expansion) • Realistic signal treatment allows for background rate studies, • Application for detector simulations including pile-up effects based on Geant4 • Macro language for detector design • Various geometry definitions: GEMC, gdml, CAD • Data card (XML) to steer application, all Geant4 macro commands supported by design

Geometry Native Input: Native, CAD, GDML. Arbitrary hierarchy, can be mixed and matched. Materials, sensitivity assigned at run-time. CAD GDML Experiments using the GEMC Framework: CLAS12 (Hall-B), EIC Beamline and detectors, HPS, Solid

Digitization, Output • Single ADC/TDC over electronic time window. • Voltage vs time signal. > BST • FADC output (4ns > True Step by Step infos (101, 0) intervals or integratal - Edep (101, 1) mode) - Pid (101, 2) • Automatic true - positions (101, 3) information > Dgtz Step by Step infos (102, 0) • All g4 steps in the - ADCL (102, 1) - ADCR (102, 2) output > True Integrated infos (103, 0) - Edep (103, 1) - Pid (103, 2) - positions (103, 3) > Dgtz Integrated infos (104, 0) - ADCL (104, 1) - ADCR (104, 2) > Voltage as a function of time (105, 0) - Identifier (105, 1) - Time (105, 2) - Voltage (105, 3) > Trigger Bank (106, 0) - Identifier (106, 1) - Time (106, 2) - Voltage (106, 3)

Graphical Interface • Geant4 • Generator OpenGL • Event time View for the window whole • Background detector. beams • Can inspect • Camera and open a views slices. view on • Axis, Scale, single Show field. volumes. • Graphical • Volumes analysis of hierarchies steps in a and hit. • Can choose properties • Output to variable to GDML display.

fun4all by Chris Pinkenburg (BNL) See talk at the EIC software meeting 07/10/2019

EicRoot by AK (BNL)

EicRoot framework building blocks Interface to GEANT, ROOT, … n FairBase PandaRoot FopiRoot TPC R&D stuff, … n “Ideal” track finder, n Interface to GenFit n EicRoot … n CbmRoot eic-smear MC generated evts import RICH stuff n n Fast smearing codes solenoid IR design n modeling configuration -> basically a yet another FairRoot software clone

End user view n No executable (steering through ROOT macro scripts) simulation digitization reconstruction PID; assembly -> MC points -> Hits -> Tracks, clusters -> Events n ROOT files for analysis available after each step n C++ class structure is well defined at each I/O stage See AK: talk at the EIC software meeting 07/10/2019

Example case studies 7 /P [%] 12 GeV pions: Hcal vs EmCal 50 GeV/c 6 25 GeV/c 10 GeV/c P σ 5 Calorimeter design optimization 1 GeV/c Momentum resolution 4 Tracker momentum slope ~1.20 resolution 3 2 1 0 -3 -2 -1 0 1 2 3 Pseudo-rapidity Neutron fluence DIS electron reconstruction

Current modeling work Possible central tracker configurations (alternatives to a TPC) n Temple University: µ RWELL µ TPC barrels LBNL: tapered all-silicon tracker

Current modeling work Outer forward GEMs Inner forward GEMs Si tracker TPC RICH Florida Tech: tracking resolution in the RICH volume BNL: far forward acceptance Events 450 100 x 10 GeV 400 350 300 first quad aperture & 250 beam pipe 200 150 Roman B0 spectrometer 100 Pots University of Birmingham: 50 vertex tracker optimization 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 DVCS proton P [GeV/c] t

PID Consortium software

Mostly RICH & ToF applications GEMC DIRC single module for JLEIC Dual radiator RICH Modular RICH All are custom GEANT4 codes n

Argonne EIC software

Full simula*on and reconstruc*on chain Event genera*on Produce the simula:on input events Detector simula*on Par:cle transport through detectors Data Model Digi*za*on Turn energy deposits in ac:ve media into detector hits Reconstruc*on of Event vertex, charged tracks, Par:cle Flow Objects (PFO) Perform analysis Collec:on of benchmark analyses

Argonne So>ware: Overview Legacy chain Evolu*on chain Adapta:on of the SiD (ILC) simula:on and Evolved from the legacy chain reconstruc:on soLware chain Geometry interface Major parts DD4HEP SLIC (wrapper around GEANT4) LCSIM (digi:za:on and event reconstruc:on) slicPandora (PFA reconstruc:on) Features Visualiza:on with JAS4pp Fully maintainable Geometry obtained from single source Limita:ons Geometry can be parametrized Geometry not constrained to be symmetric Only SiD subdetectors (e.g. no RICH) New subsystems can be easily implemented Geometry descrip:on not centralized Geometry constrained to be symmetric S:ll working on Some parts difficult to maintain Realis:c digi:za:on Full chain Generic tracking PFA reconstruc:on Available Visualiza:on Studies of F 2 reconstruc:on, :ming…

Nuclear Physics Detector Library (NPDet) Collec:on of parametrized detectors which can be developed into full concepts TOPSiDE

ProIO

Grand unification, yet another try by Dmitry Romanov, David Lawrence, Yulia Furletova & others (JLAB)

Key ingredients n (Docker) containers n Jupyter notebooks n JANA2 soLware framework n g4e GEANT-based EIC detector sandbox

Software distribution model(s) NO EFFORT AT ALL Some effort Novice Experts Efforts required axis Worksta:on Compila:on Cloud Containers A PC farm Conda EJPM Main focus at present

A side note: EIC Docker containers -> introduced in Aug,2017; went public by EICUG meeting in Nov,2017 Clear benefits for EIC user community • Allow EIC users to run the same software under standardized environment on any Linux, Mac OS or Windows machine, eventually including GRID sites, commercial cloud systems, and HPC resources • Provide consistency between software generated at different facilities • Make it easier for new users to start working on the physics program and detector design for the EIC, by minimizing the pain of “installation overhead”

Core functionality overview Generators Generators Database with various Database … Pythia Beagle Herwig MC samples (1) Fast simula:on Full simula:on g4e = Geant 4 EIC BNL & Jlab effort on g4e Geant4 Na:ve C++ GEANT4 code Fast detector prototyping JLEIC/eRHIC in Eic smear Fast mode with an EIC detector in it Geant4 (2) Reconstruc:on & analysis ejana = EIC Jana Community reference ejana – EIC JANA(2) reconstruc:on (3) (1) MC events (2) Digi:zed hits + magne:c field + material distribu:on (3) Reconstructed events -> user access (with graphics) either directly or through SSH or Web interface See D.Romanov: talk at the EIC soLware mee:ng 07/10/2019