As you can see in the output histograms, the JOs provided above contains two issues that need to be addressed.
- According to the file name, the leptoquark mass should be 1000 GeV. Modify the JOs to reflect this.
This is actually a very common problem in JOs! We often recommend parsing the JOs file name to set parameters to avoid this. Once you’ve figured out how to do it by hand, you can check here to see the recommended version:
Recommended LQ mass setting
```python # Handy function to get the name of the job options file from MCJobOptionUtils.JOsupport import get_physics_short my_physics_short = get_physics_short() # Mass of the leptoquark lq_mass = float( my_physics_short.split('_m')[1] ) # Gives 1000, units are GeV # Getting the mass this way ensures that the mass we set # actually matches the mass in the file name! ```
The physics process currently being generated by the JOs is a pair of leptoquarks that decay into an electron, muon, up quark and charm quark. For our analysis, we want the final state to have two same-flavor leptons (either two electrons or two muons). A simple solution would be to create two JOs, one for each decay mode, which is discussed in more detail in the JobOption Filters section. For this tutorial, however, we want to produce a single sample that includes both decay modes.
- Change the
generatecommand to produce all combinations of two electrons and/or muons. You will need to change the quark decay products to be inclusive as well using theqmulti-particle representation that has already been defined. For the inclusive lepton decay, you will need to define a new multi-particle representation. Call this new multi-particle representationlepton. Note that this will still include the e+mu channel, but this will be addressed in the next section.
Run the same Gen_tf.py command again as well as the derivation and validation
steps to make sure that the correct signal is now being produced.