EventMix/macro/runmb.C

0001 #pragma once
0002 #if ROOT_VERSION_CODE >= ROOT_VERSION(6,00,0)
0003 #include <fun4all/SubsysReco.h>
0004 #include <fun4all/Fun4AllServer.h>
0005 #include <fun4all/Fun4AllInputManager.h>
0006 #include <fun4all/Fun4AllDstInputManager.h>
0007 #include <fun4all/Fun4AllDstOutputManager.h>
0008 #include <trackreco/PHGenFitTrackProjection.h>
0009
0010 //#include </gpfs/mnt/gpfs02/sphenix/user/lebedev/mdc/test/analysis/EventMix/install/include/eventmix/PairMaker.h>
0011 #include </gpfs/mnt/gpfs02/sphenix/user/lebedev/mdc/test/analysis/EventMix/install/include/eventmix/PairMaker.h>
0012
0013 R__LOAD_LIBRARY(libfun4all.so)
0014 R__LOAD_LIBRARY(libeventmix.so)
0015 #endif
0016
0017 void runmb(int kkk=2000)
0018 //void runmb(const char *fname = "/sphenix/sim/sim01/sphnxpro/MDC1/embed/embedDST_sHijing_0_20fm_50kHz_bkg_0_20fm-0000000001-02990.root")
0019 {
0020   gSystem->Load("libg4dst");
0021   gSystem->Load("libeventmix");
0022
0023   char outfilename[99];
0024   char infilename1[99];
0025   char infilename2[99];
0026   sprintf(outfilename,"/sphenix/sim/sim01/sphnxpro/MDC1/embed/eePairs_mb%d.root",kkk);
0027   cout << "out: " << outfilename << endl;
0028
0029   Fun4AllServer *se = Fun4AllServer::instance();
0030   se->Verbosity(1);
0031
0032   PHGenFitTrackProjection* projection = new PHGenFitTrackProjection();
0033   projection->Verbosity(1);
0034   se->registerSubsystem(projection);
0035
0036   PairMaker *pmaker = new PairMaker("PairMaker","dummy.root");
0037   se->registerSubsystem(pmaker);
0038
0039   Fun4AllInputManager *in1 = new Fun4AllDstInputManager("in1");
0040   se->registerInputManager(in1);
0041   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("in2");
0042   se->registerInputManager(in2);
0043
0044   for(int i=0; i<10; i++) {
0045     sprintf(infilename1,"DST_TRACKS_sHijing_0_20fm_50kHz_bkg_0_20fm-0000000001-0%d.root",kkk+i);
0046     cout << "in1:  " << infilename1 << endl;
0047     in1->AddFile(infilename1);
0048     sprintf(infilename2,"DST_CALO_CLUSTER_sHijing_0_20fm_50kHz_bkg_0_20fm-0000000001-0%d.root",kkk+i);
0049     cout << "in2:  " << infilename2 << endl;
0050     in2->AddFile(infilename2);
0051   }
0052
0053   Fun4AllOutputManager *outee = new Fun4AllDstOutputManager("outee",oufilename);
0054   outee->Verbosity(1);
0055   outee->AddNode("ElectronPairs");
0056   se->registerOutputManager(outee);
0057   outee->Print();
0058
0059   se->run();
0060   se->End();
0061 }
0062