#!/usr/bin/env python

# kinvarbuilder - A library for searching kinematic variables in a systematic way
#
# Copyright 2014 University of California, San Diego
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#   http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.


from .VectorDifferenceQuantity import VectorDifferenceQuantity
from ..kinvarbuilder import IllegalArgumentTypes

import math

class DeltaR(VectorDifferenceQuantity):
    """distance in (eta,phi) plane between vectors """

    def __init__(self, vector1, vector2):

        if not hasattr(vector1, 'Eta') or not hasattr(vector2,'Eta'):
            raise IllegalArgumentTypes()

        # no requirement on the input vectors to be
        # fourvectors
        VectorDifferenceQuantity.__init__(self, vector1, vector2, False)

    def getValue(self):

        vecValues = [ vec.getValue() for vec in self.vectors ]

        for vecVal in vecValues:
            if vecVal == None:
                return None


        # calculate delta phi first
        phi1 = vecValues[0].Phi()
        phi2 = vecValues[1].Phi()

        dphi = phi1 - phi2

        while dphi > math.pi:
            dphi -= 2 * math.pi

        while dphi <= - math.pi:
            dphi += 2 * math.pi


        # now add the eta contribution
        deta = vecValues[0].Eta() - vecValues[1].Eta()

        return math.sqrt(deta * deta + dphi * dphi)

    def __str__(self):
        return "DeltaR(" + ", ".join(str(v) for v in self.vectors) +")"

#!/usr/bin/env python

# kinvarbuilder - A library for searching kinematic variables in a systematic way
#
# Copyright 2014 University of California, San Diego
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#   http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.


from .VectorDifferenceQuantity import VectorDifferenceQuantity
from ..kinvarbuilder import IllegalArgumentTypes

import math

class DeltaR(VectorDifferenceQuantity):
    """distance in (eta,phi) plane between vectors """

    def __init__(self, vector1, vector2):

        if not vector1.isFourVector() or not vector2.isFourVector():
            raise IllegalArgumentTypes()

        # no requirement on the input vectors to be
        # fourvectors
        VectorDifferenceQuantity.__init__(self, vector1, vector2, False)

    def getValue(self):

        vecValues = [ vec.getValue() for vec in self.vectors ]

        for vecVal in vecValues:
            if vecVal == None:
                return None


        # calculate delta phi first
        phi1 = vecValues[0].Phi()
        phi2 = vecValues[1].Phi()

        dphi = phi1 - phi2

        while dphi > math.pi:
            dphi -= 2 * math.pi

        while dphi <= - math.pi:
            dphi += 2 * math.pi


        # now add the eta contribution
        deta = vecValues[0].Eta() - vecValues[1].Eta()

        return math.sqrt(deta * deta + dphi * dphi)

    def __str__(self):
        return "DeltaR(" + ", ".join(str(v) for v in self.vectors) +")"

fix for bug introduced in 14fa744: now checking whether the given vector sums are Fourvectors instead of insisting that the given objects have an Eta(..) attribute