Simulation Results for CLAS12 From gemc

Similar documents
Recent Developments for the CLAS12 Simulation

Study of the Neutron Detection Efficiency of the CLAS12 Detector. Keegan Sherman

The reaction p(e,e'p)π 0 to calibrate the Forward and the Large Angle Electromagnetic Shower Calorimeters

Validation of Geant4 Physics Models Using Collision Data from the LHC

CLAS12 at Jefferson Lab

Neutrons in a Spin: Nucleon Structure at Jefferson Lab

Beyond the Born Approximation

Tagging with Roman Pots at RHIC. Proton tagging at STAR

DESY Summer Students Program 2008: Exclusive π + Production in Deep Inelastic Scattering

HASPECT in action: CLAS12 analysis of

Full-Acceptance Detector Integration at MEIC

Cascade Spectroscopy at CLAS

Identification of Central Production in the π + π π + π Channel at COMPASS

Future prospects for the measurement of direct photons at the LHC

Threshold photoproduction of J/y with the GlueX experiment. Lubomir Pentchev Jefferson Lab for the GlueX collaboration

Recent results at the -meson region from the CMD-3 detector at the VEPP-2000 collider

pp physics, RWTH, WS 2003/04, T.Hebbeker

Physics potential of ATLAS upgrades at HL-LHC

Jet Reconstruction and Energy Scale Determination in ATLAS

Muon Alignment at CMS

Particle Identification Algorithms for the Medium Energy ( GeV) MINERνA Test Beam Experiment

Muon reconstruction performance in ATLAS at Run-2

CMS Simulation Software

CLAS12 DDVCS. N. Baltzell (Jefferson Lab) ECT Workshop, October 24-28, 2016

Tracking at the LHC. Pippa Wells, CERN

Time-like Compton Scattering with transversely polarized target

The Mystery of Vus from Tau decays. Swagato Banerjee

PoS(KAON09)023. Beam Hole Photon Veto For J-PARC K O TO experiment. Yosuke Maeda Kyoto University

CALICE Test Beam Data and Hadronic Shower Models

Current Status of Luminosity Measurement with the CMD-3 Detector at the VEPP-2000 e + e collider

Status / Hadron Spectroscopy at COMPASS

1.5 TeV Muon Collider background rejection in ILCroot Si VXD and Tracker (summary report)

Event Builder Status. N. Baltzell CLAS Collaboration Meeting November 13, 2018

HERMES status and future running

The Detector Design of the Jefferson Lab EIC

12 GeV Upgrade Project DESIGN SOLUTIONS DOCUMENT. Upgrade Hall B

A New Detector for Physics at HERA - III

Physics with Tagged Forward Protons using the STAR Detector at RHIC. The Relativistic Heavy Ion Collider The pp2pp Experiment STAR 2009

Overview of validations at LHC

Some studies for ALICE

Commissioning of the ATLAS LAr Calorimeter

Preliminary Detector Design for the EIC at JLab

Absolute D Hadronic Branching Fractions at CLEO-c

Simulation study of scintillatorbased

A NEW TECHNIQUE FOR DETERMINING CHARGE AND MOMENTUM OF ELECTRONS AND POSITRONS USING CALORIMETRY AND SILICON TRACKING. Qun Fan & Arie Bodek

LHCb: Reoptimized Detector & Tracking Performance

High-energy Gamma Rays detection with the AMS-02 electromagnetic calorimeter. F. Pilo for the AMS-02 ECAL Group INFN Sezione di Pisa, Italy

Calorimetry in particle physics experiments

SciBar and future K2K physics. F.Sánchez Universitat Aútonoma de Barcelona Institut de Física d'altes Energies

Two Photon Exchange (TPE) and the Proton Form Factor Problem

Multi Lepton events at HERA

Hall B Physics Program and Upgrade Plan

Study of Baryon Form factor and Collins effect at BESIII

Particle Production Measurements at Fermilab

Particle production vs. energy: how do simulation results match experimental measurements?

High-t Meson Photoproduction: Experimental Capabilities

CALICE scintillator HCAL

PoS(PSTP 2013)034. Precession Polarimetry at JLab, 6 GeV. G.B. Franklin Carnegie Mellon University

Identifying Particle Trajectories in CMS using the Long Barrel Geometry

THE FIFTH STRUCTURE FUNCTION. Liam Murray Research Advisor: Dr. Gerard Gilfoyle

Studies of the diffractive photoproduction of isolated photons at HERA

ATLAS NOTE. August 25, Electron Identification Studies for the Level 1 Trigger Upgrade. Abstract

9/27 JUNE 2003 SUMMER STAGE PARTICLES REVELATION THROUGH CERENKOV AND SCINTILLATION COUNTER AND THE CEBAF EXPERIMENT

Measurement of charged particle spectra in pp collisions at CMS

The Electric Form Factor of the Neutron for SBS

CsI Calorimeter for KOTO experiment

Deeply Virtual Compton Scattering off 4. He: New results and future perspectives

Hall B Physics Program and Upgrade Plan

Concepts of Event Reconstruction

CLEO c. Anders Ryd Cornell University June 7, e e cc D D. D K,D K e

Cherenkov Detector Simulation

RTPC Simulation Studies for Tagged Deep Inelastic Experiment. Rachel Montgomery SBS Collaboration Meeting, Jefferson Lab, 22/07/16

MDI and detector modeling

Time-like Compton Scattering with transversely polarized target

Digital Calorimetry for Future Linear Colliders. Tony Price University of Birmingham University of Birmingham PPE Seminar 13 th November 2013

Measurement of the Neutron Magnetic Form Factor at High Q 2 Using the Ratio Method on Deuterium

LHCb Calorimetry Impact

Study of di-muon production with the ZEUS detector at HERA

MiniBooNE Progress and Little Muon Counter Overview

Luminosity determination in pp collisions using the ATLAS detector at the LHC

Proton Form Factor Puzzle and the CLAS Two-Photon Exchange Experiment

Stefan Diehl. 2nd Physics Institute, Justus-Liebig-University Giessen

Conference Report Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland

NA62: Ultra-Rare Kaon Decays

Charged Particle Identification in GLUEX

CMS Note Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland

Results From The HARP Experiment

a medium energy collider taking nucleon structure beyond the valence region

Lecture LHC How to measure cross sections...

GEANT4 simulation of the testbeam set-up for the ALFA detector

Measuring Form Factors and Structure Functions With CLAS

Recent Results of + c + X and + b + X Production Cross Sections at DØ

arxiv: v1 [hep-ex] 20 Jan 2013

Plans to measure J/ψ photoproduction on the proton with CLAS12

GMp Experiment (E ): An update

arxiv:hep-ex/ v1 31 Dec 2001

Neutron Structure Function from BoNuS

Simulation for Proton Charge Radius (PRad) Experiment at Jefferson Lab1 Li Ye Mississippi State University For the PRad Collaboration The Proton Charg

Studies of the hadronic component of the photon light-cone wave function using exclusive di-pion events at HERA

Triplet polarimeter study

Transcription:

CLAS12 Software Workshop - May 25, 20 p. 1/2 Simulation Results for CLAS12 From gemc G.P.Gilfoyle, M.Ungaro et al. CLAS12 Software Group Outline: 1. gemc Overview 2. Neutron efficiency in first TOF panel. 3. Electromagnetic calorimeter (EC) performance. 4. Hadron properties in the Silicon Vertex Tracker (SVT) 5. DVCS events in the SVT. 6. Background rates and Doses in the SVT. 7. CLAS12 Forward Tagger. 8. Conclusions.

CLAS12 Software Workshop - May 25, 20 p. 2/2 CLAS12 Simulation Essential tool for design and analysis (e.g. acceptance calculations). Quality of the results may be limited by systematic uncertainties (not statistics). Will need about four times as much Monte Carlo data as CLAS12 collects. gemc The CLAS12 detector. Modern, object-oriented, Geant4-based simulation. Still in development stage. Needed subsystems: Pre-shower calorimeter (PCal). Cerenkov counter (CC).

CLAS12 Software Workshop - May 25, 20 p. 3/2 Neutron Efficiency in First TOF Panel Neutron detection efficiency (NDE) needed for E12-07-4 (G n M experiment) and others. Study first TOF panel (OTOFP1A); we have CLAS6 NDE measurements for these scintillators. Generate elastic en events and process with gemc. Reconstruct electrons with Socrat. Found neutrons - Events where neutron is predicted to hit OTOFP1A using only electron information. 0 0 2 4 6 8 p (GeV/c) 2009-11-13 :30:05 n Reconstructed neutrons - Apply cut on angle between the predicted neutron direction and vertex-to-tof-hit vector (θ γ < ) and require E dep > 5 MeV. Counts 220 200 180 160 140 120 0 80 60 40 20 Black - Found neutrons Red - Reconstructed neutrons

CLAS12 Software Workshop - May 25, 20 p. 4/2 Neutron Efficiency in First TOF Panel Simulated NDE is about % at 6-8 GeV/c versus 9% at 3 GeV/c for CLAS6. Neutron Detection Efficiency 0.5 Preliminary CLAS12 Neutron 0.4 0.3 0.2 0.1 0 Detection Efficiency E = 11 GeV, θ e o o = 20-30 o < > 5.0 MeV, E γ E dep Elastic en events in OTOFP1A 5 6 7 8 9 (GeV/c) p n Neutron Efficiency 0.14 0.12 0.1 0.08 0.06 0.04 0.02 0 2.6 GeV 4.2 GeV CLAS6 NDE from E5 0 0.5 1 1.5 2 2.5 3 3.5 4 Neutron Momentum (GeV) To be done: timing cuts, p n dependence, simulate 1 H(e,e π + n) reaction for tagged neutrons.

CLAS12 Software Workshop - May 25, 20 p. 5/2 Electromagnetic calorimeter (EC) performance The CLAS6 EC will be re-used in CLAS12; placed 2 m downstream. EC geometry.

CLAS12 Software Workshop - May 25, 20 p. 6/2 Electromagnetic calorimeter (EC) performance The CLAS6 EC will be re-used in CLAS12; placed 2 m downstream. EC geometry.

CLAS12 Software Workshop - May 25, 20 p. 7/2 Comparison with CLAS6 - Sampling Fraction 1. CLAS6 results from EC NIM paper (NIM A460 (2001) 239-265) using GSIM. 2. Black curve in each plot based on ionization energy loss in materials preceding the EC. 3. Caveats: (1) electrons passed through different material, (2) B field off in the CLAS12 simulation and B field on in the CLAS6 simulation, (3) tracking used in CLAS6 simulation while Monte Carlo information used for CLAS12, (4) energy loss out the back,... 4. Note different horizontal ranges. EC Sampling Fraction in gemc EC Sampling Fraction 0.29 0.28 0.27 0.26 0.25 0.24 0.23 0.22 Red - No threshold > 1 MeV Blue - E dep Black curve - 0.274(E - 21 MeV)/p 500 events 0.21 2000 4000 6000 8000 000 p (MeV) 20-03-15 22:45:28 e gemc and GSIM sampling fractions are consistent to -15%.

CLAS12 Software Workshop - May 25, 20 p. 8/2 Comparison with CLAS6 - Resolution 1. CLAS6 results from EC NIM paper (NIM A460 (2001) 239-265) using GSIM. 2. Caveats: (1) electrons passed through different material, (2) B field off in the CLAS12 simulation and B field on in the CLAS6 simulation, (3) tracking used in CLAS6 simulation while Monte Carlo information used for CLAS12, (4) energy loss out the back,... 3. Note different horizontal ranges. σ (E/p) 0.06 0.05 0.04 0.03 Blue - E > 1 MeV Red - No threshold σ E = 0.08 0.02 0.01 0 2000 4000 6000 8000 000 p (MeV) 20-03-15 23:25:57 e gemc and GSIM resolutions are consistent to 5-%.

CLAS12 Software Workshop - May 25, 20 p. 9/2 Comparison with CLAS6 - Shower Size 1. CLAS6 shower size 4 cm diameter for 2.4 GeV electron (e.g. 95% of shower energy contained within 2 cm of centroid (EC NIM paper). 2. Caveats: (1) electrons passed through different material, (2) B field off in the CLAS12 simulation and B field on in the CLAS6 simulation, (3) tracking used in CLAS6 simulation while Monte Carlo information used for CLAS12, (4) energy loss out the back,... 3. Use E dep -weighted distribution in local x for θ e = 25, φ e = 0. E dep -weighted Hits 5 4 3 gemc E dep > 1 MeV p e = 4 GeV Shower Size (mm) 30 25 20 15 gemc Red - No Threshold, t < 50 ns Blue - E dep > 1 MeV, t < 50 ns N = 500 events 2 1-500-400-300-200-0 0 0 200 300 400 500 lx (mm) 5 0 20-03-16 16:00:33 2 4 6 8 (MeV) gemc and GSIM shower sizes are consistent to 5-%. p e 20-03-16 16:21:18

CLAS12 Software Workshop - May 25, 20 p. /2 Simulations of the Silicon Vertex Tracker (SVT) BST Geant4 Implementation Geant4 Design

CLAS12 Software Workshop - May 25, 20 p. 11/2 DVCS events in the SVT. Study of acceptance for DVCS events in the Forward Silicon Tracker (FST) of the SVT. Similar study made of the Barrel Silicon Tracker (BST).

CLAS12 Software Workshop - May 25, 20 p. 12/2 Event rates in the SVT. CLAS12 Luminosity: L = 35 electrons/s cm 2 and CLAS12 LH2 target: l T = 5 cm, ρ = 0.0708 gm/cm 3 Rate = 4.7 11 electrons/s SVT time window: 132 ns 62,500 electrons/event BST simulation in gemc

Event rates in the SVT. CLAS12 Software Workshop - May 25, 20 p. 13/2

CLAS12 Software Workshop - May 25, 20 p. 14/2 Event rates in the SVT. Effect of cut on deposited energy in layer 1a of the BST.

CLAS12 Software Workshop - May 25, 20 p. 15/2 Event rates in the SVT. Fluence and Dose results for Layer 1a of FST. GeV/(s cm 2 ) kilorad/year Solenoid 48.00 363.8 Off 3.56 27.042 On

Forward Tagger in CLAS12. CLAS12 Software Workshop - May 25, 20 p. 16/2

Forward Tagger in CLAS12. CLAS12 Software Workshop - May 25, 20 p. 16/2

CLAS12 Software Workshop - May 25, 20 p. 17/2 Conclusions Preliminary simulation of TOF NDE measurement consistent with CLAS6 results. EC simulation consistent with CLAS6 measurements. Simulation of DVCS events and protons and mesons in the SVT validate design. Background studies show that solenoid field will suppress background events and extend SVT lifespan. Studies of forward tagger are encouraging and a full proposal will be forthcoming. To do: PCal (Mike Wood) Cerenkov counters. Full simulation of tagged neutrons for NDE (ep e π + n).

Additional Slides CLAS12 Software Workshop - May 25, 20 p. 18/2

CLAS12 Software Workshop - May 25, 20 p. 19/2 Electromagnetic calorimeter (EC) performance 1. Sampling Fraction (f s ) (a) Ratio of energy deposited in the EC to the energy of the electron (E/p). (b) Test conditions for gemc: i. No B field (new option for gemc). ii. Fix electron momentum p e, θ e = 25, φ e = 0. iii. Convert EVIO output to Root Trees. iv. Histogram the deposited energy. v. Scan p e (0.5 GeV). 2. Integrate the distribution of E dep and divide by p e to get the sampling fraction. P ni E i dep f s = where n E e i is the number of counts in the top panel. 3. Drop observed at low p e due to threshold on E dep. Counts 3 2 1 0 0 200 300 400 500 600 700 800 E dep (MeV) 20-03-15 22:36:03 EC Sampling Fraction 0.29 0.28 0.27 0.26 0.25 0.24 0.23 0.22 0.21 Red - E = GeV Blue - E = 0.5 GeV N = 500 events EC Sampling Fraction in gemc Red - No threshold > 1 MeV Blue - E dep Black curve - 0.274(E - 21 MeV)/p 500 events 2000 4000 6000 8000 000 p (MeV) 20-03-15 22:45:28 e

CLAS12 Software Workshop - May 25, 20 p. 20/2 Electromagnetic calorimeter (EC) performance 1. Modified the identifiers in the geometry database to reflect stacking (inner/outer) and views (U, V, W) of the EC. 2. As the track is integrated determine the identifiers of the hit at each step: Sector: already known. Stack: layers 1-15 inner; layer 16-39 outer. View: use $layer % 3 to determine view from layer. Strip: 36 strips in each layer. P = (x, y) (the hit) A = (x A, y A ) = (0, pdy1) B = (x B, y B ) = (pdx2, pdy1) C = (x C, y C ) = ( pdx2, pdy1) r 1 = P C ˆr 2 = (y B y A,x A x B ) (yb y A ) 2 +(x A x B ) 2 cos θ = r 1 ˆr 2 r 1 r 2 = r 1 cos θ strip = r 2 36 + 1 B r 2 V View from Target P r 1 A O r 2 θ x strip y C

CLAS12 Software Workshop - May 25, 20 p. 21/2 Electromagnetic calorimeter (EC) performance 3. ADC signal. (a) Calculate expected number of photons from deposited energy (3.5 photons/mev). (b) Use Poisson statistics to simulate the number of photons. (c) Apply photon attenuation (λ 0 = 3760 mm). P = (x, y) A = (x A, y A ) = (0, pdy1) B = (x B, y B ) = (pdx2, pdy1) C = (x C, y C ) = ( pdx2, pdy1) Equation of BC: y = y BC = pdy1. Equation of r 1 : r 1 = r 0 + lˆr V where ˆr V = B A B A. Set y BC = y 1 and solve for l. (d) Apply PMT gain ( channels/mev) and smear with PMT resolution (use GSIM parameters). B V View from Target r 1 r r V 0 O A r V P x C strip y 4. TDC signal - convert time of hit to TDC signal using 20 ns/channel.

CLAS12 Software Workshop - May 25, 20 p. 22/2 Adding the EC to the CLAS12 Simulation 1. The CLAS6 EC geometry - defined by R. Minehart in http://www.jlab.org/ gilfoyle/clas12software/clas6ecgeometry.pdf 2. Active region defined by the following equations. The vertex of the triangle at small scattering angle and the top of the triangle at large scattering angle are y (L) = A + B (L 1) (vertex) y + (L) = A + + B + (L 1) (top) where L is the layer number, A = 1829.74 mm, B = 4.3708 mm, A + = 1899.56 mm, and B + = 4.5419 mm. The other two lines defining the active region are the following where tan θ = 1.95325. y + x tan θ = y (L ) y xtan θ = y (L) 3. EC still at θ = 25 to a perpendicular to the beamline and the same distance from the beam. 4. EC farther downstream from the target center than CLAS6. Normal distance from target center to upstream face is 7217.23 mm (JJ Roberts and D. Kashy). See http://clasweb.jlab.org/wiki/index.php/clas12 EC#Geometry

CLAS12 Software Workshop - May 25, 20 p. 23/2 Adding the EC to the CLAS12 Simulation 5. To add a new volume (the EC) to gemc: (a) Generate the parameters for a generic trapezoid (G4Trap) in Geant4. (b) Insert into a perl script (see gemc/production/database io/ clas12/geo/ec/ec build.pl). (c) Execute go tables user geometry which reads the geometry values in ec build.pl and inserts the values in the user geometry database. (d) Use the -DATABASE=user geometry option when running gemc. (e) BEWARE! Spaces versus tabs can muck up the mysql database. 6. Treating each active layer in the EC as a single trapezoid. The strips will be defined during digitization.

CLAS12 Software Workshop - May 25, 20 p. 24/2 Adding the EC to the CLAS12 Simulation 7. Geant4 G4trap geometry parameters: pdx1 Half x length of the side at y=-pdy1 of the face at -pdz pdx2 Half x length of the side at y=+pdy1 of the face at -pdz pdz Half z length ptheta Polar angle of the line joining the centres of the faces at -/+pdz pphi pdy2 pdx4 palp2 palp1 pdy1 Azimuthal angle of the line joining the centre of the face at -pdz to the centre of the face at +pdz Half y length at +pdz Half x length of the side at y=+pdy2 of the face at +pdz Angle relative to the y axis from the centre of the side (upper endcap) Angle with respect to the y axis from the centre of the side (lower endcap) Half y length at -pdz pdx3 Half x length of the side at y=-pdy2 of the face at +pdz

CLAS12 Software Workshop - May 25, 20 p. 25/2 Adding the EC to the CLAS12 Simulation 8. Add the EC banks to the simulations. (a) In gemc/production/database io/clas12/banks/ create ec.bank (see below) and edit clas12 hits def.txt. (b) Execute go tables and go hits in the same area to make entries in clas12 banks database. (c) In gemc/production/hitprocess/ create EC hitprocess.cc and EC hitprocess.h and edit Clas12 HitProcess MapRegister.cc to include the EC. (d) Recompile. (e) Talk to Mauri, fix your mistakes, and recompile...

CLAS12 Software Workshop - May 25, 20 p. 26/2 Adding the EC to the CLAS12 Simulation 12. Add the code to gemc to process an EC event; simplified version now. 13. Modify and add functions to gemc evio2root to include EC events in gemc evio2root/src/: (a) Create ec hit.cc and ec hit.h to handle Geant4 hit information. (b) Edit root tree.cc, root tree.h, and gemc evio2root.cc to include EC information. 14. First histograms from gemc gemc_evio2root for electrons. Counts 2200 2000 1800 1600 1400 1200 00 800 600 400 200 2009-11-15 20:34:09 CLAS12 7.5-GeV electrons Fixed Angle, gemc results 0 0 20 40 60 80 0 120 140 160 180 200 220 E dep in the EC (MeV) Defining an electron: t < 50 ns, N EC < 40, E dep > 5 MeV, PID = 11 Local x (mm) 400 300 200 0 0-0 -200-300 -400 CLAS12 7.5 GeV electrons 600 800 00 1200 1400 2009-11-15 20:48:14 Local y (mm) 4 3 2 1

CLAS12 Software Workshop - May 25, 20 p. 27/2 Adding the EC to the CLAS12 Simulation 15. First histograms from gemc gemc_evio2root for neutrons. E dep (MeV) hlxly Entries 4579 3 2 gemc simulation 4.3 GeV neutrons o = 30 θ n Local y (mm) 400 300 200 0 gemc simulation 4.3 GeV neutrons o = 30 θ n 2 0-0 -200 1 1 0 20 40 60 80 0 120 140 160 180 200 (MeV) 2009-11-17 22:39:02 E dep -300-400 1400 1500 1600 1700 1800 1900 2000 20 2200 Local x (mm) 2009-11-17 22:58:06 Defining a neutron: t < 50 ns, N EC < 40, E dep > 5 MeV, PID = 2112

CLAS12 Software Workshop - May 25, 20 p. 28/2 Hadron properties in the Silicon Vertex Tracker (SVT) Generated Events: 0,000 protons/kaons/pions momentum: 0.1 to 2.1 GeV θ: 0-160 uniform Longitudinal Vertex: 5 cm uniform Transverse Vertex: 1 mm uniform φ: 0-360 uniform Deposited energy (E dep ) in MeV versus Particle Momentum (p) in GeV. Protons K + π Generated θ versus primary z vertex.

2024 © sciencedocbox.com Privacy Policy | Terms of Service | Feedback