Pyramid Simulation Ralf Ehrlich University of Virginia August

Pyramid Simulation Ralf Ehrlich University of Virginia August 9, 2019

Current Status Ø Simulations with a simple pyramid model with empty voids (spheres) of different diameters and detectors around the pyramid have been run. Ø Voids with diameters of 4 m or larger could be identified in 2 D reconstructions, if the right projection plane is selected. • However, fake objects from other planes can appear. Ø A 3 D reconstruction with information from detectors on two sides is necessary. Ø No new results since May, because I was occupied with working on Mu 2 e. Ø I will restart my pyramid efforts during my vacation end of August / beginning of September. • Create new lookup tables for 2 cm wide counters. • Run new jobs on the grid which involves not only the simulation, but also the reconstruction (which was previously not done on the grid) to speed things up. • Attempt a 3 D reconstruction. Ø The following slides are the slides from my previous talk. 5/6/2019 Ralf Ehrlich - University of Virginia 2

Setup of test run Ø Currently have a basic GEANT 4 model of • A pyramid with empty chamber (spheres with diameters between 1 m and 6 m) • The correct pyramid model will be implemented later. • The detector consisting of 2 x 2 containers placed centered at 25 m away from the base of the pyramid. • Container wall: 5 mm thick aluminum • Distance between left and right bank: 2 m • Length of vertical counters: 2. 4 m • Length of horizontal counters: 4. 8 m • 16 counters per module • 6 vertical modules per bank and array, 3 horizontal modules per bank and array 5/6/2019 Ralf Ehrlich - University of Virginia 3

Setup of test run GDML drawing of the GEANT 4 model y Height 138. 75 m The empty spheres inside of the pyramid are hidden. 2 x 2 containers 0. 33 m 2 x 2 containers m 23 3 m 230. 3 25 z 25 m x Containers (L: 40’, W: 8’, H: 9. 5’ ) made of 5 mm thick aluminum not shown in drawing. 5/6/2019 Ralf Ehrlich - University of Virginia 4

Setup of test run ØTest run setup • 12 m high by 26 m wide vertical production plane at the detector. • Cosmic rate 207 Hz for an energy interval between 50 Ge. V and 1 Te. V. • Starting point of the events are projected to the surfaces of the world volume, so that GEANT can simulate the muon trajectories through the pyramid. • Production plane at +x: generated 3. 2 e 9 events, which corresponds to 178 days. • Production plane at +z: generated 1. 2 e 9 events, which corresponds to 68 days. ØBaseline run • Repeat the simulation for a pyramid without voids. • Generated 3. 9 e 9 events, which corresponds to 218 days. ØReconstruction • Use only events where only one of the 5 cm wide counters per layer has hits (vertical and horizontal layer at both detector banks). • Calculate direction and location of track (at center between both detector banks). • Projected the tracks back to a vertical plane at the center of the pyramid. • Plot the projected track locations ØResults • Subtract the plot of the baseline run from the plot of the test run. 5/6/2019 Ralf Ehrlich - University of Virginia 5

Test run for detector at +x side (before baseline subtraction) Ø Projection to the center of pyramid (x=0) • Using MC information from the end of the tracks. • No energy cut. • 1 m/bin Simulated empty spheres 5/6/2019 Ralf Ehrlich - University of Virginia 6

Test run for detector at +x side (after baseline subtraction) Ø Projection to the center of pyramid (x=0) • Using MC information from the end of the tracks. • No energy cut. • 1 m/bin Simulated empty spheres After applying TH 2 F: : Smooth() and making a cut on all bins with entries less than 70 5/6/2019 Ralf Ehrlich - University of Virginia 7

Test run for detector at +x side (before baseline subtraction) Ø Projection to the center of pyramid (x=0) • Using ”Pseudo Reco” information for 2 cm wide counters. • Based on MC information from the end of the tracks with discrete y and z positions (in steps of 2 cm) on both sides of the detector, because I don’t have the lookup tables for 2 cm wide counters, yet. • No energy cut. • 1 m/bin Simulated empty spheres 5/6/2019 Ralf Ehrlich - University of Virginia 8

Test run for detector at +x side (after baseline subtraction) Ø Projection to the center of pyramid (x=0) • Using ”Pseudo Reco” information for 2 cm wide counters. • Based on MC information from the end of the tracks with discrete y and z positions (in steps of 2 cm) on both sides of the detector, because I don’t have the lookup tables for 2 cm wide counters, yet. • No energy cut. • 1 m/bin Simulated empty spheres After applying TH 2 F: : Smooth() and making a cut on all bins with entries less than 60 5/6/2019 Ralf Ehrlich - University of Virginia 9

Test run for detector at +z side (after baseline subtraction) Ø Projection to the center of pyramid (z=+20 m) • Using MC information from the end of the tracks. • No energy cut. • 1 m/bin Simulated empty spheres After applying TH 2 F: : Smooth() and making a cut on all bins with entries less than 30 5/6/2019 Ralf Ehrlich - University of Virginia 10

Test run for detector at +z side (after baseline subtraction) Ø Projection to the center of pyramid (z=-20 m) • Using MC information from the end of the tracks. • No energy cut. • 1 m/bin Simulated empty spheres After applying TH 2 F: : Smooth() and making a cut on all bins with entries less than 30 5/6/2019 Ralf Ehrlich - University of Virginia 11