void readWtree( const Char_t *fileList="./R10081007.lis", int etaBin=7, int maxEvents=1e8, bool spinSort=true)
{
if(etaBin==7) isZ=true;
// load shared libraries
gROOT -> LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
assert( !gSystem->Load("StDetectorDbMaker"));
assert( !gSystem->Load("StTpcDb"));
assert( !gSystem->Load("StDbUtilities"));
assert( !gSystem->Load("StDbBroker"));
assert( !gSystem->Load("St_db_Maker"));
assert( !gSystem->Load("StEEmcUtil"));
assert( !gSystem->Load("StWalgo2011"));
assert( !gSystem->Load("StSpinDbMaker"));
assert( !gSystem->Load("StJets"));
// create analysis chain
chain = new StChain("chain");
TObjArray* HListZ=new TObjArray;
TObjArray* HListEta=new TObjArray;;
//define eta ranges
float etaLow[7] = {-0.9, -0.4, 0.0, 0.4, 0.9, 1.15, -1.};
float etaHigh[7] = {-0.4, 0.0, 0.4, 0.9, 1.15, 1.35, 1.};
//initiate W maker
wTreeMk = new St2011WMaker(Form("Eta%d",etaBin));
wTreeMk->setJetTreeBranch("ConeJets12_100","ConeJets12_100_noEEMC"); //select jet tree braches used
//set cuts for barrle and endcap algos
wTreeMk->setWbosonCuts(25., 0.88, 14., etaLow[etaBin-1], etaHigh[etaBin-1]); //highET, nearTotEtFrac, ptBalance, etaLow, etaHigh
wTreeMk->setE_WbosonCuts(25., 0.85, 14., etaLow[etaBin-1], etaHigh[etaBin-1]); //highET, nearTotEtFrac, ptBalance, etaLow, etaHigh
wTreeMk->setHList(HListEta);
wTreeMk->setMaxDisplayEve(1); // only first N events will get displayed
//......... spin sorting
if(spinSort){
enum {mxSM=1}; // only one for now
St2011pubSpinMaker *spinMkA[mxSM];
for(int kk=0;kk<mxSM;kk++) {
char ttx[100]; sprintf(ttx,"%cEta%dspin",'A',etaBin);
printf("add spinMaker %s %d \n",ttx,kk);
spinMkA[kk]=new St2011pubSpinMaker(ttx,Form("Eta%d",etaBin));
spinMkA[kk]->attachWalgoMaker(wTreeMk);
spinMkA[kk]->setHList(HListEta);
//assign eta bin same as W maker
spinMkA[kk]->setEta(etaLow[etaBin-1],etaHigh[etaBin-1]);
//possibly add other copies for systematic checks later
//if(kk==1) spinMkA[kk]->setEta(-1.,0.);
//if(kk==2) spinMkA[kk]->setEta(0,1.);
//if(kk==3) spinMkA[kk]->setQPT(-1);// disable Q/PT cut
//if(kk==4) spinMkA[kk]->setNoEEMC();
}
}
// Z maker
if (isZ){
ZMk=new St2011ZMaker("Z");
ZMk->attachWalgoMaker(wTreeMk);
ZMk->setHList(HListZ);
ZMk->setNearEtFrac(0.88);
ZMk->setClusterMinEt(15);
ZMk->setPhi12Min(3.1416/2.);
ZMk->setMinZMass(73.); // Zmass -20%
ZMk->setMaxZMass(114.);// Zmass +20%
}
//chain W and jet trees in W maker
chainFiles(fileList);
chain->ls(3);
chain->Init();
Int_t nevents = wTreeMk->getNumberOfEvents();
cout<<nevents<<" events in chain"<<endl;
if(nevents<1) {cout<<"No events in chain, check your list file"<<endl; return; }
Int_t stat = 0; Int_t event = 0;
while ( !stat )
{
if ( nevents>=0 )
if ( event>=nevents || event > maxEvents) break;
chain -> Clear();
stat = chain->Make();
event++;
}
char *file1=fileList;
printf("file1=%s=%s=\n",file1);
TString outFile=""; outFile+=file1;
outFile.ReplaceAll(".lis",Form("_Eta%d.wana.hist.root",etaBin));
outFile.ReplaceAll("lists/treeReader/","");
TFile *outF=new TFile(outFile,"RECREATE");
if(outF->IsOpen()){
TDirectory *Eta = outF->mkdir(Form("Eta%d",etaBin));
Eta->cd();
HListEta->Write();
if(isZ){
TDirectory *Z = outF->mkdir("Z");
Z->cd();
HListZ->Write();
}
cout<<endl<<" Histo saved -->"<<outFile<<endl;
}
else
cout<<endl<<" Couldn't open file "<<outFile<<endl;
delete outF;
}
void readWtree( const Char_t *fileList="./R13104003.lis", int etaBin=8, int maxEvents=1e8, bool spinSort=true)
{
if(etaBin==8) isZ=true;
// load shared libraries
gROOT -> LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
assert( !gSystem->Load("StDetectorDbMaker"));
assert( !gSystem->Load("StTpcDb"));
assert( !gSystem->Load("StDbUtilities"));
assert( !gSystem->Load("StDbBroker"));
assert( !gSystem->Load("St_db_Maker"));
assert( !gSystem->Load("StEEmcUtil"));
assert( !gSystem->Load("StWalgo2011"));
assert( !gSystem->Load("StSpinDbMaker"));
assert( !gSystem->Load("StJets"));
assert( !gSystem->Load("StJetEvent"));
// create analysis chain
chain = new StChain("chain");
TObjArray* HListZ=new TObjArray;
TObjArray* HListEta=new TObjArray;;
//define eta ranges (skip bin 5+6 and 9 is same as 8 except spin sorting)
float etaLow[9] = {-1.1,-0.5, 0.0, 0.5, 1.0, -2.0, 0.7,-1.5,-1.5};
float etaHigh[9] = {-0.5, 0.0, 0.5, 1.1, 2.0, -1.0, 2.5, 2.0, 2.0};
//initiate W maker
wTreeMk = new St2011WMaker(Form("Eta%d",etaBin));
wTreeMk->setJetTreeBranch("AntiKtR060NHits12","AntiKtR060NHits12_noEEMC");
//set cuts for barrle and endcap algos
wTreeMk->setWbosonCuts(25., 0.88, 14., etaLow[etaBin-1], etaHigh[etaBin-1]); //highET, nearTotEtFrac, ptBalance, etaLow, etaHigh
wTreeMk->setE_WbosonCuts(25., 0.85, 20., etaLow[etaBin-1], etaHigh[etaBin-1]); //highET, nearTotEtFrac, ptBalance, etaLow, etaHigh
wTreeMk->setHList(HListEta);
wTreeMk->setMaxDisplayEve(1e6); // only first N events will get displayed
//......... spin sorting
if(spinSort){
if(etaBin==8){
St2011WlumiMaker *lumiMk=new St2011WlumiMaker;
lumiMk->setHList(HListEta);
lumiMk->attachWalgoMaker(wTreeMk);
}
enum {mxSM=2};
St2011pubSpinMaker *spinMkA[mxSM];
for(int kk=0;kk<mxSM;kk++) {
char ttx[100]; sprintf(ttx,"%cEta%dspin",'A'+kk,etaBin);
printf("add spinMaker %s %d \n",ttx,kk);
spinMkA[kk]=new St2011pubSpinMaker(ttx,Form("Eta%d",etaBin));
spinMkA[kk]->attachWalgoMaker(wTreeMk);
spinMkA[kk]->setHList(HListEta);
//assign eta bin same as W maker
spinMkA[kk]->setEta(etaLow[etaBin-1],etaHigh[etaBin-1]);
spinMkA[kk]->setEtaE(etaLow[etaBin-1],etaHigh[etaBin-1]);
if(kk==1) spinMkA[kk]->setQPT(false); //disable Q/pT cut
//special cases
if(etaBin==7) {
spinMkA[kk]->setEta(999.,999.); //mask barrel for for etaBin 7
}
if(etaBin==8) {
spinMkA[kk]->setEta(-1.,1.); //diff eta range for spin sorting
spinMkA[kk]->setEtaE(999.,999.); //mask endcap for for etaBin 8
}
}
}
// Z maker
if (isZ){
ZMk=new St2011ZMaker("Z");
ZMk->attachWalgoMaker(wTreeMk);
ZMk->setHList(HListZ);
ZMk->setNearEtFrac(0.88);
ZMk->setClusterMinEt(14);
ZMk->setPhi12Min(3.1416/2.);
ZMk->setMinZMass(70.); // Zmass -20%
ZMk->setMaxZMass(110.);// Zmass +20%
}
//chain W and jet trees in W maker
chainFiles(fileList);
chain->ls(3);
chain->Init();
Int_t nevents = wTreeMk->getNumberOfEvents();
cout<<nevents<<" events in chain"<<endl;
if(nevents<1) {cout<<"No events in chain, check your list file"<<endl; return; }
Int_t stat = 0; Int_t event = 0;
while ( !stat )
{
if ( nevents>=0 )
if ( event>=nevents || event > maxEvents) break;
chain -> Clear();
stat = chain->Make();
event++;
}
char *file1=fileList;
printf("file1=%s=%s=\n",file1);
TString outFile=""; outFile+=file1;
outFile.ReplaceAll(".lis",Form("_Eta%d.wana.hist.root",etaBin));
outFile.ReplaceAll("lists/treeReader/","");
TFile *outF=new TFile(outFile,"RECREATE");
if(outF->IsOpen()){
TDirectory *Eta = outF->mkdir(Form("Eta%d",etaBin));
Eta->cd();
HListEta->Write();
if(isZ){
TDirectory *Z = outF->mkdir("Z");
Z->cd();
HListZ->Write();
}
cout<<endl<<" Histo saved -->"<<outFile<<endl;
}
else
cout<<endl<<" Couldn't open file "<<outFile<<endl;
delete outF;
}
