https://github.com/marcodeltutto/evtwgtanalysis
Code to study GENIE model uncertainties at MicroBooNE.
https://github.com/marcodeltutto/evtwgtanalysis
Last synced: 8 months ago
JSON representation
Code to study GENIE model uncertainties at MicroBooNE.
- Host: GitHub
- URL: https://github.com/marcodeltutto/evtwgtanalysis
- Owner: marcodeltutto
- Created: 2016-03-28T22:13:05.000Z (about 10 years ago)
- Default Branch: master
- Last Pushed: 2016-07-01T09:42:36.000Z (almost 10 years ago)
- Last Synced: 2025-03-23T02:28:20.030Z (about 1 year ago)
- Language: C++
- Homepage:
- Size: 376 KB
- Stars: 0
- Watchers: 1
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
Awesome Lists containing this project
README
# EvtWgtAnalysis
Code to study GENIE model uncertainties at MicroBooNE.
##Compile
```
cd AnaTree
root -l
> .L AnaBNB.C+
> .q
cd ..
root -l
> gSystem->Load("AnaTree/AnaBNB_C.so");
> .L EvtWgtAnalysis.cxx+
```
##Run
It works with `root v5_34_32`.
```
source SetupEvtWgtAnalysis.sh
root -l
> gSystem->Load("AnaTree/AnaBNB_C.so");
> gSystem->Load("EvtWgtAnalysis_cxx.so");
> EvtWgtAnalysis f("generic_*_to_anatree.root");
> f.MakeBackgroundPlots();
```
To make the Number of Events or efficiency plots:
```
> f.MakePlots(false, option);
```
where `option` is an `integer`:
- `option=0`: efficiency - Pmu
- `option=1`: efficiency - CosThetaMu
- `option=2`: events - Pmu
- `option=3`: events - CosThetaMu
Change `false` to `true` to have area normalized plots.
To make the plots that show the events after the selection and the background decomposition:
```
> f.MakeBackgroundPlots(option)
```
where `option` is an `integer`:
- `option=0`: as a function of Pmu
- `option=1`: as a function of CosThetaMu
To make the plots that show the cross-section percental difference:
```
> f.MakeXsecDiffPlots()
```
##Generate new AnaBNB_C.so
Change file path in makeClass.C first, then:
```
cd AnaTree
root -l makeClass.C
root -l
> .L AnaBNB.C+
```
## Production
Produced files are in `/uboone/data/users/mdeltutt/EventWeight/`.
Merged root anatree files with event weights for _all_ GENIE parameters calculated with +-1 sigma:
* standard_reco_hist_merged_all_5.root
* standard_reco_hist_merged_all_6.root
* standard_reco_hist_merged_all_7.root (with nutools fix) 1.209150e+20 POT