MSRouteProbe.cpp

Go to the documentation of this file.

/****************************************************************************/
// Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.dev/sumo
// Copyright (C) 2008-2026 German Aerospace Center (DLR) and others.
// This program and the accompanying materials are made available under the
// terms of the Eclipse Public License 2.0 which is available at
// https://www.eclipse.org/legal/epl-2.0/
// This Source Code may also be made available under the following Secondary
// Licenses when the conditions for such availability set forth in the Eclipse
// Public License 2.0 are satisfied: GNU General Public License, version 2
// or later which is available at
// https://www.gnu.org/licenses/old-licenses/gpl-2.0-standalone.html
// SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
/****************************************************************************/
// Writes route distributions at a certain edge
/****************************************************************************/
#include <config.h>
 
#include <string>
#include <microsim/MSEdge.h>
#include <microsim/MSLane.h>
#include <microsim/MSGlobals.h>
#include <microsim/MSRoute.h>
#include <microsim/MSVehicle.h>
#include <utils/common/ToString.h>
#include <utils/iodevices/OutputDevice.h>
#include <mesosim/MELoop.h>
#include <mesosim/MESegment.h>
#include "MSRouteProbe.h"
 
 
// ===========================================================================
// method definitions
// ===========================================================================
MSRouteProbe::MSRouteProbe(const std::string& id, const MSEdge* edge, const std::string& distID, const std::string& lastID,
                           const std::string& vTypes) :
    MSDetectorFileOutput(id, vTypes), MSMoveReminder(id),
    myDistID(distID),
    myLastID(lastID),
    myLastRouteDistribution(nullptr),
    myCurrentRouteDistribution(nullptr),
    myEdge(edge) {
    if (MSGlobals::gUseMesoSim) {
        MESegment* seg = MSGlobals::gMesoNet->getSegmentForEdge(*edge);
        while (seg != nullptr) {
            seg->addDetector(this);
            seg = seg->getNextSegment();
        }
        return;
    }
    for (std::vector<MSLane*>::const_iterator it = edge->getLanes().begin(); it != edge->getLanes().end(); ++it) {
        (*it)->addMoveReminder(this);
    }
}
 
 
MSRouteProbe::~MSRouteProbe() {
}
 
 
void
MSRouteProbe::clearState(SUMOTime step) {
    UNUSED_PARAMETER(step);
    myCurrentRouteDistribution = nullptr;
    myLastRouteDistribution = nullptr;
}
 
 
void
MSRouteProbe::initDistributions() {
    if (myCurrentRouteDistribution == nullptr) {
        myCurrentRouteDistribution = MSRoute::distDictionary(myDistID);
        if (myCurrentRouteDistribution == 0) {
            myCurrentRouteDistribution = new RandomDistributor<ConstMSRoutePtr>();
            MSRoute::dictionary(myDistID, myCurrentRouteDistribution, false);
        }
        myLastRouteDistribution = MSRoute::distDictionary(myLastID);
    }
}
 
 
bool
MSRouteProbe::notifyEnter(SUMOTrafficObject& veh, MSMoveReminder::Notification reason, const MSLane* /* enteredLane */) {
    if (!vehicleApplies(veh)) {
        return false;
    }
    if (reason != MSMoveReminder::NOTIFICATION_SEGMENT && reason != MSMoveReminder::NOTIFICATION_LANE_CHANGE) {
        SUMOVehicle* vehicle = dynamic_cast<SUMOVehicle*>(&veh);
        if (vehicle != nullptr) {
            initDistributions();
            myCurrentRouteDistribution->add(vehicle->getRoutePtr(), 1.);
        }
    }
    return false;
}
 
 
void
MSRouteProbe::writeXMLOutput(OutputDevice& dev,
                             SUMOTime startTime, SUMOTime stopTime) {
    if (myCurrentRouteDistribution && myCurrentRouteDistribution->getOverallProb() > 0) {
        dev.openTag(SUMO_TAG_ROUTE_DISTRIBUTION).writeAttr(SUMO_ATTR_ID, getID() + "_" + time2string(startTime, false));
        const std::vector<ConstMSRoutePtr>& routes = myCurrentRouteDistribution->getVals();
        const std::vector<double>& probs = myCurrentRouteDistribution->getProbs();
        for (int j = 0; j < (int)routes.size(); ++j) {
            dev.openTag(SUMO_TAG_ROUTE).writeAttr(SUMO_ATTR_ID, routes[j]->getID() + "_" + time2string(startTime, false))
            .writeAttr(SUMO_ATTR_EDGES, toString(routes[j]->getEdges())).writeAttr(SUMO_ATTR_PROB, probs[j]).closeTag();
        }
        dev.closeTag();
        if (myLastRouteDistribution != 0) {
            MSRoute::checkDist(myLastID);
        }
        myLastRouteDistribution = myCurrentRouteDistribution;
        myLastID = myDistID;
        myDistID = getID() + "_" + toString(stopTime);
        myCurrentRouteDistribution = new RandomDistributor<ConstMSRoutePtr>();
        MSRoute::dictionary(myDistID, myCurrentRouteDistribution, false);
    }
}
 
 
void
MSRouteProbe::writeXMLDetectorProlog(OutputDevice& dev) const {
    dev.writeXMLHeader("routes", "routes_file.xsd");
}
 
 
ConstMSRoutePtr
MSRouteProbe::sampleRoute(bool last) const {
    if (myLastRouteDistribution == 0 || !last) {
        if (myCurrentRouteDistribution && myCurrentRouteDistribution->getOverallProb() > 0) {
            return myCurrentRouteDistribution->get();
        }
        return nullptr;
    }
    return myLastRouteDistribution->get();
}