TrigMuonEFInfo.cxx File Reference

#include <sstream>
#include "GaudiKernel/MsgStream.h"
#include "TrigMuonEvent/TrigMuonEFInfo.h"
#include "TrigMuonEvent/TrigMuonEFTrack.h"
#include "TrigMuonEvent/TrigMuonEFCbTrack.h"
#include "TrigMuonEvent/TrigMuonEFInfoTrackContainer.h"

Include dependency graph for TrigMuonEFInfo.cxx:

Go to the source code of this file.

Functions
std::string	str (const TrigMuonEFInfo &d)
	Helper function for printing the object. More...
MsgStream &	operator<< (MsgStream &m, const TrigMuonEFInfo &d)
	Helper operator for printing the object. More...
bool	operator== (const TrigMuonEFInfo &a, const TrigMuonEFInfo &b)
	Operator comparing two TrigMuonEFInfo objects for equality. More...
void	diff (const TrigMuonEFInfo &a, const TrigMuonEFInfo &b, std::map< std::string, double > &variableChange)
	Comparison with feedback. More...

Function Documentation

diff()

void diff	(	const TrigMuonEFInfo &	a,
		const TrigMuonEFInfo &	b,
		std::map< std::string, double > &	variableChange
	)

Comparison with feedback.

Definition at line 499 of file TrigMuonEFInfo.cxx.

                                                    {
 
    if( a.RoINum() != b.RoINum() ) {
        variableChange[ "RoINum" ] = static_cast< double >( a.RoINum() - b.RoINum() );
    }
    if( a.NSegments() != b.NSegments() ) {
        variableChange[ "NSegments" ] = static_cast< double >( a.NSegments() -
                b.NSegments() );
    }
    if( a.NMdtHits() != b.NMdtHits() ) {
        variableChange[ "NMdtHits" ] = static_cast< double >( a.NMdtHits() -
                b.NMdtHits() );
    }
    if( a.NRpcHits() != b.NRpcHits() ) {
        variableChange[ "NRpcHits" ] = static_cast< double >( a.NRpcHits() -
                b.NRpcHits() );
    }
    if( a.NCscHits() != b.NCscHits() ) {
        variableChange[ "NCscHits" ] = static_cast< double >( a.NCscHits() -
                b.NCscHits() );
    }
    if( a.NTgcHits() != b.NTgcHits() ) {
        variableChange[ "NTgcHits" ] = static_cast< double >( a.NTgcHits() -
                b.NTgcHits() );
    }
    if( a.EtaPreviousLevel() != b.EtaPreviousLevel() ) {
        variableChange[ "EtaPreviousLevel" ] = static_cast< double >( a.EtaPreviousLevel() -
                b.EtaPreviousLevel() );
    }
    if( a.PhiPreviousLevel() != b.PhiPreviousLevel() ) {
        variableChange[ "PhiPreviousLevel" ] = static_cast< double >( a.PhiPreviousLevel() -
                b.PhiPreviousLevel() );
    }
    // simple check on number of tracks in container
    if (
            (a.hasTrack() && b.hasTrack())
            &&
            (a.TrackContainer()->size() != b.TrackContainer()->size())
    ) {
        variableChange[ "NTracks" ] = static_cast< double >(a.TrackContainer()->size() -
                b.TrackContainer()->size() );
    }
    // make a view copy of container B to rm later on elements which matched to
    // elements of container A so we don't modify container B
    TrigMuonEFInfoTrackContainer *tmpContB = new TrigMuonEFInfoTrackContainer(SG::VIEW_ELEMENTS);
    *tmpContB = *(b.TrackContainer());
 
    for (
            TrigMuonEFInfoTrackContainer::const_iterator TrkItrA = a.TrackContainer()->begin();
            TrkItrA != a.TrackContainer()->end();
            ++TrkItrA
    )
    {
        // find match in containerB
        TrigMuonEFInfoTrackContainer::iterator TrkItrB = tmpContB->begin();
        while ( TrkItrB != tmpContB->end() ) {
            
            if (*(*TrkItrB)==*(*TrkItrA)) {
                // if match rm element from tmp container, diff, and break loop
                diff( *( (*TrkItrA)->SpectrometerTrack() ), *( (*TrkItrB)->SpectrometerTrack() ), variableChange );
                diff( *( (*TrkItrA)->ExtrapolatedTrack() ), *( (*TrkItrB)->ExtrapolatedTrack() ), variableChange );
                diff( *( (*TrkItrA)->CombinedTrack() ), *( (*TrkItrB)->CombinedTrack() ), variableChange );
                tmpContB->erase(TrkItrB);
                break;
            }
            ++TrkItrB;
        }
    }
    // clean up
    delete tmpContB;
 
    if (a.hasLegacyTrack() && b.hasLegacyTrack()) {
        diff( *( a.legacySpectrometerTrack() ), *( b.legacySpectrometerTrack() ), variableChange );
        diff( *( a.legacyExtrapolatedTrack() ), *( b.legacyExtrapolatedTrack() ), variableChange );
        diff( *( a.legacyCombinedTrack() ), *( b.legacyCombinedTrack() ), variableChange );
    }
 
    return;
}

operator<<()

MsgStream& operator<<	(	MsgStream &	m,
		const TrigMuonEFInfo &	d
	)

Helper operator for printing the object.

Definition at line 434 of file TrigMuonEFInfo.cxx.

                                                                {
 
    return ( m << str( d ) );
}

operator==()

bool operator==	(	const TrigMuonEFInfo &	a,
		const TrigMuonEFInfo &	b
	)

Operator comparing two TrigMuonEFInfo objects for equality.

Definition at line 439 of file TrigMuonEFInfo.cxx.

                                                                     {
 
    if ( a.RoINum() != b.RoINum() )
        return false;
    if ( a.hasTrack() != b.hasTrack() )
        return false;
    if (a.hasTrack() && b.hasTrack() ) {
        // simple check for the number of entries in container:
        if ( a.TrackContainer()->size() != b.TrackContainer()->size() )
        {
            return false;
        }
        // sort container?
        // unfortunately sort() of DataVector just sorts by pointer value,
        // so really useless here....
        // policy: make a copy of container B (tmpContB) and loop over container A
        // for every match rm element from tmpContB.
        // at the end check whether tmpContB is empty -> condition for equality
 
        // make a view copy of container B to rm later on elements which matched to
        // elements of container A so we don't modify container B
        TrigMuonEFInfoTrackContainer *tmpContB = new TrigMuonEFInfoTrackContainer(SG::VIEW_ELEMENTS);
        *tmpContB = *(b.TrackContainer());
        for (
                TrigMuonEFInfoTrackContainer::const_iterator TrkItrA = a.TrackContainer()->begin();
                TrkItrA != a.TrackContainer()->end();
                ++TrkItrA
        )
        {
            // find match in containerB
            TrigMuonEFInfoTrackContainer::iterator TrkItrB = tmpContB->begin();
            while ( TrkItrB != tmpContB->end() ) {
                if (*(*TrkItrB)==*(*TrkItrA)) {
                    // if match rm element from tmp container and break loop
                    tmpContB->erase(TrkItrB);
                    break;
                }
                ++TrkItrB;
            }
        }
 
        // if all elements were matched tmp container must be empty
        bool result = tmpContB->empty();
        // delete tmp container
        delete tmpContB;
        // check result
        if(!result) return false;
    }
    if (a.hasLegacyTrack() && b.hasLegacyTrack()) {
        if( !(*( a.legacySpectrometerTrack() ) == *( b.legacySpectrometerTrack() ) ) )
            return false;
        if( !(*( a.legacyExtrapolatedTrack() ) == *( b.legacyExtrapolatedTrack() ) ) )
            return false;
        if( !(*( a.legacyCombinedTrack() ) == *( b.legacyCombinedTrack() ) ) )
            return false;
 
    }
    return true;
}

str()

std::string str

(

const TrigMuonEFInfo &

d

)

Helper function for printing the object.

Definition at line 399 of file TrigMuonEFInfo.cxx.

                                          {
    std::stringstream ss;
    ss << "RoI index: "   << d.RoINum()
    << "; NSegments: " << d.NSegments()
    << "; NMdtHits: " << d.NMdtHits()
    << "; NRpcHits: " << d.NRpcHits()
    << "; NTgcHits: " << d.NTgcHits()
    << "; NCscHits: " << d.NCscHits()
    << "; hasTrack: " << d.hasTrack() << std::endl;
 
    if (d.hasTrack()) {
 
        // find match in containerB
                for (const TrigMuonEFInfoTrack* dd : *d.TrackContainer()) {
            if (dd) {
 
            ss << "; muon type: " << dd->MuonType();
 
            if (dd->hasSpectrometerTrack())
                ss << "; spectrometer track: " << dd->SpectrometerTrack()
                << " : " << *(dd->SpectrometerTrack());
 
            if (dd->hasExtrapolatedTrack())
                ss << "; extrapolated track: " << dd->ExtrapolatedTrack()
                << " : " << *(dd->ExtrapolatedTrack());
 
            if (dd->hasCombinedTrack())
                ss << "; combined track: " << dd->CombinedTrack()
                << " : " << *(dd->CombinedTrack());
            }
        }
    }
    return ss.str();
}