Back to home page

EIC code displayed by LXR

 
 

    


Warning, /epic.github.io/_physics/jets_hf.md is written in an unsupported language. File is not indexed.

0001 ---
0002 title: Jets and Heavy Flavor
0003 description: Jets and Heavy Flavor
0004 name: jets_hf
0005 layout: default
0006 ---
0007 
0008 {% include layouts/wg_top.md %}
0009 
0010 __Simulations__
0011 * Repository for event generation
0012   * https://github.com/eic/HFsim-PYTHIA
0013   * https://github.com/eic/HFsim-BeAGLE 
0014 * D<sup>0</sup> samples
0015   * e+p: PYTHIA v8.306
0016     * 18x275
0017       * Each event contains at least one D<sup>0</sup> that can decay through any channels
0018       * Location: ``` epic_craterlake/SIDIS/D0_ABCONV/pythia8.306-1.0/18x275/hiDiv  ```
0019     * 10x100
0020       * Each event contains at least one D<sup>0</sup> that decays into π+K pair, and the decay daughters are within pseudorapidity of &plusmn;3.5
0021       * Location:
0022         * Q<sup>2</sup> > 1: ``` epic_craterlake/SIDIS/D0_ABCONV/pythia8.306-1.1/10x100/q2_1/hiDiv  ```
0023         * Q<sup>2</sup> > 100: ``` epic_craterlake/SIDIS/D0_ABCONV/pythia8.306-1.1/10x100/q2_100/hiDiv  ```
0024     * 5x41
0025       * Each event contains at least one D<sup>0</sup> that can decays into π+K pair
0026       * Location: ``` epic_craterlake/SIDIS/D0_ABCONV/HFsim-PYTHIA/pythia8.306-1.2/ep/5x41/q2_1to10000/hiDiv  ```
0027   * e+Au: BeAGLE v1.03.01
0028     * Each event contains at least one D<sup>0</sup> that decays into π+K pair
0029     * Q<sup>2</sup> > 1
0030     * Location:
0031       * 10x100: ``` epic_craterlake_without_zdc/SIDIS/D0_ABCONV/HFsim-BeAGLE/BeAGLE1.03.01-1.0/eAu/10x100/q2_1to10000 ``` 
0032       * 5x41: ``` epic_craterlake_without_zdc/SIDIS/D0_ABCONV/HFsim-BeAGLE/BeAGLE1.03.01-1.0/eAu/5x41/q2_1to10000 ``` 
0033 * L<sub>c</sub> samples
0034   * e+p
0035     * 18x275
0036       * Each event contains at least one L<sub>c</sub> that can decay through any channels
0037       * Location: ``` epic_craterlake/SIDIS/Lambda_ABCONV/pythia8.306-1.0/18x275/hiDiv  ```
0038     * 10x100
0039       * Each event contains at least one L<sub>c</sub> that decays to a π+K+p channel
0040       * Location: ``` epic_craterlake/SIDIS/Lc_ABCONV/HFsim-PYTHIA/pythia8.306-1.2/ep/10x100/q2_1to10000/hiDiv  ```
0041   * e+Au: BeAGLE v1.03.01
0042     * Each event contains at least one L<sub>c</sub> that decays through π+K+p channel
0043     * Q<sup>2</sup> > 1
0044     * Location:
0045       * 10x100: ``` epic_craterlake_without_zdc/SIDIS/Lc_ABCONV/HFsim-BeAGLE/BeAGLE1.03.01-1.0/eAu/10x100/q2_1to10000 ```
0046           
0047 __Analysis macros__
0048 * Jet reader: https://github.com/eic/snippets/tree/main/JetsAndHF/jetReaderExamples
0049 * D<sup>0</sup> reconstruction based on Helix swimming: https://github.com/eic/snippets/tree/main/JetsAndHF/HF/D0_helix
0050 
0051 
0052 __EICrecon Jet Reconstruction Documentation__
0053 
0054 There are currently six jet collections in the default output of EICrecon:
0055 * _GeneratedJets_, which are reconstructed using charged and neutral particles from the ```GeneratedParticles``` collection(final-state primaries) using the anti-kT algorithm with R = 1.0;
0056 * _GeneratedChargedJets_, which are reconstructed using charged particles from the ```GeneratedParticles``` collection using the anti-kT algorithm with R = 1.0;
0057 * _GeneratedCentauroJets_, which are reconstructed using charged and neutral particles from the ```GeneratedParticles``` collection using the Centauro algorithm with R = 0.8;
0058 * _ReconstructedJets_, which are reconstructed using charged and neutral particles (tracks and EMCal clusters) from the ```ReconstructedParticles``` collection using the anti-kT algorithm with R = 1.0;
0059 * _ReconstructedChargedJets_, which are reconstructed using charged particles (tracks) from the ```ReconstructedChargedParticles``` collection using the anti-kT algorithm with R = 1.0;
0060 * _ReconstructedCentauroJets_, which are reconstructed using charged and neutral particles from the ```ReconstructedParticles``` collection using the Centauro algorithm with R = 0.8;
0061 
0062 
0063 Jet reconstruction is handled by [FastJet3](https://fastjet.fr). User-configurable parameters and their default values are listed in the table below. Jets are stored as ```edm4eic::ReconstructedParticle``` objects, and since this data type was intended to describe particles rather than jets, additional information like the jet area currently is not stored..
0064 
0065 | Parameter | Name | Default Value |
0066 | --- | --- | --- |
0067 | Jet algorithm | jetAlgo | anti-kT |
0068 | Jet recombination scheme | recombScheme | E-scheme |
0069 | Jet resolution parameter | rJet | 1 |
0070 | Exponent for generalized kT algorithms | pJet | -1 |
0071 | Minimum constituent pT | minCstPt | 0.2 GeV/c |
0072 | Maximum constituent pT | maxCstPt | 100 GeV/c |
0073 | Minimum jet pT | minJetPt | 1 GeV/c |
0074 | Area type | areaType | active |
0075 | Maximum ghost rapidity | ghostMaxRap | 3.5 |
0076 | No. of repeated ghosts | numGhostRepeat | 1 |
0077 | Area per ghost | ghostArea | 0.01 |
0078 | Contributed jet algorithm | jetContribAlgo | Centauro |
0079 
0080 These parameters can be adjusted at runtime, illustrated by the code snippet below. To use a contributed fastjet algorithm, the user should set ```jetAlgo``` to ```plugin_algorithm```.
0081 ```
0082 eicrecon \
0083   -Preco:ReconstructedChargedJets:jetAlgo=ee_genkt_algorithm \
0084   -Preco:ReconstructedChargedJets:rJet=0.8 \
0085   -Preco:ReconstructedChargedJets:pJet=-1.0 \
0086   <input edm4hep file>
0087 ```