BlastWave.h

Go to the documentation of this file.

#ifndef BLASTWAVE_H
#define BLASTWAVE_H
 
//#include <cmath>
#include "TMath.h"
#include "TSpline.h"
#include "common/Acceptance.h"
#include <iostream>
//#include <omp.h>
 
 
class BlastWave {
  std::vector<double> xlag32, wlag32;
  std::vector<double> xleg32, wleg32;
  std::vector<double> xlegeta, wlegeta;
  double fMass;
  int fPDGID;
  bool fUseAcceptance;
  double fYmin, fYmax, fYcm, fEtaMax;
  //    double fNorm;
  Acceptance::AcceptanceFunction
    fAcceptance;  //fAcceptanceSTS, fAcceptanceSTSTOF;
  Acceptance::ReconstructionEfficiencyFunction fReconstructionEfficiency;
  TSpline3 fTamt;
  TSpline3 fNormT, fNormT4pi;
 
public:
  BlastWave(double mass  = 0.938,
            int PDGID    = 2212,
            bool fUseAcc = false,
            double ymin  = -3.,
            double ymax  = 3.,
            double ycm   = 2.,
            double etam  = 0.8);
 
  virtual ~BlastWave(void) {}
 
  double GetT(double amt) { return fTamt.Eval(amt); }
 
  double GetA(double multiplicity, double T) {
    if (!fUseAcceptance)
      return multiplicity;
    else
      return multiplicity * fNormT4pi.Eval(T) / fNormT.Eval(T);
  }
 
  double GetAerror(double multiplicity, double T, double multerr, double Terr) {
    if (!fUseAcceptance)
      return multerr;
    else
      return TMath::Sqrt(
        fNormT4pi.Eval(T) / fNormT.Eval(T) * fNormT4pi.Eval(T) / fNormT.Eval(T)
          * multerr * multerr
        + multiplicity * multiplicity / fNormT.Eval(T) / fNormT.Eval(T)
            * fNormT4pi.Derivative(T) * fNormT4pi.Derivative(T) * Terr * Terr
        + multiplicity * multiplicity / fNormT.Eval(T) / fNormT.Eval(T)
            / fNormT.Eval(T) / fNormT.Eval(T) * fNormT4pi.Eval(T)
            * fNormT4pi.Eval(T) * fNormT.Derivative(T) * fNormT.Derivative(T)
            * Terr * Terr);
  }
 
  double fmt(double amt, double T);
 
  double dTdmt(double amt) { return fTamt.Derivative(amt); }
 
  double mtAv(double T);
  double mtAv2(double T);
  double Normalization(double T);
  double Normalization4pi(double T);
};
 
#endif