void runPicoNpeEventMaker(const Char_t *inputFile="test.list", const Char_t *outputFile="test")
{
//Check STAR Library. Please set SL_version to the original star library used in the production from http://www.star.bnl.gov/devcgi/dbProdOptionRetrv.pl
string SL_version = "SL15c";
string env_SL = getenv ("STAR");
if(env_SL.find(SL_version)==string::npos)
{
cout<<"Environment Star Library does not match the requested library in runPicoNpeEventMaker.C. Exiting..."<<endl;
exit(1);
}
Int_t nEvents = 10000000;
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
gSystem->Load("StPicoDstMaker");
gSystem->Load("StPicoPrescales");
gSystem->Load("StPicoNpeEventMaker");
chain = new StChain();
StPicoDstMaker* picoDstMaker = new StPicoDstMaker(0,inputFile,"picoDstMaker");
StPicoNpeEventMaker* picoNpeMaker = new StPicoNpeEventMaker("picoNpeMaker",picoDstMaker,outputFile);
chain->Init();
cout<<"chain->Init();"<<endl;
int total = picoDstMaker->chain()->GetEntries();
cout << " Total entries = " << total << endl;
if(nEvents>total) nEvents = total;
for (Int_t i=0; i<nEvents; i++)
{
if(i%1000==0)
cout << "Working on eventNumber " << i << endl;
chain->Clear();
int iret = chain->Make(i);
if (iret) { cout << "Bad return code!" << iret << endl; break;}
total++;
}
cout << "****************************************** " << endl;
cout << "Work done... now its time to close up shop!"<< endl;
cout << "****************************************** " << endl;
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << nEvents << endl;
cout << "****************************************** " << endl;
delete chain;
}
void readPicoDst(const Char_t *inputFile="test.list", const Char_t *outputFile="test.root")
{
Int_t nEvents = 200000000;
// Int_t nEvents = 1000;
//Load all the System libraries
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
gSystem->Load("StPicoDstMaker");
gSystem->Load("StMyAnalysisMaker");
gSystem->Load("StMyTreeMaker");
chain = new StChain();
StPicoDstMaker *picoMaker = new StPicoDstMaker(0,inputFile,"picoDst");
StPicoAnaTreeMaker *treeMaker = new StPicoAnaTreeMaker(1,inputFile,picoMaker);
chain->Init();
cout<<"chain->Init();"<<endl;
int total = picoMaker->chain()->GetEntries();
cout << " Total entries = " << total << endl;
if(nEvents>total) nEvents = total;
for (Int_t i=0; i<nEvents; i++){
// if(i%1000==0)
cout << "Working on eventNumber " << i << endl;
chain->Clear();
int iret = chain->Make(i);
if (iret) { cout << "Bad return code!" << iret << endl; break;}
total++;
}
cout << "****************************************** " << endl;
cout << "Work done... now its time to close up shop!"<< endl;
cout << "****************************************** " << endl;
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << nEvents << endl;
cout << "****************************************** " << endl;
delete chain;
}
void runPicoHFMyAnaMaker(const Char_t *inputFile="test.list", const Char_t *outputFile="outputBaseName", unsigned int makerMode = 0 /*kAnalyze*/,
const Char_t *badRunListFileName = "picoList_bad_MB.list") {
// -- Check STAR Library. Please set SL_version to the original star library used in the production
// from http://www.star.bnl.gov/devcgi/dbProdOptionRetrv.pl
string SL_version = "SL15c";
string env_SL = getenv ("STAR");
if (env_SL.find(SL_version)==string::npos) {
cout<<"Environment Star Library does not match the requested library in runPicoHFMyAnaMaker.C. Exiting..."<<endl;
exit(1);
}
Int_t nEvents = 10000000;
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
gSystem->Load("StBTofUtil");
gSystem->Load("StPicoDstMaker");
gSystem->Load("StPicoPrescales");
gSystem->Load("StPicoHFMaker");
gSystem->Load("StPicoHFMyAnaMaker");
chain = new StChain();
// ========================================================================================
//makerMode = StPicoHFMaker::kAnalyze;
// ========================================================================================
cout << "Maker Mode " << makerMode << endl;
TString sInputFile(inputFile);
TString sInputListHF("");
if (makerMode == StPicoHFMaker::kAnalyze) {
if (!sInputFile.Contains(".list") && !sInputFile.Contains("picoDst.root")) {
cout << "No input list or picoDst root file provided! Exiting..." << endl;
exit(1);
}
}
else if (makerMode == StPicoHFMaker::kWrite) {
if (!sInputFile.Contains("picoDst.root")) {
cout << "No input picoDst root file provided! Exiting..." << endl;
exit(1);
}
}
else if (makerMode == StPicoHFMaker::kRead) {
if (!sInputFile.Contains(".list")) {
cout << "No input list provided! Exiting..." << endl;
exit(1);
}
// -- prepare filelist for picoDst from hfTrees
sInputListHF = sInputFile;
sInputFile = "tmpPico.list";
TString command = "sed 's|^.*hfTree|/project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics/P15ic|g' "
+ sInputListHF + " > " + sInputFile;
gSystem->Exec(command.Data());
command = "sed -i 's|picoHFtree|picoDst|g' " + sInputFile;
gSystem->Exec(command.Data());
}
else {
cout << "Unknown makerMode! Exiting..." << endl;
exit(1);
}
StPicoDstMaker* picoDstMaker = new StPicoDstMaker(0, sInputFile, "picoDstMaker");
StPicoHFMyAnaMaker* picoHFMyAnaMaker = new StPicoHFMyAnaMaker("picoHFMyAnaMaker", picoDstMaker, outputFile, sInputListHF);
picoHFMyAnaMaker->setMakerMode(makerMode);
StHFCuts* hfCuts = new StHFCuts("hfBaseCuts");
picoHFMyAnaMaker->setHFBaseCuts(hfCuts);
// ---------------------------------------------------
// -- Set Base cuts for HF analysis
// -- File name of bad run list
hfCuts->setBadRunListFileName(badRunListFileName);
// -- ADD USER CUTS HERE ----------------------------
// hfCuts->setCutVzMax(6.);
// hfCuts->setCutVzVpdVzMax(3.);
// hfCuts->setCutTriggerWord(0x1F);
// hfCuts->setCutNHitsFitMax(15);
// hfCuts->setCutRequireHFT(true);
// hfCuts->setCutNHitsFitnHitsMax(0.52);
// ---------------------------------------------------
// -- Channel1
picoHFMyAnaMaker->setDecayMode(StPicoHFEvent::kTwoParticleDecay);
picoHFMyAnaMaker->setDecayChannel(StPicoHFMyAnaMaker::kChannel1);
// -- ADD USER CUTS HERE ----------------------------
// ========================================================================================
// ========================================================================================
chain->Init();
cout << "chain->Init();" << endl;
int total = picoDstMaker->chain()->GetEntries();
cout << " Total entries = " << total << endl;
if(nEvents>total) nEvents = total;
for (Int_t i=0; i<nEvents; i++) {
if(i%10000==0)
cout << "Working on eventNumber " << i << endl;
chain->Clear();
int iret = chain->Make(i);
if (iret) { cout << "Bad return code!" << iret << endl; break;}
total++;
}
cout << "****************************************** " << endl;
cout << "Work done... now its time to close up shop!"<< endl;
cout << "****************************************** " << endl;
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << nEvents << endl;
cout << "****************************************** " << endl;
delete chain;
// -- clean up if in read mode
if (makerMode == StPicoHFMaker::kRead)
gSystem->Exec(Form("rm -f %s", sInputFile.Data()));
}
void runPicoMixedEvent(const Char_t *inputFile="test.list", const Char_t *outputFile="outputBaseName",
const Char_t *badRunListFileName = "picoList_bad_MB.list") {
// -- Check STAR Library. Please set SL_version to the original star library used in the production
// from http://www.star.bnl.gov/devcgi/dbProdOptionRetrv.pl
TStopwatch* stopWatch = new TStopwatch();
stopWatch->Start();
StMemStat mem;
string SL_version = "SL15c";
string env_SL = getenv ("STAR");
if (env_SL.find(SL_version)==string::npos) {
cout<<"Environment Star Library does not match the requested library in runPicoMixedEventMaker.C. Exiting..."<<endl;
exit(1);
}
// ========================================================================================
// Testing
// ========================================================================================
Int_t nEvents = 10000000;
//Int_t nEvents = 20000;
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
// gSystem->Load("StBTofUtil");
gSystem->Load("StPicoDstMaker");
gSystem->Load("StPicoPrescales");
// gSystem->Load("StPicoCutsBase");
// gSystem->Load("StPicoHFMaker");
gSystem->Load("StRefMultCorr");
gSystem->Load("StEventPlane");
gSystem->Load("StPicoMixedEventMaker");
chain = new StChain();
// ========================================================================================
// Still bone dry
// ========================================================================================
TString sInputFile(inputFile);
TString sInputListHF("");
// create list of kfVertex files
TString command = "sed 's/picodsts/hft\\\/kfVertex/g' " + sInputFile + " >correspondingkfVertex.list";
gSystem->Exec(command.Data());
command = "sed -i 's/picoDst/kfVertex/g' correspondingkfVertex.list";
gSystem->Exec(command.Data());
// ========================================================================================
StPicoDstMaker* picoDstMaker = new StPicoDstMaker(0, sInputFile, "picoDstMaker");
StRefMultCorr* grefmultCorrUtil = CentralityMaker::instance()->getgRefMultCorr();
cout<<"here"<<endl;
grefmultCorrUtil->setVzForWeight(6, -6.0, 6.0);
grefmultCorrUtil->readScaleForWeight("StRoot/StRefMultCorr/macros/weight_grefmult_vpd30_vpd5_Run14.txt");
for(Int_t i=0;i<6;i++){
cout << i << " " << grefmultCorrUtil->get(i, 0) << endl;
}
StEventPlane* eventPlaneMaker = new StEventPlane("eventPlaneMaker",picoDstMaker,grefmultCorrUtil);
StPicoMixedEventMaker* picoMixedEventMaker = new StPicoMixedEventMaker("picoMixedEventMaker", picoDstMaker, grefmultCorrUtil, eventPlaneMaker, outputFile, sInputListHF, "correspondingkfVertex.list");
// ---------------------------------------------------
// -- Set Base cuts for HF analysis
// -- File name of bad run list
//hfCuts->setBadRunListFileName(badRunListFileName);
chain->Init();
cout << "chain->Init();" << endl;
int total = picoDstMaker->chain()->GetEntries();
cout << " Total entries = " << total << endl;
if(nEvents>total) nEvents = total;
for (Int_t i=0; i<nEvents; i++) {
if(i%1000==0)
cout << "Working on eventNumber " << i << endl;
chain->Clear();
int iret = chain->Make(i);
if (iret) { cout << "Bad return code!" << iret << endl; break;}
total++;
}
cout << "****************************************** " << endl;
cout << "Work done... now its time to close up shop!"<< endl;
cout << "****************************************** " << endl;
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << nEvents << endl;
cout << "****************************************** " << endl;
delete chain;
// delete list of kfEvent list
command = "rm -f correspondingkfVertex.list";
gSystem->Exec(command.Data());
stopWatch->Stop();
stopWatch->Print();
}
void runPicoMixedEvent(const Char_t *inputFile="test.list", const Char_t *outputFile="outputBaseName",
const Char_t *badRunListFileName = "picoList_bad_MB.list") {
// -- Check STAR Library. Please set SL_version to the original star library used in the production
// from http://www.star.bnl.gov/devcgi/dbProdOptionRetrv.pl
StMemStat mem;
string SL_version = "SL15c";
string env_SL = getenv ("STAR");
if (env_SL.find(SL_version)==string::npos) {
cout<<"Environment Star Library does not match the requested library in runPicoMixedEventMaker.C. Exiting..."<<endl;
exit(1);
}
// ========================================================================================
// Testing
// ========================================================================================
Int_t nEvents = 10000000;
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
gSystem->Load("StBTofUtil");
gSystem->Load("StPicoDstMaker");
gSystem->Load("StPicoPrescales");
gSystem->Load("StPicoHFMaker");
gSystem->Load("StPicoMixedEventMaker");
chain = new StChain();
// ========================================================================================
// Still bone dry
// ========================================================================================
TString sInputFile(inputFile);
TString sInputListHF("");
// ========================================================================================
StPicoDstMaker* picoDstMaker = new StPicoDstMaker(0, sInputFile, "picoDstMaker");
StPicoMixedEventMaker* picoMixedEventMaker = new StPicoMixedEventMaker("picoMixedEventMaker", picoDstMaker, outputFile, sInputListHF);
StHFCuts* hfCuts = new StHFCuts("hfBaseCuts");
picoMixedEventMaker -> setHFBaseCuts(hfCuts);
// ---------------------------------------------------
// -- Set Base cuts for HF analysis
// -- File name of bad run list
//hfCuts->setBadRunListFileName(badRunListFileName);
// -- Event and track cuts ----------------------------
hfCuts->setCutVzMax(6.);
hfCuts->setCutVzVpdVzMax(3.);
//hfCuts->setCutTriggerWord(0x1F);
hfCuts->setCutNHitsFitMax(15);
hfCuts->setCutRequireHFT(false);
// hfCuts->setCutNHitsFitnHitsMax(0.52);
// -- Topological and PID cuts ------------------------
hfCuts->setCutSecondaryPair(0.75, 1.0, 7.0,
0.9, 0,2.0 );
//Single track pt
hfCuts->setCutPtRange(0.6,2.0,StHFCuts::kPion);
hfCuts->setCutPtRange(0.6,2.0,StHFCuts::kKaon);
//TPC setters
hfCuts->setCutTPCNSigmaPion(3.0);
hfCuts->setCutTPCNSigmaKaon(3.0);
//TOF setters, need to set pt range as well
hfCuts->setCutTOFDeltaOneOverBeta(0.05, StHFCuts::kKaon);
hfCuts->setCutPtotRangeHybridTOF(0.6,2.0,StHFCuts::kKaon);
hfCuts->setCutTOFDeltaOneOverBeta(0.05, StHFCuts::kPion);
hfCuts->setCutPtotRangeHybridTOF(0.6,2.0,StHFCuts::kPion);
// ---------------------------------------------------
chain->Init();
cout << "chain->Init();" << endl;
int total = picoDstMaker->chain()->GetEntries();
cout << " Total entries = " << total << endl;
if(nEvents>total) nEvents = total;
for (Int_t i=0; i<nEvents; i++) {
if(i%1==0)
cout << "Working on eventNumber " << i << endl;
chain->Clear();
int iret = chain->Make(i);
if (iret) { cout << "Bad return code!" << iret << endl; break;}
total++;
}
cout << "****************************************** " << endl;
cout << "Work done... now its time to close up shop!"<< endl;
cout << "****************************************** " << endl;
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << nEvents << endl;
cout << "****************************************** " << endl;
delete chain;
}
void makePicoDst(const Int_t runnumber=999999,
const Char_t *inputFile="/star/institutions/lbl_prod/mlomnitz/SSD_Simu/StEvent/hijing_**.MuDst.root",
const Char_t *inputGeant="/star/institutions/lbl_prod/mlomnitz/SSD_Simu/StEvent/hijing_**.McEvent.root",
TString outfilename="out",
const bool creatingPhiWgt = kFALSE, const int prodMod = 0, const int emcMode=0,
const bool createMcBranch = kTRUE
){
Int_t nEvents = 10000000;
//Int_t nEvents = 50;
//Load all the System libraries
cout<<inputFile<<endl;
gSystem->Load("libTable");
gSystem->Load("libPhysics");
gSystem->Load("St_base");
gSystem->Load("StChain");
gSystem->Load("St_Tables");
gSystem->Load("StUtilities"); // new addition 22jul99
gSystem->Load("StTreeMaker");
gSystem->Load("StIOMaker");
gSystem->Load("StarClassLibrary");
gSystem->Load("StTriggerDataMaker"); // new starting from April 2003
gSystem->Load("StBichsel");
gSystem->Load("StEvent");
gSystem->Load("StEventUtilities");
gSystem->Load("StDbLib");
gSystem->Load("StEmcUtil");
gSystem->Load("StTofUtil");
gSystem->Load("StPmdUtil");
gSystem->Load("StPreEclMaker");
gSystem->Load("StStrangeMuDstMaker");
gSystem->Load("StMuDSTMaker");
if( createMcBranch){
gSystem->Load("StMcEvent");
gSystem->Load("StAssociationMaker");
}
if(!creatingPhiWgt&&emcMode) {
gSystem->Load("StTpcDb");
gSystem->Load("StMcEvent");
gSystem->Load("StMcEventMaker");
gSystem->Load("StDaqLib");
gSystem->Load("libgen_Tables");
gSystem->Load("libsim_Tables");
gSystem->Load("libglobal_Tables");
gSystem->Load("StEmcTriggerMaker");
gSystem->Load("StEmcUtil");//mine
gSystem->Load("StEmcRawMaker");
gSystem->Load("StEmcADCtoEMaker");
gSystem->Load("StPreEclMaker");
gSystem->Load("StEpcMaker");
gSystem->Load("StEmcSimulatorMaker");
gSystem->Load("StEmcUtil");
gSystem->Load("StDbBroker");
gSystem->Load("StDetectorDbMaker");
gSystem->Load("StDbUtilities");
gSystem->Load("StEEmcUtil");
gSystem->Load("StEEmcDbMaker");
gSystem->Load("St_db_Maker");
gSystem->Load("StTriggerUtilities");
}
gSystem->Load("StPicoDstMaker");
chain = new StChain();
if( createMcBranch ){
StIOMaker* ioMaker = new StIOMaker("IO","r",inputGeant,"bfcTree");
ioMaker->SetDebug();
ioMaker->SetIOMode("r");
ioMaker->SetBranch("*",0,"0"); //deactivate all branches
//ioMaker->SetBranch("geantBranch",0,"r"); //activate geant Branch
ioMaker->SetBranch("McEventBranch",0,"r"); //activate geant Branch
ioMaker->SetBranch("eventBranch",0,"r"); //activate geant Branch
StAssociationMaker* associator = new StAssociationMaker;
cout<<"created new StAssociationMaker"<<endl;
associator->useInTracker();
}
StMuDstMaker *MuDstMaker = new StMuDstMaker(0,0,"",inputFile,"MuDst",100);
MuDstMaker->SetStatus("*",0);
MuDstMaker->SetStatus("MuEvent",1);
MuDstMaker->SetStatus("PrimaryVertices",1);
MuDstMaker->SetStatus("PrimaryTracks",1);
MuDstMaker->SetStatus("GlobalTracks",1);
MuDstMaker->SetStatus("CovGlobTrack",1);
//MuDstMaker->SetStatus("BTof*",1);
//MuDstMaker->SetStatus("Emc*",1);
if(!creatingPhiWgt&&emcMode) {
St_db_Maker *dbMk = new St_db_Maker("db","MySQL:StarDb","$STAR/StarDb","StarDb");
// Endcap database
StEEmcDbMaker* eemcDb = new StEEmcDbMaker;
StEmcADCtoEMaker *adc2e = new StEmcADCtoEMaker();
adc2e->setPrint(false);
//adc2e->setFillHisto(false);
//adc2e->setDebug(false); //more histograms
//adc2e->setSMDRmsCut(0,0);
adc2e->saveAllStEvent(true);
//adc2e->setRemoveGhostEvent(false);
//adc2e->setTallyHist(mTally);
//adc2e->setDbName("Calibrations/emc/y3");
StPreEclMaker *pre_ecl=new StPreEclMaker();
pre_ecl->setPrint(kFALSE);
StEpcMaker *epc=new StEpcMaker();
epc->setPrint(kFALSE);
#if 1
// Trigger simulator
StTriggerSimuMaker* trigSimu = new StTriggerSimuMaker;
trigSimu->setMC(false);
trigSimu->useBemc();
trigSimu->useEemc();
trigSimu->useOnlineDB();
trigSimu->bemc->setConfig(StBemcTriggerSimu::kOffline);
#endif
}
StPicoDstMaker *picoMaker = new StPicoDstMaker(1,Form("st_%i.picoDst.root",runnumber),"picoDst");
picoMaker->setRunNumber(runnumber);
picoMaker->setEmcMode(emcMode); // 0-No EMC, 1-EMC ON
picoMaker->setMcMode(createMcBranch);
chain->Init();
cout<<"chain->Init();"<<endl;
int total = 0;
cout<<MuDstMaker->chain()->GetEntries()<<endl;
if( nEvents > MuDstMaker->chain()->GetEntries() )
nEvents=MuDstMaker->chain()->GetEntries();
for (Int_t i=0; i<nEvents; i++){
if(i%100==0)
cout << "Working on eventNumber " << i << endl;
chain->Clear();
int iret = chain->Make(i);
if (iret) { cout << "Bad return code!" << iret << endl; break;}
total++;
}
cout << "****************************************** " << endl;
cout << "Work done... now its time to close up shop!"<< endl;
cout << "****************************************** " << endl;
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << total << endl;
cout << "****************************************** " << endl;
delete chain;
}
void runPicoDpmAnaMaker(const Char_t *inputFile="test.list", const Char_t *outputFile="outputBaseName", unsigned int makerMode = 0 /*kAnalyze*/,
const Char_t *badRunListFileName = "picoList_bad_MB.list", const Char_t *treeName = "picoHFtree",
const Char_t *productionBasePath = "/project/projectdirs/starprod/picodsts/Run14/AuAu/200GeV/physics2/P15ic",
unsigned int decayChannel = 0 /* kChannel0 */) {
// -- Check STAR Library. Please set SL_version to the original star library used in the production
// from http://www.star.bnl.gov/devcgi/dbProdOptionRetrv.pl
string SL_version = "SL15c";
string env_SL = getenv ("STAR");
if (env_SL.find(SL_version)==string::npos) {
cout<<"Environment Star Library does not match the requested library in runPicoDpmAnaMaker.C. Exiting..."<<endl;
exit(1);
}
Int_t nEvents = 10000000;
//Int_t nEvents = 2000;
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
gSystem->Load("StBTofUtil");
gSystem->Load("StPicoDstMaker");
gSystem->Load("StPicoPrescales");
gSystem->Load("StRefMultCorr");
gSystem->Load("StEventPlane");
gSystem->Load("StPicoCutsBase");
gSystem->Load("StPicoHFMaker");
gSystem->Load("StPicoHFMyAnaMaker");
gSystem->Load("StPicoDpmAnaMaker");
chain = new StChain();
// ========================================================================================
//makerMode = StPicoHFMaker::kAnalyze;
// ========================================================================================
cout << "Maker Mode " << makerMode << endl;
TString sInputFile(inputFile);
TString sInputListHF("");
TString sProductionBasePath(productionBasePath);
TString sTreeName(treeName);
if (makerMode == StPicoHFMaker::kAnalyze) {
if (!sInputFile.Contains(".list") && !sInputFile.Contains("picoDst.root")) {
cout << "No input list or picoDst root file provided! Exiting..." << endl;
exit(1);
}
}
else if (makerMode == StPicoHFMaker::kWrite) {
if (!sInputFile.Contains("picoDst.root")) {
cout << "No input picoDst root file provided! Exiting..." << endl;
exit(1);
}
}
else if (makerMode == StPicoHFMaker::kRead) {
if (!sInputFile.Contains(".list")) {
cout << "No input list provided! Exiting..." << endl;
exit(1);
}
// -- prepare filelist for picoDst from hfTrees
sInputListHF = sInputFile;
sInputFile = "tmpPico.list";
TString command = "sed 's|" + sTreeName + ".root|picoDst.root|g' " + sInputListHF + " > " + sInputFile;
cout << "COMMAND : " << command << endl;
gSystem->Exec(command.Data());
command = "sed -i 's|^.*" + sTreeName + "|" + sProductionBasePath + "|g' " + sInputFile; // + " > " + sInputFile;
cout << "COMMAND : " << command << endl;
gSystem->Exec(command.Data());
}
else {
cout << "Unknown makerMode! Exiting..." << endl;
exit(1);
}
StPicoDstMaker* picoDstMaker = new StPicoDstMaker(0, sInputFile, "picoDstMaker");
StRefMultCorr* grefmultCorrUtil = CentralityMaker::instance()->getgRefMultCorr();
StEventPlane* eventPlaneMaker = new StEventPlane("eventPlaneMaker",picoDstMaker,grefmultCorrUtil);
StPicoDpmAnaMaker* picoDpmAnaMaker = new StPicoDpmAnaMaker("picoDpmAnaMaker", picoDstMaker,
grefmultCorrUtil, eventPlaneMaker,
outputFile, sInputListHF);
picoDpmAnaMaker->setMakerMode(makerMode);
StHFCuts* hfCuts = new StHFCuts("hfBaseCuts");
picoDpmAnaMaker->setHFBaseCuts(hfCuts);
// ---------------------------------------------------
// -- Set Base cuts for HF analysis
// -- File name of bad run list
hfCuts->setBadRunListFileName(badRunListFileName);
// -- ADD USER CUTS HERE ----------------------------
hfCuts->setCutVzMax(6.);
hfCuts->setCutVzVpdVzMax(3.);
//hfCuts->setCutTriggerWord(0x1F);
hfCuts->setCutNHitsFitMax(15);
hfCuts->setCutRequireHFT(true);
// hfCuts->setCutNHitsFitnHitsMax(0.52);
// ---------------------------------------------------
// -- Channel1
picoDpmAnaMaker->setDecayMode(StPicoHFEvent::kThreeParticleDecay);
picoDpmAnaMaker->setDecayChannel(StPicoDpmAnaMaker::kChannel1);
// -- ADD USER CUTS HERE ----------------------------
// --- Lomnitz topological cuts
hfCuts->setCutSecondaryTriplet(0.01, 0.01,0.01,
0.005, 0.2,
0.995, 1.7, 2.1);
// 0.008, 0.008, 0.008);
// --- Lomnitz cuts to remove noise from ghosting
//Single track pt
hfCuts->setCutPtRange(0.6,10.0,StHFCuts::kPion);
hfCuts->setCutPtRange(0.6,10.0,StHFCuts::kKaon);
//TPC setters
hfCuts->setCutTPCNSigmaPion(3.0);
hfCuts->setCutTPCNSigmaKaon(3.0);
//TOF setters, need to set pt range as well
hfCuts->setCutTOFDeltaOneOverBeta(0.05, StHFCuts::kKaon);
hfCuts->setCutPtotRangeHybridTOF(0.6,10.0,StHFCuts::kKaon);
// ========================================================================================
// ========================================================================================
chain->Init();
cout << "chain->Init();" << endl;
int total = picoDstMaker->chain()->GetEntries();
cout << " Total entries = " << total << endl;
if(nEvents>total) nEvents = total;
for (Int_t i=0; i<nEvents; i++) {
if(i%10000==0)
cout << "Working on eventNumber " << i << endl;
chain->Clear();
int iret = chain->Make(i);
if (iret) { cout << "Bad return code!" << iret << endl; break;}
total++;
}
cout << "****************************************** " << endl;
cout << "Work done... now its time to close up shop!"<< endl;
cout << "****************************************** " << endl;
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << nEvents << endl;
cout << "****************************************** " << endl;
delete chain;
// -- clean up if in read mode
//if (makerMode == StPicoHFMaker::kRead)
//gSystem->Exec(Form("rm -f %s", sInputFile.Data()));
}
