Example of problem with .root file 1) The .zmc file - looks OK - balances energy-momentum Has 2 prompt photons. from Iurii's 'example_dir.zmc' diffractive sequence number 1407 EVENT 1407 0 46 3 11 0 0 0 0 5.19197315349E-16 -5.00603248091E-16 -2.75599994659E+01 2.75599994659E+01 5.09999983478E-04 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 3 2212 0 0 0 0 -6.45859234150E-16 7.73903510720E-16 9.20000000000E+02 9.20000488281E+02 9.38269972801E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 3 23 1 0 0 0 -6.03286456317E-03 -7.62987975031E-03 -8.65517044067E+00 8.65516757965E+00 -1.17457779124E-02 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 11 1 0 0 0 6.03286456317E-03 7.62987975031E-03 -1.89048290253E+01 1.89048328400E+01 5.09999983478E-04 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 3 990 2 0 0 0 -1.61640897393E-01 1.99645474553E-01 2.24015445709E+01 2.24014949799E+01 -2.61089861393E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 3 21 2 0 0 0 -1.03230707347E-01 1.71022355556E-01 1.33388175964E+01 1.33402509689E+01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 3 2 6 0 0 0 7.60643899441E-01 5.44139087200E-01 6.50979852676E+00 6.04990863800E+00 -2.57890486717E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 3 2 3 0 0 0 6.35681343079E+00 -2.82286548615E+00 1.64273941517E+00 7.14676427841E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 3 22 3 0 0 0 -5.60219240189E+00 3.35938715935E+00 -3.77375555038E+00 7.54395627975E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 2212 2 0 0 0 1.61640897393E-01 -1.99645474553E-01 8.97598388672E+02 8.97598937988E+02 9.38269972801E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 gamma 1 1 22 9 0 0 0 -5.60219240189E+00 3.35938715935E+00 -3.77375555038E+00 7.54395627975E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 gamma 2 1 22 0 0 0 0 -1.89648181200E-01 4.79288578033E-01 1.81301563978E-01 5.46401143074E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 2 8 0 16 16 6.35681343079E+00 -2.82286548615E+00 1.64273941517E+00 7.14676427841E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 21 2 0 16 16 -9.26906168461E-02 5.05607202649E-02 9.04816246033E+00 9.07581329346E+00 6.99999988079E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 -2 0 0 16 16 -6.39956057072E-01 -8.74355614185E-01 6.64792633057E+00 6.74372816086E+00 3.30000013113E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 92 13 15 17 21 5.62416696548E+00 -3.64666032791E+00 1.73388290405E+01 2.29663066864E+01 1.34865436554E+01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 223 16 0 28 30 5.51116323471E+00 -2.74158501625E+00 1.46881783009E+00 6.37641191483E+00 7.82299220562E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 211 16 0 0 0 5.39568483829E-01 1.50769919157E-01 7.49777138233E-01 9.46314573288E-01 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 -213 16 0 22 23 -5.25035262108E-02 -1.14756934345E-01 6.83805763721E-01 9.93309497833E-01 7.09329009056E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 331 16 0 31 33 -3.23027610779E-01 -5.68495213985E-01 6.61610364914E+00 6.71695327759E+00 9.57651197910E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 20223 16 0 24 25 -2.98226866871E-02 -3.56137037277E-01 7.87486219406E+00 7.99105310440E+00 1.30986714363E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 -211 19 0 0 0 -1.21374510229E-01 -4.06273156404E-02 7.10062265396E-01 7.34881341457E-01 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 111 19 0 34 35 8.31004306674E-02 -5.67029491067E-02 -9.16689187288E-02 1.91686451435E-01 1.34979993105E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 113 21 0 26 27 -1.86008289456E-01 -3.20338606834E-01 2.74508571625E+00 2.85877466202E+00 7.07022547722E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 22 21 0 0 0 1.56185597181E-01 -3.57984416187E-02 5.12977600098E+00 5.13227844238E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 211 24 0 0 0 -2.53643572330E-01 -4.89999175072E-01 1.67768323421E+00 1.77159106731E+00 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 -211 24 0 0 0 3.21946889162E-02 1.35777920485E-01 1.07832479477E+00 1.09623730183E+00 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 -211 17 0 0 0 1.44690990448E+00 -8.28238308430E-01 2.26269587874E-01 1.68825495243E+00 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 211 17 0 0 0 3.73508000374E+00 -1.77494359016E+00 1.10951149464E+00 4.28389310837E+00 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 111 17 0 36 37 3.29173475504E-01 -1.38403087854E-01 1.33036777377E-01 4.04263496399E-01 1.34979993105E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 111 20 0 38 39 -2.41565369070E-02 -1.96266710758E-01 9.61579203606E-01 9.90938067436E-01 1.34979993105E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 111 20 0 40 41 -2.31939479709E-01 -1.49094104767E-01 1.89597260952E+00 1.92066574097E+00 1.34979993105E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 221 20 0 42 44 -6.69316053391E-02 -2.23134398460E-01 3.75855207443E+00 3.80534934998E+00 5.47450006008E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 22 23 0 0 0 1.30199119449E-02 -1.50562394410E-02 -1.06914646924E-01 1.08751781285E-01 0.00000000000E+00 7.09068490323E-06 -4.83827489006E-06 -7.82180541137E-06 1.63559707289E-05 1 22 23 0 0 0 7.00805187225E-02 -4.16467078030E-02 1.52457281947E-02 8.29346701503E-02 0.00000000000E+00 7.09068490323E-06 -4.83827489006E-06 -7.82180541137E-06 1.63559707289E-05 1 22 30 0 0 0 4.94933500886E-02 -3.27142514288E-02 8.15997272730E-02 1.00887708366E-01 0.00000000000E+00 1.13019152195E-04 -4.75196247862E-05 4.56771413155E-05 1.38800722198E-04 1 22 30 0 0 0 2.79680103064E-01 -1.05688840151E-01 5.14370501041E-02 3.03375780582E-01 0.00000000000E+00 1.13019152195E-04 -4.75196247862E-05 4.56771413155E-05 1.38800722198E-04 1 22 31 0 0 0 -1.30872670561E-02 -8.73882919550E-02 1.67741939425E-01 1.89592584968E-01 0.00000000000E+00 -1.32877357828E-06 -1.07960022433E-05 5.28933851456E-05 5.45083239558E-05 1 22 31 0 0 0 -1.10692689195E-02 -1.08878418803E-01 7.93837249279E-01 8.01345467567E-01 0.00000000000E+00 -1.32877357828E-06 -1.07960022433E-05 5.28933851456E-05 5.45083239558E-05 1 22 32 0 0 0 1.48322805762E-02 -1.69956292957E-02 2.54664629698E-01 2.55661725998E-01 0.00000000000E+00 -1.14954718811E-05 -7.38945846024E-06 9.39689125516E-05 9.51927577262E-05 1 22 32 0 0 0 -2.46771752834E-01 -1.32098481059E-01 1.64130795002E+00 1.66500401497E+00 0.00000000000E+00 -1.14954718811E-05 -7.38945846024E-06 9.39689125516E-05 9.51927577262E-05 1 211 33 0 0 0 -8.25100839138E-02 4.33313846588E-02 9.24160957336E-01 9.39275681973E-01 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 -211 33 0 0 0 5.75364269316E-02 -1.40488669276E-01 9.46138978004E-01 9.68352377415E-01 1.39569997787E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 2 111 33 0 45 46 -4.19579483569E-02 -1.25977113843E-01 1.88825201988E+00 1.89772129059E+00 1.34979993105E-01 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 0.00000000000E+00 1 22 44 0 0 0 -2.88036167622E-02 -1.37545645237E-01 1.05890274048E+00 1.06818699837E+00 0.00000000000E+00 -2.93063258141E-05 -8.79911021912E-05 1.31888547912E-03 1.32549949922E-03 1 22 44 0 0 0 -1.31543315947E-02 1.15685323253E-02 8.29349279404E-01 8.29534292221E-01 0.00000000000E+00 -2.93063258141E-05 -8.79911021912E-05 1.31888547912E-03 1.32549949922E-03 sums -2.38418579E-07 6.05359674E-08 892.439819 947.560364 perfect energy balance for sum over isthep =1 note the 2 gammas indicated above ------------------------------------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------------------------------------ ------------------------------------------------------------------------------------------------------------------------------ .ROOT file output From Iurii's .root file event 1407 - the fmckin bank Same event after ORANGE. Has only one prompt photon. Note that examination of resolved events shows that the prompt photon that is selected is that with isthep = 3. For this event the prompt photons with isthep = 3 and = 1 have the same 4-vectors, but this is not always the case. We should select the prompt photon with isthep =1. EVENT 1407 1406 1 1407 UNBALANCEN - seq.num mass code 4-vector 1 23 0.00603286456 0.00762987975 -18.904829 18.9048328 2 29 -5.6021924 3.35938716 -3.77375555 7.54395628 <- gamma 1, gamma2 is missing 3 190 0.161640897 -0.199645475 897.598389 897.598938 4 54 0.539568484 0.150769919 0.749777138 0.946314573 5 55 -0.12137451 -0.0406273156 0.710062265 0.734881341 6 29 0.156185597 -0.0357984416 5.129776 5.13227844 7 54 -0.253643572 -0.489999175 1.67768323 1.77159107 8 55 0.0321946889 0.13577792 1.07832479 1.0962373 9 55 1.4469099 -0.828238308 0.226269588 1.68825495 10 54 3.73508 -1.77494359 1.10951149 4.28389311 11 29 0.0130199119 -0.0150562394 -0.106914647 0.108751781 12 29 0.0700805187 -0.0416467078 0.0152457282 0.0829346702 13 29 0.0494933501 -0.0327142514 0.0815997273 0.100887708 14 29 0.279680103 -0.10568884 0.0514370501 0.303375781 15 29 -0.0130872671 -0.087388292 0.167741939 0.189592585 16 29 -0.0110692689 -0.108878419 0.793837249 0.801345468 17 29 0.0148322806 -0.0169956293 0.25466463 0.255661726 18 29 -0.246771753 -0.132098481 1.64130795 1.66500401 19 54 -0.0825100839 0.0433313847 0.924160957 0.939275682 20 55 0.0575364269 -0.140488669 0.946138978 0.968352377 21 29 -0.0288036168 -0.137545645 1.05890274 1.068187 22 29 -0.0131543316 0.0115685323 0.829349279 0.829534292 sum 0.18964839 -0.479288518 892.258545 947.013977 no energy balance unbalance = -0.18964839 0.479288518 -0.18145752 -0.546020508 = missing gamma 2 (my reference: see mcdifd.exe.outdeb /analmc5/) email 28/11/13 8:39 to Janusz