This folder contains the following up-to-date analysis codes (to run after the latest version of the JPsiExtractor):
Run, for instance, the analysis with:
./run_ALL_MyAna.py --filelist filelists_date --version 1 --channel muonic --descr test
It will copy MyAnaMu.cpp as MyAna.cc, run the analysis for all the files in the filelists repository, and create rootfiles in:
Date/vX/MyAnaMu/
A log file names date_version_channel_description.log is stored in the LogMyAna repository.
If you wish to skim the extracted files, run:
./run_ALL_MyAna.py --filelist filelists_date --version X --channel skim --descr test
and then:
./createSkimFilelists.py --date 15AprYY --version X
Then, you can run:
./plotIt.py --date 15AprYY --version X --channel muonic
Using configPlotIt.template.yml, stacked histograms will be stored in:
Date/vX/PlotItMu/
This also work for the skimmed files.
If you have run the analysis for both semi-e and semi-mu channels, you can compare data with:
./xcheckElMu.py --date 15AprYY --version X
Using configXcheckElMu.yml, distributions will be stored in:
Date/vX/XcheckElMu
You can also get the number of events in:
Date/vX/NumberOfEvents.tex
after running:
./getNumberOfEvents.py --date 15AprYY --version X
hadd commands (for data, and TTJets* files) are executed when running:
./mergeChannels.py --date 15AprYY --version X
Rootfiles are stored in date/version/MyAnaAll.
Migration matrices can be obtained, as TH2F in a root file or as tables in a txt file, by running in root:
.L migration.C++
migration("date", "version", inBatch, isTrans)
This will produce plots, Chi2.txt, Matrices.txt, and Matrices.root in date/version/Migration/P or date/version/Migration/Pt.
In root:
.L sPlot.C++
sPlot("date", "version", inBatch)
Run, for instance, the analysis with:
./run_ALL_MyAna.py --filelist filelists_date --version 1 --channel muonic --descr test
It will copy MyAnaMu.cpp as MyAna.cc, run the analysis for all the files in the filelists repository, and create rootfiles in:
Date/vX/MyAnaMu/
A log file names date_version_channel_description.log is stored in the LogMyAna repository.
If you wish to skim the extracted files, run:
./run_ALL_MyAna.py --filelist filelists_date --version X --channel skim --descr test
and then:
./createSkimFilelists.py --date 15AprYY --version X
Then, you can run:
./plotIt.py --date 15AprYY --version X --channel muonic --decay semi
Using configPlotIt.template.yml, stacked histograms will be stored in:
Date/vX/PlotItMu/
To check the J/ψ enriched samples, run:
./plotItJpsi.py --date 15AprYY --version X --channel muonic --decay semi
Using configPlotItJpsi.template.yml, stacked histograms will be stored in:
Date/vX/PlotItJpsiMu/
This also work for the skimmed files.
First, create the rootfiles that contain histograms with the systematic errors as bin errors (and no bin content):
./fillSysForPlotIt.py --date 15AprYY --inf PUdown --sup PUup --out PU
Those files will be stored in date/SysForPlotIt/PU.
Then, run:
./plotItSys.py --date 15JulYY --version Ref --channel muonic --decay all
Using `configPlotItSys.template.yml stacked histograms will be store in:
Date/Ref/PlotItSysMu/
If you have run the analysis for both semi-e and semi-mu channels, you can compare data with:
./xcheckElMu.py --date 15AprYY --version X --decay semi
Using configXcheckElMu.template.yml, distributions will be stored in:
Date/vX/XcheckElMu
You can also get the number of events in:
Date/vX/NumberOfEvents.tex
after running:
./getNumberOfEvents.py --date 15AprYY --version X --decay semi
Start by morphing tW-channel in root
.L morphtW.C++
morphtW("date","version","channel")
Then, add MC bacground and signal contributions in root:
.L mergeMC.C++
mergeMC("date", "version", "electron")
mergeMC("date", "version", "muon")
In the date/version/MyAnaEl and date/version/MyAnaMu folders, rootfiles named All_mtop.root are created. They contain the binned and unbinned trilepton mass distributions, each contribution being normalized at the cross section.
hadd commands (for data, central J/ψ enriched, and All* files) are executed when running:
./mergeChannels.py --date 15AprYY --version X
Rootfiles are stored in date/version/MyAnaAll
This can be done in root:
.L calib.C++
calib("date", "version", "decay", binning)
Binning can be 0 for an unbinned fit, 2 for a binned fit, and 1 in between. Results are stored in the date/version/CalibEl, date/version/CalibMu, and date/version/CalibAll folders. Fits can be binned or unbinned.
This can be done in root:
.L simultaneousFit.C++
simultaneousFit("date", "version", "decay", blind, nEvtEl, nEvtMu)
If the blind parameter is set to false, the full stat central sample is used instead of data.
If nEvtEl or nEvtEl is not provided, the number of events in the data samples are automatically used for the toys generation.
Results are stored in the date/version/SimultaneousFitEl, date/version/SimultaneousFitMu, and date/version/SimultaneousFitAll folders.
The systematic uncertainty that is related to the fit calibration is computed. To compute other systematic uncertainties, run the analysis, mergeMCforSys.C, mergeChannels.py, and computeSys.C. Then add manually the uncertainty into sumErr.py and run it!
This can be done in root:
.L simpleFit.C++
simpleFit("date", "version", "decay", blind, nEvtEl, nEvtMu)
If the blind parameter is set to false, the full stat central sample is used instead of data.
If nEvtEl or nEvtEl is not provided, the number of events in the data samples are automatically used for the toys generation.
Results are stored in the date/version/SimpleFitEl, date/version/SimpleFitMu, and date/version/SimpleFitAll folders.
In root:
.L fitMJpsi.C++
fitMJpsi("date","version","channel","decay")
Results are store in the date/version/MJpsi directory.
To compare data and MC peak, use fitJpsiPeaks.C.
After merging channels, in root:
.L pairingStudies.C++
pairingStudies("date","version","decay")
Several J/ψ properties are considered according to their pairing (only matched events) for the central MSDecays sample. Results are stored in the date/version/PairingStudies folder.