00001 #include "TSystem.h"
00002 #include "TString.h"
00003
00004 void Run(const char *jobid="test")
00005 {
00006 Int_t nev=99999999;
00007
00008 TString input0("/star/u/russcher/gamma/analysis/data/dAu/dAuCombined_0.root");
00009 TString input1("/star/u/russcher/gamma/analysis/data/dAu/dAuCombined_1.root");
00010 TString input2("/star/u/russcher/gamma/analysis/data/dAu/UPCCombined.root");
00011 TString input3("/star/u/russcher/gamma/analysis/data/dAu/dAuUPCCombined.root");
00012
00013 TString outdir("/star/u/russcher/gamma/analysis/output/dAu/");
00014 TString psout("pi0_dAu.ps");
00015 TString rootout("pi0_dAu.root");
00016 psout.Prepend(jobid);
00017 rootout.Prepend(jobid);
00018 TString command("mkdir ");
00019 command.Append(outdir.Data());
00020 command.Append(jobid);
00021 gSystem->Exec(command.Data());
00022 cout<<endl<<"storing results in: "<<command.Data()<<endl<<endl;
00023 outdir.Append(jobid);
00024 outdir.Append("/");
00025
00026 psout.Prepend(outdir.Data());
00027 rootout.Prepend(outdir.Data());
00028
00029 gSystem->Load("$HOME/MyEvent/MyEvent.so");
00030 gSystem->Load("$HOME/gamma/analysis/lib/AnaCuts.so");
00031 gSystem->Load("$HOME/gamma/analysis/lib/EventMixer.so");
00032 gSystem->Load("$HOME/gamma/analysis/lib/Pi0Analysis.so");
00033
00034 Pi0Analysis *pi0=new Pi0Analysis(psout.Data(),"/dev/null","dAu");
00035 pi0->setMC(kFALSE);
00036 pi0->init(rootout.Data());
00037
00038 pi0->make(nev,input0.Data());
00039 pi0->make(nev,input1.Data());
00040 pi0->printPrescales();
00041
00042 pi0->make(nev,input2.Data());
00043 pi0->printPrescales();
00044 pi0->make(nev,input3.Data());
00045 pi0->printPrescales();
00046
00047 pi0->getYield();
00048 pi0->finish();
00049 }
