beam mc activity
Download
Skip this Video
Download Presentation
Beam MC activity

Loading in 2 Seconds...

play fullscreen
1 / 11

Beam MC activity - PowerPoint PPT Presentation


  • 109 Views
  • Uploaded on

Beam MC activity. A.K.Ichikawa for beam group. For more details, http ://jnusrv01.kek.jp/~ichikawa/jhfnu/nubeam/beamMCsummary.091203.ichikawa.pdf. Summary of update in 09a,b,c. 09a

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 ' Beam MC activity' - rhett


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
beam mc activity

Beam MC activity

A.K.Ichikawa for beam group

For more details,

http://jnusrv01.kek.jp/~ichikawa/jhfnu/nubeam/beamMCsummary.091203.ichikawa.pdf

summary of update in 09a b c
Summary of update in 09a,b,c
  • 09a
    • This is almost equivalent to 07a except some bugs were fixed. Those bugs would not influence the neutrino flux iteslef.
  • 09b
    • Mainly, 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 can be specified via card.
    • Neutrino tracking is discarded.
    • The density of the dump concrete was changed from 2.2 to 2.3g/cm3.
mumon

K. Matsuoka

MUMON

A few bugs introduced in 2009 were found and fixed. Mainly in horn1.

  • Magnetic field (proportional to 1/r) starts from r=2.7cm
    • Horn1 inner conductor
      • Inner radius = 2.7cm, outer radius=3cm
    • Originally it starts from 3cm
    • Actually, something between them (skin depth > 5mm)
    • For MUMON, a few% effect
muon flux at muon pit spring commissioning 1 st horn only emulsion v s mc

A.Arigaet. al.

Muon flux at muon pitSpring commissioning (1st horn only)Emulsion v.s. MC

0kA Data/MC=0.57

220kA Data/MC=0.75

273kA Data/MC=0.78

horn2 and horn3 geometry update

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 ->195.7 cm
  • Horn3
    • outer conductor radius : 70 cm  65.5 cm
    • B-field region (Z-length) : 250 cm  247.1 cm

http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=147

horn2 3 geometry update neutrino flux
Horn2&3 geometry updateneutrino flux

09b

new

less than 5% difference upto 5GeV

[email protected]

Far/Near

parents of muons in muon pit pions
Parents of muons in muon pit-pions-

Horn off

1st Horn 273kA

All horns 320kA

H.Kubo. K.Matsuoka

proton information

K. Matsuoka

Proton information
  • Store primary beam information in a common block and MUMON or neutrino ntuple.
    • libneubeam/beaminfo.inc
      • integer gipart : Primary generated particle ID
      • real gpos0(3) : Generation point
    • common /gbeam/gipart, gpos0
  • Accumulate POT w/ a flat proton beam.
  • Weight the muon profile w/ an arbitrary proton profile by using gpos0.
  • We do not need to make many MC data sets of various proton beam profile.
new nd flux calculation algorithm

H.Kubo

new ND flux calculation algorithm
  • current filling routine
    • SK : treated as a “point”.

generated neutrino is always filled

acceptance of SK for each decay is included in “norm”.

    • ND : repeat parent’s decay randomly (uniformly in CM)

by N times (default 1,000).

only neutrinos which have proper angle are filled.

  • New method : fill ND ntuple just the same method as one of SK.
    • choose a detection point (x, y) randomly in the ND plane
    • calculate weight (acceptance) for this direction

using parent 4-momentum

    • same method with SK ntuple.
  • Motivation
    • In the current version, high-energy pions are multiply used. Events are not independent. Simple error couting results in underestimate.
  • Status
    • Coding is almost done. Under check.

http://jnusrv01.kek.jp/Indico/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=162

remaining issues towards dec jan mass production
Remaining issues towards Dec.(Jan.) mass production
  • Received requests from ND280 beam group
  • Received a request from Hiro
    • To check new ND280 algorithm, want flux in both old and new algorithm
  • Need to be done
    • final(?) check of new ND algorithm
    • final(?) check of mag. field inside conductor
    • correct ND280 position
    • Optimization of proton beam area
      • two flat area?
    • check of Nicolas commitment and decide which should be stored for mass production version
    • Random seed control?
ad