HFT simulation--Hijing + HF simulation
Instructions from Yifei Zhang:
https://drupal.star.bnl.gov/STAR/blog/yfzhang/instruction-hft-simulation
1. inserted Ds+/Ds-
code position: RCF
/star/u/xlchen/work/lbl/hft/Run14/Ds_Simu/HijinSimu
first : read README
1.1 makeFZ/------------in 1 e) of Yifei's
scripts to produce fz files for hijing simulations with embedded HF decays.
If need add in new defined particles, one need check out local pams/sim/gstar/gstar_part.g to local and modify this file accordingly
And also the decay channel
| >run.sh [index] or >submit.sh [index_begin] [index_end] |
Reference class: StarPrimaryMaker.
http://www.star.bnl.gov/webdata/dox/html/classStarPrimaryMaker.html#a065785692580b6725e1a6190a40c389e
example: primary->SetSigma(0.01,0.01,2.0,0);
| void StarPrimaryMaker::SetSigma | ( | Double_t | sx, | |
| Double_t | sy, | |||
| Double_t | sz, | |||
| Double_t | rho = 0 | |||
| ) | [inline] |
Set the smearing in the x, y and z vertex positions.
- Parameters:
-
sx is the smearing in x sy is the smearing in y sz is the smearing in z rho is the correlation between x and y
particle list：vim $ROOTSYS/etc/pdg_table.txt
class: StarParticleData provides an interface to the data provided in the TDatabasePDG class from ROOT. We extend the available list of particles,
including some standard heavy ions which are lacking in the ROOT table, and provide the capability to alias names to individual particles.
1.2 Reco/
- Reco: Contains two directories:
-StRoot_sti15: Standard simulation library. Only contains StPxlSimMaker required to load Pile Up
-StRoot_detEFF_SC: Sti15 tracking library plus corrections accounting for detector (PXL+IST)
detector efficiency plus TPC space charge effect.
If you want to use StRoot_detEFF_SC, you need copy this folder to StRoot.
Then you can compile the code with command 'cons'. (Note: cons can only compile the file named StRoot)
See the BigFullChain: vim $STAR/StRoot/StBFChain/BigFullChain.h
fzin: read gstar fz-file
in: In put tree --root file
Know how many events in the input file (.fzd, .root), define the same events loop number --in runPUReco.csh -- chain->EventLoop($start,$end);
Test the code and submit jobs by
| Test: >./runPUReco.csh [FZfile_index] [pileup_index] |
| submit: >./submitPUReco.csh [FZfile_index_begin] [FZfile_index_end] |
The FZfile_index should be consistent with the file produced in 1.1
pileup_index can select one in /star/u/xlchen/work/lbl/hft/Run14/Ds_Simu/HijinSimu/pileup
Clear the useless file
| >./clear.sh |
When you want to use the new Star Library, you just need starver new star library, then compile by cons.
for example, this time I use:
| >starver SL15e >cons |
But, maybe some errors occure, because some functions may be different at different versions.
Pay attention to load correct parameters and hit errors by checking log files.
e.g.
2)To check if the pileup hits are loaded, one need find in log file:
This reasonably indicates that pileup hits are correctly added.
3)
To check if the detector geometry and are included and activated, one can check:This means IST, PXL, SSD are correctly included.
4)
hit errorsThis is a correct example that PXL hits errors are 12 micron in x and 12 micron in z. e.g. sqrt(1.44e-06).
IST hits errors are 170 micron in x and 1700 micron in z.
others
1.3 ReadMuDst-----deal with the *.MuDst.root from Reco
- ReadMuDst:
-lMuDst.C: Loads MuDSt's
-RecoSimGlobal.C: Contains sample code to load RC and MC information from MuDst. Contains example on associating mc and rc tracks,
and sample loops, etc.
1)Test
| >./Test.sh [index] [index] is the number from Reco |
2)Submit jobs:ReadMuDst/submit
| >cd submit >./getlist.sh >./SubmitJobs.sh |
3)merge the output root file by hadd.sh: Reco/submit
| >./hadd.sh [$1] [$2] [$3] [$4] [$5] the number of parameters can be 0-5 |
int myhadd(TString WorkPath,Int_t TESTMODE=1,TString SourceDir="out",TString DestinationDir="Hist",const Int_t Nfilesperjob=10,TString options="")
pwd $1 $2 $3 $4 $5
$1: 1---test mode, default value, test if the macro can be runned right; 0------run mode
$2: source root files directory
$3: out put root file directory
$4: every $4 root files to merge one root file
$5: Options : please look at hadd options. see: hadd --help
4)ReadMuDst/ana
| >./RunMySelector.sh $1 $1: dir/name.root or name.root-->current direction to produce *.mini.root in dir/ |
| >./plots.sh $1 $2 $1: dir/name.mini.root $2: decay mode 0--D_{s}#rightarrow#phi(1020)#pi D_{s}#rightarrowk#bar{k^{*0}}(892) 1--D_{s}#rightarrow#phi(1020)#pi 2--D_{s}#rightarrowk#bar{k^{*0}}(892) |
2. D0
Similar to Ds
3. Supplyment
3.1 Mc particle id, mass, life time, etc
To see the star simulation particle characters, by printing the list of particles available in STARSIM:
| >starsim starsim>make sim/gstar starsim>gprint part |
StMuMcTrack *mcTrack;
mcTrack->GePid();
3.2 Particle name used in starsim function (makefz) and decay channel, pid, etc
| >vim $ROOTSYS/etc/pdg_table.txt |
3.3 Add new particle(makefz/)
If need add in new defined particles, one need check out pams/sim/gstar/gstar_part.g to local and modify this file accordingly.
Like _phi_to_KK, code=99, pdg=333, you can also define the particle K*bar0(code=892,pdg=-313), and reconstruction the
decay channel Ds+ --> K*bar0 k+ --> k- pi+ k+
- xlchen's blog
- Login or register to post comments