study of anti-quark flavor asymmetric via e906

Study of anti-quark flavor asymmetric via E906

Shiuan-Hal Shiu 2009/7/13

1

ContentsContents Introduction

Hardware operation

Fast Monte Carlo simulation

Future work

2

Proton is composed of three valence quark,gluons and sea

By contrast with other quarks, the up and down quark are very similar.

Because of the similarity, anti-down and anti-up quark distributions in the proton are assumed to be equal.

Is this true……? 4

?

ud

Is in the proton? Is in the proton? du

=

5

Light Antiquark Flavor Light Antiquark Flavor Asymmetry: Brief HistoryAsymmetry: Brief History

Naïve Assumption:

This is a common assumption until

1991

X

6


Naïve Assumption:

Gottfried Sum Rule:

1

2 20

1

0

[( ( ) ( )) / ]

1 2( ( ) ( ))

3 3

( )1

3 p p

p nG

p p

S F x F x x dx

u x d

i

x dx

f u d

New Muon Collaboration (NMC)

, Phys. Rev. D50 (1994) R1

SG = 0.235 ± 0.026

( Significantly lower than 1/3 ! )

xdxxFxF np /1

0 22

7


Naïve Assumption:

Gottfried Sum Rule:

NA51 (Drell-Yan) NA 51 Drell-Yan confirms

d-bar(x) > u-bar(x)

8


Naïve Assumption:

Gottfried Sum Rule:

NA51 (Drell-Yan)

E866/NuSea (Drell-Yan)

E866 ExperimentE866 Experiment The cross section of Drell-Yan process is

Here q1, q2 are the beam, target quark distribution.

Detector acceptance chooses xtarget and xbeam.

9

221122112

21

2

21

2 1

9

4xqxqxqxqe

sxxdxdx

d

xtarget xbeam Xbeam

Xta

rget

u

dud

E866 ExperimentE866 Experiment

10

Cryogenic Target System

Station1

Station2

Station3

SM3 Analyzing Magnet

Hadronic Calorimeter

ElectromagneticCalorime

ter

SM12 Analyzing Magnet

Hadron absorb

er

Ring-Imaging Cherenkov

CounterMuon

Detectors

2

212

1

221

xu

xd

xx

pp

pd

E866 use hydrogen and deuterium target

E866 ExperimentE866 Experiment

11

12

Advantages of 120 GeV Main Advantages of 120 GeV Main InjectorInjector

The past: Fermilab E866/NuSeaFermilab E866/NuSea

Data in 1996-1997 1H, 2H, and nuclear targets 800 GeV proton beam

The future: Fermilab E906Fermilab E906

Data in 2009 1H, 2H, and nuclear targets 120 GeV proton Beam

Fixed Target

Beam lines

Tevatron 800 GeV

Main Injector

120 GeV

221122112

21

2

21

2 1

9

4xqxqxqxqe

sxxdxdx

d

13

Follow basic design of MEast spectrometer :

Where possible and practical, reuse elements of the E866 spectrometer. Tracking chamber electronics Hadron absorber, beam dump, muon ID walls Station 2 and 3 tracking chambers Hodoscope array PMT’s SM3 Magnet

– Two magnet spectrometer – Hadron absorber within first magnet

– Beam dump within first Magnet – Muon-ID wall before final elements

New Elements– 1st magnet (different boost)

– Sta. 1 tracking (rates)

– Scintillator (age)

– Trigger (flexibility)

E866 Meson East Spectrometer

14

256 Hodoscopes

MWPC 5500 Channels

Station 2 & 3Drift Chambers1700 ChannelsMulti-hit TDC’s

Station 4 PropTubes 400 Channels

E906 Spectrometer: Bend Plane E906 Spectrometer: Bend Plane ViewView

Target

M1

M2

Sta.1

Sta

.2

Sta

.3

Sta

.4 M

uon

ID w

all

Measurements with the Drell-Yan Measurements with the Drell-Yan processprocess

Fermilab E906/Drell-Yan will extend these measurements and reduce statistical uncertainty.

E906 expects systematic uncertainty to remain at approx. 1% in cross section ratio.

15

2

212

1

221

xu

xd

xx

pp

pd

What is CODA?What is CODA? CODA (CEBAF Online Data Acquisition) is a

software DAQ system.

17VME

EB ET ER

User proc.

Disk

cMsg

RC platform

RUN control GUI

Log messag

e

NetworkEthernet

VME

Single board computer

ROC

ROC

msqld

cMlog DataBa

se

mSQL DataBa

se

TranslateTypical system

chart

COMPUTER

The CODA control panelThe CODA control panel

18

RCPlatform

MSQL Daemon

Event Transfersystem

EBe906

ERe906

GUI

Conceptual design of DAQ based Conceptual design of DAQ based on CODAon CODA

19

Unix host running CODA

Ethernet Hub

Data Storage

Online monitoring

Data Decoder Software

MWPC, HODO, Muon

Coincidence Register System, VME

CAMAC TDC for DC

Interface to VME

Accelerator scalers, VME?

Electronics House

Counting House

ROC1ROC2ROC3

L2 Trigger

CODA toolsCODA tools

20

ROCe906roc1

Ebe906daq2

ERe906daq2

CODA file

CODA toolsCODA tools Dbedit: can edit the information of database

21

CODA toolsCODA tools Xcefdmp:

event monitoring software.

22

View File modeOffline

monitoring

Spy Event modeOnline

monitoring

Installation of CODAInstallation of CODA 1. Setup a user account with the proper environment.

2. Setup the CODA database.

3. Setup MVME6100.

4. Transfer the ROC to a DAQ crate by downloading the KERNAL and CODA_ROC program on MVME6100.

5. Add a trigger supervisor.

6. Add more VME module.

23

Hardware configuration in IPASHardware configuration in IPAS

24

MVME6100Single board computer

SIS3610Trigger

supervisor

SIS3600Multi

eventsLatch

Trigger supervisor and latch Trigger supervisor and latch modulemodule

Latch module(SIS3600): When this module received

a trigger, it will capture the signal pattern from the input channel

Will not send any interrupt to single board computer

Trigger supervisor(SIS3610): When this module were

triggered it will send an interrupt to single board computer

25

CONTROL I/O

32 Data

Channl

SIS3600

Use .crl file to control CODAUse .crl file to control CODA

Most of the components in a CODA system are pre-defined by CODA and only require configuring

We can easily write CRL code to control CODA.

26

Booted

Configured

Download

Paused Actived

Terminated

Create

Configure

Configure

Prestart

download

Pause

GO

End

Terminate

Control the SIS3600Control the SIS3600 Functions of SIS3610 is well supported in

CODA, but SIS3600 is not. In order to control the SIS3600 we need to

write the driver which compatible with CODA. I consult the driver of SIS3610 and the manual

of SIS3600 to write a driver for SIS3600. By checking the input pattern, the SIS3600

can work correctly.

27

… 1000… 00011000 …

About the fast monte carlo About the fast monte carlo simulationsimulation

Base on Fortran language. Modified version of E866/E772 “fast” Monte Carlo

Code to include E906 geometry. Only traces dimuon from Drell-Yan events, and the

decays of J/. Have 2 main parts 1. Muon pair creation2. Detector Magnetic field is simplified, but muon energy loss

and multiple scattering are included.

29

WantWantOriginal

We want to observe events of low We want to observe events of low massmass

30

Mass (GeV)

Mass (GeV)

acce

ptan

ce ac

cept

anc

e

The mass and the magnet currentThe mass and the magnet current

31

M1

M2

32


M1

M2

33


M1

M2

34


M1

M2

Configure M1 and M2Configure M1 and M2

Configuration file

35

M1M1M2M2

How to obtain the acceptanceHow to obtain the acceptance The top diagram is

the mass distribution of generated dimuon pairs.

The middle diagram is the mass distribution of reconstructed dimuon pairs.

The bottom diagram “Acceptance” as a function of dimuon mass. 36

Thrown

Reconstruction

Acceptance

Acc

epta

nce

Cou

nts

Cou

nts

GeV

GevGeV

GeV

The result of varying M1 The result of varying M1 currentcurrent

The mass region of generated dimuon pairs is from 1Gev to 15Gev

Green line is the acceptance value with the original M1 current setting.

By increasing the current we find that the peak of acceptance is shifting to high mass end.

Reducing the M1 current can increase the acceptance in the low-mass region.37

M1x0.1

M1x0.5

M1x1.0

M1x1.5

M1x2.0

GeV

Acc

epta

nce

FIX M1*1.5 varying M2FIX M1*1.5 varying M2

M2*1.5 M2*1 M2*0.5

38

Reconstruction

Thrown

M2x0.5

M2x1.0

M2x1.5

Acc

epta

nce

Cou

nts

Cou

nts

GeV GeV GeV

Acc

epta

nce

GeV

Reconstruction

Reconstruction

Thrown Thrown

Acceptance

Acceptance

Acceptance

M2*1.5 M2*1 M2*0.5

39

FIX M1*1 varying M2FIX M1*1 varying M2

M2x0.5

M2x1.0

M2x1.5

GeV GeV GeV

Acc

epta

nce

Cou

nts

Cou

nts

Acc

epta

nce

GeV

Reconstruction

Thrown

Reconstruction

Reconstruction

Thrown Thrown

Acceptance

Acceptance

Acceptance

M2*1.5 M2*1 M2*0.5

40

FIX M1*0.1 varying M2FIX M1*0.1 varying M2

M2x0.5

M2x1.0

M2x1.5

GeV GeV GeV

Acc

epta

nce

Cou

nts

Cou

nts

Acc

epta

nce

GeV

Reconstruction

Thrown

Reconstruction

Reconstruction

Thrown Thrown

Acceptance

Acceptance

Acceptance

Conclusion Reducing the M1 current can increase the

acceptance of low mass region dimuons.

Adjusting the M2 current does not change the acceptance region significantly.

Howerer, adjusting M2 current can change the acceptance.

41

Dump/Target separationDump/Target separationM1x1M1x1

Software cuts conditionsPurple: all events

Green: xF>0 and M>4.5 GeV and pz>20 GeV

Blue: Green and |ytrack|>2.25 in at z=0 (zdump)

Red: Blue and |ytrack|<10.0 in at z=-60 (zstart)

42

Thrown

Count

sC

ounts

Reconstruction

ZTarge

tTarge

tDum

pDum

p

Dump/Target separationDump/Target separation

43

Thrown

Reconstruction

Z

Count

sC

ount

s

Thrown

Reconstruction

Z

Count

sC

ount

s

M1x0.5M1x0.5 M1x0.1M1x0.1

Mass cut > 4.5

Dump/Target separationDump/Target separation

44

Thrown

Reconstruction

Z

Count

sC

ount

s

Thrown

Reconstruction

Z

Count

sC

ount

s

M1x0.5M1x0.5 M1x0.1M1x0.1

Mass cut < 4.5

Dump/Target separationDump/Target separationM1x0.5M1x0.5

45

Counts

Counts

Counts

Counts

ZZ

Z Z

1

2

3

4

Counts

Counts

Counts

Counts

No cut

Z

Z Z

Z

Green+

Y(zdump)>2.25

Xf>0Retrace M<4.5Pz>10

Blue+

Y(zstart)>10

1

2

3

4

Mass cut < 4.5

Reducing the M1 current will lead the z resolution worse.

Dump/Target separationDump/Target separationM1x0.1M1x0.1

46

Counts

Counts

Counts

Counts

No cut

Z 1

Z Z

Z

2

3

4

Green+

Y(zdump)>2.25


Blue+

Y(zstart)>10

Mass cut < 4.5

Change the target locationChange the target location

47

50 inche

s


48

130 inche

s


1

2

3

4

M1*0.5

49

Counts

Counts

Counts

Counts

Z

Z Z

Z

No cutGreen

+Y(zdump)>2.

25

Blue+

Y(zstart)>10


1

2

3

4

Mass cut < 4.5


1

2

3

4

M1*0.1

Changing target location can slightly improve the z resolution.50

Counts

Counts

Counts

Counts

Z

Z Z

Z

No cut

Green+

Y(zdump)>2.25

Blue+

Y(zstart)>10Xf>0Retrace M<4.5Pz>10

1

2

3

4

Mass cut < 4.5

Check the relations of retrace Check the relations of retrace mass and retrace zmass and retrace z After applied the

cut(mass >4.5 Gev) , we can see that the events are almost spread in the region which less than z=0.

51

18

Z(inches)

161412108

6

4

2

00

-20

-60

-80

-40

-100

20

40

60

100

80

mass

(GeV

) M1x1M1x1

Check the relations of retrace Check the relations of retrace mass and retrace zmass and retrace z

After applied the cut (mass <4.5 Gev), we can see that the events are still scattered throughout the x axis.

52

M1x0.5M1x0.518

Z(inches)

161412108

6

4

2

00

-20

-60

-80

-40

-100

20

40

60

100

80

mass

(GeV

)

Check the relations of retrace Check the relations of retrace mass and retrace zmass and retrace z

53

18

Z(inches)

161412108

6

4

2

00

-20

-60

-80

-40

-100

20

40

60

100

80

mass

(GeV

) M1x0.1M1x0.1

Conclusion Conclusion

After reducing the M1 current, more and more events can not get a correct retraced Z value.

Changing target location can slightly improve the Dump/Target separation rate.

Low mass events are affected by the multiple scattering seriously.

Low mass events generally have small opening angle, this will increase the difficulty to retrace the Z position.

54

55

Drell-Yan decay angular Drell-Yan decay angular distributionsdistributions

A general expression for Drell-Yan decay angular distributions:

2 21 31 cos sin 2 cos sin cos 2

4 2

d

d

Production plane

Decay plane

Collins-Soper frame is defined in dimuon CMS.

The Z-axis is the bisector of Pbeam and –Ptarget .

Θ is the angle between Pμ+ and Z-axis .

Φ is the angle between production plane and decay plane .

Pbea

m

Ptarge

t

Y

X

Z

Pμ+

Angular distribution reconstructAngular distribution reconstruct The top diagram

is the cosΘ distribution of generated events.

The middle diagram is the reconstructed cosΘ.

The bottom diagram “Acceptance” is obtained by

56

Counts

Counts

Counts

Counts

Acc

epta

nce

Thrown

Reconstruction

Acceptance

cosΘ (cs)

thrown

tionreconstrucacceptance

Testing the fast monte carlo Testing the fast monte carlo simulation programsimulation program Thrown a

shape angular distribution.

Fit the reconstruction data using function:

Check the fitting result λ within error bar is consist with the input value or not.

57

)(cos1 2

Counts

Thrown

cosΘ (cs)

)(cos1 2

acceptance))(cos1( 2

Monte carlo data fitting resultMonte carlo data fitting result

58

)(cos5.01 2 )(cos11 2 cosΘ (cs)

cosΘ (cs)Input function : Input function :

λ=0.687958Err=0.10103

3

λ=0.924215Err=0.0954

346

Monte carlo data fitting resultMonte carlo data fitting result

59)(cos31 2 )(cos21 2 Input function :Input function :

cosΘ (cs)

cosΘ (cs)

λ=1.50516Err=0.125

602

λ=2.34244Err=0.1655

32

ConclusionConclusion The fitting in the case λ=1 is successful, but

fails in the other case.

I will spend more time to checking the fitting program and fast monte carlo program to find out why the fitting fails.

Finally, the goal is adding the φ angular distribution to study the θ and φ acceptance and resolutions for the J/Psi and D-Y events.

60

62

Timeline of E906Timeline of E906

2002: E906 Approved by Fermilab PAC 2006: E906 funded by DOE Nuclear Physics 2008, Dec: Stage-II approval by Fermilab Director and

MOU between Fermilab and E906 Collaboration finalized. Construction and installation of spectrometer and readout

electronics to be done in 2009 and upper half of 2010. First beam expected in the fall of 2010!

Man PowerMan Power Institute of Physics, Academia Sinica (IPAS):

Wen-Chen Chang, Ping-Kun Teng, Yen-Chu Chen (stationed at FNAL), Da-Shung Su (engineer)

Shiuan-Hal Shiu (Ph.D. student, preparing his proposal and plan to go to work with Ron Gilman on level-2 trigger afterwards.)

Bo-Ru Lin (research assistant, leaving the group in July and joining Colorado group as Ph.D. student.)

Jia-Ye Chen (Expect to work as postdoc next summer.) Ling-Tung University (LTU):

Ting-Hua Chang Summer-study students: will help on the mass production of

electronics. There are people from “National Kaohsiung Normal

University” (NKNU) who wish to join E906 experiment: Rurng-Sheng Guo (member of FNAL E665 and CDF Collaboration) Su-Yin Wang (Ph.D. student in NKNU, under the supervision of

Rurng-Sheng Guo)

63

Responsibility of Taiwan groupResponsibility of Taiwan group

Build 400 preamplifier-discriminators cards (6400 channels in total) for 5500 tracking channels in Station 1 MWPC and for 400 channels Proportional tubes in Station 4,

Build the readout system for Coincidence Registers (CR), which consists of 110 CR modules (7040 channels in total) for 5500 tracking channels in Station 1 MWPC, for 400 channels of Proportional tubes in Station 4, and for 320 channels of hodoscope planes in all stations,

Participate in the LVL2-Trigger project – two CEAN V1495 FPGA logic units purchased and one Ph.D. student will be available in fall to work on this project.

64

preamplifier-discriminators cardspreamplifier-discriminators cards

65

16

ch

EC

L i

np

uts

16

ch

EC

L i

np

uts

VME Backplane

NIM/ECL ↔ TTL

A24

~A

31

D16

~D

31

Eth

ern

et

10/1

00M

bp

sARM processorAT91SAM9260

64

ch

EC

L t

o T

TL

[2

]

Front panel

TransceiverDM9161

Deb

ug

U

AR

T

COM port

JTAG-ICE

16

ch

EC

L i

np

uts

FPGAActel A3P-1000

VM

E b

us

P1

VM

E b

us

P2

2ch

NIM

I/O

64MB SDRAM

256MB x 16 NAND Flash [5]

Bu

s b

uff

er

x5

[6

]B

us

bu

ffe

r x

3 [

6]

A01

~A

23,

D0~

D15

, D

S, A

S, W

rite

, DT

AC

K

Trigger, Fast clear

Reset

64MB x 32 SDRAM [4]

BusyReady/

- +

16

1

CLKDATAx xx x

6ch

EC

L c

on

tro

l I/O

8ch TTL to ECL [1]JTAG

16

ch

EC

L i

np

uts

64MB SDRAM4Kbit x 32 DP-SRAM [3]

Coincidence Register (CR) modulesCoincidence Register (CR) modules

[1] ON Semiconductor, MC10H124[2] ON Semiconductor, MC10H125[3] IDT, IDT70V24, 4Kbit x16 Dual-Port SRAM[4] Winbond, W9825G6, 256Mbit x16 SDRAM[5] Micron, MT29F2G16AACWP, 2Gb x16 NAND Flash[6] TI, SN74VMEH22501A, VME bus transceiver

A32

, Ad

dre

ssin

g

rota

ry s

wit

ches

Th

e d

efin

itio

n o

f E

CL

co

nn

ecto

r

66

E906 detector readout and DAQ systemE906 detector readout and DAQ system

67

Some information and estimationSome information and estimation Main Injector RF clock frequency: 53 MHz.

Beam structure: 1013 protons in a 5 s slow extraction spill every minute. Beam intensity: 2*1012/sec.

Level-2 trigger latency: Master Trigger OR decision time= 91ns.

Level-1 trigger rate (X): MWPC designed Singles rates:53MHz.(?) The total rate of single muons traversing the detector and

passing the trigger matrix tracking will be approximately 100 kHz with the LH2 target and 150 kHz with the LD2 target (both cases include tracks originating in the beam dump).

Event size of MWPC W/O data reduction: 5500bit=0.8kB

Depth of memory buffer: 53MHz*91ns~ 5 events. Size of total memory to buffer 5 events=0.8kB*5=4kB.

VME transfer throughput:0.8kB*1kHz=0.8MB/per sec << Optic fiber ~100 MB/sec and VME 160 MB/sec. Deadtime-free is possible.68

What I have learned from recent What I have learned from recent works. works. How to write a driver for VME module.

How to install and use CODA DAQ system.

Simulation and analysis methods in high energy physics.

More knowledge about nuclear physics.

69

backup

72

The physics motivation Why we choose measuring Drell-Yan process. What we can learn from Drell-Yan process.

*The past experiments *The results of E866 Etc.

E906 experiments Introduce the experiment’s framework. Comparing with the past experiments.

73

Why we want to study the low mass region. Introduce the goal we want to reach.

What the simulation program do. Introduce the program.

The detector geometry and configure file. Introduce the detector geometry and the fast

monte carlo simulation flow paht. The results of this simulation.

74

Why we want to study the angle distribution

How I define the axis

The result of this study

75

What is CODA? CODA flow chart.

Introduce sis3610 and sis3600 What is trigger supervisor and latch.

The flow chart of controlling the vme module Make a description of how the data acquisition

work.

76

Introduce the works Taiwan group is going to do.

What parts I will do in the future. What is Level2 trigger.

What I have learned from recent works.

77

What will we learn?

– d-bar/u-bar in the proton– Nuclear effects in the sea quark distributions– High-x valence distributions– Partonic energy loss in cold nuclear matter

78

Momentum Difference(X) (M1*0.5)

79

Momentum Difference(Y) (M1*0.5)

80

Momentum Difference(Z) (M1*0.5)

81

Pt/Pz(positive) (M1*0.5)

82

Pt/Pz(negative) (M1*0.5)

83

Opening angle of muon pairs

fe198v5m101mrtr.le.4.5zrtr.le.0

fe198v5m101mrtr.le.4.5zrtr.ge.0

fe198v5m101_dumpmrtr.le.4.5zrtr.le.0

fe198v5m101_dumpmrtr.le.4.5zrtr.ge.0

84

Momentum Difference(X) (M1*0.1)

85

Momentum Difference(Y) (M1*0.1)

86

Momentum Difference(Z) (M1*0.1)

87

Pt/Pz(positive) (M1*0.1)

88

Pt/Pz(negative) (M1*0.1)

89

Opening angle of muon pairs

fe198v5m101mrtr.le.4.5zrtr.le.0

fe198v5m101mrtr.le.4.5zrtr.ge.0

fe198v5m101_dumpmrtr.le.4.5zrtr.le.0

fe198v5m101_dumpmrtr.le.4.5zrtr.ge.0

90

Theta y distribution of positive muon at zstart

92

Theta y distribution of positive muon at zstart

93

Theta y distribution of negative muon at zstart

94

Theta y distribution of negative muon at zstart

95

Momentum distribution of positive muon at zstart

96

Momentum distribution of positive muon at zstart

97

Momentum distribution of negative muon at zstart

98

Momentum distribution of negative muon at zstart

99

X distribution of positive muon at zstart

100

X distribution of positive muon at zstart

101

X distribution of negative muon at zstart

102

X distribution of negative muon at zstart

103

Theta x distribution of positive muon at zstart

104

Theta x distribution of positive muon at zstart

105

Theta x distribution of negative muon at zstart

106

Theta x distribution of negative muon at zstart

107

Xf distribution of muon at zstart

108

Xf distribution of muon at zstart

109

X1 distribution of muon at zstart

110


111


112


113

Total momentum distribution of muon at zstart

114

Total momentum distribution of muon at zstart

115

116

Station 1 Chamber RatesOccasionally a muon showers in the

absorber If this happens in the center of the

absorber, no effect is seen as shower is also absorbed

If this happens in the last few inches of the absorber, shower can create extremely large rates in Station 1 (of low momentum particles)

Solution is to have an absorber-free region at the end of the field volume and use field as a sweeper

In Solid Iron magnet, there is no absorber-free sweeper region! (Can we find a wide gap sweeper magnet?)

Requires GEANT MC to see magnitude of effect

Absorber and B Field

Sta. 1



B Field only

119

Data Sources

256 Hodoscopes—bit latches

MWPC 5500 Channels

bit latches

Station 2 & 3Drift Chambers1700 ChannelsMulti-hit TDC’s

Station 4 PropTubes 400 Channels

Bit latches

The problem we meetThe problem we meet There can be thousands to millions of

events occurring per second.

Detectors are very large - containing many thousands of individual channels.

Events are different sizes.

Events occur at random.

Only a few events are interesting.

120

cut

121

Use the E906 Fast Monte Carlo simulation for study configuration data.

The configuration file is “fe198v5.dat”.

By changing the entry “current and step to scale” we can adjust the M1 or M2 current.

Fix M2 current and vary M1 currentFix M2 current and vary M1 current

The input to the simulation is decided by Ykick*input/2000 and

the tracking plane from #2 to #13 will affected

by this factor.

The input to the simulation is decided by Ykick*input/2000 and

the tracking plane from #2 to #13 will affected

by this factor.122

Fix the M1 current, and change the M2 current (Ykick).

Fix M1 current and vary M2 currentFix M1 current and vary M2 current

123

Cuts conditionsPurple: all events

Green: xF>0 and M>4.5 GeV and pz>20 GeV

Blue: Green and |ytrack|>2.25 in at z=0 (zdump)

Red: Blue and |ytrack|<10.0 in at z=-60 (zstart)

Moving the cut condition at zdump to Station 1 and change the value.

124

Original Changed

125

Original Changed

Try |ytrack|>8 in at z=238 (before station1)126

M1*0.5

1

2

3

127

4

M1*0.1

1

2

3

4

Changing cut condition can not improve the z resolution.128

Gottfried-Jackson frame(GJ)Gottfried-Jackson frame(GJ) Gottfried-Jackson

frame is defined in dimuon CMS .

The Z-axis is the direction of Pbeam .

130

Production plane

Decay plane

Helicity frame(HX)Helicity frame(HX) Helicity frame is

defined in dimuon CMS .

The Z-axis is Pbeam + Ptarget.

131

Production plane

Decay plane

TriggerSIS3610

VME flow chartVME flow chart

135

Reach the event number?

No

Yes

EBCreateevent

CODAStart

Initialize3600/3610

Clear/Reset 3600/3610

End

Events into SIS3600

Setup phase

SIS3610Send

Interrupt toMVME6100

Latch phase

Exit phase

Data to storagephase

MVME6100Take DataFrom 6100

ERET

DISK

No

Yes

EXIT?

Why we need DAQ?Why we need DAQ? The goal of a nuclear physics experiment is to get data

about nuclear interactions.

Particles pass through detectors which generate electrical signals contain information about the particles type, energy, trajectory .

The complete set of signals which describe a single nuclear interaction is called an Event.

The data acquisition system digitizes, formats and stores this information in a way which can be retrieved for later analysis.

136

The structure of DAQThe structure of DAQ Triggering (choosing events we want)

Readout (digitizing detector signals)

Event formatting (standardize what we’re saving)

Event building (putting fragments together)

Event transport (make events available to all)

Event storage (save data for analysis)

Run Control (configure-start-stop experiments)

Monitoring (tell me what’s going on)

137

A DAQ system exampleA DAQ system example

138

Trigger

ReadoutEvent

formattingEvent

building Event

transport

Event storage

Run control

Monitoring

Why using CODA?Why using CODA? CODA is a software toolkit with some specialized

hardware support.

Modular software components use the network for inter-process communication and event transport.

Use open standards and minimize the use of commercial software while maximizing use of commercial hardware.

DAQ systems for each experimental Hall can be “built-up” from common components to fit their needs.

139

VME VME (Versa Module Europa)(Versa Module Europa) VMEbus is a computer bus

standard, originally developed for the Motorola 68000 line of CPUs, but later widely used for many applications.

8/16/24/32/64 bit bus,24/32/64 support BLT(block transfer).

The data is transfering in the back plane bus.

140

Angular distribution Angular distribution The angular distribution may be changed by

the detector acceptance.

141

?!!

What is Level 2 triggerWhat is Level 2 trigger

Level 1 trigger

Level 2 trigger

142

Backplane bus

SIS3600

Computingmodule

SIS3610

interrupt data

Hardware configuration in IPASHardware configuration in IPAS

143

MVME6100Single board computer

SIS3610Trigger

supervisor

SIS3600Multi

eventsLatch

BackPlane Bus


The target original location is at -70 to -50. We change it to -150 to -130.

144

147

is more suppressed than in the proton since

(Field and Feynman 197

(pQCD calculati

Paul

ons b

Origins of ( ) ( )?

y Ross and Sachrajda)

(Bag model calculation by Signal, Thomas, Schreib

7

i blocking by the valenc

)

e quarks

u x d x

g uu g dd

p uud

Quark spectrum includes a bound state plus the

polarized negative

er)

(Diakonov

Chiral quark-soli

, Pobylitsa, Poly

and positive D

akov, Wakamats

irac continuu

u, Kubota)

ton model

Instanton m l

m

ode

,

Sta

(Dorokhov, Kochelev)

(Bourrely, Butistical model ccella, Sof

,

f er)

etc. L R R L L R R Lu u d d d d u u

The valence quarks affect the Dirac vacuum and the quark-antiquark sea

質子 p (uud). 中子 n (udd)

148

What is the origin of ?What is the origin of ?du,

In general : 1, 0, 0

Fermilab E866 Measurements

149

800 ( ) / ( )GeV p d X p p X

150

Proton Economics Total of 5.2 X 1018 protons (over 2 years)

Maximum instantaneous rate of 2 X 1012 proton/sec Based on E866 experience with target related rate

dependence—balance systematic and statistical uncertainties Station 1 chamber rates.

Possible delivery scenario: 5 sec spill of 1 X 1013 protons each minute Longer spill (5 sec) desirable over 5-1 sec spills

Leading order structure function

151

],,[, 22222 QxqQxqxeQxF ii

ii

The leading-order structure function F2 is given by quark-momentum distributions in the nucleon

How the fast monte carlo generate muon 1. random create virtual photon mass. 2. random create feynman X 3. determine Pt^2(max), X1, X2 4. determine Pt 5. random phi angle in ? Frame. 6. from dimuon CMS ,the virtual photon average separate it’s mass

to two muons. 7.by the difference mass condition to determine cos(theta). 8.calculate Pl, Pt. 9.boost along Pt and create phimu, phimu is the mu+ with respect

to the pt of the gamma 10.boost along Pl to overall CMS frame, and check Pz grate or less

pmumin. 11.boost to lab.

152

study of anti-quark flavor asymmetric via e906

Documents

brief historynave assumption

gottfried sum rule

na51 drellyanna

common assumption

valence quark

quark distributions

target quark distribution

muon id wallsstation