Cite as: Geiser, Achim; Dutta, Irene; Hirvonsalo, Harri; Sheeran, Bridget; (2017). Example code to produce the di-muon spectrum from a CMS 2011 or 2012 primary dataset. CERN Open Data Portal. DOI:10.7483/OPENDATA.CMS.D00J.UVB1
This simple analysis example is set up at Research level, i.e. it requires university-student-level programming experience. Minimal acquaintance with Linux and the ROOT analysis package (https://root.cern.ch/) as well as a basic text editor is needed.
Use this with 2011 or 2012 primary datasets containing muons.
Only the list of validated runs are accepted:
If you do not have the CERN Virtual Machine for 2011 CMS data installed, follow the instructions in step 1 at How to install a CERN Virtual Machine.
You can setup the example with the detailed instructions below, if you want to get insight on how to build your own analysis and what are the files needed. Install and run the Demo (demo analyzer) program as instructed in step 2 at 'How to install a CERN Virtual Machine'. Then
Demo/DemoAnalyzer/datasetsdirectory. Only the first index file is actually needed for initial running. You can also use other primary datasets containing muons from 2011 or 2012 as listed under "Use With" above.
src/DemoAnalyzer.ccwith the ones from this record and read the comments in
DemoAnalyzer.ccfor more details.
You should get an output file DoubleMu.root, which contains histograms for 10000 input events (a small subset of the data). These can be looked at using a ROOT Browser. The most interesting histogram is GM_mass_log. In order to compare it with the invariant-mass spectrum of di-muons in MUO-10-004, it should be viewed with the logy option. With 10000 events the comparison is poor, but the J/ψ (at log10(mass)=0.5), Υ (at log10(mass=0.98)) and Z (at log10(mass)=1.95, one event only) peaks should be visible. Furthermore, since the triggers for 2010 and 2011, as well as SingleMu and DoubleMu are different, the general shape of the distributions will differ somewhat depending on input.
To run over more or even the full statistics, edit the relevant parts of the Python file demoanalyzer_cfg.py (see comments therein, change max. number of events and/or insert the proper index and JSON files for 2011/2012 single/dimuon) and rerun. Add up the output histograms from different (nonoverlapping, i.e. use either DoubleMu or SingleMu, not both) input index files using Root tools.
In addition to these instructions, which guide you through the example in detail, a github repository (https://github.com/cms-opendata-analyses/DimuonSpectrum2011) based on this original code with minor modifications for direct download is available.
GNU General Public License (GPL) version 3