{
// For getting pdf uncertainty

//setup cdfsoft2 5.3.3
//setenv  MRST2003_DATA=/home/miller/analysis/ewk/Electroweak/mrst2003
//setenv CTEQ6_DATA=/home/miller/analysis/ewk/Electroweak/cteq6/data
//stnfit.exe
// .x run_pdf_example.C


gROOT->Reset();

gROOT->ProcessLine(".! cp TopTree.h.mc.pdf2 TopTree.h");

TTree *tree;

chain = new TChain("TopTree"); 
// this is your MC ntuple file
chain->Add("/data/nglas08/a/stdenis/topnt_topnt_top_533_nt5_ttopel_1-Inclusive.root");

//TopTree->MakeClass();

gSystem->SetIncludePath("-I./ -I./include -I$CDFSOFT2_DIR/include");
//gSystem->SetIncludePath(" -I$MYDIR/JetCorrRel/include -I$MYDIR " );
//gROOT->ProcessLine(".L $ROOTSYS/lib/libPhysics.so");  // load Physics lib
gROOT->ProcessLine(".! cp Pdfweight_example.C TopTree.C");
gROOT->ProcessLine(".L TopTree.C+");  // Cint compiler
gROOT->ProcessLine("TopTree t(chain); t.Loop(); ");



}