00001 #include "TSystem.h"
00002 #include "TString.h"
00003
00004 void RunPythia(const char *jobid="test",const char *bin)
00005 {
00006 Int_t nev=99999999;
00007
00008 TString input("/star/u/russcher/gamma/analysis/data/pythia/");
00009 input+=TString(bin);
00010 input+=TString(".root");
00011 TString outdir("/star/u/russcher/gamma/analysis/output/pythia/");
00012 outdir=outdir+TString(bin)+"_"+TString(jobid)+"/";
00013 TString psout("pi0_pythia.ps");
00014 TString psout2("pi0_pythia2.ps");
00015 TString rootout("pi0_pythia.root");
00016 psout.Prepend(bin);
00017 psout2.Prepend(bin);
00018 rootout.Prepend(bin);
00019 psout.Prepend(jobid);
00020 psout2.Prepend(jobid);
00021 rootout.Prepend(jobid);
00022 TString command("mkdir ");
00023 command.Append(outdir.Data());
00024 gSystem->Exec(command.Data());
00025 cout<<endl<<"storing results in: "<<command.Data()<<endl<<endl;
00026
00027 psout.Prepend(outdir.Data());
00028 psout2.Prepend(outdir.Data());
00029 rootout.Prepend(outdir.Data());
00030
00031 gSystem->Load("$HOME/MyEvent/MyEvent.so");
00032 gSystem->Load("$HOME/gamma/analysis/lib/AnaCuts.so");
00033 gSystem->Load("$HOME/gamma/analysis/lib/EventMixer.so");
00034 gSystem->Load("$HOME/gamma/analysis/lib/Pi0Analysis.so");
00035
00036 Pi0Analysis *pi0=new Pi0Analysis(psout.Data(),psout2.Data(),"pp05");
00037 pi0->setMC(kTRUE);
00038 pi0->setPythia(kTRUE);
00039 pi0->init(rootout.Data());
00040 pi0->make(nev,input.Data());
00041 pi0->getYield();
00042 pi0->finish();
00043 }
