Page History

 bottom*
----
h4. 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()

h4. 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 |W-event tree ver 4.1]

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.

h4.

 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

.
----
h3. Version 4.42, by

 Justin


h4. Code is located at /star/institutions/iucf/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.

----
{anchor:Walgo-in-cvs

}

h4. W-algo in CVS on Nov 23, 2009

{code

}
cvs co StRoot/StSpinPool/StWalgoB2009

cvs checkout: Updating StRoot/StSpinPool/StWalgoB2009
U StRoot/StSpinPool/StWalgoB2009/St2009WMaker.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009WMaker.h
U StRoot/StSpinPool/StWalgoB2009/St2009W_accessMuDst.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009W_algo.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009W_histo.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009W_trigger.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009WlumiMaker.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009WlumiMaker.h
U StRoot/StSpinPool/StWalgoB2009/St2009Wlumi_histo.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubJSMaker.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubJSMaker.h
U StRoot/StSpinPool/StWalgoB2009/St2009pubJS_histo.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubMcMaker.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubMcMaker.h
U StRoot/StSpinPool/StWalgoB2009/St2009pubMc_histo.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubSpinMaker.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubSpinMaker.h
U StRoot/StSpinPool/StWalgoB2009/St2009pubSpin_histo.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubWanaMaker.cxx
U StRoot/StSpinPool/StWalgoB2009/St2009pubWanaMaker.h
U StRoot/StSpinPool/StWalgoB2009/St2009pubWana_histo.cxx
U StRoot/StSpinPool/StWalgoB2009/WanaConst.h
U StRoot/StSpinPool/StWalgoB2009/Wevent2009.h
U StRoot/StSpinPool/StWalgoB2009/WeventDisplay.cxx
U StRoot/StSpinPool/StWalgoB2009/WeventDisplay.h

cvs checkout: Updating StRoot/StSpinPool/StWalgoB2009/macros
U StRoot/StSpinPool/StWalgoB2009/macros/haddWana.C
U StRoot/StSpinPool/StWalgoB2009/macros/plPubJS.C
U StRoot/StSpinPool/StWalgoB2009/macros/plTransMass.C
U StRoot/StSpinPool/StWalgoB2009/macros/plWana.C
U StRoot/StSpinPool/StWalgoB2009/macros/rdMuWana.C

{code}

h4. Code used for DNP2009 results  was tagged in CVS as 'DNP2009'

To retrieve the tagged version (w/o any changes appled after December 30, 2009) do

{code

}cvs co -r DNP2009 StRoot/StSpinPool/StWalgoB2009 {code}
To retrieve the latest dev version do

{code

}cvs co  StRoot/StSpinPool/StWalgoB2009 {code}
----
h4. W-algo Update on Jan 6, 2010

\-Remove Justin's "homemade jet finder" (real jet finder will be included shortly)

\-Change all track cuts to use primary component instead of global component

.
----
h4. *W-algo Upadate on Jan 8, 2010

*

\-Set useETOW=3 as default (use ETOW in near and away side pT sums)

\-Add jet finder to algo: rdMuWana.C now reads in muDst and jet trees simultaneously (set option useJetFinder=2 to read trees and

 muDst)


                                        jet trees for use at APS are located at /star/institutions/iucf/stevens4/wAnalysis/aps2010/jetTree

                                        define and set cut on ptBalance variable