Version 4.1
Code is located at
/star/institutions/mit/balewski/freezer/2009-W-algoVer4.1
It is set up to run on 300 Pythia W-events - takes ~1 minute.
execute
root4star rdMuWana.C
It will display the first 2 identified Ws and save collection of histograms in the local file:
mcSetD1_ppWprod.wana.hist.root
To display all of them use macro
root.exe plWana.C'(1)'
change page index in the range 1-15. Macro can display all & save as PS files execute
.L plWana.C
doAllMC()
Key components of the W-reco code.
The St2009WMaker::Make() does the job of finding all potential Ws & Zs in given event consisting of many primary vertices.
The result is a Wevent-tree defined as follows
User can access this tree from another maker:
St2009pubWanaMaker::Make()
which is 'friend' of St2009WMaker and has access to all private variables.
This way one can study W-algo code w/o affecting the code itself. This way future changes of W-algo code are better isolated from the evaluation code.
A short example of W-evaluation code is in St2009pubWanaMaker::evalWeleTrackSign()
St2009pubWanaMaker saves its histograms to the common output root file and those can be viewed by the macro
plPub1.C - the same concept of plotting pages as in plWana.C.
Run9 data processing
Two scripts:
multiBfcJobTempl.xml , subBfcTempl-multi.sh
allow execution of the rdMuWana.C. Sending all 423 jobs to scheduler takes about 15 minutes, results are available 10-20 minutes later.
The directory: fillListA/ contains sorted muDst file lists for SL09b production.
All works together . Instruction on demand.
Version 4.42, by Justin
Code is located at /star/u/stevens4/wAnalysis/ver4.42/
-"Away Cone" calculation: There's a new function called findAwayCone() that searches for seeds for a "jet" on the away side and calculates it's pt.
-Access to MC record in geant.root files: All of this is contained in a new maker St2009pubMcMaker, so it can be switched on and off in the macro. It also requires submitting separate jobs for each MC file (for StIOMaker). There is a file list there (mcSetD1_ppWprod_fileList.lis) to do this.
-Transverse Mass Calculation: This function to calculate the "hadronic recoil" is in St2009W_accessMuDst.cxx, but calculating the transverse mass is done in St2009pubJSMaker.cxx and St2009pubMcMaker.cxx for data and mc respectively.
I haven't added any cuts from any of these new features so the yields should be the same as before. I added a new macro plTransMass.C to display figures related to the transverse mass calculation.