[Rivet] Fwd: FourMomentum mass

Sercan Sen Sercan.Sen at cern.ch
Thu Nov 11 00:26:34 GMT 2010


And, please can somebody explain why we have those values (NaN) ?

Thanks in advance
Sercan




Begin forwarded message:

> From: Sercan Sen <Sercan.Sen at cern.ch>
> Date: November 11, 2010 1:10:28 AM GMT+01:00
> To: Rivet <rivet at projects.hepforge.org>
> Subject: FourMomentum mass
>
>
> 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/247d32b4/attachment.html>


More information about the Rivet mailing list