Implement a new energy loss module in JETSCAPE - TianyuDai/JETSCAPE-rhic-ags GitHub Wiki

Adding a new module to energy loss directory

Find CustomModule.cc, CustomModule.h, Lorenz.h in analysis/new_modules, and copy them to src/jet/

cp /rhic-ags-school/JETSCAPE-rhic-ags/analysis/* ../src/jet

Re-compile JETSCAPE

cd  /rhic-ags-school/JETSCAPE-rhic-ags/build
cmake –DUSE_MUSIC=ON ..
make -j

Create a new XML file, config/jetscape_user_custom.xml, add the new module into it

<jetscape>
  
  <nEvents> 10000 </nEvents>
  
  <!--  JetScape Writer Settings -->
  <outputFilename>brick_out</outputFilename>
  <JetScapeWriterAscii> on </JetScapeWriterAscii>

  <!-- Hard Process -->
  <Hard>
    <PGun>
      <pT>120</pT>
      <parID>21</parID>
    </PGun>
  </Hard>
  
  <!-- Hydro  Module  -->
  <Hydro>
    <Brick>
      <T>0.3</T>
    </Brick>
  </Hydro>
  
  <!--Eloss Modules -->
  <Eloss>

    <deltaT>0.01</deltaT>
    <maxT>1</maxT>

    <CustomModule>
      <name> CustomModule </name>
      <alpha_s> 0.3 </alpha_s>
      <pcut> 1.0 </pcut>
      <cutoff> 1.0 </cutoff>
    </CustomModule>

  </Eloss>
  
</jetscape>

Run JETSCAPE with the new energy loss module

cd  /rhic-ags-school/JETSCAPE-rhic-ags/build
./runJetscape ../config/jetscape_user_custom.xml

Run FinalStatePartons to process Ascii output file, and plot the final parton energy distribution. Remember to change the file name in plot_energy_distribution.py!

./FinalStatePartons custom_brick_out.dat ../../JETSCAPE-output/custom_final_state_partons.txt
cd ../analysis
python3 plot_energy_distribution.py
⚠️ **GitHub.com Fallback** ⚠️