[Rivet] FourMomentum mass

[Sercan Sen]

Hi,

The mass member in the FourMomentum class sometimes returns NaN, in  
particularly for the photons.
Now it's going to be a question for me that how should I consider  
those values in my analysis ?


00368                         {
00369       assert(Rivet::isZero(mass2()) || mass2() > 0);
00370       return sqrt(mass2());
00371     }


Rivet.Analysis.Handler: DEBUG  Event #100 weight = 1
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[0.93957, 0.882792, 2112, 938.54, (0.426217, 0.418899, 938.539)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[0.93957, 0.882792, 2112, 211.161, (0.214592, 0.0281939, 211.159)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[0.49765, 0.247656, 130, 11.0824, (0.516175, -0.0481894, 11.0591)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[0.49765, 0.247656, 130, 5.22922, (-0.123332, -0.202639, 5.20008)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[3.33818e-08, 1.11434e-15, 22, 1.65032, (-0.08561, -0.0331697, 1.64776)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[4.68262e-08, 2.19269e-15, 22, 4.94035, (-0.145798, 0.106414, 4.93705)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[1.60095e-08, 2.56305e-16, 22, 1.45722, (-0.246436, 0.0124048, 1.43617)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[2.52661e-08, 6.38378e-16, 22, 3.71863, (-0.441186, -0.0710671,  
3.69168)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[nan, -3.09561e-15, 22, 7.35588, (-0.145765, 0.0699382, 7.3541)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[1.13449e-07, 1.28707e-14, 22, 7.48494, (-0.0504706, -0.0235015,  
7.48474)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[nan, -5.82867e-16, 22, 0.710786, (-0.189912, 0.135686, 0.671371)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[1.20236e-05, 1.44567e-10, 22, 487.876, (0.134386, 0.0574909, 487.876)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[nan, -2.1452e-12, 22, 468.457, (-0.00154584, 0.0341641, 468.457)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[8.57925e-06, 7.36036e-11, 22, 411.51, (0.134967, -0.0703667, 411.51)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[1.04885e-07, 1.10008e-14, 22, 1.17373, (-0.00163007, 0.00672073,  
1.17371)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[2.12389e-06, 4.51093e-12, 22, 68.9775, (-0.0401148, 0.0518824,  
68.9775)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[7.09915e-06, 5.03979e-11, 22, 283.375, (-0.00347369, -0.00781676,  
283.375)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[nan, -4.27735e-14, 22, 10.7512, (-0.129995, 0.0453852, 10.7503)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[nan, -5.52943e-17, 22, 1.50234, (-0.0526513, 0.0431477, 1.5008)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[9.04115e-08, 8.17423e-15, 22, 4.85834, (-0.141602, 0.0259915, 4.85621)]
Rivet.Analysis.LHC_2010_EFWDDIFFGAP: DEBUG  [mass, mass2, id, E, p] =  
[1.29497e-07, 1.67696e-14, 22, 7.48867, (-0.262596, -0.121533, 7.48308)]
Rivet.Analysis.Handler: DEBUG  Event #101 weight = 1


Cheers
Sercan

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.hepforge.org/lists-archive/rivet/attachments/20101111/ace095a5/attachment.html>