Merge pull request #97 from JETSCAPE/3DGlauber-MUSIC-FOR_RC2

3d glauber music for rc2
JETSCAPE · Feb 28, 2023 · 4c81992 · 4c81992
2 parents 568be55 + be73737
commit 4c81992
Show file tree

Hide file tree

Showing 11 changed files with 516 additions and 48 deletions.
diff --git a/.github/workflows/test-events-Pythia-Isr-Dat.yaml b/.github/workflows/test-events-Pythia-Isr-Dat.yaml
@@ -33,17 +33,19 @@ jobs:
       with:
         path: ${{ github.event.repository.name }}
 
-    - name: Download 3DGlouber
+    - name: Download 3DGlauber
       run: |
         cd ${GITHUB_WORKSPACE}/${REPO_NAME}/external_packages
         ./get_3dglauber.sh
+        ./get_music.sh
+        ./get_iSS.sh
 
     - name: Build Application
       run: |
         cd ${GITHUB_WORKSPACE}/${REPO_NAME}
         mkdir build
         cd build
-        cmake .. -DCMAKE_CXX_STANDARD=14 -DUSE_3DGlauber=ON
+        cmake .. -DCMAKE_CXX_STANDARD=14 -DUSE_3DGlauber=ON -DUSE_MUSIC=ON -DUSE_ISS=ON
         make -j2
 
     - name: Checkout TEST-EXAMPLES

diff --git a/.github/workflows/test-events-Pythia-Isr-Hadron.yaml b/.github/workflows/test-events-Pythia-Isr-Hadron.yaml
@@ -33,17 +33,19 @@ jobs:
       with:
         path: ${{ github.event.repository.name }}
 
-    - name: Download 3DGlouber
+    - name: Download 3DGlauber
       run: |
         cd ${GITHUB_WORKSPACE}/${REPO_NAME}/external_packages
         ./get_3dglauber.sh
+        ./get_music.sh
+        ./get_iSS.sh
 
     - name: Build Application
       run: |
         cd ${GITHUB_WORKSPACE}/${REPO_NAME}
         mkdir build
         cd build
-        cmake .. -DCMAKE_CXX_STANDARD=14 -DUSE_3DGlauber=ON
+        cmake .. -DCMAKE_CXX_STANDARD=14 -DUSE_3DGlauber=ON -DUSE_MUSIC=ON -DUSE_ISS=ON
         make -j2
 
     - name: Checkout TEST-EXAMPLES

diff --git a/.github/workflows/test-events-Pythia-Isr-Parton.yaml b/.github/workflows/test-events-Pythia-Isr-Parton.yaml
@@ -33,17 +33,19 @@ jobs:
       with:
         path: ${{ github.event.repository.name }}
 
-    - name: Download 3DGlouber
+    - name: Download 3DGlauber
       run: |
         cd ${GITHUB_WORKSPACE}/${REPO_NAME}/external_packages
         ./get_3dglauber.sh
+        ./get_music.sh
+        ./get_iSS.sh
 
     - name: Build Application
       run: |
         cd ${GITHUB_WORKSPACE}/${REPO_NAME}
         mkdir build
         cd build
-        cmake .. -DCMAKE_CXX_STANDARD=14 -DUSE_3DGlauber=ON
+        cmake .. -DCMAKE_CXX_STANDARD=14 -DUSE_3DGlauber=ON -DUSE_MUSIC=ON -DUSE_ISS=ON
         make -j2
 
     - name: Checkout TEST-EXAMPLES

diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -104,10 +104,10 @@ endif (USE_SMASH)
 message("Compiler and Linker flags ...")
 set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -fPIC -pipe -Wall -std=c++11")
 ## can turn off some warnings
-set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wno-reorder -Wno-unused-variable -Wno-deprecated-copy -Wno-unused-value")
+set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wno-reorder -Wno-unused-variable -Wno-inconsistent-missing-override -Wno-unused-private-field")
 ## Then set the build type specific options. These will be automatically appended.
 set(CMAKE_CXX_FLAGS_DEBUG "-g -O0")
-set(CMAKE_CXX_FLAGS_RELEASE "-O3 -w")
+set(CMAKE_CXX_FLAGS_RELEASE "-O3")
 
 ## can turn on debugging information
 # set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -g")
@@ -117,7 +117,7 @@ set(CMAKE_CXX_FLAGS_RELEASE "-O3 -w")
 if (APPLE)
     message( STATUS "Apple : " ${CMAKE_HOST_SYSTEM})
     if ("${CMAKE_CXX_COMPILER_ID}" STREQUAL "AppleClang")
-        set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -stdlib=libc++ -Wno-inconsistent-missing-override -Wno-unused-private-field -Wno-delete-non-abstract-non-virtual-dtor")
+        set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -stdlib=libc++" )
     elseif("${CMAKE_CXX_COMPILER_ID}" STREQUAL "GNU")
         set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wno-sign-compare -Wno-unused-but-set-variable -Wno-parentheses -fext-numeric-literals")
     endif()
@@ -477,20 +477,25 @@ target_link_libraries(FinalStateHadrons JetScape )
 add_executable(FinalStatePartons ./examples/FinalStatePartons.cc)
 target_link_libraries(FinalStatePartons JetScape )
 
+add_executable(PythiaBrickTest ./examples/custom_examples/PythiaBrickTest.cc)
+target_link_libraries(PythiaBrickTest JetScape )
+
 if(USE_3DGlauber)
   add_executable(PythiaIsrTest ./examples/PythiaIsrTest.cc)
   target_link_libraries(PythiaIsrTest JetScape )
+  add_executable(PythiaIsrMUSIC ./examples/custom_examples/PythiaIsrMUSIC.cc)
+  target_link_libraries(PythiaIsrMUSIC JetScape )
 endif (USE_3DGlauber )
 
-#add_executable(PythiaBrickTest ./examples/custom_examples/PythiaBrickTest.cc)
-#target_link_libraries(PythiaBrickTest JetScape )
-
 #add_executable(PythiaBDMTest ./examples/custom_examples/PythiaBDMTest.cc)
 #target_link_libraries(PythiaBDMTest JetScape )
 
 #add_executable(PythiaNoBDMTest ./examples/custom_examples/PythiaNoBDMTest.cc)
 #target_link_libraries(PythiaNoBDMTest JetScape )
 
+#add_executable(PythiaBDMTest_Hydro_Cascade ./examples/custom_examples/PythiaBDMTest_Hydro_Cascade.cc)
+#target_link_libraries(PythiaBDMTest_Hydro_Cascade JetScape )
+
 if ( USE_SMASH )
     add_executable(SMASHTest ./examples/custom_examples/SMASHTest.cc)
     target_link_libraries(SMASHTest JetScape )

diff --git a/config/jetscape_main.xml b/config/jetscape_main.xml
@@ -84,6 +84,17 @@
 
     <!-- Options to read initial conditions from saved file  -->
     <initial_profile_path>../examples/test_hydro_files</initial_profile_path>
+
+    <MCGlauber> 
+        <ran_seed_MCGla> -1 </ran_seed_MCGla>
+        <rapidity_loss_method> 4 </rapidity_loss_method>
+        <ylossParam4At2> 1.666 </ylossParam4At2>
+        <ylossParam4At4> 2.10 </ylossParam4At4>
+        <ylossParam4At6> 21.20 </ylossParam4At6>
+        <ylossParam4At10> 2.30 </ylossParam4At10>
+        <ylossParam4var> 0.555 </ylossParam4var>
+        <remnant_energy_loss_fraction> 0.5 </remnant_energy_loss_fraction>
+    </MCGlauber>
   </IS>
 
   <!-- Hard Process -->
@@ -260,22 +271,27 @@
     <MUSIC>
       <name>MUSIC</name>
       <MUSIC_input_file>music_input</MUSIC_input_file>
+      <InitialProfile>13</InitialProfile>
+      <beastMode>0</beastMode>
       <Initial_time_tau_0>0.5</Initial_time_tau_0>
-      <output_evolution_to_file>1</output_evolution_to_file>
+      <output_evolution_to_file>0</output_evolution_to_file>
+      <store_hydro_info_in_memory>0</store_hydro_info_in_memory>
       <shear_viscosity_eta_over_s>0.08</shear_viscosity_eta_over_s>
-      <T_dependent_Shear_to_S_ratio>0</T_dependent_Shear_to_S_ratio>
-      <eta_over_s_T_kink_in_GeV>0.16</eta_over_s_T_kink_in_GeV>
-      <eta_over_s_low_T_slope_in_GeV>0.</eta_over_s_low_T_slope_in_GeV>
-      <eta_over_s_high_T_slope_in_GeV>0.</eta_over_s_high_T_slope_in_GeV>
-      <eta_over_s_at_kink>0.08</eta_over_s_at_kink>
-      <temperature_dependent_bulk_viscosity>0</temperature_dependent_bulk_viscosity>
-      <zeta_over_s_max>0</zeta_over_s_max>
-      <zeta_over_s_T_peak_in_GeV>0.18</zeta_over_s_T_peak_in_GeV>
+      <T_dependent_Shear_to_S_ratio>3</T_dependent_Shear_to_S_ratio>
+      <shear_viscosity_3_at_kink>0.16</shear_viscosity_3_at_kink>
+      <shear_viscosity_3_low_T_slope_in_GeV>0</shear_viscosity_3_low_T_slope_in_GeV>
+      <shear_viscosity_3_high_T_slope_in_GeV>0</shear_viscosity_3_high_T_slope_in_GeV>
+      <temperature_dependent_bulk_viscosity>3</temperature_dependent_bulk_viscosity>
+      <bulk_viscosity_3_max>0.06</bulk_viscosity_3_max>
+      <zeta_over_s_T_peak_in_GeV>0.17</zeta_over_s_T_peak_in_GeV>
       <zeta_over_s_width_in_GeV>0.02</zeta_over_s_width_in_GeV>
-      <zeta_over_s_lambda_asymm>0.0</zeta_over_s_lambda_asymm>
+      <zeta_over_s_lambda_asymm>0.1</zeta_over_s_lambda_asymm>
       <Include_second_order_terms>1</Include_second_order_terms>
-      <freezeout_temperature>-1</freezeout_temperature>
+      <freezeout_temperature>0.3</freezeout_temperature>
       <Perform_CooperFrye_Feezeout>0</Perform_CooperFrye_Feezeout>
+      <string_source_sigma_x>0.5</string_source_sigma_x>
+      <string_source_sigma_eta>0.5</string_source_sigma_eta>
+      <stringPreEqFlowFactor>0.13</stringPreEqFlowFactor>
     </MUSIC>
 
     <!-- CLVisc  -->
@@ -366,3 +382,4 @@
     </SMASH>
   </Afterburner>
 </jetscape>
+
diff --git a/config/jetscape_user_iMATTERMCGlauberMUSIC.xml b/config/jetscape_user_iMATTERMCGlauberMUSIC.xml
@@ -0,0 +1,100 @@
+<?xml version="1.0"?>
+
+<jetscape>
+
+  <enableAutomaticTaskListDetermination> false </enableAutomaticTaskListDetermination>
+
+  <setReuseHydro> false </setReuseHydro>
+  <!-- <nReuseHydro> 10 </nReuseHydro> -->
+
+  <vlevel> 4 </vlevel>
+
+  <nEvents> 10 </nEvents>
+
+  <Random>
+    <seed>1</seed>
+  </Random>
+
+  <JetScapeWriterAscii> on </JetScapeWriterAscii>
+
+  <!-- Hard Process -->
+  <Hard>
+    <PythiaGun>
+      <pTHatMin>4</pTHatMin>
+      <pTHatMax>-1</pTHatMax>
+      <eCM>5020</eCM>
+      <LinesToRead>
+        PhaseSpace:bias2Selection = on
+        PhaseSpace:bias2SelectionPow = 4
+        PhaseSpace:bias2SelectionRef = 10
+        SigmaProcess:alphaSvalue = 0.16
+        MultipartonInteractions:alphaSvalue = 0.16
+        MultipartonInteractions:pTmin = 7.0
+      </LinesToRead>
+      <useHybridHad>1</useHybridHad>
+    </PythiaGun>
+  </Hard>
+
+  <!--Preequilibrium Dynamics Module -->
+  <Preequilibrium>
+    <NullPreDynamics> </NullPreDynamics>
+  </Preequilibrium>
+
+  <!-- Initial condition  Module  -->
+  <IS>
+    <MCGlauber> 
+        <ran_seed_MCGla> -1 </ran_seed_MCGla>
+        <ylossParam4At2> 1.666 </ylossParam4At2>
+        <ylossParam4At4> 2.10 </ylossParam4At4>
+        <ylossParam4At6> 2.20 </ylossParam4At6>
+        <ylossParam4At10> 2.30 </ylossParam4At10>
+        <ylossParam4var> 0.555 </ylossParam4var>
+    </MCGlauber>
+  </IS>
+
+  <!-- Hydro  Module  -->
+  <Hydro>
+        <MUSIC>
+            <InitialProfile>13</InitialProfile>
+            <string_source_sigma_x>0.5</string_source_sigma_x>
+            <string_source_sigma_eta>0.5</string_source_sigma_eta>
+            <stringPreEqFlowFactor>0.13</stringPreEqFlowFactor>
+            <output_evolution_to_file>0</output_evolution_to_file>
+            <store_hydro_info_in_memory>0</store_hydro_info_in_memory>
+            <T_dependent_Shear_to_S_ratio>3</T_dependent_Shear_to_S_ratio>
+            <temperature_dependent_bulk_viscosity>3</temperature_dependent_bulk_viscosity>
+            <bulk_viscosity_3_max>0.1</bulk_viscosity_3_max>
+            <shear_viscosity_3_at_kink>0.13</shear_viscosity_3_at_kink>
+            <freezeout_temperature>0.3</freezeout_temperature>
+        </MUSIC>
+  </Hydro>
+
+  <!--Eloss Modules -->
+  <Eloss>
+    <maxT>50</maxT>
+    <Matter>
+      <Q0> 2.0 </Q0>
+      <in_vac> 1 </in_vac>
+      <useHybridHad>1</useHybridHad>
+      <vir_factor> 0.25 </vir_factor>
+      <recoil_on> 0 </recoil_on>
+      <broadening_on> 0 </broadening_on>
+      <brick_med> 0 </brick_med>
+    </Matter>
+  </Eloss>
+
+  <!-- Jet Hadronization Module -->
+  <JetHadronization>
+    <name>colorless</name>
+  </JetHadronization>
+
+  <!-- Particlization Module  -->
+  <SoftParticlization>
+    <iSS>
+        <number_of_repeated_sampling>100</number_of_repeated_sampling>
+        <Perform_resonance_decays>1</Perform_resonance_decays>
+    </iSS>
+  </SoftParticlization>
+
+</jetscape>
+