Page History

...

   1. If you want to use Pythia8, you can load the config file using this function. pythia8->set_config_file("phpythia8_DY.cfg") or pythia8->set_config_file("phpythia8_Jpsi.cfg"); 
depending on your work.  
   2. In the configuration file, you can change the configuration based on your need.  For example, you can set the mass cut at Pythia8 by putting this line: PhaseSpace:mHatMin = 4.0. 
   3. In addition to the pythia configuration file, you can add more variables using PHPy8ParticleTrigger. A basic setup has shwn below, but you can find more function here.

// pythia8
  if(gen_pythia8) {    
    PHPythia8 *pythia8 = new PHPythia8();
    //pythia8->Verbosity(99);
    pythia8->set_config_file("phpythia8_DY.cfg");
    if(SQ_vtx_gen) pythia8->enableLegacyVtxGen();
    else{
      pythia8->set_vertex_distribution_mean(0, 0, target_coil_pos_z, 0);
    } 
    pythia8->set_embedding_id(1);
    se->registerSubsystem(pythia8);

    pythia8->set_trigger_AND();

    PHPy8ParticleTrigger* trigger_mup = new PHPy8ParticleTrigger();
    trigger_mup->AddParticles("-13");
    //trigger_mup->SetPxHighLow(7, 0.5);
    //trigger_mup->SetPyHighLow(6, -6);
    trigger_mup->SetPzHighLow(120, 10);
    pythia8->register_trigger(trigger_mup);

    PHPy8ParticleTrigger* trigger_mum = new PHPy8ParticleTrigger();
    trigger_mum->AddParticles("13");
    //trigger_mum->SetPxHighLow(-0.5, 7);
    //trigger_mum->SetPyHighLow(6, -6);
    trigger_mum->SetPzHighLow(120, 10);
    pythia8->register_trigger(trigger_mum);
  }
  }

gen_e906legacy:

This generator is adopted from the SeaQuest experiment. If you interested in generating dimuon mass, especially in the low mass regions, you can use the script, because generating 
more statistics using gen_pythia8 is computationally much more expensive.  And If you use gen_e906legacy generator, you will need to save the event weight information, since 
each event has a generator weight.

...