L1GtCorrelationCondition Class Reference

#include <L1GtCorrelationCondition.h>

Inheritance diagram for L1GtCorrelationCondition:

Public Member Functions
const bool	evaluateCondition () const override
	the core function to check if the condition matches More...
const L1GtCorrelationTemplate *	gtCorrelationTemplate () const
	get / set the pointer to a L1GtCondition More...
const L1GlobalTriggerGTL *	gtGTL () const
	get / set the pointer to GTL More...
const unsigned int	gtNrBinsPhi () const
	get / set the number of phi bins More...
const L1GlobalTriggerPSB *	gtPSB () const
	get / set the pointer to PSB More...
	L1GtCorrelationCondition ()
	L1GtCorrelationCondition (const L1GtCondition *, const L1GtCondition *, const L1GtCondition *, const int, const int, const int, const int, const L1GlobalTriggerGTL *, const L1GlobalTriggerPSB *, const L1GtEtaPhiConversions *)
	from base template condition (from event setup usually) More...
	L1GtCorrelationCondition (const L1GtCorrelationCondition &)
L1GtCorrelationCondition &	operator= (const L1GtCorrelationCondition &)
void	print (std::ostream &myCout) const override
	print condition More...
void	setGtCorrelationTemplate (const L1GtCorrelationTemplate *)
void	setGtGTL (const L1GlobalTriggerGTL *)
void	setGtNrBinsPhi (const unsigned int)
void	setGtPSB (const L1GlobalTriggerPSB *)
	~L1GtCorrelationCondition () override
Public Member Functions inherited from L1GtConditionEvaluation
bool	condLastResult () const
	get the latest result for the condition More...
int	condMaxNumberObjects () const
void	evaluateConditionStoreResult ()
	call evaluateCondition and save last result More...
CombinationsInCond const &	getCombinationsInCond () const
	get all the object combinations evaluated to true in the condition More...
virtual std::string	getNumericExpression () const
	get numeric expression More...
	L1GtConditionEvaluation ()
	constructor More...
void	setCondMaxNumberObjects (int condMaxNumberObjectsValue)
void	setVerbosity (const int verbosity)
virtual	~L1GtConditionEvaluation ()
	destructor More...

Private Member Functions
void	copy (const L1GtCorrelationCondition &cp)
	copy function for copy constructor and operator= More...

Private Attributes
int	m_cond0EtaBits
int	m_cond0NrL1Objects
int	m_cond1EtaBits
int	m_cond1NrL1Objects
const L1GtCondition *	m_gtCond0
	pointer to first sub-condition More...
const L1GtCondition *	m_gtCond1
	pointer to second sub-condition More...
const L1GtCorrelationTemplate *	m_gtCorrelationTemplate
	pointer to a L1GtCorrelationTemplate More...
const L1GtEtaPhiConversions *	m_gtEtaPhiConversions
	pointer to eta and phi conversion class More...
const L1GlobalTriggerGTL *	m_gtGTL
	pointer to GTL, to be able to get the trigger objects More...
const L1GlobalTriggerPSB *	m_gtPSB
	pointer to PSB, to be able to get the trigger objects More...
bool	m_isDebugEnabled
unsigned int	m_nrBinsPhi
	number of bins for delta phi More...

Additional Inherited Members
Protected Member Functions inherited from L1GtConditionEvaluation
template<class Type1 >
const bool	checkBit (const Type1 &mask, const unsigned int bitNumber) const
	check if a bit with a given number is set in a mask More...
template<class Type1 , class Type2 >
const bool	checkThreshold (const Type1 &threshold, const Type2 &value, const bool condGEqValue) const
CombinationsInCond &	combinationsInCond () const
	get all the object combinations (to fill it...) More...
Protected Attributes inherited from L1GtConditionEvaluation
CombinationsInCond	m_combinationsInCond
	store all the object combinations evaluated to true in the condition More...
bool	m_condLastResult
	the last result of evaluateCondition() More...
int	m_condMaxNumberObjects
int	m_verbosity
	verbosity level More...

Detailed Description

Description: evaluation of a CondCorrelation condition.

Implementation: <TODO: enter implementation details>

Author

: Vasile Mihai Ghete - HEPHY Vienna

Definition at line 34 of file L1GtCorrelationCondition.h.

Constructor & Destructor Documentation

L1GtCorrelationCondition::L1GtCorrelationCondition

(

)

constructors default

Definition at line 58 of file L1GtCorrelationCondition.cc.

                                                   :
    L1GtConditionEvaluation(),
    m_isDebugEnabled(edm::isDebugEnabled())
{

    // empty

}

L1GtCorrelationCondition::L1GtCorrelationCondition	(	const L1GtCondition *	corrTemplate,
		const L1GtCondition *	cond0Condition,
		const L1GtCondition *	cond1Condition,
		const int	cond0NrL1Objects,
		const int	cond1NrL1Objects,
		const int	cond0EtaBits,
		const int	cond1EtaBits,
		const L1GlobalTriggerGTL *	ptrGTL,
		const L1GlobalTriggerPSB *	ptrPSB,
		const L1GtEtaPhiConversions *	etaPhiConversions
	)

from base template condition (from event setup usually)

Definition at line 68 of file L1GtCorrelationCondition.cc.

References L1GtConditionEvaluation::m_condMaxNumberObjects, and m_nrBinsPhi.

                                                        :
            L1GtConditionEvaluation(),
            m_gtCorrelationTemplate(static_cast<const L1GtCorrelationTemplate*>(corrTemplate)),
            m_gtCond0(cond0Condition), m_gtCond1(cond1Condition),
            m_cond0NrL1Objects(cond0NrL1Objects),
            m_cond1NrL1Objects(cond1NrL1Objects), m_cond0EtaBits(cond0EtaBits),
            m_cond1EtaBits(cond1EtaBits), m_gtGTL(ptrGTL), m_gtPSB(ptrPSB),
            m_gtEtaPhiConversions(etaPhiConversions),
            m_isDebugEnabled(edm::isDebugEnabled())
            {

    m_condMaxNumberObjects = 2; // irrelevant for correlation conditions
    m_nrBinsPhi = 0;
}

L1GtCorrelationCondition::L1GtCorrelationCondition

(

const L1GtCorrelationCondition &

cp

)

Definition at line 120 of file L1GtCorrelationCondition.cc.

References copy().

                                            :
    L1GtConditionEvaluation() {
    copy(cp);
}

L1GtCorrelationCondition::~L1GtCorrelationCondition ( )

override

Definition at line 127 of file L1GtCorrelationCondition.cc.

                                                    {

    // empty

}

Member Function Documentation

void L1GtCorrelationCondition::copy ( const L1GtCorrelationCondition &  cp )

private

copy function for copy constructor and operator=

Definition at line 90 of file L1GtCorrelationCondition.cc.

References L1GtConditionEvaluation::m_combinationsInCond, m_cond0EtaBits, m_cond0NrL1Objects, m_cond1EtaBits, m_cond1NrL1Objects, L1GtConditionEvaluation::m_condLastResult, L1GtConditionEvaluation::m_condMaxNumberObjects, m_gtCond0, m_gtCond1, m_gtCorrelationTemplate, m_gtEtaPhiConversions, m_gtGTL, m_gtPSB, m_isDebugEnabled, m_nrBinsPhi, and L1GtConditionEvaluation::m_verbosity.

Referenced by gtPSB(), L1GtCorrelationCondition(), and operator=().

                                                                      {

    m_gtCorrelationTemplate = cp.m_gtCorrelationTemplate;

    m_gtCond0 = cp.m_gtCond0;
    m_gtCond1 = cp.m_gtCond1;

    m_cond0NrL1Objects = cp.m_cond0NrL1Objects;
    m_cond1NrL1Objects = cp.m_cond1NrL1Objects;
    m_cond0EtaBits = cp.m_cond0EtaBits;
    m_cond1EtaBits = cp.m_cond1EtaBits;

    m_nrBinsPhi = cp.m_nrBinsPhi;

    m_gtCorrelationTemplate = cp.m_gtCorrelationTemplate;
    m_gtGTL = cp.m_gtGTL;
    m_gtPSB = cp.m_gtPSB;

    m_gtEtaPhiConversions = cp.m_gtEtaPhiConversions;

    m_condMaxNumberObjects = cp.m_condMaxNumberObjects;
    m_condLastResult = cp.m_condLastResult;
    m_combinationsInCond = cp.m_combinationsInCond;

    m_verbosity = cp.m_verbosity;
    m_isDebugEnabled = cp.m_isDebugEnabled;


}

const bool L1GtCorrelationCondition::evaluateCondition ( ) const

overridevirtual

the core function to check if the condition matches

Implements L1GtConditionEvaluation.

Definition at line 172 of file L1GtCorrelationCondition.cc.

References CenJet, L1GtConditionEvaluation::checkBit(), hitfit::clear(), L1GtConditionEvaluation::combinationsInCond(), L1GtCorrelationTemplate::cond0Category(), L1GtCorrelationTemplate::cond1Category(), CondCalo, CondEnergySum, L1GtConditionEvaluation::condLastResult(), CondMuon, L1GtEtaPhiConversions::convertEtaIndex(), L1GtEtaPhiConversions::convertPhiIndex(), L1GtCorrelationTemplate::correlationParameter(), L1GtCorrelationTemplate::CorrelationParameter::deltaEtaRange, L1GtCorrelationTemplate::CorrelationParameter::deltaPhiRange, ETM, L1GtConditionEvaluation::evaluateConditionStoreResult(), ForJet, L1GlobalTriggerPSB::getCandL1CenJet(), L1GlobalTriggerPSB::getCandL1ETM(), L1GlobalTriggerPSB::getCandL1ForJet(), L1GlobalTriggerPSB::getCandL1HTM(), L1GlobalTriggerPSB::getCandL1IsoEG(), L1GlobalTriggerGTL::getCandL1Mu(), L1GlobalTriggerPSB::getCandL1NoIsoEG(), L1GlobalTriggerPSB::getCandL1TauJet(), L1GtConditionEvaluation::getCombinationsInCond(), L1GtEtaPhiConversions::gtObjectNrBinsPhi(), L1GtEtaPhiConversions::gtObjectPairIndex(), hexStringToInt64(), HTM, IsoEG, l1GtObjectEnumToString(), LogTrace, m_cond0EtaBits, m_cond0NrL1Objects, m_cond1EtaBits, m_cond1NrL1Objects, m_gtCond0, m_gtCond1, m_gtCorrelationTemplate, m_gtEtaPhiConversions, m_gtGTL, m_gtPSB, m_isDebugEnabled, L1GtConditionEvaluation::m_verbosity, Mu, NoIsoEG, L1GtCondition::objectType(), L1GctEtMiss::phi(), L1GctHtMiss::phi(), L1GtEnergySumCondition::print(), L1GtMuonCondition::print(), L1GtCaloCondition::print(), findQualityFiles::size, and TauJet.

                                                             {

    // std::cout << "m_isDebugEnabled = " << m_isDebugEnabled << std::endl;
    // std::cout << "m_verbosity = " << m_verbosity << std::endl;

    bool condResult = false;
    bool reqObjResult = false;

    // number of objects in condition (it is 2, no need to retrieve from
    // condition template) and their type
    int nObjInCond = 2;
    std::vector<L1GtObject> cndObjTypeVec(nObjInCond);

    // evaluate first the two sub-conditions (Type1s)

    const L1GtConditionCategory cond0Categ = m_gtCorrelationTemplate->cond0Category();
    const L1GtConditionCategory cond1Categ = m_gtCorrelationTemplate->cond1Category();

    const L1GtMuonTemplate* corrMuon = nullptr;
    const L1GtCaloTemplate* corrCalo = nullptr;
    const L1GtEnergySumTemplate* corrEnergySum = nullptr;

    // FIXME copying is slow...
    CombinationsInCond cond0Comb;
    CombinationsInCond cond1Comb;

    switch (cond0Categ) {
        case CondMuon: {
            corrMuon = static_cast<const L1GtMuonTemplate*>(m_gtCond0);
            L1GtMuonCondition muCondition(corrMuon, m_gtGTL,
                    m_cond0NrL1Objects, m_cond0EtaBits);

            muCondition.evaluateConditionStoreResult();
            reqObjResult = muCondition.condLastResult();

            cond0Comb = (muCondition.getCombinationsInCond());
            cndObjTypeVec[0] = (corrMuon->objectType())[0];

            if (m_verbosity && m_isDebugEnabled ) {
                std::ostringstream myCout;
                muCondition.print(myCout);

                LogTrace("L1GlobalTrigger") << myCout.str() << std::endl;
            }
        }
            break;
        case CondCalo: {
            corrCalo = static_cast<const L1GtCaloTemplate*>(m_gtCond0);

            L1GtCaloCondition caloCondition(corrCalo, m_gtPSB,
                    m_cond0NrL1Objects, m_cond0NrL1Objects, m_cond0NrL1Objects,
                    m_cond0NrL1Objects, m_cond0NrL1Objects, m_cond0EtaBits);

            caloCondition.evaluateConditionStoreResult();
            reqObjResult = caloCondition.condLastResult();

            cond0Comb = (caloCondition.getCombinationsInCond());
            cndObjTypeVec[0] = (corrCalo->objectType())[0];

            if (m_verbosity && m_isDebugEnabled) {
                std::ostringstream myCout;
                caloCondition.print(myCout);

                LogTrace("L1GlobalTrigger") << myCout.str() << std::endl;
            }
        }
            break;
        case CondEnergySum: {
            corrEnergySum = static_cast<const L1GtEnergySumTemplate*>(m_gtCond0);
            L1GtEnergySumCondition eSumCondition(corrEnergySum, m_gtPSB);

            eSumCondition.evaluateConditionStoreResult();
            reqObjResult = eSumCondition.condLastResult();

            cond0Comb = (eSumCondition.getCombinationsInCond());
            cndObjTypeVec[0] = (corrEnergySum->objectType())[0];

            if (m_verbosity && m_isDebugEnabled ) {
                std::ostringstream myCout;
                eSumCondition.print(myCout);

                LogTrace("L1GlobalTrigger") << myCout.str() << std::endl;
            }
        }
            break;
        default: {
            // should not arrive here, there are no correlation conditions defined for this object
            return false;
        }
            break;
    }

    // return if first subcondition is false
    if (!reqObjResult) {
        if (m_verbosity && m_isDebugEnabled) {
            LogTrace("L1GlobalTrigger")
                    << "\n  First sub-condition false, second sub-condition not evaluated and not printed."
                    << std::endl;
        }
        return false;
    }

    // second object
    reqObjResult = false;

    switch (cond1Categ) {
        case CondMuon: {
            corrMuon = static_cast<const L1GtMuonTemplate*>(m_gtCond1);
            L1GtMuonCondition muCondition(corrMuon, m_gtGTL,
                    m_cond1NrL1Objects, m_cond1EtaBits);

            muCondition.evaluateConditionStoreResult();
            reqObjResult = muCondition.condLastResult();

            cond1Comb = (muCondition.getCombinationsInCond());
            cndObjTypeVec[1] = (corrMuon->objectType())[0];

            if (m_verbosity && m_isDebugEnabled ) {
                std::ostringstream myCout;
                muCondition.print(myCout);

                LogTrace("L1GlobalTrigger") << myCout.str() << std::endl;
            }
        }
            break;
        case CondCalo: {
            corrCalo = static_cast<const L1GtCaloTemplate*>(m_gtCond1);
            L1GtCaloCondition caloCondition(corrCalo, m_gtPSB,
                    m_cond1NrL1Objects, m_cond1NrL1Objects, m_cond1NrL1Objects,
                    m_cond1NrL1Objects, m_cond1NrL1Objects, m_cond1EtaBits);

            caloCondition.evaluateConditionStoreResult();
            reqObjResult = caloCondition.condLastResult();

            cond1Comb = (caloCondition.getCombinationsInCond());
            cndObjTypeVec[1] = (corrCalo->objectType())[0];

            if (m_verbosity && m_isDebugEnabled) {
                std::ostringstream myCout;
                caloCondition.print(myCout);

                LogTrace("L1GlobalTrigger") << myCout.str() << std::endl;
            }

        }
            break;
        case CondEnergySum: {
            corrEnergySum = static_cast<const L1GtEnergySumTemplate*>(m_gtCond1);
            L1GtEnergySumCondition eSumCondition(corrEnergySum, m_gtPSB);
            eSumCondition.evaluateConditionStoreResult();
            reqObjResult = eSumCondition.condLastResult();

            cond1Comb = (eSumCondition.getCombinationsInCond());
            cndObjTypeVec[1] = (corrEnergySum->objectType())[0];

            if (m_verbosity && m_isDebugEnabled ) {
                std::ostringstream myCout;
                eSumCondition.print(myCout);

                LogTrace("L1GlobalTrigger") << myCout.str() << std::endl;
            }
        }
            break;
        default: {
            // should not arrive here, there are no correlation conditions defined for this object
            return false;
        }
            break;
    }

    // return if second sub-condition is false
    if (!reqObjResult) {
        return false;
    } else {
        LogTrace("L1GlobalTrigger") << "\n"
                << "    Both sub-conditions true for object requirements."
                << "    Evaluate correlation requirements.\n" << std::endl;

    }

    //
    // evaluate the delta_eta and delta_phi correlations
    // if here, the object requirements are satisfied for both sub-conditions
    //

    // get the correlation parameters

    L1GtCorrelationTemplate::CorrelationParameter corrPar =
        *(m_gtCorrelationTemplate->correlationParameter());

    // convert the eta template requirements from string to 64-bit integers
    // number of 64-bit integers: string length / 16
    size_t deltaEtaRangeConvSize = (corrPar.deltaEtaRange).size() / 16
            + 1;
    std::vector<unsigned long long> deltaEtaRangeConv(
            deltaEtaRangeConvSize);

    if (!(hexStringToInt64(corrPar.deltaEtaRange, deltaEtaRangeConv))) {
        return false;
    }

    // convert the phi template requirements from string to 64-bit integers
    // number of 64-bit integers: string length / 16
    size_t deltaPhiRangeConvSize = (corrPar.deltaPhiRange).size() / 16
            + 1;
    std::vector<unsigned long long> deltaPhiRangeConv(
            deltaPhiRangeConvSize);

    if (!(hexStringToInt64(corrPar.deltaPhiRange, deltaPhiRangeConv))) {
        return false;
    }

    // get the index of L1 GT object pair
    unsigned int objPairIndex = (m_gtEtaPhiConversions->gtObjectPairIndex(
            cndObjTypeVec[0], cndObjTypeVec[1]));

    // get the maximum number of bins for the delta phi scales
    unsigned int corrParDeltaPhiNrBins =
            (m_gtEtaPhiConversions->gtObjectNrBinsPhi(objPairIndex)) / 2 + 1;


    // vector to store the indices of the calorimeter objects
    // from the combination evaluated in the condition
    SingleCombInCond objectsInComb;
    objectsInComb.reserve(nObjInCond);

    // clear the m_combinationsInCond vector
    (combinationsInCond()).clear();

    // pointers to objects
    const std::vector<const L1MuGMTCand*>* candMuVec = nullptr;
    const std::vector<const L1GctCand*>* candCaloVec = nullptr;
    //    only ETM and HTM  can appear in correlation conditions
    const L1GctEtMiss* candETM = nullptr;
    const L1GctHtMiss* candHTM = nullptr;

    // make the conversions of the indices, depending on the combination of objects involved
    // (via pair index)

    unsigned int phiIndex0 = 0;
    unsigned int phiIndex1 = 0;

    unsigned int phiIndex0Converted = 0;
    unsigned int phiIndex1Converted = 0;

    unsigned int etaIndex0 = 0;
    unsigned int etaIndex1 = 0;

    unsigned int etaIndex0Converted = 0;
    unsigned int etaIndex1Converted = 0;

    LogTrace("L1GlobalTrigger")
            << "  Sub-condition 0: std::vector<SingleCombInCond> size: "
            << (cond0Comb.size()) << std::endl;
    LogTrace("L1GlobalTrigger")
            << "  Sub-condition 1: std::vector<SingleCombInCond> size: "
            << (cond1Comb.size()) << std::endl;


    // loop over all combinations which produced individually "true" as Type1s
    for (std::vector<SingleCombInCond>::const_iterator it0Comb =
            cond0Comb.begin(); it0Comb != cond0Comb.end(); it0Comb++) {

        // Type1s: there is 1 object only, no need for a loop, index 0 should be OK in (*it0Comb)[0]
        // ... but add protection to not crash
        int obj0Index = -1;

        if (!(*it0Comb).empty()) {
            obj0Index = (*it0Comb)[0];
        } else {
            LogTrace("L1GlobalTrigger")
                    << "\n  SingleCombInCond (*it0Comb).size() "
                    << ((*it0Comb).size()) << std::endl;
            return false;
        }

        bool convStatus = false;

        switch (cond0Categ) {
            case CondMuon: {
                candMuVec = m_gtGTL->getCandL1Mu();
                phiIndex0 = (*candMuVec)[obj0Index]->phiIndex();
                etaIndex0 = (*candMuVec)[obj0Index]->etaIndex();

                convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                        objPairIndex,  0, phiIndex0, phiIndex0Converted);

                if (!convStatus) {
                    // conversion failed, message written in L1GtEtaPhiConversions
                    return false;
                }

                convStatus = m_gtEtaPhiConversions->convertEtaIndex(
                        Mu, etaIndex0, etaIndex0Converted);

                if (!convStatus) {
                    // conversion failed, message written in L1GtEtaPhiConversions
                    return false;
                }

            }
                break;
            case CondCalo: {
                switch (cndObjTypeVec[0]) {
                    case NoIsoEG: {
                        candCaloVec = m_gtPSB->getCandL1NoIsoEG();
                    }
                        break;
                    case IsoEG: {
                        candCaloVec = m_gtPSB->getCandL1IsoEG();
                    }
                        break;
                    case CenJet: {
                        candCaloVec = m_gtPSB->getCandL1CenJet();
                    }
                        break;
                    case ForJet: {
                        candCaloVec = m_gtPSB->getCandL1ForJet();
                    }
                        break;
                    case TauJet: {
                        candCaloVec = m_gtPSB->getCandL1TauJet();
                    }
                        break;
                    default: {
                        // do nothing
                    }
                        break;
                }

                phiIndex0 = (*candCaloVec)[obj0Index]->phiIndex();
                etaIndex0 = (*candCaloVec)[obj0Index]->etaIndex();

                convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                        objPairIndex,  0, phiIndex0, phiIndex0Converted);

                if (!convStatus) {
                    // conversion failed, message written in L1GtEtaPhiConversions
                    return false;
                }

                convStatus = m_gtEtaPhiConversions->convertEtaIndex(
                        cndObjTypeVec[0], etaIndex0, etaIndex0Converted);

                if (!convStatus) {
                    // conversion failed, message written in L1GtEtaPhiConversions
                    return false;
                }
            }
                break;
            case CondEnergySum: {
                switch (cndObjTypeVec[0]) {
                    case ETM: {
                        candETM = m_gtPSB->getCandL1ETM();
                        phiIndex0 = candETM->phi();

                        convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                                objPairIndex,  0, phiIndex0, phiIndex0Converted);

                        if (!convStatus) {
                            // conversion failed, message written in L1GtEtaPhiConversions
                            return false;
                        }

                    }
                        break;
                    case HTM: {
                        candHTM = m_gtPSB->getCandL1HTM();
                        phiIndex0 = candHTM->phi();

                        convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                                objPairIndex,  0, phiIndex0, phiIndex0Converted);

                        if (!convStatus) {
                            // conversion failed, message written in L1GtEtaPhiConversions
                            return false;
                        }

                    }
                        break;
                    default:
                        // do nothing
                        break;
                }
            }
                break;
            default: {
                // should not arrive here, there are no correlation conditions defined for this object
                return false;
            }
                break;
        }

        for (std::vector<SingleCombInCond>::const_iterator it1Comb =
                cond1Comb.begin(); it1Comb != cond1Comb.end(); it1Comb++) {

            // Type1s: there is 1 object only, no need for a loop (*it1Comb)[0]
            // ... but add protection to not crash
            int obj1Index = -1;

            if (!(*it1Comb).empty()) {
                obj1Index = (*it1Comb)[0];
            } else {
                LogTrace("L1GlobalTrigger")
                        << "\n  SingleCombInCond (*it1Comb).size() "
                        << ((*it1Comb).size()) << std::endl;
                return false;
            }

            switch (cond1Categ) {
                case CondMuon: {
                    candMuVec = m_gtGTL->getCandL1Mu();
                    phiIndex1 = (*candMuVec)[obj1Index]->phiIndex();
                    etaIndex1 = (*candMuVec)[obj1Index]->etaIndex();

                    convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                            objPairIndex,  1, phiIndex1, phiIndex1Converted);

                    if (!convStatus) {
                        // conversion failed, message written in L1GtEtaPhiConversions
                        return false;
                    }

                    convStatus = m_gtEtaPhiConversions->convertEtaIndex(
                            cndObjTypeVec[1], etaIndex1, etaIndex1Converted);

                    if (!convStatus) {
                        // conversion failed, message written in L1GtEtaPhiConversions
                        return false;
                    }
                }
                    break;
                case CondCalo: {
                    switch (cndObjTypeVec[1]) {
                        case NoIsoEG:
                            candCaloVec = m_gtPSB->getCandL1NoIsoEG();
                            break;
                        case IsoEG:
                            candCaloVec = m_gtPSB->getCandL1IsoEG();
                            break;
                        case CenJet:
                            candCaloVec = m_gtPSB->getCandL1CenJet();
                            break;
                        case ForJet:
                            candCaloVec = m_gtPSB->getCandL1ForJet();
                            break;
                        case TauJet:
                            candCaloVec = m_gtPSB->getCandL1TauJet();
                            break;
                        default:
                            // do nothing
                            break;
                    }

                    phiIndex1 = (*candCaloVec)[obj1Index]->phiIndex();
                    etaIndex1 = (*candCaloVec)[obj1Index]->etaIndex();

                    convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                            objPairIndex,  1, phiIndex1, phiIndex1Converted);

                    if (!convStatus) {
                        // conversion failed, message written in L1GtEtaPhiConversions
                        return false;
                    }

                    convStatus = m_gtEtaPhiConversions->convertEtaIndex(
                            cndObjTypeVec[1], etaIndex1, etaIndex1Converted);

                    if (!convStatus) {
                        // conversion failed, message written in L1GtEtaPhiConversions
                        return false;
                    }

                }
                    break;
                case CondEnergySum: {
                    switch (cndObjTypeVec[1]) {
                        case ETM: {
                            candETM = m_gtPSB->getCandL1ETM();
                            phiIndex1 = candETM->phi();

                            convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                                    objPairIndex,  1, phiIndex1, phiIndex1Converted);

                            if (!convStatus) {
                                // conversion failed, message written in L1GtEtaPhiConversions
                                return false;
                            }
                        }
                            break;
                        case HTM: {
                            candHTM = m_gtPSB->getCandL1HTM();
                            phiIndex1 = candHTM->phi();

                            convStatus = m_gtEtaPhiConversions->convertPhiIndex(
                                    objPairIndex,  1, phiIndex1, phiIndex1Converted);

                            if (!convStatus) {
                                // conversion failed, message written in L1GtEtaPhiConversions
                                return false;
                            }
                        }
                            break;
                        default:
                            // do nothing
                            break;
                    }
                }
                    break;
                default: {
                    // should not arrive here, there are no correlation conditions defined for this object
                    return false;
                }
                    break;
            }

            if (m_verbosity && m_isDebugEnabled ) {
                LogTrace("L1GlobalTrigger") << "    Correlation pair ["
                        << l1GtObjectEnumToString(cndObjTypeVec[0]) << ", "
                        << l1GtObjectEnumToString(cndObjTypeVec[1])
                        << "] with collection indices [" << obj0Index << ", "
                        << obj1Index << "] " << " has: \n      phi indices = ["
                        << phiIndex0 << ", " << phiIndex1 << "] converted to ["
                        << phiIndex0Converted << ", " << phiIndex1Converted
                        << "] \n      eta indices [" << etaIndex0 << ", "
                        << etaIndex1 << "] converted to ["
                        << etaIndex0Converted << ", " << etaIndex1Converted
                        << "] \n" << std::endl;
            }

            // the conversions were successful, need to evaluate requirements now

            bool reqEtaPhiResult = false;

            // evaluate candDeltaPhi requirements

            unsigned int candDeltaPhi;

            // calculate absolute value of candDeltaPhi
            if (phiIndex0Converted > phiIndex1Converted) {
                candDeltaPhi = phiIndex0Converted - phiIndex1Converted;
            } else {
                candDeltaPhi = phiIndex1Converted - phiIndex0Converted;
            }

            // check if candDeltaPhi > 180 (via delta_phi_maxbits)
            // delta_phi contains bits for 0..180 (0 and 180 included)
            // protect also against infinite loop...

            int nMaxLoop = 10;
            int iLoop = 0;

            while (candDeltaPhi >= corrParDeltaPhiNrBins) {

                unsigned int candDeltaPhiInitial = candDeltaPhi;

                // candDeltaPhi > 180 ==> take 360 - candDeltaPhi
                candDeltaPhi = (corrParDeltaPhiNrBins - 1) * 2 - candDeltaPhi;
                if (m_verbosity) {
                    LogTrace("L1GlobalTrigger")
                            << "    Initial candDeltaPhi = "
                            << candDeltaPhiInitial
                            << " > corrParDeltaPhiNrBins = "
                            << corrParDeltaPhiNrBins
                            << "  ==> candDeltaPhi rescaled to: "
                            << candDeltaPhi << " [ loop index " << iLoop
                            << "; breaks after " << nMaxLoop << " loops ]\n"
                            << std::endl;
                }

                iLoop++;
                if (iLoop > nMaxLoop) {
                    return false;
                }
            }


            // template requirements already converted from string to 64-bit integers
            // ...now check for each 64-bit integer against template requirements
            bool indResult = true;

            for (size_t iDeltaPhi = 0; iDeltaPhi < deltaPhiRangeConv.size(); ++iDeltaPhi) {
                if (!checkBit(deltaPhiRangeConv[iDeltaPhi], candDeltaPhi)) {
                    indResult = false;
                }
            }

            if (!indResult) {

                if (m_verbosity && m_isDebugEnabled ) {
                    LogTrace("L1GlobalTrigger") << "      object delta phi = "
                            << candDeltaPhi << " fails delta phi requirements."
                            << "\n      Pair fails correlation condition.\n"
                            << std::endl;
                }

                reqEtaPhiResult = false;
                continue;

            } else {
                if (m_verbosity && m_isDebugEnabled ) {
                    LogTrace("L1GlobalTrigger") << "      object delta phi = "
                            << candDeltaPhi
                            << " passes delta phi requirements."
                            << std::endl;
                }
            }


            // evaluate candDeltaEta requirements

            unsigned int candDeltaEta;

            // calculate absolute value of candDeltaEta
            if (etaIndex0Converted > etaIndex1Converted) {
                candDeltaEta = etaIndex0Converted - etaIndex1Converted;
            } else {
                candDeltaEta = etaIndex1Converted - etaIndex0Converted;
            }

            // template requirements already converted from string to 64-bit integers
            // ...now check for each 64-bit integer against template requirements
            indResult = true;

            for (size_t iDeltaEta = 0; iDeltaEta < deltaEtaRangeConv.size(); ++iDeltaEta) {
                if (!checkBit(deltaEtaRangeConv[iDeltaEta], candDeltaEta)) {
                    indResult = false;
                }
            }

            if (!indResult) {

                if (m_verbosity && m_isDebugEnabled ) {
                    LogTrace("L1GlobalTrigger") << "      object delta eta = "
                            << candDeltaEta << " fails delta eta requirements."
                            << "\n      Pair fails correlation condition.\n"
                            << std::endl;
                }

                reqEtaPhiResult = false;
                continue;

            } else {
                if (m_verbosity && m_isDebugEnabled ) {
                    LogTrace("L1GlobalTrigger") << "      object delta eta = "
                            << candDeltaEta
                            << " passes delta eta requirements."
                            << "\n      Pair passes correlation condition.\n"
                            << std::endl;
                }

                reqEtaPhiResult = true;

            }


            // clear the indices in the combination
            objectsInComb.clear();

            objectsInComb.push_back(obj0Index);
            objectsInComb.push_back(obj1Index);

            // if we get here all checks were successful for this combination
            // set the general result for evaluateCondition to "true"

            if (reqEtaPhiResult) {

                condResult = true;
                (combinationsInCond()).push_back(objectsInComb);

            }

        }

    }

    if (m_verbosity && m_isDebugEnabled && condResult) {
        LogTrace("L1GlobalTrigger") << (combinationsInCond()).size()
                << " correlation pair(s) [" << l1GtObjectEnumToString(
                cndObjTypeVec[0]) << ", " << l1GtObjectEnumToString(
                cndObjTypeVec[1]) << "] pass(es) the correlation condition.\n"
                << std::endl;
    }

    return condResult;

}

const L1GtCorrelationTemplate* L1GtCorrelationCondition::gtCorrelationTemplate ( ) const

inline

get / set the pointer to a L1GtCondition

Definition at line 76 of file L1GtCorrelationCondition.h.

References m_gtCorrelationTemplate, and setGtCorrelationTemplate().

                                                                        {
        return m_gtCorrelationTemplate;
    }

const L1GlobalTriggerGTL* L1GtCorrelationCondition::gtGTL ( ) const

inline

get / set the pointer to GTL

Definition at line 83 of file L1GtCorrelationCondition.h.

References m_gtGTL, and setGtGTL().

                                                   {
        return m_gtGTL;
    }

const unsigned int L1GtCorrelationCondition::gtNrBinsPhi ( ) const

inline

get / set the number of phi bins

Definition at line 69 of file L1GtCorrelationCondition.h.

References m_nrBinsPhi, and setGtNrBinsPhi().

                                                  {
        return m_nrBinsPhi;
    }

const L1GlobalTriggerPSB* L1GtCorrelationCondition::gtPSB ( ) const

inline

get / set the pointer to PSB

Definition at line 90 of file L1GtCorrelationCondition.h.

References copy(), SimDataFormats::CaloAnalysis::cp, m_gtPSB, and setGtPSB().

                                                   {
        return m_gtPSB;
    }

L1GtCorrelationCondition & L1GtCorrelationCondition::operator=

(

const L1GtCorrelationCondition &

cp

)

Definition at line 134 of file L1GtCorrelationCondition.cc.

References copy().

{
    copy(cp);
    return *this;
}

void L1GtCorrelationCondition::print ( std::ostream &  myCout ) const

overridevirtual

print condition

Reimplemented from L1GtConditionEvaluation.

Definition at line 861 of file L1GtCorrelationCondition.cc.

References m_gtCorrelationTemplate, L1GtConditionEvaluation::print(), and L1GtCorrelationTemplate::print().

Referenced by L1GlobalTriggerGTL::run().

                                                             {

    m_gtCorrelationTemplate->print(myCout);
    L1GtConditionEvaluation::print(myCout);

}

void L1GtCorrelationCondition::setGtCorrelationTemplate

(

const L1GtCorrelationTemplate *

corrTempl

)

Definition at line 150 of file L1GtCorrelationCondition.cc.

References m_gtCorrelationTemplate.

Referenced by gtCorrelationTemplate().

                                                  {

    m_gtCorrelationTemplate = corrTempl;

}

void L1GtCorrelationCondition::setGtGTL

(

const L1GlobalTriggerGTL *

ptrGTL

)

Definition at line 158 of file L1GtCorrelationCondition.cc.

References m_gtGTL.

Referenced by gtGTL().

                                                                        {

    m_gtGTL = ptrGTL;

}

void L1GtCorrelationCondition::setGtNrBinsPhi

(

const unsigned int

nrBins

)

Definition at line 143 of file L1GtCorrelationCondition.cc.

References m_nrBinsPhi.

Referenced by gtNrBinsPhi().

                                                                       {

    m_nrBinsPhi = nrBins;

}

void L1GtCorrelationCondition::setGtPSB

(

const L1GlobalTriggerPSB *

ptrPSB

)

Definition at line 165 of file L1GtCorrelationCondition.cc.

References m_gtPSB.

Referenced by gtPSB().

                                                                        {

    m_gtPSB = ptrPSB;

}

Member Data Documentation

int L1GtCorrelationCondition::m_cond0EtaBits

private

Definition at line 115 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

int L1GtCorrelationCondition::m_cond0NrL1Objects

private

Definition at line 113 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

int L1GtCorrelationCondition::m_cond1EtaBits

private

Definition at line 116 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

int L1GtCorrelationCondition::m_cond1NrL1Objects

private

Definition at line 114 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

const L1GtCondition* L1GtCorrelationCondition::m_gtCond0

private

pointer to first sub-condition

Definition at line 107 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

const L1GtCondition* L1GtCorrelationCondition::m_gtCond1

private

pointer to second sub-condition

Definition at line 110 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

const L1GtCorrelationTemplate* L1GtCorrelationCondition::m_gtCorrelationTemplate

private

pointer to a L1GtCorrelationTemplate

Definition at line 104 of file L1GtCorrelationCondition.h.

Referenced by copy(), evaluateCondition(), gtCorrelationTemplate(), print(), and setGtCorrelationTemplate().

const L1GtEtaPhiConversions* L1GtCorrelationCondition::m_gtEtaPhiConversions

private

pointer to eta and phi conversion class

Definition at line 128 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

const L1GlobalTriggerGTL* L1GtCorrelationCondition::m_gtGTL

private

pointer to GTL, to be able to get the trigger objects

Definition at line 122 of file L1GtCorrelationCondition.h.

Referenced by copy(), evaluateCondition(), gtGTL(), and setGtGTL().

const L1GlobalTriggerPSB* L1GtCorrelationCondition::m_gtPSB

private

pointer to PSB, to be able to get the trigger objects

Definition at line 125 of file L1GtCorrelationCondition.h.

Referenced by copy(), evaluateCondition(), gtPSB(), and setGtPSB().

bool L1GtCorrelationCondition::m_isDebugEnabled

private

Definition at line 132 of file L1GtCorrelationCondition.h.

Referenced by copy(), and evaluateCondition().

unsigned int L1GtCorrelationCondition::m_nrBinsPhi

private

number of bins for delta phi

Definition at line 119 of file L1GtCorrelationCondition.h.

Referenced by copy(), gtNrBinsPhi(), L1GtCorrelationCondition(), and setGtNrBinsPhi().