Documentation for the photon-jet reconstruction code

Documentation for the photon-jet reconstruction code (Ilya Selyuzhenkov)

 

Analysis flow chart

  • compile everything with cons
     
  • To run makers (expect starsim and bfc) use macros from here and run:

    root4star -b -q 'RunXXXMaker.C("inputFileName.extension.root")'

    Note: You can use files from iucf disk
     

Anylysis flow chart:

starsim (Kumac) -> fzd.gz

    bfc.C (fzd) -> MuDst.root / geant.root

      JetFinder (MuDst) -> jet.root / skim.root

      EEmcDstMaker (MuDst) -> eemc.root

            GammaJetMaker (jet/skim) -> dijet.root

                  GammaJetAnaMaker (dijet/eemc) -> ana.root

                       GammaJetDrawMaker (ana) -> draw.root / mlp.root

EEMC related MuDst/jet/skim/gamma trees location

Main directory for Ilya's private directory files at RCF IUCF disk:

/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC

CDF Tune A simulations (private production)

fzd.gz, geant.root, and MuDst.root are in files/ subdirectory
logs (MuDst.log.gz, sim.log.gz) in logs/ subdirectory

Prompt photons (partonic pt range 3-25, single bin):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/GammaJet_pt3_25_pytune100

QCD (partonic pt range 4-35, bins: 4-6, 6-9, 9-15, 15-35):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt4_6_pytune100
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt6_9_pytune100
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt9_15_pytune100
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt15_35_pytune100

Perugia0 simulations (private production)

fzd.gz, geant.root, and MuDst.root are in files/ subdirectory
logs (MuDst.log.gz, sim.log.gz) in logs/ subdirectory

Prompt photons (partonic pt range 3-25, single bin):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/GammaJet_pt3_25_pytune320

QCD (partonic pt range 4-35, bins: 4-6, 6-9, 9-15, 15-35):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt15_35_pytune320
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt4_6_pytune320
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt6_9_pytune320_1
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt6_9_pytune320_2
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100630/QCD_pt9_15_pytune320

Pro-pT0 simulations (private production)

fzd.gz, geant.root, and MuDst.root are in files/ subdirectory
logs (MuDst.log.gz, sim.log.gz) in logs/ subdirectory

Prompt photons (partonic pt range 3-25, single bin):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100727/GammaJet_pt3_25_pytune329

QCD (partonic pt range 4-35, bins: 4-6, 6-9, 9-15, 15-35):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100727/QCD_pt4_6_pytune329
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100727/QCD_pt6_9_pytune329
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100727/QCD_pt9_15_pytune329
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/20100727/QCD_pt15_35_pytune329

Pro-pT0 (official production)

trees (ana, dijet, draw, eemc, gamma, geant, jet, mlp, MuDst, skim) are in:
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/

logs for:
ana, dijet, draw, eemc, gamma, jet, mlp, skim:
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/prodlog/trees/
MuDst and geant:
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/prodlog/P10ii/log/trs/

Prompt photons (partonic pt range 2-35, bins: 2-3, 3-4, 4-6, 6-9, 9-15, 15-35):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_2-3gev/eemcgammafilt100_gamma
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_3-4gev/eemcgammafilt100_gamma
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_4-6gev/eemcgammafilt100_gamma
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_6-9gev/eemcgammafilt100_gamma
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_9-15gev/eemcgammafilt100_gamma
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_15-35gev/eemcgammafilt100_gamma

QCD (partonic pt range 4-35, bins: 4-6, 6-9, 9-15, 15-35):
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_4-6gev/eemcgammafilt100_qcd
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_6-9gev/eemcgammafilt100_qcd
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_9-15gev/eemcgammafilt100_qcd
/star/institutions/iucf/IlyaSelyuzhenkov/gammaFilterMC/official/pp200/pythia6_423/pt_15-35gev/eemcgammafilt100_qcd

Run 6 Jet/skim trees

/star/institutions/iucf/IlyaSelyuzhenkov/jetTrees/2006/ppProductionLong/

Run 6 gamma trees

/star/institutions/anl/Run6GammaTrees/log/
/star/institutions/anl/Run6GammaTrees/root/

 

 

 

Kumac file examples

Examples of different kumac files:

  • Single particle Monte-Carlo:

    Combine singleParticle_begin.kumac with singleParticle_end.kumac
    using needed geometry tag (example: detp geom y2006h)
     
  • Prompt photon Pythia Monte-Carlo testGammaJet.kumac

    MSUB (14)=1
    MSUB (18)=1
    MSUB (29)=1
    MSUB (114)=1
    MSUB (115)=1
     
  • QCD 2->2 processes testQCD.kumac

    MSUB (11) = 1
    MSUB (12) = 1      
    MSUB (13) = 1      
    MSUB (28) = 1
    MSUB (53) = 1      
    MSUB (68) = 1

 

L2Egamma trigger emulator howto

Running Run 6 L2 gamma trigger emulator

  TObjArray* HList=new TObjArray;
  StTriggerSimuMaker *simuTrig = new StTriggerSimuMaker("StarTrigSimu");
  simuTrig->setHList(HList);
  simuTrig->setMC(true); // must be before individual detectors, to be passed
  simuTrig->useBbc();
  simuTrig->useBemc();
  simuTrig->useEemc(0);
  simuTrig->bemc->setConfig(StBemcTriggerSimu::kOffline);
  StGenericL2Emulator* simL2Mk = new StL2_2006EmulatorMaker;
  assert(simL2Mk);
  simL2Mk->setSetupPath("/afs/rhic.bnl.gov/star/users/kocolosk/public/StarTrigSimuSetup/");
  simL2Mk->setOutPath("/star/institutions/iucf/IlyaSelyuzhenkov/data/MCFilter/StGenericL2Emulator_log/");
  simuTrig->useL2(simL2Mk);
 

Run 6 and Run 9 bfc chain examples

Run 9 bfc options with EEMC slow and fast simulators (EEfs EEss):

"trs,fss,Idst,IAna,l0,tpcI,fcf,ftpc,Tree,logger,ITTF,Sti,MakeEvent,McEvent,
geant,evout,IdTruth,tags,bbcSim,tofsim,emcY2,EEfs,EEss,
GeantOut,big,-dstout,fzin,-MiniMcMk,beamLine,clearmem,eemcDB,VFPPVnoCTB"

Run 6 bfc options with EEMC slow and fast simulators (EEfs EEss):

"trs fss y2006h Idst IAna l0 tpcI fcf ftpc Tree logger
ITTF Sti VFPPV bbcSim tofsim tags emcY2 EEfs EEss evout
-dstout IdTruth geantout big fzin MiniMcMk clearmem eemcDb beamLine sdt20060523"

Run 6 bfc options with EEMC gamma filter in the chain (FiltEemcGamma):

"FiltEemcGamma trs fss y2006h Idst IAna l0 tpcI fcf ftpc Tree logger
ITTF Sti VFPPV bbcSim tofsim tags emcY2 EEfs EEss evout
-dstout IdTruth geantout big fzin MiniMcMk clearmem eemcDb beamLine sdt20060523"

Note:
need to use a special macro RunEemcGammaFilterBfc.C
to run FiltEemcGamma with bfc

Scheduler xml template file examples

Official scheduler documentation

Schedule template file examples:

Catalog request to get official Gamma filtered Monte-Carlo files:

 get_file_list.pl -distinct -keys 'path,filename' -cond 'production=P10ii,path~pt_4-6gev/eemcgammafilt100_qcd/y2006h,filetype=MC_reco_MuDst,storage=nfs'

 or

 get_file_list.pl -distinct -keys 'path,filename' -cond 'production=P10ii,runnumber=2000010060,filetype=MC_reco_MuDst,storage=nfs

 where run number = 2000000000 + 10060

 

StEemcDstMaker (Emc dst event container - similar to the gamma maker code structure)

StEemcDstMaker is Eemc dst event container (similar to the gamma maker code structure).

Creates root tree from MuDst which stores the relevant information
for the photon-jet analysis

 

StEemcGammaFilter (Pythia level EEMC gamma filter)

StEemcGammaFilter is the Pythia level EEMC gamma filter.

The code is available at STAR/cvs: StEemcGammaFilter.h / StEemcGammaFilter.cxx
Basic algo description:

  • Loop over particles and search for the ones with energy higher than threshold and falls into the fiducial (in rapidity) volume
  • Search for clusters around each seed include tracks in eta and phi (detector Eta and Phi) space within the cone radious

Algo parameters:

mConeRadius - eta-phi cluster cone radius
mSeedThreshold - seed track threshold
mClusterThreshold - track cluster threshold
mEtaLow - lowerst rapidity cut
mEtaHigh - highest rapidity cut
mMaxVertex - vertex cut

Other parameters:

mCalDepth - calorimeter depth at which tracks are extrapolated
mMinPartEnergy - minimum particle energyto be included in the cluster
mHadronScale - down scale factor for hadrons (be careful when playing with this) No scaling by default
mFilterMode - filter mode: 0 - test mode (no event rejection), 1 - filter reject events
mPrintLevel - print level (0 - no output, 1 or 2 print some logs)

StEemcGammaFilterMaker (BFC level Endcap gamma filter)

StEemcGammaFilterMaker big full chain (BFC) Endcap gamma filter

Code is accessible in CVS/STRoot: StEemcGammaFilterMaker.h / StEemcGammaFilterMaker.cxx

Parameter can be stored in the data base (see eemcGammaFilterMakerParams.idl file),
but this is not enabled in the current implementation.

Basic idea of StEemcGammaFilterMaker algo:

  • Search for the 3x3 tower cluster with high tower and cluster Et above thresholds

Available parameters to vary:
Seed energy threshold (mSeedEnergyThreshold, GeV)
Cluster eT threshold (mClusterEtThreshold GeV)
Maximum z vertex (mMaxVertex, cm)

Macro to run StEemcGammaFilterMaker with BFC
(fixes the problem of loading StJetSkimEvent library):

StRoot/StFilterMaker/macros/RunEemcGammaFilterBfc.C

StGammaJetAnaEvent (event container which stores the gamma-jet sided residual info)

StGammaJetAnaEvent is an event container which stores the information
on gamma-jet candidates, output from shower shape fits and sided residual analysis
Also it has an information from the crude pi0 (multi-photon event) finder

 

StGammaJetDraw (applying analsyis cutsand generate histograms)

StGammaJetDraw - applying analsyis cuts
(such as photon isolation, photon and jet pt cuts, etc)
and generate pre-shower sorted histograms

StRoot/StGammaJetDraw/StGammaJetDraw:

StGammaJetEvent (simple container for the gamma-jet events)

StGammaJetEvent - container for the gamma-jet candidate events
(stores the information on selected di-jets from the jet finder and Eemdc Dst trees).

 

StGammaJetMaker (read jet, skim, and EemcDst trees and generate gamma-jet tree)

StGammaJetMaker read jet, skim, and EemcDst trees and select/write gamma-jet tree.