Corroboration of Neutron Capture Rates in Rock for Physics Study

Investigating neutron capture rates in rock for physics simulations to expand the understanding of low-energy physics. Utilizing specialized software and grid computing, the study aims to validate capture rates and spectra to enhance the physics scope and fidelity of the simulations.

• Physics Study
• Neutron Capture
• Rock Simulation
• Low-Energy Physics
• Data Analysis

tan_bel Follow

Uploaded on Feb 21, 2025 | 0 Views

Download Presentation

Please find below an Image/Link to download the presentation.

The content on the website is provided AS IS for your information and personal use only. It may not be sold, licensed, or shared on other websites without obtaining consent from the author. Download presentation by click this link. If you encounter any issues during the download, it is possible that the publisher has removed the file from their server.

Presentation Transcript

1. Validation or Corroboration of external n s Eric Church July 22, 2020

2. Quick check of n production from Rock See Aran s presentation from Monday PNNL runs these simulations, as well. Our interests are in altering the input spectra for different alpha,n models and generally checking n production from various detector elements Investigating backgrounds to low energy physics which may be open to DUNE to expand/strengthen its physics scope These form the input into our BackgroundExplorer (Chris J s talk) Before we do that we want to show we can corroborate Sussex and SDSMT capture rates from rock neutrons. We use the new full 10 kT geometry and decay neutrons from as deep as 2m into the cavern rock 2

3. Corroboration details We have a larsoft installation at PNNL Currently dunetpc v08_37_00 With updated gdml and fcl files from Sussex Want to use full geometry and proper input spectra Looking at Th_tot and U_late_total here. Look at fiducialization We run on the PNNL grid True information from G4 decays and propagation No reco here, as is true with Aran s work too. I use my own ana module to produce an analysis TTree Then run a python script to chain up the 100 output Ttrees w their 200 evts each and plot the output capture spectra. 3

4. Spectra of Th n-capture gammas in Argon fid. Vol. 1000 1200 200 400 600 800 0 0 MCPdg==22&&MCProcess ==\"nCapture\"&&MCMater ial==\"LAr\"&&abs(MCOrigi n[0])<695.&&abs(MCOrigin [1])<600.&&MCOrigin[2]>5 0&&MCOrigin[2]<5950" MCMomentum[3] {MCPdg==22&&MCProcess=="nCapture"&&MCMaterial=="LAr"&&abs(MCOrigin[0])<950.&&abs(MCOrigin[1])<600.&&MCOrigin[2]>50&&MCOrigin[2]<5950} 0.002 0.004 0.006 0.008 MCMomentum[3] 0.01 Std Dev 0.002704 Std Dev 0.002704 Mean Mean Entries Entries htemp htemp 0.002995 0.002995 8898 8898 ~1500 evts summed above 6.1 MeV GeV 4

5. Thorium simulation tsim = Nneutrons/m_rock/Rate / IntSpec m_rock = 7.79E10 gm Rate = 7.11 ppm IntSpec = 1.70E-11 [n/sec/gm/ppb] Ncap = 1500 tsim = 20291296 / 7.79E10 /1000 / 7.11 / 1.70E-11 = 2155.0 sec Capture Rate = 1500/tsim = 0.69 Hz 5

6. Late Uranium simulation tsim = Nneutrons/m_rock/Rate / IntSpec m_rock = 7.79E10 gm Rate = 3.44 ppm IntSpec = 3.23E-11 [n/sec/gm/ppb] Ncap = 1757 {5.0,10.0 MeV} tsim = 20293590 / 7.79E10 /1000 / 3.44 / 3.23E-11 = 2344 sec Capture Rate = 2159/tsim = 0.92 Hz Compare to Aran s 1.4 Hz. 6

7. X coordinate for captures 1000 1200 1400 1600 1800 2000 2200 2400 200 400 600 800 0 800 - MCOrigin[0] {MCPdg==22&&MCProcess=="nCapture"&&MCMaterial=="LAr"&&abs(MCOrigin[0])<1200&&abs(MCOrigin[1])<600.&&MCOrigin[2]>50&&MCOrigin[2]<5950} neutron capture gammas may clearly be fiducialized away. 600 - 400 - 200 - 0 200 400 600 MCOrigin[0] Std Dev 680.5 Std Dev 680.5 Mean Mean Entries Entries 800 htemp cm htemp 19.64 19.64 8748 8748 7

8. Input for BackgroundExplorer db Need to form a json file that holds the gamma spectrum and metadata about this part of the simulated background model. Then upload the record into the BackgroundExplorer mongo db. We use our j4bgexp.py to do this. chur558@constance03 ana]$ python j4bgexp.py chur558@constance03 ana]$ more bgdexp_Rock-neutrons.json {"nentries": 2000000, "nprimaries": 151556, "ncaptures": 10, "hits": [[0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 1, 0, 2, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0.0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6000000000000001, 0.700000 0000000001, 0.8, 0.9, 1.0, 1.1, 1.2000000000000002, 1.3, 1.4000000000000001, 1.5, 1.6, 1.7000000000000002, 1.8, 1.9000000000000001, 2.0, 2.1, 2.2, 2. 3000000000000003, 2.4000000000000004, 2.5, 2.6, 2.7, 2.8000000000000003, 2.9000000000000004, 3.0, 3.1, 3.2, 3.3000000000000003, 3.4000000000000004, 3 .5, 3.6, 3.7, 3.8000000000000003, 3.9000000000000004, 4.0, 4.1000000000000005, 4.2, 4.3, 4.4, 4.5, 4.6000000000000005, 4.7, 4.800000000000001, 4.9, 5 .0, 5.1000000000000005, 5.2, 5.300000000000001, 5.4, 5.5, 5.6000000000000005, 5.7, 5.800000000000001, 5.9, 6.0, 6.1000000000000005, 6.2, 6.3000000000 00001, 6.4, 6.5, 6.6000000000000005, 6.7, 6.800000000000001, 6.9, 7.0, 7.1000000000000005, 7.2, 7.300000000000001, 7.4, 7.5, 7.6000000000000005, 7.7, 7.800000000000001, 7.9, 8.0, 8.1, 8.200000000000001, 8.3, 8.4, 8.5, 8.6, 8.700000000000001, 8.8, 8.9, 9.0, 9.1, 9.200000000000001, 9.3, 9.4, 9.5, 9. 600000000000001, 9.700000000000001, 9.8, 9.9, 10.0]], "primary": "neutrons", "volume": "volworld_1", "sam_metadata": {"prodradiological_neutron_larfd _gen-dot-root": {"last_event": 200, "art-file_format_era": "ART_2011a", "application-version": "development", "group": "dune", "application-family": "art", "file_format": "artroot", "file_type": "mc", "art-first_event": [1, 2, 1], "start_time": "2020-07-20T19:56:52", "event_count": 200, "runs": [[ 1, 2, "physics"]], "art-process_name": "MARLEYGen", "first_event": 1, "end_time": "2020-07-20T22:21:22", "art-file_format_version": 13, "data_tier": "generated", "art-last_event": [1, 2, 200], "data_stream": "out1", "art-run_type": "physics"}}, "config_metadata": {"largeant": {"DumpParticleList": false, "GeantCommandFile": "LArG4.mac", "SmartStacking": 0, "module_type": "LArG4", "VisualizeEvents": false, "DebugVoxelAccumulation": 0, "KeepParti clesInVolumes": [], "CheckOverlaps": false, "DumpLArVoxelList": false}, "out1": {"compressionLevel": 1, "module_type": "RootOutput", "streamName": "o ut1", "dataTier": "generated", "fileName": "prodradiological_neutron_larfd_gen.root"}, "source": {"firstRun": 1, "firstSubRun": 2, "timestampPlugin": {"plugin_type": "GeneratedEventTimestamp"}, "maxEvents": 200, "module_type": "EmptyEvent", "firstEvent": 1}, "largana": {"module_type": "LArG4Ana", "GeantModuleLabel": "largeant", "TruthModuleLabel": "cNeutronGen", "Ndaughters": 20}, "cNeutronGen": {"BqPercc": [7.6e-06], "Nuclide": ["Rock_Late_U_ Total"], "T0": [-2246000.0], "T1": [2246000.0], "module_type": "RadioGen", "Material": ["RadioRockAverage"], "Y1": [2036], "Y0": [-1078], "X0": [-119 0], "X1": [1190], "Z0": [-736], "Z1": [14724]}, "rns": {"module_type": "RandomNumberSaver"}}, "spectrum": "Rock_Late_U_Total"} 8