SUMO - Simulation of Urban MObility
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
MSMeanData_HBEFA.cpp
Go to the documentation of this file.
1 /****************************************************************************/
8 // Emission data collector for edges/lanes
9 /****************************************************************************/
10 // SUMO, Simulation of Urban MObility; see http://sumo-sim.org/
11 // Copyright (C) 2001-2013 DLR (http://www.dlr.de/) and contributors
12 /****************************************************************************/
13 //
14 // This file is part of SUMO.
15 // SUMO is free software: you can redistribute it and/or modify
16 // it under the terms of the GNU General Public License as published by
17 // the Free Software Foundation, either version 3 of the License, or
18 // (at your option) any later version.
19 //
20 /****************************************************************************/
21 
22 
23 // ===========================================================================
24 // included modules
25 // ===========================================================================
26 #ifdef _MSC_VER
27 #include <windows_config.h>
28 #else
29 #include <config.h>
30 #endif
31 
32 #include <microsim/MSNet.h>
33 #include <microsim/MSLane.h>
34 #include <microsim/MSVehicle.h>
35 #include <utils/common/SUMOTime.h>
36 #include <utils/common/ToString.h>
37 #include <utils/iodevices/OutputDevice.h>
38 #include "MSMeanData_HBEFA.h"
39 #include <utils/common/HelpersHBEFA.h>
40 #include <limits>
41 
42 #ifdef CHECK_MEMORY_LEAKS
43 #include <foreign/nvwa/debug_new.h>
44 #endif // CHECK_MEMORY_LEAKS
45 
46 
47 // ===========================================================================
48 // method definitions
49 // ===========================================================================
50 // ---------------------------------------------------------------------------
51 // MSMeanData_HBEFA::MSLaneMeanDataValues - methods
52 // ---------------------------------------------------------------------------
53 MSMeanData_HBEFA::MSLaneMeanDataValues::MSLaneMeanDataValues(MSLane* const lane,
54  const SUMOReal length, const bool doAdd,
55  const std::set<std::string>* const vTypes,
56  const MSMeanData_HBEFA* parent)
57  : MSMeanData::MeanDataValues(lane, length, doAdd, vTypes),
58  CO2(0), CO(0), HC(0), NOx(0), PMx(0), fuel(0), myParent(parent) {}
59 
60 
61 MSMeanData_HBEFA::MSLaneMeanDataValues::~MSLaneMeanDataValues() {
62 }
63 
64 
65 void
66 MSMeanData_HBEFA::MSLaneMeanDataValues::reset(bool) {
67  sampleSeconds = 0.;
68  travelledDistance = 0.;
69  CO2 = 0;
70  CO = 0;
71  HC = 0;
72  NOx = 0;
73  PMx = 0;
74  fuel = 0;
75 }
76 
77 
78 void
79 MSMeanData_HBEFA::MSLaneMeanDataValues::addTo(MSMeanData::MeanDataValues& val) const {
80  MSLaneMeanDataValues& v = (MSLaneMeanDataValues&) val;
81  v.sampleSeconds += sampleSeconds;
82  v.travelledDistance += travelledDistance;
83  v.CO2 += CO2;
84  v.CO += CO;
85  v.HC += HC;
86  v.NOx += NOx;
87  v.PMx += PMx;
88  v.fuel += fuel;
89 }
90 
91 
92 void
93 MSMeanData_HBEFA::MSLaneMeanDataValues::notifyMoveInternal(SUMOVehicle& veh, SUMOReal timeOnLane, SUMOReal speed) {
94  sampleSeconds += timeOnLane;
95  travelledDistance += speed * timeOnLane;
96  const double a = veh.getAcceleration();
97  CO += (timeOnLane * HelpersHBEFA::computeCO(veh.getVehicleType().getEmissionClass(), (double) speed, a));
98  CO2 += (timeOnLane * HelpersHBEFA::computeCO2(veh.getVehicleType().getEmissionClass(), (double) speed, a));
99  HC += (timeOnLane * HelpersHBEFA::computeHC(veh.getVehicleType().getEmissionClass(), (double) speed, a));
100  NOx += (timeOnLane * HelpersHBEFA::computeNOx(veh.getVehicleType().getEmissionClass(), (double) speed, a));
101  PMx += (timeOnLane * HelpersHBEFA::computePMx(veh.getVehicleType().getEmissionClass(), (double) speed, a));
102  fuel += (timeOnLane * HelpersHBEFA::computeFuel(veh.getVehicleType().getEmissionClass(), (double) speed, a));
103 }
104 
105 
106 void
107 MSMeanData_HBEFA::MSLaneMeanDataValues::write(OutputDevice& dev, const SUMOTime period,
108  const SUMOReal /*numLanes*/, const SUMOReal defaultTravelTime, const int /*numVehicles*/) const {
109  const SUMOReal normFactor = SUMOReal(3600. / STEPS2TIME(period) / myLaneLength);
110  dev.writeAttr("CO_abs", OutputDevice::realString(CO, 6)).writeAttr("CO2_abs", OutputDevice::realString(CO2, 6))
111  .writeAttr("HC_abs", OutputDevice::realString(HC, 6)).writeAttr("PMx_abs", OutputDevice::realString(PMx, 6))
112  .writeAttr("NOx_abs", OutputDevice::realString(NOx, 6)).writeAttr("fuel_abs", OutputDevice::realString(fuel, 6));
113  dev << "\n ";
114  dev.writeAttr("CO_normed", OutputDevice::realString(normFactor * CO, 6))
115  .writeAttr("CO2_normed", OutputDevice::realString(normFactor * CO2, 6))
116  .writeAttr("HC_normed", OutputDevice::realString(normFactor * HC, 6))
117  .writeAttr("PMx_normed", OutputDevice::realString(normFactor * PMx, 6))
118  .writeAttr("NOx_normed", OutputDevice::realString(normFactor * NOx, 6))
119  .writeAttr("fuel_normed", OutputDevice::realString(normFactor * fuel, 6));
120  if (sampleSeconds > myParent->myMinSamples) {
121  dev << "\n ";
122  SUMOReal vehFactor = myParent->myMaxTravelTime / sampleSeconds;
123  SUMOReal traveltime = myParent->myMaxTravelTime;
124  if (travelledDistance > 0.f) {
125  vehFactor = MIN2(vehFactor, myLaneLength / travelledDistance);
126  traveltime = MIN2(traveltime, myLaneLength * sampleSeconds / travelledDistance);
127  }
128  dev.writeAttr("traveltime", OutputDevice::realString(traveltime))
129  .writeAttr("CO_perVeh", OutputDevice::realString(CO * vehFactor, 6))
130  .writeAttr("CO2_perVeh", OutputDevice::realString(CO2 * vehFactor, 6))
131  .writeAttr("HC_perVeh", OutputDevice::realString(HC * vehFactor, 6))
132  .writeAttr("PMx_perVeh", OutputDevice::realString(PMx * vehFactor, 6))
133  .writeAttr("NOx_perVeh", OutputDevice::realString(NOx * vehFactor, 6))
134  .writeAttr("fuel_perVeh", OutputDevice::realString(fuel * vehFactor, 6));
135  } else if (defaultTravelTime >= 0.) {
136  dev << "\n ";
137  const MSVehicleType* t = MSNet::getInstance()->getVehicleControl().getVType();
138  const SUMOReal speed = MIN2(myLaneLength / defaultTravelTime, t->getMaxSpeed());
139  dev.writeAttr("traveltime", OutputDevice::realString(defaultTravelTime))
140  .writeAttr("CO_perVeh", OutputDevice::realString(HelpersHBEFA::computeDefaultCO(t->getEmissionClass(), speed, t->getCarFollowModel().getMaxAccel(), defaultTravelTime), 6))
141  .writeAttr("CO2_perVeh", OutputDevice::realString(HelpersHBEFA::computeDefaultCO2(t->getEmissionClass(), speed, t->getCarFollowModel().getMaxAccel(), defaultTravelTime), 6))
142  .writeAttr("HC_perVeh", OutputDevice::realString(HelpersHBEFA::computeDefaultHC(t->getEmissionClass(), speed, t->getCarFollowModel().getMaxAccel(), defaultTravelTime), 6))
143  .writeAttr("PMx_perVeh", OutputDevice::realString(HelpersHBEFA::computeDefaultPMx(t->getEmissionClass(), speed, t->getCarFollowModel().getMaxAccel(), defaultTravelTime), 6))
144  .writeAttr("NOx_perVeh", OutputDevice::realString(HelpersHBEFA::computeDefaultNOx(t->getEmissionClass(), speed, t->getCarFollowModel().getMaxAccel(), defaultTravelTime), 6))
145  .writeAttr("fuel_perVeh", OutputDevice::realString(HelpersHBEFA::computeDefaultFuel(t->getEmissionClass(), speed, t->getCarFollowModel().getMaxAccel(), defaultTravelTime), 6));
146  }
147  dev.closeTag();
148 }
149 
150 
151 
152 // ---------------------------------------------------------------------------
153 // MSMeanData_HBEFA - methods
154 // ---------------------------------------------------------------------------
155 MSMeanData_HBEFA::MSMeanData_HBEFA(const std::string& id,
156  const SUMOTime dumpBegin,
157  const SUMOTime dumpEnd,
158  const bool useLanes, const bool withEmpty,
159  const bool printDefaults,
160  const bool withInternal,
161  const bool trackVehicles,
162  const SUMOReal maxTravelTime,
163  const SUMOReal minSamples,
164  const std::set<std::string> vTypes)
165  : MSMeanData(id, dumpBegin, dumpEnd, useLanes, withEmpty, printDefaults,
166  withInternal, trackVehicles, maxTravelTime, minSamples, vTypes) {
167 }
168 
169 
170 MSMeanData_HBEFA::~MSMeanData_HBEFA() {}
171 
172 
173 MSMeanData::MeanDataValues*
174 MSMeanData_HBEFA::createValues(MSLane* const lane, const SUMOReal length, const bool doAdd) const {
175  return new MSLaneMeanDataValues(lane, length, doAdd, &myVehicleTypes, this);
176 }
177 
178 
179 /****************************************************************************/