Eclipse SUMO - Simulation of Urban MObility
MSMeanData_Emissions.cpp
Go to the documentation of this file.
1/****************************************************************************/
2// Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.org/sumo
3// Copyright (C) 2001-2022 German Aerospace Center (DLR) and others.
4// This program and the accompanying materials are made available under the
5// terms of the Eclipse Public License 2.0 which is available at
6// https://www.eclipse.org/legal/epl-2.0/
7// This Source Code may also be made available under the following Secondary
8// Licenses when the conditions for such availability set forth in the Eclipse
9// Public License 2.0 are satisfied: GNU General Public License, version 2
10// or later which is available at
11// https://www.gnu.org/licenses/old-licenses/gpl-2.0-standalone.html
12// SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
13/****************************************************************************/
19// Emission data collector for edges/lanes that
20/****************************************************************************/
21#include <config.h>
22
23#include <limits>
24#include <microsim/MSNet.h>
25#include <microsim/MSLane.h>
26#include <microsim/MSVehicle.h>
33
34
35// ===========================================================================
36// method definitions
37// ===========================================================================
38// ---------------------------------------------------------------------------
39// MSMeanData_Emissions::MSLaneMeanDataValues - methods
40// ---------------------------------------------------------------------------
42 const double length, const bool doAdd,
43 const MSMeanData_Emissions* parent)
44 : MSMeanData::MeanDataValues(lane, length, doAdd, parent),
45 myEmissions() {}
46
47
49}
50
51
52void
54 sampleSeconds = 0.;
55 travelledDistance = 0.;
56 myEmissions = PollutantsInterface::Emissions();
57}
58
59
60void
63 v.sampleSeconds += sampleSeconds;
64 v.travelledDistance += travelledDistance;
65 v.myEmissions.addScaled(myEmissions);
66}
67
68
69void
70MSMeanData_Emissions::MSLaneMeanDataValues::notifyMoveInternal(const SUMOTrafficObject& veh, const double /* frontOnLane */, const double timeOnLane, const double /*meanSpeedFrontOnLane*/, const double meanSpeedVehicleOnLane, const double /*travelledDistanceFrontOnLane*/, const double travelledDistanceVehicleOnLane, const double /* meanLengthOnLane */) {
71 if (myParent != nullptr && !myParent->vehicleApplies(veh)) {
72 return;
73 }
74 if (veh.isVehicle()) {
75 sampleSeconds += timeOnLane;
76 travelledDistance += travelledDistanceVehicleOnLane;
77 const double a = veh.getAcceleration();
79 // XXX: recheck, which value to use here for the speed. (Leo) Refs. #2579
80 meanSpeedVehicleOnLane, a, veh.getSlope(),
81 static_cast<const SUMOVehicle&>(veh).getEmissionParameters()), timeOnLane);
82 }
83}
84
85bool
87 if (veh.isVehicle()) {
89 0., 0., 0.,
90 static_cast<const SUMOVehicle&>(veh).getEmissionParameters()), TS);
91 }
92 return true;
93}
94
95
96void
97MSMeanData_Emissions::MSLaneMeanDataValues::write(OutputDevice& dev, long long int attributeMask, const SUMOTime period,
98 const double /*numLanes*/, const double /*speedLimit*/, const double defaultTravelTime, const int /*numVehicles*/) const {
99 const double normFactor = double(3600. / STEPS2TIME(period) / myLaneLength);
100 dev.writeOptionalAttr(SUMO_ATTR_CO_ABS, OutputDevice::realString(myEmissions.CO, 6), attributeMask);
101 dev.writeOptionalAttr(SUMO_ATTR_CO2_ABS, OutputDevice::realString(myEmissions.CO2, 6), attributeMask);
102 dev.writeOptionalAttr(SUMO_ATTR_HC_ABS, OutputDevice::realString(myEmissions.HC, 6), attributeMask);
103 dev.writeOptionalAttr(SUMO_ATTR_PMX_ABS, OutputDevice::realString(myEmissions.PMx, 6), attributeMask);
104 dev.writeOptionalAttr(SUMO_ATTR_NOX_ABS, OutputDevice::realString(myEmissions.NOx, 6), attributeMask);
105 dev.writeOptionalAttr(SUMO_ATTR_FUEL_ABS, OutputDevice::realString(myEmissions.fuel, 6), attributeMask);
106 dev.writeOptionalAttr(SUMO_ATTR_ELECTRICITY_ABS, OutputDevice::realString(myEmissions.electricity, 6), attributeMask);
107 if (attributeMask == 0) {
108 dev << "\n ";
109 }
110 dev.writeOptionalAttr(SUMO_ATTR_CO_NORMED, OutputDevice::realString(normFactor * myEmissions.CO, 6), attributeMask);
111 dev.writeOptionalAttr(SUMO_ATTR_CO2_NORMED, OutputDevice::realString(normFactor * myEmissions.CO2, 6), attributeMask);
112 dev.writeOptionalAttr(SUMO_ATTR_HC_NORMED, OutputDevice::realString(normFactor * myEmissions.HC, 6), attributeMask);
113 dev.writeOptionalAttr(SUMO_ATTR_PMX_NORMED, OutputDevice::realString(normFactor * myEmissions.PMx, 6), attributeMask);
114 dev.writeOptionalAttr(SUMO_ATTR_NOX_NORMED, OutputDevice::realString(normFactor * myEmissions.NOx, 6), attributeMask);
115 dev.writeOptionalAttr(SUMO_ATTR_FUEL_NORMED, OutputDevice::realString(normFactor * myEmissions.fuel, 6), attributeMask);
116 dev.writeOptionalAttr(SUMO_ATTR_ELECTRICITY_NORMED, OutputDevice::realString(normFactor * myEmissions.electricity, 6), attributeMask);
117
118 if (sampleSeconds > myParent->getMinSamples()) {
119 double vehFactor = myParent->getMaxTravelTime() / sampleSeconds;
120 double traveltime = myParent->getMaxTravelTime();
121 if (travelledDistance > 0.f) {
122 vehFactor = MIN2(vehFactor, myLaneLength / travelledDistance);
123 traveltime = MIN2(traveltime, myLaneLength * sampleSeconds / travelledDistance);
124 }
125 if (attributeMask == 0) {
126 dev << "\n ";
127 }
128 dev.writeOptionalAttr(SUMO_ATTR_TRAVELTIME, OutputDevice::realString(traveltime), attributeMask);
129 dev.writeOptionalAttr(SUMO_ATTR_CO_PERVEH, OutputDevice::realString(vehFactor * myEmissions.CO, 6), attributeMask);
130 dev.writeOptionalAttr(SUMO_ATTR_CO2_PERVEH, OutputDevice::realString(vehFactor * myEmissions.CO2, 6), attributeMask);
131 dev.writeOptionalAttr(SUMO_ATTR_HC_PERVEH, OutputDevice::realString(vehFactor * myEmissions.HC, 6), attributeMask);
132 dev.writeOptionalAttr(SUMO_ATTR_PMX_PERVEH, OutputDevice::realString(vehFactor * myEmissions.PMx, 6), attributeMask);
133 dev.writeOptionalAttr(SUMO_ATTR_NOX_PERVEH, OutputDevice::realString(vehFactor * myEmissions.NOx, 6), attributeMask);
134 dev.writeOptionalAttr(SUMO_ATTR_FUEL_PERVEH, OutputDevice::realString(vehFactor * myEmissions.fuel, 6), attributeMask);
135 dev.writeOptionalAttr(SUMO_ATTR_ELECTRICITY_PERVEH, OutputDevice::realString(vehFactor * myEmissions.electricity, 6), attributeMask);
136 } else if (defaultTravelTime >= 0.) {
138 const double speed = MIN2(myLaneLength / defaultTravelTime, t->getMaxSpeed());
139
140 if (attributeMask == 0) {
141 dev << "\n ";
142 }
143 dev.writeOptionalAttr(SUMO_ATTR_TRAVELTIME, OutputDevice::realString(defaultTravelTime), attributeMask);
151 }
152 dev.closeTag();
153}
154
155
156
157// ---------------------------------------------------------------------------
158// MSMeanData_Emissions - methods
159// ---------------------------------------------------------------------------
161 const SUMOTime dumpBegin,
162 const SUMOTime dumpEnd,
163 const bool useLanes, const bool withEmpty,
164 const bool printDefaults,
165 const bool withInternal,
166 const bool trackVehicles,
167 const double maxTravelTime,
168 const double minSamples,
169 const std::string& vTypes,
170 const std::string& writeAttributes,
171 const std::vector<MSEdge*>& edges,
172 bool aggregate) :
173 MSMeanData(id, dumpBegin, dumpEnd, useLanes, withEmpty, printDefaults,
174 withInternal, trackVehicles, 0, maxTravelTime, minSamples, vTypes, writeAttributes, edges, aggregate)
175{ }
176
177
179
180
182MSMeanData_Emissions::createValues(MSLane* const lane, const double length, const bool doAdd) const {
183 return new MSLaneMeanDataValues(lane, length, doAdd, this);
184}
185
186
187/****************************************************************************/
long long int SUMOTime
Definition: GUI.h:36
#define STEPS2TIME(x)
Definition: SUMOTime.h:54
#define TS
Definition: SUMOTime.h:41
@ SUMO_ATTR_NOX_ABS
@ SUMO_ATTR_FUEL_PERVEH
@ SUMO_ATTR_HC_NORMED
@ SUMO_ATTR_FUEL_NORMED
@ SUMO_ATTR_CO_PERVEH
@ SUMO_ATTR_ELECTRICITY_ABS
@ SUMO_ATTR_PMX_PERVEH
@ SUMO_ATTR_HC_PERVEH
@ SUMO_ATTR_CO2_NORMED
@ SUMO_ATTR_NOX_NORMED
@ SUMO_ATTR_TRAVELTIME
@ SUMO_ATTR_PMX_ABS
@ SUMO_ATTR_ELECTRICITY_PERVEH
@ SUMO_ATTR_CO_ABS
MSMeanData_Emissions.
@ SUMO_ATTR_FUEL_ABS
@ SUMO_ATTR_CO_NORMED
@ SUMO_ATTR_NOX_PERVEH
@ SUMO_ATTR_HC_ABS
@ SUMO_ATTR_CO2_PERVEH
@ SUMO_ATTR_PMX_NORMED
@ SUMO_ATTR_ELECTRICITY_NORMED
@ SUMO_ATTR_CO2_ABS
T MIN2(T a, T b)
Definition: StdDefs.h:71
double getMaxAccel() const
Get the vehicle type's maximum acceleration [m/s^2].
Definition: MSCFModel.h:254
Representation of a lane in the micro simulation.
Definition: MSLane.h:84
Data structure for mean (aggregated) edge/lane values.
Definition: MSMeanData.h:66
double travelledDistance
The sum of the distances the vehicles travelled.
Definition: MSMeanData.h:182
Data structure for mean (aggregated) edge/lane values.
void write(OutputDevice &dev, long long int attributeMask, const SUMOTime period, const double numLanes, const double speedLimit, const double defaultTravelTime, const int numVehicles=-1) const
Writes output values into the given stream.
MSLaneMeanDataValues(MSLane *const lane, const double length, const bool doAdd, const MSMeanData_Emissions *parent)
Constructor.
bool notifyIdle(SUMOTrafficObject &veh)
Computes idling emission values and adds them to the aggregate emission sums.
PollutantsInterface::Emissions myEmissions
Collected values.
void reset(bool afterWrite=false)
Resets values so they may be used for the next interval.
void addTo(MSMeanData::MeanDataValues &val) const
Add the values of this to the given one and store them there.
void notifyMoveInternal(const SUMOTrafficObject &veh, const double, const double timeOnLane, const double, const double meanSpeedVehicleOnLane, const double travelledDistanceFrontOnLane, const double travelledDistanceVehicleOnLane, const double)
Internal notification about the vehicle moves.
Emission data collector for edges/lanes.
MSMeanData::MeanDataValues * createValues(MSLane *const lane, const double length, const bool doAdd) const
Create an instance of MeanDataValues.
virtual ~MSMeanData_Emissions()
Destructor.
MSMeanData_Emissions(const std::string &id, const SUMOTime dumpBegin, const SUMOTime dumpEnd, const bool useLanes, const bool withEmpty, const bool printDefaults, const bool withInternal, const bool trackVehicles, const double minSamples, const double maxTravelTime, const std::string &vTypes, const std::string &writeAttributes, const std::vector< MSEdge * > &edges, bool aggregate)
Constructor.
Data collector for edges/lanes.
Definition: MSMeanData.h:57
static MSNet * getInstance()
Returns the pointer to the unique instance of MSNet (singleton).
Definition: MSNet.cpp:183
MSVehicleControl & getVehicleControl()
Returns the vehicle control.
Definition: MSNet.h:379
MSVehicleType * getVType(const std::string &id=DEFAULT_VTYPE_ID, SumoRNG *rng=nullptr, bool readOnly=false)
Returns the named vehicle type or a sample from the named distribution.
The car-following model and parameter.
Definition: MSVehicleType.h:63
const EnergyParams * getEmissionParameters() const
retrieve parameters for the energy consumption model
double getMaxSpeed() const
Get vehicle's (technical) maximum speed [m/s].
const MSCFModel & getCarFollowModel() const
Returns the vehicle type's car following model definition (const version)
SUMOEmissionClass getEmissionClass() const
Get this vehicle type's emission class.
Static storage of an output device and its base (abstract) implementation.
Definition: OutputDevice.h:61
static std::string realString(const double v, const int precision=gPrecision)
Helper method for string formatting.
bool closeTag(const std::string &comment="")
Closes the most recently opened tag and optionally adds a comment.
OutputDevice & writeOptionalAttr(const SumoXMLAttr attr, const T &val, long long int attributeMask)
writes a named attribute unless filtered
Definition: OutputDevice.h:268
static double computeDefault(const SUMOEmissionClass c, const EmissionType e, const double v, const double a, const double slope, const double tt, const EnergyParams *param)
Returns the amount of emitted pollutant given the vehicle type and default values for the state (in m...
static Emissions computeAll(const SUMOEmissionClass c, const double v, const double a, const double slope, const EnergyParams *param)
Returns the amount of all emitted pollutants given the vehicle type and state (in mg/s or ml/s for fu...
Representation of a vehicle, person, or container.
virtual bool isVehicle() const
Whether it is a vehicle.
virtual double getAcceleration() const =0
Returns the object's acceleration.
virtual const MSVehicleType & getVehicleType() const =0
Returns the object's "vehicle" type.
virtual double getSlope() const =0
Returns the slope of the road at object's position in degrees.
Representation of a vehicle.
Definition: SUMOVehicle.h:60
virtual EnergyParams * getEmissionParameters() const =0
Returns the vehicle's emission model parameter.
Storage for collected values of all emission types.
void addScaled(const Emissions &a, const double scale=1.)
Add the values of the other struct to this one, scaling the values if needed.