Beam mc activity
This presentation is the property of its rightful owner.
Sponsored Links
1 / 11

Beam MC activity PowerPoint PPT Presentation


  • 86 Views
  • Uploaded on
  • Presentation posted in: General

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

Download Presentation

Beam MC activity

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.


Comparison among different versions

K. Matsuoka

Comparison among different versions


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

  • [email protected] decreased by 3%, profile shape is same

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?


  • Login