Dubna1mHBC dp
Deuterónové ožiarenia 1m HBC
DODAT POPIS STUKTURY DST PRE FULL A MINI FORMAT
Expozície číslo: 33, 34, 35 hybnosť d zväzku 3,35 GeV/c,
v rokoch 1971 a 1972. Zväzok bol produkovaný strippingom a obsahoval určitú prímes protónov. V expozícii č. 35 bola použitá separácia deuterónov. (Odvolávame sa na ne ako staré alebo tiež nepolarizované expozície.)
Expozície čísla: 49, 51, 59, 62, 63, 64, 66, 68 boli uskutočnené vo vektorovo-polarizovaných ( P1-88-6, Dubna 1988 https://lib-extopc.kek.jp/preprints/PDF/1988/8806/8806011.pdf ; 2003_vvg_trudy/139.doc: " Byli predusmotreny tri metki: "1" i "2" oznacali poljarizovannyj pucok dejtronov s protivopoloznym napravleniem spina jadra, a metka "0" sootvetstvuet nepoljarizovannomu pucku.") zväzkoch deuterónov priamo urýchlených so stupňom polarizácie 41% (ibid Zakljucenie ... Pz = -0,41 +- 0,09 ) a s hybnosťou 3,34 GeV/c (ibid). (V expozíciách 49 a 51 zväzok vstupoval zo strany mechanizmu rozšírenia – piest, pre expozície 59, 62,63,64,68 komora bola otočená o 180°. )
Formáty DST plný (FULL) formát a MINI formát.
Vysvetlenie rozdielu medzi hypotezou a kanalom reakcie
Hypoteza pod danym cislom, ako je uvedena v zozname, moze zahrnut — v dosledku zakona zachovania elektrickeho naboja a baryonoveho cisla — niekolko kanalov reakcie :
- kanal bez neutralnej castice (4C-, 3C-, 2C-fit)
- kanal s 1 (jednou) neutralnou casticou (1C-fit)
- kanal s 2 (dvoma) a viac neutralnymi casticami (0C-fit tiez nazyvany ako NOFIT).
- Ak zakon zachovania naboja automaticky zabezpeci aj splnenie zakona zachovania baryonoveho cisla, hypoteza obsahuje vsetky tri vyssie uvedene typy kanalov reakcie, priklad:
- [math]\displaystyle{ dp \rightarrow dp \pi^+ \pi^- }[/math],
- [math]\displaystyle{ dp \rightarrow dp \pi^+ \pi^- \pi^0 }[/math],
- [math]\displaystyle{ dp \rightarrow dp \pi^+ \pi^- X^0 }[/math], kde [math]\displaystyle{ X^0 = k \pi^0, k\ge2. }[/math]
- Ak pri splnenom zakone zachovania naboja, zakon zachovania baryonoveho cisla vyzaduje doplenie vystupneho kanala reakcie o jeden neutron , hypoteza obsahuje dva posledne typy (2. a 3.) vyssie uvedenych kanalov reakcie, priklad:
- [math]\displaystyle{ dp \rightarrow ppn }[/math],
- [math]\displaystyle{ dp \rightarrow ppX^0 }[/math], kde [math]\displaystyle{ X^0 = n + k \pi^0, k\ge1. }[/math]
- Ak zakon zachovania baryonoveho cisla vyzaduje doplenie vystupneho kanala reakcie o dva alebo viac neutronov , hypoteza obsahuje iba posledny typ (3.) vyssie uvedenych kanalov reakcie, priklady:
- [math]\displaystyle{ dp \rightarrow p \pi^+ X^0 }[/math], kde [math]\displaystyle{ X^0 = n n + k \pi^0, k\ge0, }[/math]
- [math]\displaystyle{ dp \rightarrow \pi^+ \pi^+ X^0 }[/math], kde [math]\displaystyle{ X^0 = n n n + k \pi^0, k\ge0. }[/math]
Pripad (event) v novych expoziciach (49, 51, 59, 62, 63, 64, 66, 68) moze obsahovat maximalne 2 (dve) hypotezy v starych expoziciach (33, 34, 35) boli akceptovatelne 3 (tri) hypotezy.
Plný format (FULL) binárny file
| File | Počet prípadov | Inkonzistentná štruktúra | Expozícia |
|---|---|---|---|
| dp6668.dvo | 59612 | 0 | 66, 68 |
| dp62.dvo | 16870 | 34 | 59 ? |
| dp59.dvo | 14573 | 30 | 45 |
| dp34.dvo | 18777 | 35 | 34 |
| Total | 109832 | 99 | all |
Štruktúra zápisu
Nbytes, Nwords, b[0], b[1], …, b[Nwords-1], Nbytes,
kde
- Nbytes je štvorbyteový integer a Nbytes = 4 * (Nwords + 1) dlžka zápisu v štvorbyteových slovách, samotné slová Nbytes sa do dlžky nezapočítavajú,
- Nwords je štvorbyteový integer a vyjadruje dlžku poľa
- b[0], b[1], … b[Nwords-1], samotné slovo Nwords sa do dl´žky nezapočítavá,
- b[0], b[1], … b[Nwords-1] sú štvorbyteový integer alebo float,
- b[0], b[1] sú integer (pre lab Dubna)
- b[2], … b[Nwords-1] sú float
- Full format DST pre dp obsahuje iba prípady s LAB = 1 (Dubna)
- NOFIT hypotesis number > 100.5
Inkozistentná štruktúra znamená, že Nbytes = 4, Nwords = -1 (integer), Nbytes = 4 . Takéto prípady ignorujem.
Bezpečné čítanie (ROOT, C++):
union u1 { int ib [1000]; float b[1000]; } equiv;
union u2 { int N; float AN; } len, lenw, lentr;
int *pN = &len.N;
int *pNtr = &lentr.N;
int *pNw = &lenw.N;
...
f.read((char*)pN, 4); // record length in bytes+4 , header length !!!
if (f.eof()) { break ; }
f.read((char*)equiv.b, len.N);
f.read((char*)pNtr, 4); // record length in bytes+4 , trailor length !!!
...
if (len.N != lentr.N) {
cout << "Error header length " << len.N << " not equal trailor length " << lentr.N << endl;
break;
}
if (len.N != 4*(equiv.ib[0]+1 )) {
inconsinstency++;
continue; // skip this record
}
Výsledky
Flagy polarizacie:
- pre Expo <= 34 treba dodať ručne flag -1, non-polarized beam, DST ho neobsahuje
- pre Expo => 59 integer ib[1] = polflag*10000 + frame*10 + ordre
Číslo expozície na DST nie vždy odpovedá názvu súboru, správne číslo expozície je v názve súboru :
dp6668.dvo LAB-EXP-INDEX = LAB * 10000 + EXPONUMBER = 10066., 10068. dp62.dvo LAB-EXP-INDEX = LAB * 10000 + EXPONUMBER = 10059. dp59.dvo LAB-EXP-INDEX = LAB * 10000 + EXPONUMBER = 10045. dp34.dvo LAB-EXP-INDEX = LAB * 10000 + nieco *100 + EXPONUMBER = 10334.
Nedôslednosti:
Exposition 34, “dp34.dvo”
- počet hypotéz=1 , ale w =0.333, 0.5, 1.0
- počet hypotéz=1 pre cely file, w=0.333, ale v hlavicke su pritomne 2 hypotezy. ktoru brat?
- hmota primarky PROTON preskocit ? Co s nimi ? PRESKOCIL SOM
- hyp number > 100 a ND = 4
- v oblasti trackov je (ScanTop/100 +2) drah, co je v poriadku !
- SCANTOP=200, two hyps HYP1=6 ND1=4; HYP2=8 ND2=4 su CISTE NOFIT
- SCANTOP = 700 !!!
- po testoch na hmotnost beamu a topologiu z 18777 pripadov zostane 854 !!!
- osetrit FITNOFIT = 0
Exposition 59, “dp59.dvo”
- exp 59(?) v oblasti trackov je (ScanTop/100 +2) drah, co je v poriadku !
- po testoch na hmotnost beamu a topologiu zo 14573 pripadov zostane 14543, 30 pripadov vypadlo na inkonzistenciu
- cislo expozicie v samotnom file je ale 45. Ovsem to je 3He pri 13 GeV/c! Takze drzat sa cislovania exp co je v nazve file-u
Exposition 62, “dp62.dvo”
- exp 62(?) v oblasti trackov je (ScanTop/100 +2) drah, co je v poriadku !
- po testoch na hmotnost beamu a topologiu zo 16870 pripadov zostane 16836, 34 pripadov vypadlo na inkonzistenciu
- cislo expozicie je ale 59. Dohady: nespravne zadali cislo expozicie pri pisani DST. Exp 59 podla zoznamu exp. ma 31105 snimkov, exp 62 podla zoznamu exp. ma 65110 snimkov. Predpoklad ze ide tu o exp 62 nie je v rozpore s poctom pripadov.
Exposition 66, 68 “dp6668.dvo”
- v oblasti trackov je (ScanTop/100 +2) drah, co je v poriadku !
- po testoch na hmotnost beamu a topologiu zo 59612 pripadov zostane 59612, 0 pripadov vypadlo na inkonzistenciu
Medzi dvoma track regionmi alebo aj za poslednym mozu byt urcite neuzitocne slova.
Asi kvoli tomu, aby track region pre druhu hypotezu zacal od okruhleho cisla.
| Input file | Events read in | Accepted events | Inconsinstency | Proton beam in hypo |
|---|---|---|---|---|
| dp6668.dvo | 59612 | 59612 | 0 | 0 |
| dp62.dvo | 16870 | 16836 | 34 | 0 |
| dp59.dvo | 14573 | 14543 | 30 | 0 |
| dp34.dvo | 18777 | 854 | 35 | 17888 |
| Total statistics | 109832 | 91845 | 99 | 17888 |
Suma vah (WEIGHT) = 91532.8, pri entries 94649 (hypotez) (rozdelenie dlzky primarky )
190 hypos 3C fit, test pre 4C ma byt 1.5 < ND < 4.5 Expos zostali : 34, 45, 59, 66, 68
Poznamka:
Proton beam :
NUMofHYPS = 1 vzdy a vahy su rozne:
WEIGHTS = 0.33 (2174.4 hyp)
WEIGHTS = 0.5 (1504.5 hyp)
WEIGHTS = 1.0 (8349 hyp)
--------------
Total = 12024 hyp
Deuteron beam: =
NUMofHYPS = 2 WEIGHTS = 0.5 (1402.0 hyp)
NUMofHYPS = 1 WEIGHTS = 1 (88485.0 hyp)
NUMOfHYP = 1 WEIGHTS = 0.5 (190.0 hyp)
NUMOfHYP = 1 WEIGHTS = 0.33 (63.999203 hyp)
---------------
Total 90140 hyp
Rolls from dp FULL DST pre exp 59 a 62
Exp 59
18 19 20
21 23 24 25 26 27 28 29 36 37 38 39 40
41
Exp 62
40
45 46 47 48 49 50
52 54 55 57 59 61 62 65 87 88 89 90
92 98
109
137
Flagy polarizacie po expoziciach -- dp FULL DST iba deuteron beam
Exp Pol flag
|
4 | 32 | 33 | 34 | 59(45) | 62(59) | 63 | 64 | 66 | 68 | Sum |
|---|---|---|---|---|---|---|---|---|---|---|---|
| -1 | 0 | 0 | 0 | 854 | 0 | 0 | 0 | 0 | 0 | 0 | 854 |
| 0 | 0 | 0 | 0 | 0 | 6080 | 8389 | 0 | 0 | 4 | 9083 | 23556 |
| 1 | 0 | 0 | 0 | 0 | 7310 | 7096 | 0 | 0 | 20530 | 459 | 35395 |
| 2 | 0 | 0 | 0 | 0 | 1152 | 1350 | 0 | 0 | 20515 | 9012 | 32029 |
| 3 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 5 | 0 | 6 |
| 4 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 5 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 4 | 0 | 4 |
| 6 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| Sum | 0 | 0 | 0 | 854 | 14542 | 16836 | 0 | 0 | 41058 | 18554 | 91844 |
| Footnote | Overflow Polarization flag: 1, weight = 1 | ||||||||||
Redukovaný format (MINI, MIKRO) binárny file
| File | Number of hypo | Exposition |
|---|---|---|
| allppn.dvo | 102778 | 34, 59, 62, 63, 64, 66, 68 (all “dp→ppn” ) |
| all2pr.dvo | 125673 | 32, 33, 34, 59, 62, 63, 64, 66, 68 (without “dp→ppn” ) |
| all4pr.dvo | 8962 | 34, 59, 62, 63, 64, 66, 68 |
| Total | 237413 | all |
Čítanie (ROOT, C++):
union u1 { int ib [1000]; float b[1000]; } equiv ;
union u2 { int N; float AN;} len, lenw, lentr;
int *pN = &len.N ;<br>
int *pNtr = &lentr.N ;<br>
int *pNw = &lenw.N ;<br>
...
f.read((char*)pN, 4); // record length in bytes !!!
if (f.eof()) { break; }
f.read((char*)pNw, 4); // record length in 4-byte words !!!
f.read((char*)equiv.b, 4*lenw.N);
f.read((char*)pNtr, 4);
Výsledky
Vsetky tri MINI formaty fileov citane v jednom chode: “testRead2Wdp.C”
| File | Events read in | Accepted d beam | Inconsistency | Proton beam |
|---|---|---|---|---|
| all2pr.dvo | 125673 | 121085 | 0 | 4588 |
| all4pr.dvo | 8962 | 8962 | 0 | 0 |
| allppn.dvo | 102778 | 102778 | 0 | 0 |
| Total | 237413 | 232825 | 0 | 4588 |
Žiadne prípady s inkozistentnou štruktúrou zápisu.
Entries = 232825, suma vah = 220190.80 (rozdelenie Beam Lab Momentum)
Pre MINI format pod “event” alebo “pripad” treba chapat hypotezu. Totiz na DST kazdy zaznam obsahuje iba JEDNU hypotezu.
| Hodnota vahy | 0.33 | 0.5 | 1.0 | Suma | Pocet hypotez |
| Suma vah | 2289.81 | 7979 | 209928 | 220196.81 | 232825 |
| Exp | 4 | 32 | 33 | 34 | 59 | 62 | 63 | 64 | 66 | 68 | Suma | Pocet hypotez |
| Suma vah | 18.5 | 218.38 | 5415.87 | 32230.4 | 19273.5 | 39338 | 35873.5 | 28221 | 41055 | 18553 | 220197 | 232825 |
| # of frames | 51400 | 31105 | 65110 | 56735 | 35610 | 52594 | 31068 |
|---|
Podla zoznamu expozicii este existuju:
35 1972-12-15 / 1972-12-28 d(p) сепар. 49 1981-09-17 / 1981-09-19 d(↑) 10810 + 6525 Добавка от 1981-10-20 51 1982-07-03 / 1982-07-12 d(↑) 55800
Hyp Expo
|
-108 | -107 | -106 | -105 | -104 | -103 | -102 | -101 | 1 | 2 | 101 | 102 | 103 | 104 | 105 | Sum Run |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 4 | 18.5 | 18.5 | ||||||||||||||
| 32 | 135.12 | 51.78 | 1 | 0.83 | 20.33 | 2.99 | 3 | 3.33 | 218.38 | |||||||
| 33 | 3495.34 | 612.528 | 6.66 | 4.65 | 8.3 | 1046.5 | 93.05 | 54.59 | 94.2 | 5415.82 | ||||||
| 34 | 7933.71 | 1739.27 | 6.15 | 3.65 | 19.93 | 1794 | 188.64 | 145.54 | 235.49 | 18683.7 | 30750 | |||||
| 59 | 27.5 | 2420 | 2220.5 | 2446.5 | 102.5 | 97 | 82.5 | 83.5 | 673 | 776 | 220 | 196 | 227 | 248.5 | 8533 | 18353.5 |
| 62 | 52 | 5067.5 | 4308 | 5050.5 | 245 | 192 | 210.5 | 219 | 1019 | 1681 | 467 | 477 | 471.5 | 546 | 17661 | 37667 |
| 63 | 55 | 4553 | 4213 | 5212.5 | 122 | 98 | 106.5 | 96 | 1360 | 1741.5 | 213.5 | 208 | 303.5 | 357 | 15669.5 | 34309 |
| 64 | 58 | 4051 | 3133 | 3881.5 | 102 | 79.5 | 73.5 | 89.5 | 930.5 | 1082.5 | 168 | 165.5 | 222 | 300 | 12715.5 | 27052 |
| 66 | 77 | 6334 | 5182.5 | 6466 | 153.5 | 122 | 130.5 | 137.5 | 436 | 544.5 | 272.5 | 263 | 322.5 | 461 | 18780.5 | 39683 |
| 68 | 37 | 2727.5 | 2350.5 | 2587.5 | 89.5 | 48 | 58.5 | 45 | 571.5 | 600 | 126 | 118.5 | 162.5 | 243.5 | 8132.5 | 17898 |
| Sum Hyp | 306.5 | 25153 | 32971.7 | 28048.1 | 827.31 | 645.8 | 662 | 699.56 | 7869.33 | 6425.5 | 1751.68 | 1428 | 1912.13 | 2489.02 | 100176 | 211360 |
Pre top = 200 grand total sum of weights = 211360
Hyp Expo
|
-110 | -109 | -108 | -107 | -106 | -105 | -104 | -103 | -102 | -101 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 101 | 102 | 103 | 104 | 105 | 106 | 107 | 108 | 109 | 110 | Sum Run |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 34 | 2 | 1.83 | 2 | 0.5 | 1 | 1.5 | 37.5 | 29 | 6 | 7 | 28 | 19.5 | 1171.83 | 29.33 | 68.82 | 15.5 | 58.98 | 1480.29 | ||||||||||
| 59 | 3.5 | 21.5 | 13 | 8.5 | 13.5 | 8 | 28 | 17 | 5 | 9 | 11 | 4 | 14 | 590 | 2 | 1 | 1 | 21 | 56 | 33 | 60 | 920 | ||||||
| 62 | 2 | 1 | 1.5 | 3 | 30.5 | 25 | 18.5 | 29.5 | 18.5 | 26 | 23 | 6 | 17 | 20 | 5 | 28 | 1177.5 | 2 | 29 | 86 | 18 | 104 | 1671 | |||||
| 63 | 4.5 | 2 | 10 | 41 | 29.5 | 23 | 23.5 | 12.5 | 10.5 | 29.5 | 27 | 12 | 19 | 21 | 7 | 33 | 1010.5 | 2 | 2 | 2 | 2 | 3 | 30 | 84 | 17 | 107 | 1564.5 | |
| 64 | 16 | 3 | 4 | 7 | 30 | 16 | 26 | 15 | 20 | 32 | 11 | 6 | 16 | 13 | 5 | 26.5 | 767.5 | 1 | 3 | 1 | 1 | 24 | 48 | 14 | 63 | 1169 | ||
| 66 | 2 | 5 | 175.5 | 87.5 | 74 | 83 | 80.5 | 169 | 7 | 4 | 9 | 7 | 3 | 16 | 401 | 2 | 2 | 3 | 6 | 1 | 23.5 | 70 | 35 | 106 | 1372 | |||
| 68 | 3 | 54 | 22 | 20.5 | 9 | 24.5 | 53 | 14 | 2 | 12 | 2 | 2 | 13 | 339 | 1 | 1 | 8 | 24 | 8 | 43 | 655 | |||||||
| Sum Hyp | 24.5 | 6 | 17.5 | 64.33 | 343 | 186.5 | 171.5 | 163.5 | 163.5 | 337.5 | 136.5 | 64 | 88 | 81 | 54 | 150 | 5457.33 | 6 | 5 | 10 | 12 | 1 | 5 | 164.83 | 436.82 | 140.5 | 541.98 | 8824 |
Pre top = 400 grand total sum of weights = 8824
| Hodnota vahy | 0.33 | 0.5 | 1.0 | Suma | Pocet hypotez |
| Suma vah | 844.148 | 453.5 | 1123 | 2420.6480 | 4588 |
| Exp | 32 | 33 | 34 | Suma | Pocet hypotez |
| Suma vah | 54.24 | 679.219 | 1687.16 | 2420.6190 | 4588 |
| Exp | 32 | 33 | 34 | Suma | Pocet hypotez |
| Cislo hyp | -106 | -106 | -106 | ||
| Suma vah | 54.24 | 679.219 | 1687.16 | 2420.6190 | 4588 |
Preverit sekundarne castice, pretoze neviem co hyp -106 znamena pre proton beam. Tlac hmot castic hypotezy -106 pre proton beam. Masses of: beam = 0.93826 S1 = 0.93955 S2 = 0.93826 S3 = 0.1396; ? pp → npπ+ ? je to 1C fit neviem preco ma zaporne cislo.
Nedôslednosti:
- Vahy pripadov su 0.33 (miesto 1/3), 0.5 (OK) a 1.0 (OK) pre exp <= 34
- Vahy pripadov su 0.5 (OK) a 1.0 (OK) pre exp >= 59
- Roll number (cislo filmu) je -99.00 pre expozice 32, 33, 34 (alias 10332, 10333, 10334)!
- Exposition 4; deuteron beam; 24 hypotez cislo 1, 24 = 5.5(2) + 13 (1) vo file "all2pr.dvo", Roll number -99.0; hypotezy patria roznym eventom - rozdelenie pxlab pre beam, sec1, sec2. Hypoteza c 1 je pruzny rozptyl dp → dp 4C fit a mala by byt jednoznacna, aspon podla kriterii pre exp >= 40. Ci su tieto pripady dublovane ostatnymi expoziciami - nezistoval som, tiez neviem k akej expozicii mozu patrit. FULL DST neobsahuje ziadny exposition 4 !
Vypis : all2pr.dvo Pol flag = -1 Lab-Exp = 10204 Roll = -99 Hypo = 1 Masses: 1.87549 1.87549 0.93826
| Poradie | File | Pol fl | LabEXp | Vaha | Roll | Hyp | Beam mass | Sec1 mass | Sec2 mass |
|---|---|---|---|---|---|---|---|---|---|
| 57539 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57540 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57541 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57542 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57543 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57544 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57545 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57546 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57547 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57548 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57549 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57550 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57551 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57552 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57553 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57554 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57555 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57556 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57557 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57558 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57559 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57560 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57561 | all2pr.dvo | -1 | 10204 | 0.5 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
| 57562 | all2pr.dvo | -1 | 10204 | 1 | -99 | 1 | 1.87549 | 1.87549 | 0.93826 |
Eventy(hypotezy) idu za sebou vo file all2pr.dvo maju poradove cisla 57539, ..., 57562. Nie su duplikaty ani jedneho eventu z dp MINI DST (compare4.C).
Rolls from dp MINI DST pre exp 59 a 62
Exp 59
18 19 20 21 23 24 25 26 27 28 29 30 32 33 34 35 36 37 38 39 40 41
Exp 62
40 45 46 47 48 49
50 52 54 55 57 58 59
61 62 65 66 67 68 69
70 71 87 88 89
90 92 98 99
102 103 104 105 106 107 108 109
110 111 112 113 114 115 117 118 119
120 121 122 125 126 127 128 129
130 131 137
Filmy 36, 37, 38 3Hep exp 45 a dp exp 59 sa prekryvaju. Grind checking tychto expozicii uz 12.10.1982 bol rozbehnuty, Pokyny od G.D. Pestovej. Exp 59 bola robena 1985. july 15 - 23, takze filmy sa nijak premiesat nemohli okrem toho 3He a d beam sa vizualne krasne da odlisit .Nechapem ako vznikol problem, ale mozem zalezitost povazovat za bezpredmetnu.
Hodnoty polarizacneho flagu v tabulke (proton beam vyluceny):
| Flag | -1 | 0 | +1 | +2 | +3 | +4 | >+4 | Total weighted |
|---|---|---|---|---|---|---|---|---|
| Počet hypo | 46889 | 61989 | 59763 | 64126 | 30 | 21 | 7 | 232818 |
Exp Pol flag
|
4 | 32 | 33 | 34 | 59 | 62 | 63 | 64 | 66 | 68 | Sum |
|---|---|---|---|---|---|---|---|---|---|---|---|
| -1 | 24 | 303 | 7580 | 38982 | 0 | 0 | 0 | 0 | 0 | 0 | 46889 |
| 0 | 0 | 0 | 0 | 0 | 8897 | 20586 | 13668 | 9710 | 4 | 9124 | 61989 |
| 1 | 0 | 0 | 0 | 0 | 10144 | 17285 | 2017 | 9256 | 20599 | 462 | 59763 |
| 2 | 0 | 0 | 0 | 0 | 1190 | 3162 | 20861 | 9294 | 20571 | 9048 | 64126 |
| 3 | 0 | 0 | 0 | 0 | 0 | 1 | 5 | 19 | 5 | 0 | 30 |
| 4 | 0 | 0 | 0 | 0 | 0 | 0 | 20 | 1 | 0 | 0 | 21 |
| 5 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 4 | 0 | 4 |
| 6 | 0 | 0 | 0 | 0 | 0 | 2 | 0 | 0 | 0 | 0 | 2 |
| Sum | 24 | 303 | 7580 | 38982 | 20231 | 41036 | 36571 | 28280 | 41183 | 18634 | 232824 |
| # of frames | 51400 | 31105 | 65110 | 56735 | 35610 | 52594 | 31068 | ||||
| Footnote | Overflow Polarization flag: 1, weight = 1 | ||||||||||
- Flagy po expozíciach
- Expo < 35 flag -1, non-polarized beam, OK
- Expo < 30 Co robi na DST?
- Expo 66 only 2 flags: 1 and 2
- Expo 59, 62, 63, 64 a 68 maju flagy 0, 1, 2
- Flagy po laboratoriách
- Lab=1 pol flags: -1, 0, 1, 2
- Lab=2 pol flags: 0, 1
- Lab=3 pol flags: 0, 1, 2
- Lab=6 (Kosice) nie je na MINI DST ?
- π- mass < 0, OK
- NOFIT hypo number < -100.5, OK
Aktualne MINI DST a povodne FULL DST na CDC
Prace (2003_vvg_trudy: 125.doc):
- Z. Phys. A - Atomic Nuclei 325, 391-397(1986) resp.
- Preprint JINR E1-86-58 https://lib-extopc.kek.jp/preprints/PDF/1986/8606/8606032.pdf
"In the following analysis we use the same experimental sample as in [5], namely, 20507 events of the dp → ppn reaction."
...
[5] Siemiarczuk. T. et al., Phys.Lett. 128B, 367(1983).
"A sample of 18593 deuteron break-up dp → pnp events has been collected. The break-up events can be distinctly separated into two processes: the charge exchange reaction dp → (pp)n corresponding to the configuration when the neutron momentum is higher than that of any proton in the deuteron rest system, and the remaining charge retention events proceeding without charge exchange between the projectile and the target: dp → (pn)p [14]."
...
[14] B.S. Aladashvili et al., Nucl. Instrum. Methods 129(1975) 109. (2003_vvg_trudy: 072.doc cela procedura odhadu primesi od proton beamu)
Aktualne MINI DST tiez obsahuje 20507 = 1229(0.33) + 2000(0.5) + 17278(1) pripadov dp → ppn (presnejsie je to hypoteza c. 105) a suma vah je 18683.565 = 405.565(1/0.33) + 1000(2) + 17278(1) vo file="allppn.dvo" pre test "exposition <=34". Co sa tyka poctu teda mozeme doverovat. Este preverit proton beam ! Precitat "B.S. Aladashvili et al., Nucl. Instrum. Methods 129(1975) 109." pre klud duse. Clanok stiahnut z INSPIRE, fotokopia z 2003_vvg_trudy 072.doc sa zle cita. Viac asi neurobime, musime zit s primesou a riadit sa vahami.
Kontrola:
- sumy vah z poctu entries 1229.*0.33 + 2000.*0.5 + 17278. = 18683.570 vs 18683.565
- poctu hypotez zo sumy vah 405.565*(1./0.33) + 1000*2 + 17278 = 20506.985 vs 20507
Porovnanie štatistík na FULL a MINI DST
| Exposition | Unpolarized beam | Polarized beam | Σ | |||||||||
| 32 | 33 | 34 | 59 | 62 | 63 | 64 | 66 | 68 | ||||
| Events | FULL DST | Lab = 1 | 854 | 14543 | 16836 | 41058 | 18554 | 91845 | ||||
| MINI DST | Lab = 1 | 218.38 | 5415.87 | 32230.39 | 14081.5 | 34767 | 28692.0 | 28221.0 | 41055.0 | 18553.0 | 203234.14 | |
| Lab = 2 | 5192.0 | 5192.0 | ||||||||||
| Lab = 3 | 4571.0 | 7181.0 | 11752.0 | |||||||||
Full DST obsahuje iba prípady Lab=1 (Dubna) a naviac neobsahuje celú štatistiku prípadov. Expozície 32, 33, 63, 64 celkom chýbajú na FULL DST.
Prítomné expozície:
- 34, 62 na FULL DST je oveľa menej prípadov Lab = 1, ako na MINI,
- 59 na FULL DST je viac prípadov Lab = 1, ako na MINI,
- 66, 68 možno považovať za rovnaké na oboch DST.Rozdelenie flagov polarizácie možno považovať za totožné ako ukazuje nasledujúca tabuľka.
| Exposition | 59 | 62 | 66 | 68 | |||||
| DST | FULL | MINI | FULL | MINI | FULL | MINI | FULL | MINI | |
| Polflag | 0 | 6080 | 5893.5 | 8389 | 17534.5 | 4 | 4 | 9083 | 9082.5 |
| 1 | 7310 | 7076.5 | 7096 | 14506 | 20530 | 20529 | 459 | 459 | |
| 2 | 1152 | 1110.5 | 1350 | 2723.5 | 20515 | 20513 | 9012 | 9011.5 | |
| 3 | 0 | 0 | 1 | 1 | 5 | 5 | 0 | 0 | |
| ovflw | 0 | 1 | 0 | 2 | 0 | 4 | 0 | 0 | |
- Cisla filmov sa podla expozicii zretelne oddeluju pre polarizovane zvazky 63, 64, 66, 68 okrem expozicie 59 a 62, ked obe obsahuju film c. 40 s nasledujucimi poctami pripadov:
| Expozicia | Film | MINI DST format | FULL DST format |
|---|---|---|---|
| 59 | 40 | 936.5 | 970 |
| 62 | 40 | 504.5 | 518 |
- Su to dublovane pripady ?
- FULL DST roll 40 z exp. 59 a 62 zapisane do roznych fileov: FULLexp59roll40.dvo a FULLexp62roll40.dvo. Testy na duplikaty:
- rovnost roll, frame, ordre, meas da 362 doubles
- rovnost roll, frame, ordre, meas, polflag da 168 doubles
- rovnost roll, frame, ordre, meas, polflag, scantop da 162 doubles
- v najdenych "duplikatoch" trojice suradnic vertexu (Vx, Vy, Vz), podla vypisu, sa vyrazne lisia
- rovnost Vx, Vy, Vz neda ZIADNE doubles
- Rozne suradnice vertexov pripadov pre roll=40 s rovnakymi frame, ordre, meas, polflag a scantop v exp 59 a 62 hovoria, ze pripady nie su duplikaty
Ako k tomu doslo? Cisla expozicii sa nastavovali rucne. Pri starte exp 62 sa omylom nastavilo cislo prveho filmu 40, ked sa to zistilo, potom bolo nastavene nove cislo filmu 45. Vidiet z prilozenej tabulky a histogramu, osahujucich pocty pripadov na jednotlivych filmoch, ze 41 je posledne cislo filmu v exp 59 a po oprave exp 62 pokracuje s cislom filmu 45.
Roll Expo
|
35 | 36 | 37 | 38 | 39 | 40 | 41 | 42 | 43 | 44 | 45 | 46 | 47 | 48 | 49 | 50 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 59 | 0 | 974 | 1134 | 1013 | 1014 | 970 | 1084 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 62 | 0 | 0 | 0 | 0 | 0 | 518 | 0 | 0 | 0 | 0 | 391 | 708 | 18 | 489 | 642 | 682 |
Medzery v cislovani filmov expozicie 62 na FULL DST: 66-71,72-86,93-97,99,100-101,102-108,110-115,116,117-122,123-124,125-131,132-136 - FULL DST (iba Lab=1), ---->,72-86,93-97,->,100-101,------>,------>,116,------>,123-124,------>,132-136 - MINI DST (vsetky Lab), su ciastocne vysvetlitelne rozdelenim filmov (--->) do jednotlivych laboratorii.
MINI DST analyza kanalov [math]\displaystyle{ dp \rightarrow dp }[/math] a [math]\displaystyle{ dp \rightarrow ppn }[/math]
Pruzny rozptyl [math]\displaystyle{ dp \rightarrow dp }[/math] hyp 1 alebo 2
- MINI DST
- vsetky laboratoria a expozicie
- proton beam vyluceny
- vsetky polarizacne flagy
Velmi pomale rozptylene deuterony lab hybnost az do 0.7 GeV/c, velmi rychle odrazene protony lab hybnost az nad 2.0 GeV/c, pricom hybnost deuteronoveho zvazku je v intervale (3.0, 3.6) GeV/c ma uzke rozdelenie sirka std. dev. = 0.06 GeV/c. Kvadraty prenesenej stvorhybnosti t su tiez velmi velke, ich sklon je okolo 20 GeV**(-2) podla Glagoleva VV. Fit a*exp(b*t) v intervale -0.5<= t <= -0.25 je na obrazku. Porovnat s 4Hep pruznym rozptylom, ake hodnoty t su rozumne a kde bolo difrakcne minimum (alebo pik ?). Mal som hypotezu, ze rozptylene deuterony a odrazene protony na DST su prehodene - poradie hmotnosti na DST j d a p, cize moja hypoteza neplati. Tieto nepravidelnosti su pritomne vo vsetkych expoziciach a laboratoriach. Pravdepodobne ide o zle identifikovane protony:
- - bola horna hranica pre hybnost odrazenych protonov?
- - su este niekde grind checking kriteria pre dp oziarenia ?
Skumat dalej.
Azimutalne rozdelenia pre rozne polarizacne flagy dp elastic
Flag polarizacie -1 a 0 znamena nepolarizovany zvazok, azimutalne rozdelenia deuteronov a protonov s tymito flagmi zosumujem pre kazdu casticu a bude sluzit ako pozadie pre flagy 1 a 2. Azimutalne rozdelenia pre deuteron a proton s flagom 1 a 2 vydelim so svojim pozadim, pricom deleny a deliaci histogram bude odnormovany na 1 a k pomeru bude fitovana funkcia [math]\displaystyle{ dN/d\phi = p_{0} (1 + p_{1}cos(\phi) + p_{2}cos(\phi)) }[/math]. Vysledky su uvedne na nasledujucich obrazkoch:
Vysledky pre rozptyleny deuteron flag 2 () p1= -0.16 ± 0.03 sa zhoduju s vysledkami P1-88-6 (Glagolev trudy 139) p1 = -0.13 ± 0.032 p2 = 0.02 ± 0.03 pre Flag = 2. Pre Flag 1 p1 = 0.08 ± 0.03 polarizaciu nevidim? Kde je chyba?
Bezpionova dezintegracia deuteronu [math]\displaystyle{ dp \rightarrow ppn }[/math] hyp 105 charge retention
Cast [math]\displaystyle{ dp \rightarrow ppn }[/math] priamy kanal (charge retention) rozdelenie laboratorneho polarneho uhola vs modul hybnosti pre rychly (scattered a spectator) a pomaly (odrazeny, recoiled) proton. Nechapem v rozdeleni pomalsieho protona pasmo male uhly (0., 0.1) rad. Zeby to boli odrazene protony z [math]\displaystyle{ np \rightarrow np }[/math] charge retention, kde rychly proton je spectator? Pasmo p Scat/Spec (1.6, 1.8) GeV/c a p Rec (0.1, 1.1) GeV/c ? SVG file Rec p vs Scat/Spect p (all Lab) nemohol som uploadnut - velkost presahuje hranicu 2 MB.
Flag polarizacie -1 a 0 znamena nepolarizovany zvazok, azimutalne rozdelenia nukleonov F, M, S s tymito flagmi zosumujem pre kazdy nukleon a bude sluzit ako pozadie pre flagy 1 a 2. Azimutalne rozdelenia pre F, M, S s flagom 1 a 2 vydelim so svojim pozadim, pricom deleny a deliaci histogram bude odnormovany na 1 a k pomeru bude fitovana funkcia [math]\displaystyle{ dN/d\phi = p_{0} (1 + p_{1}cos(\phi) + p_{2}cos(\phi)) }[/math]. Vysledky su uvedne na nasledujucich obrazkoch:
Vysledky pre flag 2 (M - Recoiled v ALAB, scattered v LAB) sa zhoduju s vysledkami P1-88-6 (Glagolev trudy 139) p1 = -0.13 ± 0.032 p2 = 0.02 ± 0.03 pre Flag = 2. Pre Flag 1 polarizaciu nevidim. Kde je chyba?
Pre kontrolu som urobil vypocet azimutalnych uhlov [math]\displaystyle{ \phi = atan2(PZ*, PY*) }[/math]. Vysledok sa nemenil.
Nove fakty a okolnosti okolo polarizacie zvazku
Mail od: gabriela.martinska@upjs.sk dna 25/Apr/2024 15:39
V starom maili VV sa piše: Teper na sčjot flagov. V raznych expoziciach byla raznaja nastrojka istočnika (píšem rýchle foneticky) poljarizovannych dejtronov. Eto lechko vidno jesli postrojit raspredelenia otdeľno po expoziciam dlja flagov 0, 1, 2. My berjom lučšij nabor po poljarizacii a eto flag 1 dlja 59 exp. a flag2 dlja expozicij 63 i 64. I konečno phi=phi + pi dlja odnogo iz flagov. (potom mám aj inú informáciu že nič netreba otáčat alebo meniť. V podstate je jasné, že to nefungovalo optimálne, keď sa vyberá pre danú expozíciu určitý flag. Potom z toho vyplýva, že sa údaje flagov líšili. Berú sa zalivky 59, 63, 64, 66,68. G.
Poprosil som o original mailu od VVG
Brat:
- exp. 59 flag 1,
- exp. 63 a 64 flag 2,
- exp. 66 a 68 co???
Skusim rozdelenie azimutalneho uhla oddelene pre kazdu expoziciu a flag pre M casticu (odrazena v ALAB, rozptylena v LAB sustave) z priameho kanalu reakcie [math]\displaystyle{ dp \rightarrow ppn }[/math].Pozadie bude vytvorene z suctu azimutalneho rozdelenia M castic s flagom polarizacie = 0 zo vsetkych expozicii (59, 62, 63,64,66,68) dohromady. Rozdelenia azimutalneho uhla s oddelene pre expozicie a flagy (1,2) budu normovane na 1 a budu vydelene pozadim tiez normovanym na 1. K tymto rozdeleniam bude fitovana funkcia [math]\displaystyle{ dN/d\phi = p_{0} (1 + p_{1}cos(\phi) + p_{2}cos(\phi)) }[/math]. Ucelom je urobit pokus pochopit kvalitu polarizacie deuteronoveho zvazku.
Vysledky su uvedene na nasledujucich obrazkoch.
Predbezne :
- Exp 59 flag=1 OK, flag=2 mala statistika p1=? prijatelna
- Exp 62 vyhodit nespolahlivy zvazok p1 pre oba flagy [math]\displaystyle{ 0.117 \pm 0.019 }[/math], [math]\displaystyle{ 0.14 \pm 0.05 }[/math] rovnaky. Alebo brat obe casti ako flag=1 ?
- Exp 63 flag=2 , flag=1 mala statistika, ale p1 OK.
- Exp 64 flag=1, flag=2 rozdelenie nepochopitelne
- Exp 66 flag=1 a flag=2, oba OK
- Exp 68 flag=2, flag=1 mala statistika, ale p1 OK.
Zopakovat to aj pre pruzny kanal?
Záver
Navrhujem zverejniť dp DST MINI format s dorazom na dolezitost testovat hmotnost zvazku. Neuskodila by aj sucasna previerka identifikacnych (grind checking) kriterii:
- ale pravdepodobnost fit hypotez nie je k dispozicii,
- missing mass by sa dala dopocitat, ale komu sa chce
- moduly hybnosti sa daju testovat
Pre metodické štúdie je k dispozícii aj plná verzia DST, ktora obsahuje, iba pripady z Lab=1, Dubna.
- Pre kanaly bez neutralnej alebo s jednou neutralnou casticou :
- Pocet nabitych sekundarnych castic * 1000 + (pocet neutralnych castic=0/1)*100 + poradove cislo skupiny vystupneho kanalu ako je to v popise struktury DST.
- Pre kanaly obsahujuce dve a viac neutralnych castic :
- –(Pocet nabitych sekundarnych castic * 1000 + (pocet neutralnych systemov=1)*100 + poradove cislo skupiny vystupneho kanalu), ako je to v popise struktury DST.Pre tieto kanaly neutralny system oznacujeme [math]\displaystyle{ X^{0} }[/math] a podla zakona zachovania baryonoveho cisla a energie reakcie moze obsahovat bud :
- dva a viac [math]\displaystyle{ \pi^{0} }[/math] alebo
- neutron a niekolko [math]\displaystyle{ \pi^{0} }[/math] alebo
- dva neutrony alebo
- dva neutrony a niekolko [math]\displaystyle{ \pi^{0} }[/math]
- tri neutrony alebo
- tri neutrony a niekolko [math]\displaystyle{ \pi^{0} }[/math]
- Prahy niektorych reakcii
- dva neutrony, t.j. [math]\displaystyle{ dp\rightarrow p\pi^{+}nn }[/math] (prah pd(Lab) = 1.37 GeV/c, sme nad prahom)
- dva neutrony a tri π0 mezony t.j. [math]\displaystyle{ dp\rightarrow p\pi^{+}nn\pi^{0}\pi^{0}\pi^{0} }[/math] (prah pd(Lab) = 3.17 GeV/c, sme nad prahom)
- dva neutrony a styri π0 mezony t.j. [math]\displaystyle{ dp\rightarrow p\pi^{+}nn\pi^{0}\pi^{0}\pi^{0}\pi^{0} }[/math] (prah pd(Lab) = 3.73 GeV/c, sme pod prahom)
- tri neutrony, t.j. [math]\displaystyle{ dp\rightarrow\pi^{+}\pi^{+}nnn }[/math] (prah pd(Lab) = 2.05 GeV/c, sme nad prahom)
- tri neutrony a dva π0 mezony, t.j. [math]\displaystyle{ dp\rightarrow\pi^{+}\pi^{+}nnn\pi^{0}\pi^{0} }[/math] (prah pd(Lab) = 3.2 GeV/c, sme tesne nad prahom),
- tri neutrony a tri π0 mezony, t.j. [math]\displaystyle{ dp\rightarrow\pi^{+}\pi^{+}nnn\pi^{0}\pi^{0}\pi^{0} }[/math] (prah pd(Lab) = 3.76 GeV/c, sme pod prahom).
Počet neutrónov sa dá určiť zo zákona zachovania baryónového čísla. Otázka:
• zadať počet v kode hypotézy - bude nejednoznacne, problematicke
• alebo kod nechať nezmenený ale neutrálny systém rozpísať explicitne v kanále reakcie?
Príklad:
kod= -2101 a kanal [math]\displaystyle{ dp \rightarrow dp(x\pi^{0}) }[/math], kde x = 2, 3, ... ;
kod= -2102 a kanal [math]\displaystyle{ dp \rightarrow d\pi^{+}(n+x\pi^{0}) }[/math], kde x =1, 2, ...
kod= -2105 a kanal
[math]\displaystyle{ dp \rightarrow \pi^{+}\pi^{+} (3n+x\pi^{0}) }[/math], kde x =0, 1,...;
Ak by niečo bol by som pre tento druhý spôsob.
kodovanie
| No neutral | One neutral | Two or more neutrals | ||||
| Kanal reakcie | Kod | Kanal reakcie | Kod | Kanal reakcie | [math]\displaystyle{ \textbf X^0 }[/math] | Kod |
| [math]\displaystyle{ dp \rightarrow dp }[/math] | 2001 | [math]\displaystyle{ dp \rightarrow dp\pi^{0} }[/math] | 2101 | [math]\displaystyle{ dp \rightarrow dpX^{0} }[/math] | [math]\displaystyle{ k\pi^0, k\ge2 }[/math] | -2101 |
| [math]\displaystyle{ dp \rightarrow d\pi^{+}n }[/math] | 2102 | [math]\displaystyle{ dp \rightarrow d\pi^{+}X^{0} }[/math] | [math]\displaystyle{ n+k\pi^0, k\ge1 }[/math] | -2102 | ||
| [math]\displaystyle{ dp \rightarrow ppn }[/math] | 2103 | [math]\displaystyle{ dp \rightarrow ppX^{0} }[/math] | [math]\displaystyle{ n+k\pi^0, k\ge1 }[/math] | -2103 | ||
| [math]\displaystyle{ dp \rightarrow p\pi^{+}X^{0} }[/math] | [math]\displaystyle{ 2n+k\pi^0, k\ge0 }[/math] | -2104 | ||||
| [math]\displaystyle{ dp \rightarrow \pi^{+}\pi^{+}X^{0} }[/math] | [math]\displaystyle{ 3n+k\pi^0, k\ge0 }[/math] | -2105 | ||||
| [math]\displaystyle{ dp \rightarrow dp\pi^{+}\pi^{-} }[/math] | 4001 | [math]\displaystyle{ dp \rightarrow dp\pi^{+}\pi^{-}\pi^{0} }[/math] | 4101 | [math]\displaystyle{ dp \rightarrow dp\pi^{+}\pi^{-}X^{0} }[/math] | [math]\displaystyle{ k\pi^0, k\ge2 }[/math] | -4101 |
| [math]\displaystyle{ dp \rightarrow ppp\pi^{-} }[/math] | 4002 | [math]\displaystyle{ dp \rightarrow ppp\pi^{-}\pi^{0} }[/math] | 4102 | [math]\displaystyle{ dp \rightarrow ppp\pi^{-}X^{0} }[/math] | [math]\displaystyle{ k\pi^0, k\ge2 }[/math] | -4102 |
| [math]\displaystyle{ dp \rightarrow pp\pi^{+}\pi^{-}n }[/math] | 4103 | [math]\displaystyle{ dp \rightarrow pp\pi^{+}\pi^{-}X^{0} }[/math] | [math]\displaystyle{ n+k\pi^0, k\ge1 }[/math] | -4103 | ||
Prechod ku kodom DP MINI format interna zalezitost
- kanal [math]\displaystyle{ dp \rightarrow ppn }[/math], file: allppn
kod 2103
- 2 prong bez [math]\displaystyle{ dp \rightarrow ppn }[/math], file: all2pr
[5] = NO OF HYP = NHYP
| NHYP | Kod | Kanal reakcie | |
| 1≤ NHYP ≤2 | 2001 | [math]\displaystyle{ dp \rightarrow dp }[/math] | |
| 101≤ NHYP ≤102 | 2101 | [math]\displaystyle{ dp \rightarrow dp\pi^{0} }[/math] | |
| -102≤ NHYP ≤-101 | -2101 | [math]\displaystyle{ dp \rightarrow dpX^{0} }[/math] | [math]\displaystyle{ X^0=k\pi^{0}, k\ge2 }[/math] |
| 103≤ NHYP ≤104 | 2102 | [math]\displaystyle{ dp \rightarrow d\pi^{+}n }[/math] | |
| -104 ≤ NHYP ≤-103 | -2102 | [math]\displaystyle{ dp \rightarrow d\pi^{+}X^{0} }[/math] | [math]\displaystyle{ X^0=n+k\pi^{0}, k\ge 1 }[/math] |
| NHYP = 105 | 2103 | [math]\displaystyle{ dp \rightarrow ppn }[/math] | |
| NHYP = -105 | -2103 | [math]\displaystyle{ dp \rightarrow ppX^{0} }[/math] | [math]\displaystyle{ X^0=n+k\pi^{0}, k\ge 1 }[/math] |
| -107≤ NHYP ≤-106 | -2104 | [math]\displaystyle{ dp \rightarrow p\pi^{+}X^0 }[/math] | [math]\displaystyle{ X^0=2n+k\pi^{0}, k\ge0 }[/math] |
| NHYP = -108 | -2105 | [math]\displaystyle{ dp \rightarrow \pi^{+}\pi^{+}X^{0} }[/math] | [math]\displaystyle{ X^0=3n+k\ge0 }[/math] |
- 4 prong, file: all4pr
[5] = NO OF HYP = NHYP
| NHYP | Kod | Kanal reakcie | |
| 1≤ NHYP ≤6 | 4001 | [math]\displaystyle{ dp \rightarrow dp\pi^{+}\pi^{-} }[/math] | |
| 101≤ NHYP ≤ 106 | 4101 | [math]\displaystyle{ dp \rightarrow dp\pi^{+}\pi^{-}\pi^{0} }[/math] | |
| -106≤ NHYP ≤ -101 | -4101 | [math]\displaystyle{ dp \rightarrow dp\pi^{+}\pi^{-}X^{0} }[/math] | [math]\displaystyle{ X^0=k\pi^0, k\ge2 }[/math] |
| NHYP = 7 | 4002 | [math]\displaystyle{ dp \rightarrow ppp\pi^{-} }[/math] | |
| NHYP = 107 | 4102 | [math]\displaystyle{ dp \rightarrow ppp\pi^{-}\pi^{0} }[/math] | |
| NHYP = -107 | -4102 | [math]\displaystyle{ dp \rightarrow ppp\pi^{-}X^0 }[/math] | [math]\displaystyle{ X^0=k\pi^0, k\ge2 }[/math] |
| 108 ≤ NHYP ≤110 | 4103 | [math]\displaystyle{ dp \rightarrow pp\pi^{+}\pi^{-}n }[/math] | |
| -110 ≤ NHYP ≤-108 | -4103 | [math]\displaystyle{ dp \rightarrow pp\pi^{+}\pi^{-}X^0 }[/math] | [math]\displaystyle{ X^0=n+k\pi^{0}, k\ge1 }[/math] |
Vysvetlenie kodu:
- kod > 0
- 1. cislica pocet sekundarnych nabitych castic
- 2. cislica (jedna neutralna castica v track region na DST) :
- 0 reakcia neobsahuje neutralnu casticu
- 1 reakcia obsahuje prave jednu neutralnu casticu
- kod < 0
- 1. cislica pocet sekundarnych nabitych castic
- 2. cislica :(jeden neutralny system v track region na DST), t.j. reakcia obsahuje dve a/alebo viac neutralnych castic
Posledne 2 cislice sluzia na oznacenie kanalov reakcie.