1 / 33

Beam MC progresses

Beam MC progresses. for beam MC sub-group. Summary of update in 09b,c,10a. 09b Geometry of baffle, target, 1st horn, dump and MUMON is updated. 09c MUMON structures were included. Energy deposit in the MUMON detector can be stored. Emittance and Twiss parameters via card.

zona
Download Presentation

Beam MC progresses

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. Beam MC progresses for beam MC sub-group

  2. Summary of update in 09b,c,10a 09b • Geometry of baffle, target, 1st horn, dump and MUMON is updated. 09c • MUMON structures were included. • Energy deposit in the MUMON detector can be stored. • Emittance and Twiss parameters via card. 10a to be released soon • Horn2&3 geometry update • Mag. field inside inner conductor • Store primary proton vector information to enable weighting method • New ND280 flux algorithm • Store particle interaction history • K±µ3 and K0µ3 decay for neutrino • Random number generation seeds control

  3. K. Matsuoka Comparison among different versions

  4. A.K.Ichikawa T2K horn effect (jnubeam 09c) X9.4 x16 On-axis Off-axis

  5. K. Matsuoka Horn Magnetic Field How should we treat the magnetic field inside the inner conductor? Horn1 inner conductor Inner radius = 2.7cm, outer radius=3cm Assuming that the elec. current is uniform in the conductor, (skin depth > 5mm) Significant effect on MUMON signal was found Modest effect on neutrino flux Realistic. B-field B r Inner conductor

  6. mag. field inside inner conductor-MUMON, all horns on- (* per 3.3 x 1014 POT) 8~9% effect. (difference from realistic B field) • Max. B-field • 1/r-shape field starts from inner surface of the inner conductor • Min. B-field • starts from the outer surface of the inner condcutor

  7. mag. field inside inner conductor En spectra (SK) Red: Max./realistic Blue: Min./realistic All nm Red: Max. Black: Realistic Blue: Min. <3% effect (need more statistics)

  8. mag. field inside inner conductor En spectra (ND on-axis) Red: Max./realistic Blue: Min./realistic All nm Red: Max. Black: Realistic Blue: Min. (Error bars may be under-estimated.See later slides.)

  9. mag. field inside inner conductor INGRID profile ND2 All nm ND2 All nm Red: Max. Black: Realistic Blue: Min. Red: Max. Black: Realistic Blue: Min. abs(x) ≤ 5.5 m && abs(y) ≤ 0.5 m abs(x) ≤ 0.5 m && abs(y) ≤ 5.5 m • Due to the magnetic field in the inner conductors, n flux at ND on the axis gets sharper than that of 09c (min B-field). • Peak value: (min) 5.66 x 1017 (realistic) 5.79 x 1017 /m2/1021 POT

  10. H.Kubo Horn2 and Horn3 Geometry update update from conceptual shape to real shape • Horn2 • outer conductor radius : 40 cm -> 49.048 cm • B-field region (Z-length) : 200 cm ->199.7 cm • Horn3 • outer conductor radius : 70 cm  65.5 cm http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=147

  11. Horn2&3 geometry updatemuon flux @mumon (Si-plane, horn 320kA) • flux decreased by 3%, profile shape is same 3horns, 320kA @3.4x1011POT

  12. Horn2&3 geometry updateneutrino flux 09b new less than 5% difference upto 5GeV nm@SK Far/Near

  13. Parents of muons in muon pit-pions- Horn off 1st Horn 273kA All horns 320kA H.Kubo. K.Matsuoka

  14. Parents of muons in muon pit-kaons- Horn off 1st Horn 273kA All horns 320kA H.Kubo. K.Matsuoka

  15. Parents of muons in muon pit-K/pi ratio-

  16. K. Matsuoka Proton information • Store primary beam information • Accumulate POT w/ a flat proton beam and weight t w/ an arbitrary proton profile to simulate that profile • No need to make many MC data sets of various proton beam profile. Demonstration w/ MUMON profile (* m+–/3.4 x 1011 POT) sy: 1.7 mm

  17. N. Abgrall Production history • Fill ntuples with neutrino history, taking decay chains into account. • information of primary, secondary,...interactions) • Include additional decay modes for pions and kaons, updated branching ratios • π± → e± νe • K±µ3 and K0µ3 • neglect K0S semileptonic decays ? (e.g. K2K case) http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=4&resId=0&materialId=slides&confId=133

  18. H.Kubo new ND flux calculation algorithm • current filling routine • SK : treated as a “point”. for every decay of p/K/m, neutrino is forced to go towards SK probability is calculated and stored as “norm”. • ND : repeat parent’s decay randomly (uniformly in CM) by 1,000 times only neutrinos which have proper angle are filled. • New method : same method as SK case. • choose a detection point (x, y) randomly in the ND plane • calculate weight (acceptance) for this direction • Motivation • In the current version, high-energy pions are multiply used. Events are not independent. Simple error couting results in underestimate. http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=162

  19. Enu • spectrum seems to be consistent • χ2 = 19.7 / 39 • to small • due to using same set of parents ?

  20. error histogram (Enu) • low energy (< 1GeV) : same or smaller error • Original method had been giving underestimated error original new

  21. on-axis xnu (fitting) • fit with Gaussian • large chi2 & mean offset (10 sigma) in original algorism indicates under estimation of error original new

  22. Other activity • Detailed check of dimensions by P.Perio • Treatment of Random numbers • M.Hartz, K.Sakashita • code is modified to select 215 good seed-pairs for GRNDM by K.Sakashita http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=2&resId=0&materialId=slides&confId=155 • OTR simulation by OTR group • Target scan simulation • K.Matsuoka, M.Hartz • CPU saving effort • Review on gcalor (secondary interaction model) • Review on INGRID study A.Minamino

  23. Prospect Flux Mass production • Received requests from ND280 beam group • Need to be done • Implementation of the correct ND280 position • Optimization of proton beam area • two flat area? • Release 10a • In Next week at 250kA horn current Remaining update • Striplines • Transfer matrix with new ND280 algorithm • Inclusion of NA61 results w/ NA61-T2K group • And studies.

  24. Other geometry update • MUMON structure has been added. • Geometry of the collimator at the entrance of DV has been update based on the measurement. • The size of the DV entrance has been changed based on the measurement. • Density of dump material • concrete from 2.2 to 2.3 g/cm3 • ~1% effect on MUMON • concrete rebars 2.3g/cm3 -> 2.377 g/cm3 • <1% effect on MUMON

  25. mag. field inside inner conductor-MUMON,1st horn only- Primary proton profile sx= 0.36, sy= 0.17 (mm) Only Horn1 same as April ’09 commissioning (*per 3.4x1011 POT, 2D fit peak) A few % effect http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=152

  26. Horn2&3 geometry updatepion production point (mumon) horn1 09b new horn3 horn2 horn1 horn2 dump horn3 103entries / 5.0 x 107 POT

  27. Horn2&3 geometry updateeffect of horn2&3 material absorption ~ 10% (14%@ peak) less effect than horn1

  28. Horn2&3 geometry update horn 320kA pions whose daughter muons goes through mumon Si-plane particle flux [104 / cm2 / 3.4 x 1011POT] @MUMON Si-plane ( peak of X projection-fit )

  29. A.Ariga et. al. Muon flux at muon pitEmulsion v.s. MC

  30. Comaprison Emulsion, Si, MC Proton : from beam summary (result_run24.root) Emulsion : cutoff 0.05GeV/c, q<0.3rad Si : using only 1 line (7 sensors) which corresponds emulsion modules. MC : muon, position at emulsion, momentum>0.05GeV/c, q<0.3rad normalized by Proton(CT05). Normalized at POT=3.4e11

  31. M.Hartz Target scan Data for y scan

  32. H.Kubo Review on gcalor Consists of : • NMTC : nucleons < 3.5GeV, π±< 2.5GeV • SCALE : Scaling Model (3GeV to 10GeV) • MICAP : neutron < 20MeV • (FLUKA) : >10GeV & other particles NMTC & Scaling • For nucleons below 3.5GeV and π±below 2.5GeV, NMTC is used. • Above 10GeV, FLUKA is used. • Scaling energy range (3-10GeV) • FLUKA or (scaled)NMTC is called for each interaction with a linear probability function for smooth transition For more details, http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=2&resId=0&materialId=slides&confId=133

More Related