vhlle/src/main.eostest.cpp

0001 #include <iostream>
0002 #include <fstream>
0003 #include <iomanip>
0004 #include <cstring>
0005 #include <ctime>
0006 #include <sstream>
0007 #include <TApplication.h>
0008 #include <TCanvas.h>
0009 #include <TGraph.h>
0010
0011 #include "eos.h"
0012 #include "eo3.h"
0013 #include "eo1.h"
0014 #include "eoChiral.h"
0015 #include "eoHadron.h"
0016 #include "rmn.h"
0017
0018 using namespace std ;
0019
0020
0021 int main(int argc, char **argv)
0022 {
0023   TApplication theApp("App", &argc, argv);
0024   // pointers to all the main objects
0025   EoS *eos ;
0026
0027   char * eosfile = "eos/Laine_nf3.dat" ;
0028   int ncols = 3, nrows = 286 ;
0029   //eos = new EoSs(eosfile,ncols) ;
0030   eos = new EoSChiral() ;
0031   EoS* eosH = new EoSHadron("eos/eosHadron3D.dat") ;
0032
0033   double e=0.5, p, nb=0.2, nq=0., ns=0., vx=0.8, vy=0., vz=0. ;
0034   double Q [7] ;
0035   transformCV(e, eos->p(e,nb,nq,ns), nb, nq, ns, vx, vy, vz, Q) ;
0036   transformPV(eosH, Q, e, p, nb, nq, ns, vx, vy, vz) ;
0037   cout<<"HelloWorld;\n" ;
0038   transformCV(e, eosH->p(e,nb,nq,ns), nb, nq, ns, vx, vy, vz, Q) ;
0039   cout<<"HelloWorld;\n" ;
0040
0041   const int N = 500 ;
0042   double x[N], y[N] ;
0043   for(int i=0; i<N; i++){
0044    double Tt, mubt, muqt, must, pt ;
0045    double e = i*1.0/N ;
0046    double nb = i*0.2/N ;
0047    double nq = i*0.1/N ;
0048    eosH->eos(0.5, nb, 0., 0., Tt, mubt, muqt, must, pt) ;
0049    x[i] = nb;
0050    y[i] = mubt ;
0051   }
0052   TGraph *g = new TGraph(N,x,y) ;
0053   g->SetMarkerStyle(22) ;
0054   g->SetMarkerSize(0.8) ;
0055   g->Draw("AP") ;
0056   theApp.Run() ;
0057
0058   delete eos ;
0059   delete eosH ;
0060 }