DeltaPhi(phi1 = { 0, 0, 0, 1, 1, -0.5 },
         phi2 = { 1, -0.5, 4.14159, -0.5, 4.14159, 4.14159 })
 = { 1, -0.5, -2.14159, -1.5, 3.14159, -1.64159 }

DeltaR(eta1 = { 2.4, 2.4, 2.4, -1.5, -1.5, 1 },
       eta2 = { -1.5, 1, 0, 1, 0, 0 },
       phi1 = { 0, 0, 0, 1, 1, -0.5 },
       phi2 = { 1, -0.5, 4.14159, -0.5, 4.14159, 4.14159 })
 = { 4.02616, 1.48661, 3.21659, 2.91548, 3.48132, 1.92219 }

InvariantMass(pt1 = { 40, 20, 30 },
              eta1 = { 2.5, 0.5, -1 },
              phi1 = { -0.5, 0, 1 },
              mass1 = { 10, 10, 10 },
              pt2 = { 20, 10, 40 },
              eta2 = { 0.5, -0.5, 1 },
              phi2 = { 0, 1, -1 },
              mass2 = { 2, 2, 2 })
 = { 69.0799, 23.6971, 101.326 }

InvariantMass(pt = { 40, 20, 30 },
              eta = { 2.5, 0.5, -1 },
              phi = { -0.5, 0, 1 },
              mass = { 10, 10, 10 })
 = 220.308