Page History

...

• This is the slurm script for Rivanna (https://www.rc.virginia.edu/userinfo/hpc/slurm/)
• sbatch -A gives the Rivanna allocation.  Both spinquest and spinquest_standard should be available.
• You can use "squeue -u user", to check the status of your jobs (or use the "Active Jobs" tab on your UVA OpenOnDemand web page: https://ood.hpc.virginia.edu/pun/sys/dashboard/activejobs?jobcluster=all&jobfilter=user).
• Navigate to "/project/ptgroup/script"  and then implement the following command (this script by-default will skip bad-files/corrupted-files/etc.)
  $ ./merge_mc_prod.sh /scratch/<your_MC_output_file_location>
  
  This will create a folder inside /project/ptgrpup/spinquest/MC_merge_files" with the same name as your MC_output_folder
• Copy the "merged_trackQA_v2.root" file from the folder that you've newly generated by ./merge_mc_prod.sh and copy it to "/project/ptgroup/spinquest/MC_storage" location. Rename the file into <Channel>_<Vertex Origin>_<Number of accepted event>.root. for example: JPsi_Dump_300K.root

Additional Information about changing the Setup

The script that you need to change for any particular configuration is the Fun4Sim.C

Most of this configuration file you will not need to change, but depending on what you want to do here are some pointers

This is how you would select what physics generator you want to use (for normal DY or J/Psi use pythia):

  const bool gen_pythia8  = true;
  const bool gen_cosmic   = false;
  const bool gen_gun      = false;
  const bool gen_particle = false;
  const bool read_hepmc   = false;
  const bool gen_e906legacy = false;

Then you will need to switch to the appropriate cfg

   pythia8->set_config_file("phpythia8_DY.cfg");
   pythia8->set_config_file("phpythia8_Jpsi.cfg");

For making a single muons you can use multi particle gun

if(gen_particle) {

    PHG4SimpleEventGenerator *genp = new PHG4SimpleEventGenerator("MUP");

    //genp->set_seed(123);
    genp->add_particles("mu+", nmu);  // mu+,e+,proton,pi+,Upsilon
    if (legacyVtxGen) genp->enableLegacyVtxGen();
    else{
    genp->set_vertex_distribution_function(PHG4SimpleEventGenerator::Uniform,
        PHG4SimpleEventGenerator::Uniform,
        PHG4SimpleEventGenerator::Uniform);
    genp->set_vertex_distribution_mean(0.0, 0.0, -300);
    genp->set_vertex_distribution_width(100.0, 100.0, 400.0);
    genp->set_vertex_size_function(PHG4SimpleEventGenerator::Uniform);
    genp->set_vertex_size_parameters(0.0, 0.0);

    }

Here you can change the muon to mu+ or mu-, or change the or change the vertex interaction volume, for example this one creates a cylinder 1m radius around the beamline from 700 cm upstream of the dump to 100 cm inside the dump.

• Beam profile
  • The profile here means the distribution shape of beam protons in X and Y (= R). 
  • When "legacyVtxGen = true"
    • The shape is "Gaussian" at "R < 5*sigma" and "1 / R" at "R >= 5*sigma" (cf. https://github.com/E1039-Collaboration/e1039-core/blob/master/generators/E906LegacyVtxGen/SQPrimaryVertexGen.C#L290 )
    • The Gaussian width ("sigma" in cm) is defined by "SIGX_BEAM" and "SIGY_BEAM" in Fun4Sim.C.
  • When "legacyVtxGen = false"
    • The shape is defined by the following functions of event generators.
    • "set_vertex_distribution_function()" sets the shape to "Uniform" or "Gaussian".
    • "set_vertex_distribution_mean()" sets the mean of the distribution.
    • "set_vertex_distribution_width()" sets the half width in case of "Uniform" or "sigma" in case of "Gaussian".