Particle identification by energy loss measurement in the na61 shine experiment
Download
1 / 44

Particle identification by energy loss measurement in the NA61 (SHINE) experiment - PowerPoint PPT Presentation


  • 77 Views
  • Uploaded on

Particle identification by energy loss measurement in the NA61 (SHINE) experiment. Magdalena Posiadala University of Warsaw. Outline:. Energy loss measurements in NA61 experiment at CERN Energy loss parametrisation

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about ' Particle identification by energy loss measurement in the NA61 (SHINE) experiment' - sydney


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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Particle identification by energy loss measurement in the na61 shine experiment

Particle identification by energy loss measurement in the NA61 (SHINE)experiment

Magdalena Posiadala

University of Warsaw


Outline
Outline:

  • Energy loss measurements in NA61 experiment at CERN

  • Energy loss parametrisation

  • Methods of Particle Identification in NA61 in the region of interest for the T2K.

  • Energy loss measurement in the region below 1GeV/c

    • Two methods of approach to analysis:

    • Method 1 – fits done using TMinuit package

    • Method 2- using probability functions

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Na61 shine detector
NA61 (SHINE) detector

Max 72 points

Max 90 points

TPCs are filled with mixtures of Ar+ CO2(90:10) for VTPCs and (95:5) for MTPCs

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Energy loss measurement
Energy loss measurement

  • The optimization of the parameters required for the determination of energy loss in the TPCs was performed using the method developed by NA49. (For details see Eur. Phys. J , C45 (2006) 343)

  • Corrections for the following effects were applied during thecalibration of the NA61 data:

    • signal loss due to threshold cuts; corrections were obtained using Monte Carlo calculations for gas mixture used in NA61,

    • time dependence of the TPC gas pressure,

    • residual time dependence of the measurements (day/night),

    • charge absorption during the drift,

    • differences in the TPC sector gain factors,

    • differences in the amplification of the preamplifiers and edge effects at sectors boundaries.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Energy loss measurement1
Energy loss measurement

Preliminary distributions of corrected dE/dx values forpositively (white) and negatively (yellow) charged particles with momentum chosen around 5 and 15 GeV/c.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Energy loss parametrisation
Energy loss parametrisation

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Energy loss parametrisation1
Energy loss parametrisation

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


De dx vs number of measured points on the track
dE/dx vs number of measured points on the track

1 island-> vtpc

2 island –> mtpc+vtpc

Estimated resolution of dE/dx measurements of 4-5% can be achieved for particles passing through both the vertex and the main TPC chambers.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Total momentum
Total momentum

N>10

N>30

N>40

N>50

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


De dx for total momentum p 2 gev c
dE/dx for total momentum p< 2 GeV/c

tracks with high dE/dx errors are rejected

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Methods of particle identification in na61 in the region of interest for the t2k
Methods of Particle Identification in NA61 in the region of interest for the T2K

  • In NA61 (SHINE) experiment particle identification is possible using ionization measurements in active volume of the TPCs and time-of-flight measurements in ToF counters.

  • For ToF-F measurements see next presentation of Sebastien Murphy.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Comparison of interest for the T2K

NA61 data with T2K simulation

  • JNUBEAM

JNUBEAM –simulation of the T2K beam

NA61, 2C data, N>30

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Tracks with and without tof f n 30
Tracks with and without TOF-F N>30 interest for the T2K

dE/dx info

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Energy loss measurement in the region below 1gev c

Energy loss measurement interest for the T2Kin the region below 1GeV/c


Statistics below 1gev c na61 data
Statistics below 1GeV/c interest for the T2K – NA61 Data

in %

In T2K bins

ptot: 200MeV/bin x 50 bins till10GeV/c

:

20 rad/bin x 20bins

till 0.4 rad

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Statistics for n 30 p 1gev c
Statistics for N>30, p<1GeV/c interest for the T2K

p (GeV/c)

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Statistics below 1GeV/c interest for the T2K -2C Data

in %

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


De dx region below 1gev c
dE/dx region below 1GeV/c interest for the T2K

  • Preliminary plot ofthe reconstructed dE/dx values versus momenta (p<1.2 GeV) for positively charged particles, together with the Bethe-Bloch curves for positrons,pions, kaons, protons and deuterons.

  •  of dE/dx vs dE/dx

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Smearing with gaussian theoretical distribution s of energy loss
Smearing with Gaussian theoretical distribution interest for the T2Ks of energy loss

  • For whole tracks we can calculate theoretical energy loss dE/dxteor using track total momentum and applied parametrisation functions.

  • In this way we get theoretical distribution of energy loss in momentum bins (see example plot on the right).

  • Smearing with Gaussian function

Example ->

dE/dx data

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Smearing with gaussian theoretical distribution s of energy loss 2
Smearing with Gaussian theoretical distribution interest for the T2Ks of energy loss (2)

  • The theoretical distributions of energy loss are smeared with Gaussian function following this rules:

  • Mean = dE/dxtheor calculated for each track

  • Sigma = dE/dx for the each track; from NA61 data

After Gaussian smearing

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Methods of analysis
Methods of analysis interest for the T2K

  • Having theoretical energy loss distributions smeared for each of the particle hypotheses I start to extract raw particle yields using 2 methods:

  • Method 1-> Making fits using TMinuit Package in Root

  • Method 2 ->Calculating probability functions for 4 particles hypotheses

  • For more details look at the example in the next few slides.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 16 0 3 gev c
Momentum range (0.16-0.3) GeV/c interest for the T2K

  • Preliminary plot of the reconstructed dE/dx values versus momentum (p<1.2 GeV) for positively charged particles, together with the Bethe-Bloch curves for positrons, pions, and only kaons in selectedmomentum range.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 16 0 3 gev c1
Momentum range (0.16-0.3) GeV/c interest for the T2K

  • We have theoretical smeared with Gaussian energy loss distributions (right) and NA61 data (left)

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 16 0 3 gev c2
Momentum range (0.16-0.3) GeV/c interest for the T2K

  • Method 1 -> uses 2 function for minimization

  • Theoretical distribution of energy loss is fitted to data

  • Calculation of the fit parameters ,,,

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 16 0 3 gev c3
Momentum range (0.16-0.3) GeV/c interest for the T2K

  • Method 1

  • TMinuit fits

  • Results

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 16 0 3 gev c4
Momentum range (0.16-0.3) GeV/c interest for the T2K

  • Method 2->Probability functions; for each of the bin i assuming particle hypotheses I calculate probability function defined as:

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 16 0 3 gev c5
Momentum range (0.16-0.3) GeV/c interest for the T2K

  • Method 2-> Multiply probability functions by distribution of the energy loss for the data in selected momentum slice to obtain colourful plots at the bottom picture;corresponding to particle species.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 3 0 4 gev c
Momentum range (0.3-0.4) GeV/c interest for the T2K

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 3 0 4 gev c1
Momentum range (0.3-0.4) GeV/c interest for the T2K

  • Method 1

  • TMinuit fits

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 3 0 4 gev c2
Momentum range (0.3-0.4) GeV/c interest for the T2K

  • Method 2

  • Probability functions

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 5 0 6 gev c
Momentum range (0.5-0.6) GeV/c interest for the T2K

protons are seen

(yellow line)

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 5 0 6 gev c1
Momentum range (0.5-0.6) GeV/c interest for the T2K

  • Method 1

  • TMinuit fits

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 5 0 6 gev c2
Momentum range (0.5-0.6) GeV/c interest for the T2K

  • Method 2

  • Probability functions

Almost the same number of positrons and kaons, which is

of course not true.

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 6 0 7 gev c
Momentum range (0.6-0.7) GeV/c interest for the T2K

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 6 0 7 gev c1
Momentum range (0.6-0.7) GeV/c interest for the T2K

  • Method 1

  • TMinuit fits

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 7 0 8 gev c
Momentum range (0.7-0.8) GeV/c interest for the T2K

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 7 0 8 gev c1
Momentum range (0.7-0.8) GeV/c interest for the T2K

  • Method 1

  • TMinuit fits

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 8 0 9 gev c
Momentum range (0.8-0.9) GeV/c interest for the T2K

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 8 0 9 gev c1
Momentum range (0.8-0.9) GeV/c interest for the T2K

  • Method 1

  • TMinuit fits

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 9 1 1 gev c
Momentum range (0.9-1.1) GeV/c interest for the T2K

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 9 1 1 gev c1
Momentum range (0.9-1.1) GeV/c interest for the T2K

  • Method 1

  • TMinuit fits

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Momentum range 0 9 1 1 gev c2
Momentum range (0.9-1.1) GeV/c interest for the T2K

  • Method 2

  • Probability functions

Method 2 at this stage of studies works well only for momentum slices where theoretical lines do not cross completely !!!

Needs further improvements (/K ratio parameterised by function)

"Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


Summary
Summary interest for the T2K

  • The resolution dE/dx strongly depends on number of points per track. For tracks passing through all TPCs the resolution at the level of 4 - 5% .

  • NA61 data cover acceptance in p,  variables for the T2K experiment.

    • Important dE/dx analysis below 1 GeV/c region.

  • Presented two methods towards raw particle yields

    • Method 1 –>making fits using TMinuit package

    • Method 2 with probability functions still needs improvement->approaching towards max log-likelihood.

  • "Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


    Further steps
    Further steps interest for the T2K

    • On the right relative error vs energy loss distribution is presented.

    • Next step is to smear theoretical lines with Gaussians that have error equal to relative error of dedx and see how it influences the presented results

    "Neutrino physics at accelerators", JINR, Dubna 27-30.01.2009


    ad