2009 Endcap geometry update (ecalgeo6.g)

2009/2010 GEANT geometry file update summary (ecalgeo v6.0):

• Endcap geometry update and QA tests

Final updates (ecalgeo v6.2):

• Ecalgeo-v6.2: embedded LOW_EM cuts after TPC/EEMC overlap fix
• Ecalgeo-v6.2: embedded LOW_EM cuts in the calorimeter geometry
• Low EM timing study: 10KeV vs. 100KeV cut settings
• Low EM study: LOW_EM option, 100KeV cuts, and DCUTE=100KeV

See also:

• Endcap geometry description
• Endcap EMC project page (1998-2008 archives: MIT, IUCF)

Analysis flow chart

• Copy the StRoot directory from here (for example $my_dir)
   
• Change to the parent directory for StRoot:
  
  cd $my_dir

• 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);
 

Macros to run TMVA analysis over the EEMC photon-jet events

TMVAnalysis/StRoot/macros

• TMVAnalysis.C
• TMVAplots.C
• TMVArun.C
• TMVAcombinePlots.C
• TMVApplication.C
• TMVAcorrelations.C
• tmvaglob.C
• TMVApplication2.C
• TMVApplication1.C
• MultiLayerPerceptron.C
• TMVAcombinePreShowerPlots.C

Macros to run various makers

macros to run various makers for the gamma-jet analysis:

StRoot/macros/L2GammaTests

• DrawL2gammaTestPlots.C
• ReadMuDst.C
• dumpL2info.C

StRoot/macros/draw

• DrawAllPtPionCuts.C
• DrawHistFromRootFile.C
• DrawHistRatios.C
• DrawImagesForHistsInFile.C
• DrawLoopHist.C
• DrawMCtoDataRatios.C
• DrawTuneComparison.C
• DrawscaleGammaJetBGsimu.C
• compareShowerShapesWithPibero.C
• getchain.C

StRoot/macros/drawFilterQA

• ComparePyTunes.C
• ComparePyTunesXsection.C
• DrawTimingPlots.C
• drawAll.C
• draw_candidatePt_vsFilter.C
• draw_filterQA_plots.C
• draw_reject_plots.C

StRoot/macros/drawL2perRun

• DrawOutputHist.C
• drawL2eventsPerRun.C
• jetSkimTreeTest.C

Macros to run various makers:

• RunCreateFilterTree.C
• RunEemcMakerMC.C
• RunGammaJetAnaMC.C
• RunGammaJetDrawMC.C
• RunGammaJetMakerMC.C
• RunGammaMaker2006.C
• RunJetFinder2006mc.C

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:

• Schedule-template_FullChain.xml -
  advanced template to submit jobs for all (or selected) makers in one job (run over private Monte-Carlo files)
• Schedule-template_FullChain_official.xml -
  advanced template to submit jobs for all (or selected) makers in one job (for official Monte-Carlo)
• Schedule-template_GammaJetDraw.xml - submit jobs for RunGammaJetDrawMC.C for private Monte-Carlo
• Schedule-template_GammaJetDraw_official.xml - submit jobs for RunGammaJetDrawMC.C for private Monte-Carlo
• Schedule-template_candidate.xml - submit jobs for RunGammaMaker.C / RunCreateFilterTree.C

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

• StEemcDstEventMaker.h / StEemcDstEventMaker.cxx - gets pointer to the MuDst information
• StEemcDstMaker.h / StEemcDstMaker.cxx - init EemcDstRawMaker
• StEemcDstRawMaker.h / StEemcDstRawMaker.cxx - copy the information from MuDst with GetEndcap() and GetTracks() methods
• StEemcDstStrip.h / StEemcDstStrip.cxx - container for the EEMC SMD strip
• StEemcDstTower.h / StEemcDstTower.cxx - container for the  EEMC tower and pre/post shower layers
• StEemcDstTrack.h / StEemcDstTrack.cxx - container for the TPC track
• StGeantRecord.h / StGeantRecord.cxx - copy the relevant GEANT/Pythia information from geant.root
• StGeantRecordMaker.h / StGeantRecordMaker.cxx - TBD

 

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

StGammaJetAna (maker which does the shower shape analysis)

StGammaJetAna is the main maker which actually do the shower shape analysis.

StGammaJetAna.h / StGammaJetAna.cxx
 

StRoot/StGammaJetAna/StGammaJetAna:

• CheckNeighbors.h - short method to get the adjacent towers
• Cluster.h -
• Constants.h - parameters
• CreateHistLists.h -
• Cuts.h
• FillAnaTree.h
• FillHists.h
• FindMaxEnergySector.h
• FindSmdCluster.h - redundant
• FindSmdGammaCluster.h
• FitSMDstripsGaus.h - redundant
• FitSMDstripsSidedResidual.h
• GetEemcEnergiesForMaxTower.h
• GetPolarization.h
• Init.h
• Make.h
• SetPreshowerBin.h
• getNextEemcTowerId.h

 

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

• StEemcSector.h / StEemcSector.cxx
• StGammaJetAnaEvent.h / StGammaJetAnaEvent.cxx
• StGammaJetGamma.h / StGammaJetGamma.cxx
• StGammaJetPion.h / StGammaJetPion.cxx

 

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

• StGammaJetDraw.h / StGammaJetDraw.cxx

StRoot/StGammaJetDraw/StGammaJetDraw:

• Constants.h
• Constants_1Dhists.h
• Constants_1Dhists_Geant.h
• Constants_1Dprofiles.h
• Constants_2Dhists.h
• Constants_2Dhists_Geant.h
• Constants_2Dprofiles.h
• Constants_partonicPtBins.h
• FillCandidateTxtFile.h - dumps gamma-jet info in txt file in format request by Hal Spinka
• GetJetDetEta.h - short code calculates jet detector rapidity
• GetPolarization.h
• Init.h
• Make.h
• ReadSpinInfo.h
• SetDrawOptions.h
• SetPartonicWeight.h - sets partonic pt weigts for Monte-Carlo
• createResidualHists.h
• drawEemc.h
• drawResidual.h - applying analysis cuts and fill histograms

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).

• StEmcJet.h / StEmcJet.cxx
• StEmcJetTower.h / StEmcJetTower.cxx
• StGammaJetEvent.h / StGammaJetEvent.cxx
• diJetTreeConstants.h

 

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.

• StGammaJetMaker.h / StGammaJetMaker.cxx
   
• StRoot/StGammaJetMaker/StGammaJetMaker subdirectory has:
  
  CheckNeighbors.h - simple code to check if the towers are adjacent
  Constants.h - list of constant used by the StGammaJetMaker
  getNextEemcTowerId.h - get the id of the adjacent tower
  Make.h - pice of code which reads the trees and fills the EemcDst event