24.12.2012 Views

References - Bogoliubov Laboratory of Theoretical Physics - JINR

References - Bogoliubov Laboratory of Theoretical Physics - JINR

References - Bogoliubov Laboratory of Theoretical Physics - JINR

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

The MC production<br />

comprises three<br />

steps: first the events<br />

are generated, then<br />

the particles pass<br />

through a simulated<br />

spectrometer using<br />

a program based on<br />

GEANT3 [6] and finally<br />

the events are<br />

reconstructed using<br />

the same procedure<br />

applied to real data.<br />

For the first step<br />

the LEPTO 6.5 [7]<br />

event generator is<br />

used together with a<br />

leading order parametrisation<br />

<strong>of</strong> the<br />

unpolarised parton<br />

distributions. The<br />

MRST04LO set <strong>of</strong><br />

parton distributions<br />

is used in a fixedflavour<br />

scheme generation.<br />

This set <strong>of</strong><br />

parton distributions<br />

has a good description<br />

<strong>of</strong> F2 in the<br />

COMPASS kinematic<br />

region.<br />

NLO corrections<br />

are simulated partially<br />

by including<br />

Entries<br />

Data / MC<br />

Entries<br />

Data / MC<br />

3<br />

× 10<br />

Data<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c) MC<br />

T<br />

10<br />

5<br />

0 -3<br />

10<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

5<br />

10<br />

4<br />

10<br />

3<br />

10<br />

2<br />

10<br />

Entries<br />

Data / MC<br />

10<br />

-2<br />

10<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c)<br />

T<br />

Preliminary<br />

-1<br />

10 1<br />

xBj<br />

Data<br />

MC<br />

Preliminary<br />

1<br />

0.5 1 1.5 2 2.5 3<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

5<br />

10<br />

4<br />

10<br />

3<br />

10<br />

2<br />

10<br />

10<br />

p [GeV/c]<br />

T1<br />

1<br />

0 20 40 60 80 100<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c)<br />

T<br />

Data<br />

MC<br />

Preliminary<br />

p [GeV/c]<br />

1<br />

Entries<br />

Data / MC<br />

Entries<br />

Data / MC<br />

3<br />

× 10<br />

Data<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c) MC<br />

T<br />

6<br />

4<br />

2<br />

Preliminary<br />

0<br />

0 0.2 0.4 0.6 0.8 1<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

5<br />

10<br />

4<br />

10<br />

3<br />

10<br />

2<br />

10<br />

Entries<br />

Data / MC<br />

10<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c)<br />

T<br />

1<br />

0.5 1 1.5 2 2.5<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

5<br />

10<br />

4<br />

10<br />

3<br />

10<br />

2<br />

10<br />

10<br />

y<br />

Data<br />

MC<br />

Preliminary<br />

p [GeV/c]<br />

T2<br />

1<br />

0 20 40 60 80 100<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c)<br />

T<br />

Data<br />

MC<br />

Preliminary<br />

p [GeV/c]<br />

2<br />

Entries<br />

Data / MC<br />

Entries<br />

Data / MC<br />

3<br />

× 10<br />

Data<br />

COMPASS 2004<br />

2<br />

2 MC<br />

15<br />

High-p , Q >1 (GeV/c)<br />

T<br />

10<br />

5<br />

0<br />

1 10<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

5<br />

10<br />

4<br />

10<br />

3<br />

10<br />

2<br />

10<br />

Entries<br />

Data / MC<br />

10<br />

1<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

5<br />

10<br />

4<br />

10<br />

3<br />

10<br />

2<br />

10<br />

10<br />

1<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c)<br />

T<br />

5 10<br />

COMPASS 2004<br />

2<br />

2<br />

High-p , Q >1 (GeV/c)<br />

T<br />

5 10<br />

Preliminary<br />

2<br />

10<br />

2<br />

2<br />

Q [(GeV/c) ]<br />

Data<br />

MC<br />

Preliminary<br />

2<br />

2<br />

Σ(p<br />

) [(GeV/c) ]<br />

T<br />

Data<br />

MC:Compass<br />

MC:Default<br />

Preliminary<br />

2<br />

2<br />

Σ(p<br />

) [(GeV/c) ]<br />

T<br />

Figure 2: Comparison between data and MC simulations – The distributions<br />

and ratios <strong>of</strong> Data/MC for: inclusive variables: xBj, Q 2 ,y (1st row). For<br />

hadrons pT (2nd raw). For hadron momenta (3rd row), and also the comparison<br />

<strong>of</strong> MC with LEPTO default tuning.<br />

gluon radiation in the initial and final states (parton shower – PS).<br />

The fragmentation is based on the Lund string model [9] implemented in JETSET<br />

[10]. In this model the probability that a fraction z <strong>of</strong> the available energy will be<br />

carried by a newly created hadron is expressed by the Lund symmetric function f(z) =<br />

z−1 (1 − z) ae−bm2 ⊥ /z ,withm2 ⊥ = m2 + p2 ⊥ ,wheremis the hadron mass.<br />

To improve the agreement between MC and data, the parameters (a,b) in the fragmentation<br />

function are modified from their default values (0.3,0.58) to (0.6, 0.1).<br />

The transverse momentum <strong>of</strong> the hadrons, pT , at the fragmentation level is given by<br />

the sum <strong>of</strong> the pT <strong>of</strong> each hadron quarks. Then the pT <strong>of</strong> the newly created hadrons is<br />

described by three steering paramrters JETSET parameters: PARJ(21), PARJ(23) and<br />

PARJ(24). The default values <strong>of</strong> these three parameters are (0.36, 0.01, 2.0), and were<br />

modified to (0.30, 0.02, 3.5).<br />

354

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!