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