HelpersHBEFA4.h

Go to the documentation of this file.

/****************************************************************************/
// Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.dev/sumo
// Copyright (C) 2001-2025 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
/****************************************************************************/
// Helper methods for HBEFA4-based emission computation
/****************************************************************************/
#pragma once
#include <config.h>
 
#include <vector>
#include <limits>
#include <cmath>
#include <utils/common/StdDefs.h>
#include <utils/geom/GeomHelper.h>
#include <utils/common/SUMOVehicleClass.h>
#include "EnergyParams.h"
#include "PollutantsInterface.h"
 
 
// ===========================================================================
// class definitions
// ===========================================================================
class HelpersHBEFA4 : public PollutantsInterface::Helper {
private:
    static const int HBEFA4_BASE = 7 << 16;
 
public:
    HelpersHBEFA4();

    SUMOEmissionClass getClass(const SUMOEmissionClass base, const std::string& vClass, const std::string& fuel, const std::string& eClass, const double weight) const;

    std::string getAmitranVehicleClass(const SUMOEmissionClass c) const;

    std::string getFuel(const SUMOEmissionClass c) const;

    int getEuroClass(const SUMOEmissionClass c) const;
 

    inline double compute(const SUMOEmissionClass c, const PollutantsInterface::EmissionType e, const double v, const double a, const double slope, const EnergyParams* param) const {
        if (param != nullptr && param->isEngineOff()) {
            return 0.;
        }
        const std::string& fuel = getFuel(c);
        if (fuel != "Electricity" && v > ZERO_SPEED_ACCURACY && a < getCoastingDecel(c, v, a, slope, param)) {
            return 0.;
        }
        const int index = (c & ~PollutantsInterface::HEAVY_BIT) - HBEFA4_BASE;
        double scale = 1.;
        if (e == PollutantsInterface::FUEL && myVolumetricFuel) {
            if (fuel == "Diesel") {
                scale *= 836.;
            } else if (fuel == "Gasoline") {
                scale *= 742.;
            }
        }
        const double a2 = a + GRAVITY * sin(DEG2RAD(slope));
        const double* f = myFunctionParameter[index][e];
        double result = (f[0] + f[1] * v + f[2] * a2 + f[3] * v * v + f[4] * v * v * v + f[5] * a2 * v + f[6] * a2 * v * v) / scale;
        if (e != PollutantsInterface::ELEC) {
            // no negative emissions
            result = MAX2(0.0, result);
        }
        return result;
    }
 
 
private:
    static double myFunctionParameter[833][7][7];
 
};