StRoot: StJetMaker/macros/RunJetFinder2009.C Source File

00001 //
00002 // Pibero Djawotho <pibero@tamu.edu>
00003 // Texas A&M University
00004 // 10 March 2010
00005 //
00006 
00007 void RunJetFinder2009(int nevents = 1e6,
00008                       const char* mudstfile = "/star/data60/reco/production2009_200Gev_Single/ReversedFullField/P10ic/2009/143/10143008/st_physics_10143008_raw_6020001.MuDst.root",
00009                       const char* jetfile   = "st_physics_10143008_raw_6020001.jets.root",
00010                       const char* skimfile  = "st_physics_10143008_raw_6020001.skim.root",
00011                       bool useL2 = false)
00012 {
00013   cout << "Read MuDst file:\t" << mudstfile << endl;
00014   cout << "Write jet file:\t" << jetfile << endl;
00015   cout << "Write skim file:\t" << skimfile << endl;
00016 
00017   gROOT->Macro("loadMuDst.C");
00018   gROOT->Macro("LoadLogger.C");
00019 
00020   gSystem->Load("StTpcDb");
00021   gSystem->Load("StDetectorDbMaker");
00022   gSystem->Load("StDbUtilities");
00023   gSystem->Load("StMcEvent");
00024   gSystem->Load("StMcEventMaker");
00025   gSystem->Load("StDaqLib");
00026   gSystem->Load("StEmcRawMaker");
00027   gSystem->Load("StEmcADCtoEMaker");
00028   gSystem->Load("StEpcMaker");
00029   gSystem->Load("StEmcSimulatorMaker");
00030   gSystem->Load("StDbBroker");
00031   gSystem->Load("St_db_Maker");
00032   gSystem->Load("StEEmcUtil");
00033   gSystem->Load("StEEmcDbMaker");
00034   gSystem->Load("StSpinDbMaker");
00035   gSystem->Load("StEmcTriggerMaker");
00036   gSystem->Load("StTriggerUtilities");
00037   gSystem->Load("StMCAsymMaker");
00038   gSystem->Load("StRandomSelector");
00039   gSystem->Load("libfastjet.so");
00040   gSystem->Load("StJetFinder");
00041   gSystem->Load("StJetSkimEvent");
00042   gSystem->Load("StJets");
00043   gSystem->Load("StJetEvent");
00044   gSystem->Load("StJetMaker");
00045   gSystem->Load("StTriggerFilterMaker");
00046 
00047   StChain* chain = new StChain; 
00048 
00049   // MuDst reader
00050   StMuDstMaker* muDstMaker = new StMuDstMaker(0,0,"",mudstfile,"",100000,"MuDst");
00051 
00052   // StMuDbReader...
00053   StMuDbReader* db = StMuDbReader::instance();
00054 
00055   // Trigger filter
00056   StTriggerFilterMaker* filterMaker = new StTriggerFilterMaker;
00057 
00058   // 2009 pp500
00059   filterMaker->addTrigger(230410); // JP1
00060   filterMaker->addTrigger(230411); // JP2
00061   filterMaker->addTrigger(230420); // AJP
00062 
00063   // 2009 pp200
00064   // http://www.star.bnl.gov/protected/common/common2009/trigger2009/triggers2009.html
00065   // BHT3
00066   filterMaker->addTrigger(240530);
00067   // L2BGamma
00068   filterMaker->addTrigger(240620);
00069   // L2EGamma
00070   filterMaker->addTrigger(240630);
00071   filterMaker->addTrigger(240631);
00072   // L2JetHigh
00073   filterMaker->addTrigger(240650);
00074   filterMaker->addTrigger(240651);
00075   filterMaker->addTrigger(240652);
00076   // JP1
00077   filterMaker->addTrigger(240410);
00078   filterMaker->addTrigger(240411);
00079   // Upsilon
00080   filterMaker->addTrigger(240641);
00081   // Luminosity monitors
00082   // BBCMB
00083   filterMaker->addTrigger(240010);
00084   filterMaker->addTrigger(240110);
00085   filterMaker->addTrigger(240120);
00086   filterMaker->addTrigger(240220);
00087   // BBCMB-Cat0
00088   filterMaker->addTrigger(240011);
00089   filterMaker->addTrigger(240111);
00090   filterMaker->addTrigger(240121);
00091   filterMaker->addTrigger(240221);
00092   // BBCMB-Cat1
00093   filterMaker->addTrigger(240012);
00094   filterMaker->addTrigger(240112);
00095   filterMaker->addTrigger(240122);
00096   filterMaker->addTrigger(240222);
00097   // BBCMB-Cat2
00098   filterMaker->addTrigger(240013);
00099   filterMaker->addTrigger(240113);
00100   filterMaker->addTrigger(240123);
00101   filterMaker->addTrigger(240223);
00102   // BBCMB-Cat3
00103   filterMaker->addTrigger(240014);
00104   filterMaker->addTrigger(240114);
00105   filterMaker->addTrigger(240124);
00106   filterMaker->addTrigger(240224);
00107   // BBCMB-Cat4
00108   filterMaker->addTrigger(240015);
00109   filterMaker->addTrigger(240115);
00110   filterMaker->addTrigger(240125);
00111   filterMaker->addTrigger(240225);
00112   // VPDMB
00113   filterMaker->addTrigger(240020);
00114   filterMaker->addTrigger(240025);
00115 
00116   // star database
00117   St_db_Maker* starDb = new St_db_Maker("StarDb","MySQL:StarDb");
00118 
00119   // Endcap database
00120   StEEmcDbMaker* eemcDb = new StEEmcDbMaker;
00121 
00122   // Spin database
00123   StSpinDbMaker* spinDb = new StSpinDbMaker;
00124 
00125   // Barrel ADC to energy maker
00126   StEmcADCtoEMaker* adc = new StEmcADCtoEMaker;
00127 
00128   // Trigger simulator
00129   StTriggerSimuMaker* simuTrig = new StTriggerSimuMaker;
00130   simuTrig->setMC(false); // Must be before individual detectors, to be passed
00131   // BBC was not used in Run 9
00132   //simuTrig->useBbc();
00133   simuTrig->useBemc();
00134   simuTrig->useEemc();
00135   simuTrig->bemc->setConfig(StBemcTriggerSimu::kOffline);
00136 
00137   // L2 (only L2btowCalib, L2etowCalib, L2ped, L2jet in CVS as of 17 April 2010)
00138   if (useL2) {
00139     StL2_2009EmulatorMaker* simL2Mk = new StL2_2009EmulatorMaker;
00140     assert(simL2Mk);
00141     simL2Mk->setSetupPath("/star/u/pibero/public/StarTrigSimuSetup/");
00142     simL2Mk->setOutPath("./");
00143     simuTrig->useL2(simL2Mk);
00144   }
00145 
00146   // Add Mike's 4p maker.
00147   // Here we also tag whether or not to do the swap.
00148   // The classes available for correcting tower energy for tracks are:
00149   // 1. StjTowerEnergyCorrectionForTracksMip
00150   // 2. StjTowerEnergyCorrectionForTracksFraction
00151   bool doTowerSwapFix = true;
00152   StBET4pMaker* bet4pMaker = new StBET4pMaker("BET4pMaker",muDstMaker,doTowerSwapFix,new StjTowerEnergyCorrectionForTracksFraction(1.00));
00153   bet4pMaker->setUseTPC(true);
00154   bet4pMaker->setUseBEMC(true);
00155   bet4pMaker->setUseEndcap(true);
00156   bet4pMaker->setUse2003Cuts(false);
00157   bet4pMaker->setUse2005Cuts(false);
00158   bet4pMaker->setUse2006Cuts(true);
00159 
00160   // Jet maker
00161   StJetMaker* jetMaker = new StJetMaker("jetMaker",muDstMaker,jetfile);
00162   // Use StJetEvent format. Comment out the line below to use old StJets format.
00163   jetMaker->SetTreeWriter(new StjeJetEventTreeWriter(jetfile));
00164 
00165   // Skim event maker
00166   StJetSkimEventMaker* skimEventMaker = new StJetSkimEventMaker("StJetSkimEventMaker",muDstMaker,skimfile);
00167 
00168   // Set the analysis cuts (See StJetMaker/StppJetAnalyzer.h -> class StppAnaPars)
00169   StppAnaPars* anapars = new StppAnaPars;
00170   anapars->setFlagMin(0); // track->flag() > 0
00171   anapars->setCutPtMin(0.2); // track->pt() > 0.2
00172   anapars->setAbsEtaMax(2.5); // abs(track->eta()) < 2.5
00173   anapars->setJetPtMin(5.0);
00174   anapars->setJetEtaMax(100.0);
00175   anapars->setJetEtaMin(0);
00176   anapars->setJetNmin(0);
00177 
00178   // Setup the cone finder (See StJetFinder/StConeJetFinder.h -> class StConePars)
00179   StConePars* cpars = new StConePars;
00180   cpars->setGridSpacing(105,-3.0,3.0,120,-TMath::Pi(),TMath::Pi());
00181   cpars->setSeedEtMin(0.5);
00182   cpars->setAssocEtMin(0.1);
00183   cpars->setSplitFraction(0.5);
00184   cpars->setPerformMinimization(true);
00185   cpars->setAddMidpoints(true);
00186   cpars->setRequireStableMidpoints(true);
00187   cpars->setDoSplitMerge(true);
00188   cpars->setDebug(false);
00189 
00190   jetMaker->addAnalyzer(anapars,cpars,bet4pMaker,"ConeJets12");
00191 
00192   anapars->setNhits(5);
00193   jetMaker->addAnalyzer(anapars,cpars,bet4pMaker,"ConeJets5");
00194 
00195   anapars->setNhits(1000000);
00196   jetMaker->addAnalyzer(anapars,cpars,bet4pMaker,"ConeJetsEMC");
00197 
00198   chain->Init();
00199   chain->EventLoop(nevents);
00200 }