chorus physics workshop, 7 june 2004, rome luca scotto lavina kinematical selection and fc analysis...

CHORUS Physics Workshop, 7 June 2004, Rome

Luca Scotto Lavina

Kinematical selection and FC analysis for oscillations

Montecarlo chain used for 0mu channel

Events generated in emulsion

Loca

tion

pro

gra

m

Selected eventsSato

Mu

rat

an

d

Gold

en

sele

ctio

n

Stopped in bulk

• PRDC recontruction,• CS, SS and bulk scanback• [4,35] fiducial volume cut• 1.0 cm fiducial volume border

cut

No spec muon in PRTK

First 0mu definitionFor each track in PRVTPRTK:1) read MUON-ID = q(JSMTR+3) =

= q(JMUTR+3)2) If MUON-ID>0, the track is a

spectrometer muonIf there is no muon, the event is

0mu

Production of .man files and charm.kin file

1.5 cm fiducial volume border cut

Events in 0mu.listSecond 0mu definition

Based on MURECO flg_mu<4, IP cut respect to vertex, 2 cut

This definition includes the first one

SatoMurat: SatoMurat_mc-ver04.cpp Selection: goldenMC-0mu-04-02-02.awk

EMUL 7131

In bulk 2777

Spec 0mu 209

.man files 201

0mu.list 117

Selected 23

7131 events generated in emulsion with

Location efficiency + 0mu request =0mu.lis

t

EMUL= (1.65 ± 0.15) %

Selection efficiency =Selected

0mu.list= (20.2 ± 3.7) %

D+, Ds, c 1 prong

MC charm 0mu, 1 prong

0mu request = 0mu.lis

t

In bulk

· 1.043 = 4.4%

Spec 0mu

In bulk

= 7.5%

If I consider only the first 0mu definition, based only on spectrometer muons:

1.043 = 1.0 cm f.v. border

cut

1.5 cm f.v. border cut

(Full statistics)

EMUL 9078

In bulk 3585

Spec 0mu 410

.man files 384

0mu.list 197

Selected 100



t

EMUL= (2.18 ± 0.15) %


0mu.list= (50.8 ± 3.5) %

D+, Ds, c 3 prong

MC charm 0mu, 3 prong


t

In bulk

· 1.043 = 5.7%

Spec 0mu

In bulk

= 11.4%


1.043 = 1.0 cm f.v. border

cut


(Full statistics)

EMUL 6831

In bulk 1479

Spec 0mu 1393

.man files 1322

0mu.list 875

Selected 209



t

EMUL= (12.8 ± 0.4) %


0mu.list= (23.9 ± 1.4) %

e,h

MC e,h


t

In bulk

· 1.043 = 61.7%

Spec 0mu

In bulk

= 94.2%


1.043 = 1.0 cm f.v. border

cut


EMUL 1552

In bulk 418

Spec 0mu 394

.man files 367

0mu.list 272

Selected 143



t

EMUL= (17.6 ± 1.0) %


0mu.list= (52.6 ± 3.0) %

3h

MC 3h


t

In bulk

· 1.043 = 67.9%

Spec 0mu

In bulk

= 94.3%


1.043 = 1.0 cm f.v. border

cut


EMUL 1765

In bulk 636

Spec 0mu 107

.man files 99

0mu.list 74

Selected 18



t

EMUL= (4.2 ± 0.5) %


0mu.list= (25.0 ± 4.9) %

MC


t

In bulk

· 1.043 = 12.1%

Spec 0mu

In bulk

= 16.8%


1.043 = 1.0 cm f.v. border

cut


EMUL 2233

In bulk 938

Spec 1mu 845

.man files 840

CC.list 771


Location efficiency + 1mu request =CC.list

EMUL= (34.5 ± 1.0) %

CC

MC CC

1mu request = CC.list

In bulk· 1.043 = 85.7%

Spec 1mu

In bulk

= 90.1%


1.043 = 1.0 cm f.v. border

cut


Summary of efficiencies

Charm 0mu 1 prong

Charm 0mu 3 prong

e,h 3h

Location eff.+ 0mu request

1.65 ± 0.15

2.18 ± 0.15

12.8 ± 0.4 17.6 ± 1.0 4.2 ± 0.5

Selection efficiency

20.2 ± 3.7 50.8 ± 3.5 23.9 ± 1.4 52.6 ± 3.0 25.0 ± 4.9

numuCC

Location eff.+ 1mu request

34.5 ± 1.0

Nmax(e,h) expected events

Nmax(e,h) = N0 loc · r · rA · BR(e,h) · sel

N0 loc = 23678

r =

CC

NC

CC

CC

CC

NC

= · = 0.53/0.31 = 1.71

rA = A

0

ANC 0

= 1.06 ± 0.07

BR(e,h) = 0.673

Nmax(e,h) = 6903sel = (23.9 ± 1.4)%

Nmax(3h) expected events

Nmax(3h) = N0 loc · r · rA · BR(3h) · sel

N0 loc = 23678

r =

CC

NC

CC

CC

CC

NC

= · = 0.53/0.31 = 1.71

rA = A

0

ANC 0

= 1.20 ± 0.12

BR(3h) = 0.152

Nmax(3h) = 3884sel = (52.6 ± 3.0)%

Nmax() expected events

Nmax() = N0 loc · r · rA · BR() · sel

N0 loc = 23678

r =

CC

NC

CC

CC

CC

NC

= · = 0.53/0.31 = 1.71

rA = A

0

ANC 0

= 0.35 ± 0.09

BR() = 0.174

Nmax() = 616sel = (25.0 ± 4.9)%

Background: NcharmC1 0 expected events

NcharmC1 0 = N1 loc · r · rA · BR(DC1) · sel

N1 loc = 95450

r =

charm+

CC

= 0.03

rA = A

charm 0

ACC 1

= 0.048 ± 0.006

BR(DC1) = 0.65

NcharmC1 0 = 18sel = (20.2 ± 3.7)%

Background: NcharmC3 0 expected events

NcharmC3 0 = N1 loc · r · rA · BR(DC1) · sel

N1 loc = 95450

r =

charm+

CC

= 0.03

rA = A

charm 0

ACC 1

= 0.063 ± 0.006

BR(DC1) = 0.35

sel = (50.8 ± 3.5)%

NcharmC1 0 = 32

Background: NWSK 1prong expected events

NWSK 1prong =

N0 loc = 23678

sel =

NWSK 1prong = 6.0

N0 loc

ANC 0

· sel

ANC 0 = 0.142

= 3.60 · 10-5wi

selected

NWSK generated in emulsion

Using WHINTER=0.720

19973

Signal and background expected without post-scanning cuts

BG: charm BG: WSK BG: Total Nmax

Old 1 0.1 - 0.1 5014

0 + 1prong 18 6 24 7519

0 + 3prong 32 0 32 3884

Using Feldman & Cousins (NOMAD) approach we obtain:

Sensitivity (zero events observed): P < 2.8 · 10-4

of the signal

= angle between parent particle and the mean of all primary tracks

c = angle between parent particle and the mean of all primary tracks, but the mean is done without the most far track from the parent

c

D

of the BG

= angle between parent particle and the mean of all primary tracks

c = angle between parent particle and the mean of all primary tracks, but the mean is done without the most far track from the parent

D

c cut

c > 1.8 rad

BG reduced to 15%

Tau reduced to 70%

Signal and background expected applying c cut

BG: TotalObserved (assumed)

Nmax

Old 1 0.1 0 5014

0 + 1prong 3.6 4 5263

0 + 3prong 4.8 5 2719

Using Feldman & Cousins (NOMAD) approach we obtain:

Sensitivity (zero events observed): P < 2.7 · 10-4

Oscillation limit: P < 2.9 · 10-4

chorus physics workshop, 7 june 2004, rome luca scotto lavina kinematical selection and fc analysis...

Documents