我有一个对有经验的用户来说可能天真的问题，但我有点卡住了。我试图让ApplGrid代码正常工作（https://applgrid.hepforge.org/index.htm），所以我希望运行提供的示例。
据我所知，依赖项是ROOT（root.cern.ch）和LHAPDF库（lhapdf.hepforge.org/index.html），它们在我的机器上运行良好。 上述示例是为了使用旧版本的LHAPDF而构建的，因此我不得不更改调用方法的名称。 当我尝试制作时，我得到关于该函数的未定义引用的错误，但如果我隔离该指令并执行
g++ extract.cxx `lhapdf-config --cflags --ldflags`
它就像一个魅力。
config命令的输出为：
-I / home / matteo / software / gcc-4.7 / lhapdf / install / include -I / home / matteo / software / gcc-4.7 / boost / install / include -L / home / matteo / software / gcc- 4.7 / lhapdf / install / lib -lLHAPDF
如果我删除＆＃34; LHAPDF相关命令＆＃34;从原来的cxx文件我可以运行＆＃34; make＆＃34;没有任何问题。
作为最后的测试，我试图找出＆＃34; applgrid相关代码＆＃34;并尝试编译：
g++ extract_2.cxx `applgrid-config --cxxflags --ldfflags`
但是我得到了一个未定义的引用错误。
config命令的输出为：
-I / home / matteo / software / gcc-4.7 / applgrid / install / include -pthread -m64 -I / home / matteo / software / gcc-4.7 / root / install / include -L / home / matteo / software / gcc-4.7 / applgrid / install / lib -lfAPPLgrid -lAPPLgrid -m64 -L / home / matteo / software / gcc-4.7 / root / install / lib -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d - lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic -L / home / matteo / software / gcc-4.7 / hoppet / install / lib -lhoppet_v1
makefile的内容如下：
CXX = g++
F77 = gfortran
FFLAGS += -O3 -fPIC
CXXFLAGS += -O3 -fPIC
# root
ROOTINCS = $(shell root-config --cflags)
ROOTLIBS = $(shell root-config --glibs)
ROOTARCH = $(findstring -m64, $(ROOTINCS) )
#LHAPDF
LHAPDFINCS = $(shell lhapdf-config --cflags)
LHAPDFDIR = $(shell lhapdf-config --libdir)
LHAPDFLIBS = $(shell lhapdf-config --libs)
# applgrid
APPLCXXFLAGS = $(shell applgrid-config --cxxflags)
APPLCLIBS = $(shell applgrid-config --ldcflags)
APPLFLIBS = $(shell applgrid-config --ldflags)
# hoppet
HOPPETLIBS = $(shell hoppet-config --libs)
# get the fotran runtime library for linking fortran
FRTLLIB = $(shell gfortran $(CXXFLAGS) -print-file-name=libgfortran.a)
FRTLIB = -L$(subst /libgfortran.a, ,$(FRTLLIB) ) -lgfortran
# now set up the compile and link flags and libs
CXXFLAGS += $(ROOTARCH) $(ROOTINCS) $(APPLCXXFLAGS) $(LHAPDFINCS)
LDFLAGS += $(ROOTARCH)
FFLAGS += $(ROOTARCH)
CLIBS += $(ROOTLIBS) $(LHAPDFLIBS) $(HOPPETLIBS) $(APPLCLIBS)
FLIBS += $(ROOTLIBS) $(LHAPDFLIBS) $(HOPPETLIBS) $(APPLFLIBS) $(APPLCLIBS) $(FRTLIB)
install : all
all : stand fnlo fstand getdata
stand: stand.o
$(CXX) $(LDFLAGS) -o $@ $< $(CLIBS)
fstand: fstand.o fmain.o
$(CXX) $(LDFLAGS) -o $@ fstand.o fmain.o $(FLIBS)
% : %.o
$(CXX) $(LDFLAGS) -o $@ $< $(CLIBS)
fnlo: fnmain.o
$(CXX) $(LDFLAGS) -o $@ $< $(CLIBS)
.SUFFIXES : .cxx .o .f .c
.f.o :
$(F77) $(FFLAGS) -c $<
.cxx.o:
$(CXX) $(CXXFLAGS) -c $<
#########################################################
# just type make getdata and it will download all
# the required data for these examples
#########################################################
.PRECIOUS : fnlodata appldata
getdata : getfnlodata getappldata
getfnlodata: fnl0004.tab
fnl0004.tab:
@ echo "\n************************\ndownloading fastnlo data\n************************\n"
(curl http://fastnlo.hepforge.org/code/src/fnl0004.tab.gz | gunzip > fnl0004.tab )
getappldata: atlas-incljets-arxiv-1009.5908v2
atlas-incljets-arxiv-1009.5908v2:
@ echo "\n*************************************\n downloading atlas inclusive jet data\n***********************************\n"
(curl http://www.hepforge.org/archive/applgrid/atlas-incljets-arxiv-1009.5908v2.tgz | tar -xzf - )
clean:
rm -rf ./.libs ./.obj *.lo *.o *.la stand fnlo fstand *~
而目录中的其他文件如下：
fmain.cxx
extern "C" void fstand_();
int main() {
fstand_();
return 0
fnmain.cxx
#include <iostream>
#include <stdio.h>
#include "appl_grid/fastnlo.h"
#include "TFile.h"
#include "TH1D.h"
// lhapdf routines
#include "LHAPDF/LHAPDF.h"
extern "C" void evolvepdf_(const double& , const double& , double* );
extern "C" double alphaspdf_(const double& Q);
int main(int argc, char** argv) {
if ( argc<2 ) {
std::cerr << "usage: fnlo fastnlogrid.tab" << std::endl;
return -1;
}
std::cout << "set up lhapdf..." << std::endl;
const std::string _pdfname = "cteq6mE.LHgrid";
int iset = 0;
LHAPDF::initPDFSet( _pdfname, iset );
// don't need to hard code, can use a runtime parameter...
std::string gridname = argv[1];
// std::string gridname = "fnt1007midp.tab";
// std::string gridname = "fnt1008midp.tab";
// std::string gridname = "fnt2004.tab";
// std::string gridname = "fnl2004.tab";
fastnlo f( gridname );
std::vector<appl::grid*> g = f.grids();
// g.push_back( new appl::grid("atlas-incljets04-eta1.root") );
/// histograms and file for seeing the results
std::vector<TH1D*> hc(g.size());
std::string foutname = "appl.root";
TFile fout( foutname.c_str(),"recreate");
for ( int i=0 ; i<g.size() ; i++ ) {
// trim the grids (not actually needed, done in
// the fastnlo constructor
g[i]->trim();
char hname[64];
sprintf(hname, "hist%02d", i);
/// optionally print out the grid documentation
std::cout << "\n" << g[i]->getDocumentation() << std::endl;
/// perform the convolution
hc[i] = g[i]->convolute( evolvepdf_, alphaspdf_ );
hc[i]->SetName(hname);
hc[i]->SetDirectory(&fout);
hc[i]->Write();
/// print out the results
for ( int j=1 ; j<=hc[i]->GetNbinsX() ; j++ ) {
std::cout << "xsec(" << j-1 << ")=" << hc[i]->GetBinContent(j) << std::endl;
}
}
std::cout << "writing file " << foutname << std::endl;
fout.Close();
return 0;
}
}
fstand.f
C----------------------------------------------------------
C dummy routines that call the pdf and alphas routines
C this example, they just call the lhapdf routine
C for apdf fit, call your own routines
C----------------------------------------------------------
double precision function fnalphas(Q)
double precision Q
double precision alphaspdf
fnalphas = alphaspdf(Q)
return
end
subroutine fnpdf(x, Q, xf)
double precision x, Q
double precision xf(13)
call evolvePDF(x, Q, xf)
return
end
C----------------------------------------------------------
stand.cxx
#include <iostream>
#include <string>
#include <vector>
#include "appl_grid/appl_grid.h"
#include "appl_grid/appl_timer.h"
#include "TH1D.h"
#include "TPad.h"
// lhapdf routines
#include "LHAPDF/LHAPDF.h"
extern "C" void evolvepdf_(const double& , const double& , double* );
extern "C" double alphaspdf_(const double& Q);
int main(int argc, char** argv) {
// use a default atlas inclusive grid
std::string gridname = "atlas-incljets04-eta1.root";
if ( argc>1 ) gridname = argv[1];
std::cout << "reading grid " << gridname << std::endl;
// get name of grid from user and create from grid file
appl::grid g(gridname);
g.trim(); // trim away uneeded memory
/// print the grid documentation
std::cout << g.getDocumentation() << std::endl;
// initialise lhapdf
std::cout << "setting up lhapdf" << std::endl;
// const std::string _pdfname = "cteq6mE.LHgrid";
const std::string _pdfname = "CT10.LHgrid";
int iset = 0;
LHAPDF::initPDFSet( _pdfname, iset );
// initpdfset_(_pdfname.c_str());
// initpdf_( iset );
// do the convolution into a vector
std::cout << "doing standalone convolution" << std::endl;
// struct timeval atimer = appl_timer_start();
std::vector<double> xsec = g.vconvolute( evolvepdf_, alphaspdf_ );
// double atime = appl_timer_stop(atimer);
// std::cout << "time " << atime << " ms" << std::endl;
for ( int i=0 ; i<xsec.size() ; i++ ) std::cout << "xsec[" << i << "]\t= " << xsec[i] << std::endl;
// do the convolution into a TH1D
// atimer = appl_timer_start();
TH1D* hxsec = g.convolute( evolvepdf_, alphaspdf_ );
// atime = appl_timer_stop(atimer);
// std::cout << "time " << atime << " ms" << std::endl;
hxsec->SetLineColor(kRed);
hxsec->SetMarkerColor(kRed);
hxsec->SetMarkerStyle(20);
hxsec->DrawCopy();
gPad->SetLogy(true);
gPad->Print("xsec.pdf");
return 0;
}
待机full.cxx
#include <iostream>
#include <string>
#include <vector>
#include "appl_grid/appl_grid.h"
#include "appl_grid/appl_timer.h"
#include "TH1D.h"
#include "TPad.h"
// lhapdf routines
#include "LHAPDF/LHAPDF.h"
extern "C" void evolvepdf_(const double& , const double& , double* );
extern "C" double alphaspdf_(const double& Q);
int main(int argc, char** argv) {
// use a default atlas inclusive grid
std::string gridname = "atlas-incljets04-eta1.root";
if ( argc>1 ) gridname = argv[1];
std::cout << "reading grid " << gridname << std::endl;
// get name of grid from user and create from grid file
appl::grid g(gridname);
g.trim(); // trim away uneeded memory
std::cout << "setting up lhapdf" << std::endl;
// initialise lhapdf
// const std::string _pdfname = "cteq6mE.LHgrid";
const std::string _pdfname = "CT10.LHgrid";
int iset = 0;
LHAPDF::initPDFSet( _pdfname, iset );
// initpdfset_(_pdfname.c_str());
// initpdf_( iset );
// do the convolution into a vector
std::cout << "doing standalone convolution" << std::endl;
std::cout << g.getDocumentation() << std::endl;
struct timeval atimer = appl_timer_start();
std::vector<double> xsec = g.vconvolute( evolvepdf_, alphaspdf_ );
double atime = appl_timer_stop(atimer);
const int Nsubproc = g.subProcesses();
std::cout << "time : " << atime << " ms\tNsubproc : " << Nsubproc << std::endl;
std::vector<std::vector<double> > xsb(Nsubproc);
for ( int i=0 ; i<Nsubproc ; i++ ) {
xsb[i] = g.vconvolute( i, evolvepdf_, alphaspdf_ );
}
// print results
for ( int i=0 ; i<xsec.size() ; i++ ) {
std::cout << "xsec(" << i << ")=" << xsec[i] << std::endl;
}
for ( int i=0 ; i<Nsubproc ; i++ ) {
std::cout << "subproc: " << i << std::endl;
for ( int j=0 ; j<xsec.size() ; j++ ) {
std::cout << "xsec(" << i << ")=" << xsb[i][j] << std::endl;
}
}
TH1D* hxsec = g.convolute( evolvepdf_, alphaspdf_ );
std::vector<TH1D*> hxsb(Nsubproc);
for ( int i=0 ; i<Nsubproc ; i++ ) {
hxsb[i] = g.convolute( i, evolvepdf_, alphaspdf_ );
}
hxsec->SetLineColor(kRed);
hxsec->SetMarkerColor(kRed);
hxsec->SetMarkerStyle(20);
hxsec->SetMinimum(0.000001);
hxsec->DrawCopy();
gPad->SetLogy(true);
for ( int i=0 ; i<Nsubproc ; i++ ) {
hxsb[i]->SetLineColor(kRed+i%5);
hxsb[i]->SetMarkerColor(kRed+i%5);
hxsb[i]->SetLineStyle(1+i%5);
hxsb[i]->SetMarkerStyle(20+i%5);
hxsb[i]->DrawCopy("same");
}
gPad->Print("xsec.pdf");
double escale = 7e3/8e3;
TH1D* hxsec8 = g.convolute( escale, evolvepdf_, alphaspdf_ );
escale = 7e3/2.76e3;
TH1D* hxsec276 = g.convolute( escale, evolvepdf_, alphaspdf_ );
hxsec->DrawCopy();
hxsec8->SetLineColor(kRed+1);
hxsec8->SetLineStyle(2);
hxsec8->DrawCopy("same");
gPad->Print("xsec8.pdf");
hxsec276->SetLineColor(kRed+2);
hxsec276->SetLineStyle(3);
hxsec276->DrawCopy("same");
gPad->Print("xsec276.pdf");
// or get into a histogram
// TH1D* hxsec = g.convolute( evolvepdf_, alphaspdf_ );
// hxsec->SetName("xsec");
return 0;
}
待机save.cxx
#include <iostream>
#include <string>
#include <vector>
#include "appl_grid/appl_grid.h"
#include "appl_grid/appl_timer.h"
#include "TH1D.h"
#include "TPad.h"
// lhapdf routines
#include "LHAPDF/LHAPDF.h"
extern "C" void evolvepdf_(const double& , const double& , double* );
extern "C" double alphaspdf_(const double& Q);
int main(int argc, char** argv) {
// use a default atlas inclusive grid
std::string gridname = "atlas-incljets04-eta1.root";
if ( argc>1 ) gridname = argv[1];
std::cout << "reading grid " << gridname << std::endl;
// get name of grid from user and create from grid file
appl::grid g(gridname);
g.trim(); // trim away uneeded memory
/// print the grid documentation
std::cout << g.getDocumentation() << std::endl;
// initialise lhapdf
std::cout << "setting up lhapdf" << std::endl;
// const string _pdfname = "cteq6mE.LHgrid";
const string _pdfname = "CT10.LHgrid";
int iset = 0;
LHAPDF::initPDFSet( _pdfname, iset );
// initpdfset_(_pdfname.c_str());
// initpdf_( iset );
// do the convolution into a vector
std::cout << "doing standalone convolution" << std::endl;
// struct timeval atimer = appl_timer_start();
std::vector<double> xsec = g.vconvolute( evolvepdf_, alphaspdf_ );
// double atime = appl_timer_stop(atimer);
// std::cout << "time " << atime << " ms" << std::endl;
for ( int i=0 ; i<xsec.size() ; i++ ) std::cout << "xsec[" << i << "]\t= " << xsec[i] << std::endl;
// do the convolution into a TH1D
// atimer = appl_timer_start();
TH1D* hxsec = g.convolute( evolvepdf_, alphaspdf_ );
// atime = appl_timer_stop(atimer);
// std::cout << "time " << atime << " ms" << std::endl;
hxsec->SetLineColor(kRed);
hxsec->SetMarkerColor(kRed);
hxsec->SetMarkerStyle(20);
hxsec->DrawCopy();
gPad->SetLogy(true);
gPad->Print("xsec.pdf");
return 0;
}
待机tutorial.cxx
#include <iostream>
#include <string>
#include <vector>
#include "appl_grid/appl_grid.h"
#include "appl_grid/appl_timer.h"
#include "TH1D.h"
#include "TLatex.h"
#include "TFile.h"
#include "TH2D.h"
#include "TStyle.h"
#include "TCanvas.h"
#include "TPad.h"
#include "TStyle.h"
#include "TColor.h"
#include "TPaveStats.h"
// lhapdf routines
#include "LHAPDF/LHAPDF.h"
extern "C" void evolvepdf_(const double& x, const double& Q, double* xf);
extern "C" double alphaspdf_(const double& Q);
void setpalette()
{
const Int_t NRGBs = 7;
const Int_t NCont = 98;
Double_t stops[NRGBs] = { 0.00, 0.20, 0.40, 0.60, 0.70, 0.90, 1.00 };
Double_t red[NRGBs] = { 0.00, 0.00, 0.00, 0.60, 0.90, 1.00, 0.50 };
Double_t green[NRGBs] = { 0.00, 0.00, 0.40, 0.95, 1.00, 0.20, 0.00 };
Double_t blue[NRGBs] = { 0.00, 1.00, 1.00, 0.10, 0.00, 0.00, 0.00 };
TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
gStyle->SetNumberContours(NCont);
}
TH2D* smooth(TH2D* h) {
int Nnbins = 10;
TH2D* hi = new TH2D("h2","",
h->GetXaxis()->GetNbins()*Nnbins,
h->GetXaxis()->GetBinLowEdge(1),
h->GetXaxis()->GetBinLowEdge(h->GetXaxis()->GetNbins()+1),
h->GetYaxis()->GetNbins()*Nnbins,
h->GetYaxis()->GetBinLowEdge(1),
h->GetYaxis()->GetBinLowEdge(h->GetYaxis()->GetNbins()+1) );
hi->GetXaxis()->SetTitle( h->GetXaxis()->GetTitle() );
hi->GetYaxis()->SetTitle( h->GetYaxis()->GetTitle() );
hi->GetYaxis()->SetTitleOffset(1.1);
for ( int i=1 ; i<=h->GetXaxis()->GetNbins() ; i++ ) {
for ( int j=1 ; j<=h->GetYaxis()->GetNbins() ; j++ ) {
double z1 = h->GetBinContent(i,j);
double z2 = h->GetBinContent(i+1,j);
double z3 = h->GetBinContent(i,j+1);
double z4 = h->GetBinContent(i+1,j+1);
if ( z2==0 ) z2 = z1;
if ( z3==0 ) z3 = z1;
if ( z4==0 ) z4 = z3;
for ( int ii=0 ; ii<Nnbins ; ii++ ) {
double fx = 1.0*ii/Nnbins;
for ( int jj=0 ; jj<Nnbins ; jj++ ) {
double fy = 1.0*jj/Nnbins;
double z = ( (z4-z3-z2+z1)*fx + z3 - z1 )*fy + (z2-z1)*fx + z1;
// std::cout << (i-1)*Nnbins+ii+1 << " " << (j-1)*Nnbins+jj+1 << "\t z : " << z << "\tz1-4: " << z1 << " " << z2 << " " << z3 << " " << z4 << std::endl;
hi->SetBinContent( (i-1)*Nnbins+ii+1, (j-1)*Nnbins+jj+1, z);
// hi->SetBinContent( (i-1)*Nnbins+ii+1, (j-1)*Nnbins+jj+1, h->GetBinContent(i,j) );
}
}
}
}
return hi;
}
int main(int argc, char** argv) {
gStyle->SetOptStat(0);
setpalette();
// use a default atlas inclusive grid
std::string gridname = "atlas-incljets-arxiv-1009.5908v2/r04/atlas-incljets-eta1.root";
if ( argc>1 ) gridname = argv[1];
std::cout << "reading grid " << gridname << std::endl;
appl::grid g(gridname);
// g.trim(); // trim away uneeded memory
/// print the grid documentation
std::cout << g.getDocumentation() << std::endl;
// initialise lhapdf
std::cout << "setting up lhapdf" << std::endl;
// const std::string _pdfname = "cteq6mE.LHgrid";
const std::string _pdfname = "CT10.LHgrid";
int iset = 0;
LHAPDF::initPDFSet( _pdfname, iset );
// initpdfset_(_pdfname.c_str());
// initpdf_( iset );
/// ----------------- example 1 --------------------
if ( true ) {
// do the convolution into a vector
std::cout << "example 1: standalone convolution" << std::endl;
struct timeval atimer = appl_timer_start();
std::vector<double> xsec = g.vconvolute( evolvepdf_, alphaspdf_ );
double atime = appl_timer_stop(atimer);
std::cout << "time : " << atime << " ms" << std::endl;
for ( int i=0 ; i<xsec.size() ; i++ ) std::cout << "xsec[" << i << "]\t= " << xsec[i] << std::endl;
TH1D* hxsec = g.convolute( evolvepdf_, alphaspdf_ );
hxsec->SetTitle("");
hxsec->DrawCopy();
gPad->SetLogy(true);
gPad->Print("example-1.pdf");
}
/// ----------------- example 2 (corrections) --------------------
if ( false ) {
std::cout << "example 2: multiplicative corrections" << std::endl;
std::cout << g.getApplyCorrections() << std::endl;
g.setApplyCorrections(false);
TH1D* huncor = g.convolute( evolvepdf_, alphaspdf_ );
g.setApplyCorrections(true);
TH1D* hcor = g.convolute( evolvepdf_, alphaspdf_ );
huncor->SetLineStyle(2);
huncor->SetTitle("");
huncor->DrawCopy();
hcor->DrawCopy("same");
gPad->Print("example-2.pdf");
}
/// ----------------- example 3 (cms scaling) --------------------
if ( false ) {
std::cout << "example 3: centre-of-mass energy scaling " << std::endl;
TH1D* h7 = g.convolute( evolvepdf_, alphaspdf_ );
double Escale = 7/2.46;
TH1D* h246 = g.convolute( Escale, evolvepdf_, alphaspdf_ );
std::cout << "cms scale " << g.getCMSScale() << std::endl;
h7->SetTitle("");
h7->SetMinimum(0.001);
gStyle->SetOptStat(0);
h7->DrawCopy();
h246->SetLineColor(kRed);
h246->DrawCopy("same");
gPad->SetLogy(true);
gPad->Print("example-3.pdf");
}
/// ------------------ example 4 (scale variation) -----------------
if ( false ) {
std::cout << "example 4: scale variation " << std::endl;
int nloops = 1; // loop order
int Nbins = 30;
TH2D* hscale = new TH2D("scale", "", Nbins, 0.45, 2.05, Nbins, 0.45, 2.05 );
struct timeval atimer = appl_timer_start();
for ( int i=1 ; i<=hscale->GetNbinsX() ; i++ ) {
double fscale = hscale->GetXaxis()->GetBinCenter(i);
for ( int j=1 ; j<=hscale->GetNbinsY() ; j++ ) {
double rscale = hscale->GetYaxis()->GetBinCenter(j);
std::vector<double> xs = g.vconvolute( evolvepdf_, alphaspdf_, nloops, rscale, fscale );
double total = 0;
for ( int k=xs.size() ; k-- ; ) total += xs[k];
hscale->Fill( fscale, rscale, total );
}
}
double atime = appl_timer_stop(atimer);
std::cout << "time " << atime << " ms" << std::endl;
/// plot ....
gStyle->SetPalette(52);
gStyle->SetPadRightMargin(0.17);
TCanvas* tc = new TCanvas( "","", 650, 600 );
hscale->GetXaxis()->SetTitle("factorisation scale factor");
hscale->GetYaxis()->SetTitle("renormalisation scale factor");
hscale->GetYaxis()->SetTitleOffset(1.1);
TH2D* _hscale = smooth(hscale);
_hscale->DrawCopy("colz");
gPad->Print("example-4.png");
TFile fs("scale.root","recreate");
fs.cd();
hscale->Write();
fs.Write();
fs.Close();
}
/// ------------------ example 5 (subprocesses) -----------------
if ( false ) {
std::cout << "example 4: subprocess contributions " << std::endl;
std::vector<double> xsec = g.vconvolute( evolvepdf_, alphaspdf_);
std::vector<TH1D*> xsec_sub(g.subProcesses());
for ( int i=0 ; i<g.subProcesses() ; i++ ) {
xsec_sub[i] = g.convolute_subproc( i, evolvepdf_, alphaspdf_);
for ( int j=0 ; j<xsec_sub[i]->GetNbinsX() ; j++ ) {
xsec_sub[i]->SetBinContent( j+1, xsec_sub[i]->GetBinContent(j+1)/xsec[j] );
}
}
/// plot ...
int col[7] = { kRed, kBlue, kMagenta, kCyan, kGreen, kBlack, kViolet };
xsec_sub[0]->SetMaximum(0.6);
xsec_sub[0]->SetMinimum(0);
xsec_sub[0]->GetXaxis()->SetTitle("P_{T} [GeV]");
xsec_sub[0]->GetYaxis()->SetTitle("subprocess fraction");
xsec_sub[0]->SetTitle("");
xsec_sub[0]->DrawCopy();
for ( int i=0 ; i<g.subProcesses() ; i++ ) {
xsec_sub[i]->SetLineColor(col[i]);
xsec_sub[i]->DrawCopy("same");
TLatex* tl = new TLatex(0,0,"");
tl->SetTextColor(col[i]);
tl->SetTextSize(0.04);
char label[64];
sprintf(label, "subprocess %d", i);
tl->DrawLatex(430, 0.55-0.035*i, label );
}
gPad->SetLogy(false);
gPad->Print("example-5.pdf");
}
/// -------------------- example 6 (pdf sets) -----------------------
if ( false ) {
// const std::string _pdfname = "cteq6mE.LHgrid";
const std::string _pdfname = "CT10.LHgrid";
int iset = 0;
LHAPDF::initPDFSet( _pdfname, iset );
std::vector<double> _xs = g.vconvolute( evolvepdf_, alphaspdf_ );
int Nset = LHAPDF::numberPDF();
std::vector<TH1D*> xs(Nset);
for ( int i=0 ; i<Nset ; i++ ) {
struct timeval atimer = appl_timer_start();
LHAPDF::initPDFSet( _pdfname, i );
double atime = appl_timer_stop(atimer);
std::cout << "set " << i << "\tinitPDF time " << atime << " ms" << std::endl;
atimer = appl_timer_start();
xs[i] = g.convolute( evolvepdf_, alphaspdf_ );
atime = appl_timer_stop(atimer);
std::cout << "convolution time " << atime << " ms" << std::endl;
for ( int j=_xs.size() ; j-- ; ) xs[i]->SetBinContent( j+1, xs[i]->GetBinContent(j+1)/_xs[j]);
}
xs[0]->SetMaximum(1.1);
xs[0]->SetMinimum(0.9);
gPad->SetLogy(false);
gPad->SetLogx(true);
xs[0]->DrawCopy("l");
for ( int i=0 ; i<Nset ; i++ ) xs[i]->DrawCopy("lsame");
gPad->Print("example-6.pdf");
gPad->SetLogx(false);
}
/// ------------------- example 7 (timers) ---------------------
if ( false ) {
gStyle->SetOptStat(2220);
std::cout << "example 7 : timings" << std::endl;
TH1D* htrim = new TH1D("htrim", "", 100, 38, 45);
// TH1D* htrim2 = new TH1D("htrim2", "", 100, 37, 45);
TH1D* huntrim = new TH1D("huntrim", "", 100, 38, 45);
g.trim();
std::cout << "grid size : " << g.size() << std::endl;
for ( int i=400 ; i-- ; ) {
struct timeval atimer = appl_timer_start();
std::vector<double> xsec = g.vconvolute( evolvepdf_, alphaspdf_ );
double atime = appl_timer_stop( atimer );
htrim->Fill(atime);
}
g.untrim();
std::cout << "grid size : " << g.size() << std::endl;
for ( int i=400 ; i-- ; ) {
struct timeval atimer = appl_timer_start();
std::vector<double> xsec = g.vconvolute( evolvepdf_, alphaspdf_ );
double atime = appl_timer_stop( atimer );
huntrim->Fill(atime);
}
if ( huntrim->GetMaximum()<htrim->GetMaximum() ) huntrim->SetMaximum(1.1*htrim->GetMaximum());
else htrim->SetMaximum(1.1*huntrim->GetMaximum());
huntrim->SetLineStyle(2);
huntrim->SetLineColor(kRed);
htrim->GetXaxis()->SetTitle("time [ms]");
huntrim->GetXaxis()->SetTitle("time [ms]");
huntrim->Draw();
htrim->Draw("same");
// htrim2->DrawCopy("same");
huntrim->Draw("same");
htrim->Draw("sames+");
huntrim->Draw("sames+");
gPad->Update();
TPaveStats* tp1 = (TPaveStats*)htrim->GetListOfFunctions()->FindObject("stats");
TPaveStats* tp2 = (TPaveStats*)huntrim->GetListOfFunctions()->FindObject("stats");
std::cout << tp1 << " " << tp2 << std::endl;
tp2->SetTextColor(kRed);
tp1->SetY1NDC(0.93-0.15);
tp1->SetY2NDC(0.93);
tp2->SetY1NDC(0.93-0.16-0.15);
tp2->SetY2NDC(0.93-0.16);
tp1->Draw();
tp2->Draw();
gPad->Update();
gPad->Print("example-7.pdf");
}
return 0;
}
任何帮助将不胜感激，谢谢。