ThermalModelFitQuantities.h

Go to the documentation of this file.

/*
 * Thermal-FIST package
 * 
 * Copyright (c) 2019 Volodymyr Vovchenko
 *
 * GNU General Public License (GPLv3 or later)
 */
#ifndef THERMALMODELFITQUANTITIES_H
#define THERMALMODELFITQUANTITIES_H

#include <string>

#include "HRGBase/ParticleDecay.h"

namespace thermalfist {
  
  struct ExperimentMultiplicity {
    long long fPDGID;

    double fValue;
    
    double fError;

    Feeddown::Type fFeedDown; 

    ExperimentMultiplicity(long long PDGID = -211, double value = 300., double error = 20., Feeddown::Type fd = Feeddown::StabilityFlag) :
      fPDGID(PDGID), fValue(value), fError(error), fFeedDown(fd) { }

    void addSystematicError(double percent) {
      fError = sqrt(fError*fError + percent * percent * fValue * fValue);
    }
  };

  struct ExperimentRatio {
    long long fPDGID1;

    long long fPDGID2;

    double fValue;

    double fError;

    Feeddown::Type fFeedDown1;

    Feeddown::Type fFeedDown2; 

    ExperimentRatio(long long PDGID1 = 211, long long PDGID2 = -211, double value = 1., double error = 0.1, Feeddown::Type fd1 = Feeddown::StabilityFlag, Feeddown::Type fd2 = Feeddown::StabilityFlag) :
      fPDGID1(PDGID1), fPDGID2(PDGID2), fValue(value), fError(error), fFeedDown1(fd1), fFeedDown2(fd2) { }
    
    ExperimentRatio(long long PDGID1, long long PDGID2, double value1, double error1, double value2, double error2, Feeddown::Type fd1 = Feeddown::StabilityFlag, Feeddown::Type fd2 = Feeddown::StabilityFlag) :
      fPDGID1(PDGID1), fPDGID2(PDGID2), fFeedDown1(fd1), fFeedDown2(fd2) {
      fValue = value1 / value2;
      fError = sqrt(error1*error1 / value2 / value2 + value1 * value1 / value2 / value2 / value2 / value2 * error2 * error2);
    }
  };

  struct FittedQuantity {
    enum FittedQuantityType { 
        Multiplicity = 0, 
        Ratio = 1 
    };

    FittedQuantityType type;

    bool toFit;

    ExperimentMultiplicity mult;

    ExperimentRatio ratio;

    FittedQuantity() {
      toFit = true;
      type = FittedQuantity::Multiplicity;
      mult = ExperimentMultiplicity(-211, 10., 1.);
    }

    FittedQuantity(const ExperimentMultiplicity & op) {
      toFit = true;
      type = FittedQuantity::Multiplicity;
      mult = op;
    }

    FittedQuantity(const ExperimentRatio & op) {
      toFit = true;
      type = FittedQuantity::Ratio;
      ratio = op;
    }

    double Value() const {
      if (type == Multiplicity)
        return mult.fValue;
      else
        return ratio.fValue;
    }

    double ValueError() const {
      if (type == Multiplicity)
        return mult.fError;
      else
        return ratio.fError;
    }
  };

} // namespace thermalfist

#endif