1 / 33

Particle Identification with TRD Prashant Shukla Institute of Physics University of Heidelberg

Particle Identification with TRD Prashant Shukla Institute of Physics University of Heidelberg. Presentation at International Workshop on TRD 27 th September 2005 at Cheile Gradistei, Romania. Summary of 1. Implementation of TRD particle ID in reconstruction.

tatum
Download Presentation

Particle Identification with TRD Prashant Shukla Institute of Physics University of Heidelberg

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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Particle Identification with TRD Prashant Shukla Institute of Physics University of Heidelberg Presentation at International Workshop on TRD 27th September 2005 at Cheile Gradistei, Romania.

  2. Summary of 1. Implementation of TRD particle ID in reconstruction. 2. Identification of protons and kaons. 3. Particle ID from test beam 2004. 4. New methods for improvements in PID. 5. Preparation for analysis of TRD events in Physics Data Challenge Outline

  3. Implementaion of TRD PID: 1. Each particle track is assigned with dE/dx signals and position of maximum charge from all 6 planes of TRD. 2. Class AliTRDprobdist contains Detector response or Distributions of dE/dx and MaxTimeBin for electron, pions, muons, kaons and protons between 0.5 GeV to 6 GeV momentum. 3. Probability for a track to be electron or pion or any other particle is assigned to each ESDtrack.

  4. Modification Scheme of ESD and TRD track AliRootV4-02-09 and v4-02-Release

  5. The dE/dx distribution and the Time Bin of max. Cluster distribution AliRoot V4-02-Release (Modified) P =2 GeV electrons(red) and pions (blue)

  6. Particle probabilities with ESDtrack The combined probability P over all layers is obtained as With ESD track the probabilities are as follows: P[0] for e = Pe/(Pe+Ppi+Pmu+Pk+Pp) p[1] for mu p[2] for pi p[3] for K p[4] for p

  7. Liklihood and efficiency AliRoot V4-02-Release electrons(red) and pions (blue) L- QX P=2 GeV

  8. V4-02-Release (Modified)Pion efficiency vs Momentum

  9. dE/dx dist for Pions,Electrons,kaonsandprotons p=560 MeV/c Minimum Ionizing Pions (p=560 MeV/c) Truncated Mean (Cluster aplitudes between 5 % and 70%) over all six layers

  10. PID for pions, kaons and protons generated with equal probabilities Submitted for PPR Volume II Chapter 5

  11. V4-02-Release (Modified)TRD signals (p=2 GeV) with beam data 2002 Reconstructed Vs. Data With Tail Cancellation Digits Vs. Data

  12. MACROS for test and generation of Detector response Generating response functions: 1. Run job.sh: Runs AliRoot (sem.C+rec.C) 50 times for Config.C, for 200 pions and 200 electrons. Repeat it all momenta. 2. Run Set5.C to generate files for response functions for electrons and pions. Repeat it at all momenta. 3. Run handata.C to create the function AliTRDprobdist::FillData(). 4. Cut and paste it to the AliTRDprobdist.cxx file at the end.

  13. MACROS for test and generation of Detector response Comparing Reconstructed and beam data: 1. Run job.sh: Runs AliRoot (sem.C+rec.C) 50 times for Config.C, for 200 pions and 200 electrons. 2. Run Set5.C to get energy loss distribution histograms for electrons and pions. 3. Run rplot.C to plot Digits and the experimental Data. Comparing Digits and beam data: 1. Run digits.sh: Runs Aliroot (sem.C+rec.C) 5 times each for Configpi.C and Config.C. 2. Run TRDDigits00.C to get energy loss distribution histograms for electrons and pions. 3. Run plotdigits1.C to plot Digits and the experimental Data. 4. Run plotdigits2.C to plot Digits and the Reconstructed data.

  14. HEAD (25 July)TRD signals (p=2 GeV) with beam data 2002 Reconstructed signals TRD Digits

  15. HEAD (25 July) Number of Clusters per track

  16. Test Beam: Old (2002) and New Data (2004) (4, 6 GeV)

  17. Beam test 2004 : Time Bin of Maximum Charge (p=6 GeV)

  18. Pulse Shape Descrimination Test Beam Data 2004 (p=6 GeV)

  19. Pulse Shape Test Beam Data 2004 (p=6 GeV) Electrons

  20. Pulse Shape Test Beam Data 2004 (p=6 GeV) pions

  21. Pulse slope distribution Test Beam Data 2004 (p=6 GeV)

  22. Pulse shape discriminationTest Beam Data 2004 (p=6 GeV)

  23. Pulse shape discrimination Test Beam Data 2004 (p=6 GeV)

  24. Principal Component Analysis Test Beam Data 2004 (p=6 GeV)

  25. Physics Data Challenge : TRD Signal events 10 ( J/’) ---->e- + e+ 10 ( Y, Y’, Y”) ----> e- + e+ 10 ( Ds and Dbars) ----> e- + e+ 10 ( Bs and Bbars) ----> e- + e+ 10 ( Bs and Bbars) -->J/ --> e- + e+ 10 f ---->e- + e+ 10 Omega ---->e- + e+ These are the parameterizations of Ramona Vogt results Given in AliGenMUONlib Underlying Events: Hijingcent1 -->b (0, 5) Hijingper1 --> b (5, 8.6)

  26. Physics Data Challenge : TRD Signal events 1. Total 5000 X 20 = 100000 events will be produced 2. With the aim to Study the efficiencies and resolutions for physics signals after full reconstruction. Aliroot v4-02-Release 1. Ramona Vogt parameterization for J/psi does not work any more ---> Use CDF scaled. 2. J/psi, Upsilon Pt values does not differ from event to event 3. There is new tracking code (impoved) after this version

  27. dir/nucleon pair (b) PbPb 5.5 GeV Calculated by NLO J/’ dir/nucleon pair (b) 11.7 2.65 Y/evt (5% Centrality) 0.31 0.07 e pairs/evt 0.02 0.005*10^(-3) Y Y’ Y” dir/nucleon pair (b) 0.15 0.094 0.057 Y/evt (5% Centrality) 0.004 0.0025 0.0015 e pairs/evt 0.0001 0.000034 0.00003 Cross sections and yields of Quarkonia: L = 1027 cm-2 s-1 = 1 mb-1 s-1

  28. D+ D0 Ds+ Lc+ dir/nucleon pair (b) Total= 45 barn Y/evt (5% Centrality) 22.3 70.4 13.4 9.0 e pairs/evt (10.18) 3.8 4.9 1.07 0.405 Yields for open charm mesons Pb–Pb collisions at √sNN = 5.5 TeV centrality selection of 5% These values correspond to the average of the result obtained with MRST HO and CTEQ 5M1 parton distribution functions. EKS98 parameterization of nuclear shadowing.

  29. B0 B+ Bs0 Lb0 dir/nucleon pair (b) Total=1.79 barn Y/evt (5% Centrality) 1.83 1.82 0.53 0.34 e pairs/evt (0.47) 0.19 0.19 0.056 0.034 Yields for open beauty mesons Pb–Pb collisions at √sNN = 5.5 TeV centrality selection of 5% B/evt 4.6 B->J/psi / evt 0.046 e pairs/evt 0.0030 These values correspond to the average of the result obtained with MRST HO and CTEQ 5M1 parton distribution functions. EKS98 parameterization of nuclear shadowing.

  30. CDF scaled J/ and Upsilon Pt Distribution: Upsilon Pt Distribution: J/ Pt Distribution:

  31. Invariant Mass of electron pairs: HEAD 25 July Ccbar Bbbar

  32. Invariant Mass of electron pairs: HEAD July 25 pt Cut = 1 GeV No pt Cut Ccbar Bbbar

  33. 1. Implementation of TRD particle ID in reconstruction. 2. Study of TRD Identification of other particles protons and kaons submitted to PPR Volume II 3. Particle ID from test beam 2004. 4. Pulse Shape method gives upto 4 times improvement in pion efficiency. Next: Implement it in Aliroot TRD Reconstruction. 5. TRD events in Physics Data Challenge: v4-02-Release Next: Tag a new version with improved tracking Summary and Future Outlook

More Related