The** HadronicRecoilPFO2016 ** package is developed originally for the W mass analysis at 13 TeV ; however it is also intended for other measurements that may find interest in having a proper recoil definition (e.g. pTW measurement). It runs only at analysis level making use of derivations containing the relevant information. In the next release of STDM derivations, STDM4 should contain all that is needed. Please email to your derivation contact and this page's responsible if you want it included in your derivation (NB : the size increase due to the additional containers can be non-negligible).

At the moment, there is no calibration or systematic derived ; it is planned to do it the same way as it was done for the 7 TeV W mass analysis, i.e. using Z events.

This Twiki provides technical details about the use of the package.

svn co svn+ssh://svn.cern.ch/reps/atlasinst/Institutes/Saclay/Run2AnalysisW/HadRecoil4Winners/trunk HadRecoil4WinnersThen include the package in the dependencies of your code (AnalysisCode/cmt/Makefile).

Exemple code on how to use it ; in the header file :

#include "HadRecoil4Winners/hadrecoil.h" hadrecoil *m_hadrecoil;In the code at initialization :

m_hadrecoil = new hadrecoil("recoilTool");Then in the event loop you can access the HR 4-vectors like this (several definitions are available, see below) ; I put the type of the argument to be passed to each function (of course to be removed in the code). Be careful to take the final definition of your lepton.

TLorentzVector PFOHR = m_hadrecoil -> getHadRecoilPFO(xAOD::TEvent& m_event, xAOD::IParticleContainer* leptons); TLorentzVector PFOVoronoiHR = m_hadrecoil -> getHadRecoilPFOVoronoi(xAOD::TEvent& m_event, xAOD::IParticleContainer* leptons); TLorentzVector PFOVoronoiSpreadingHR = m_hadrecoil -> getHadRecoilPFOVoronoiSpreading(xAOD::TEvent& m_event, xAOD::IParticleContainer* leptons); TLorentzVector PFOSKVoronoiHR = m_hadrecoil -> getHadRecoilPFOSKVoronoi(xAOD::TEvent& m_event, xAOD::IParticleContainer* leptons);Eventually don't forget to delete the tool in finalization:

if(m_hadrecoil){ delete m_hadrecoil; m_hadrecoil = 0;}

- getHadRecoilPFO : uses the sum of all particle flow objects (at LCW scale) ; charged PFOs (cPFOs) have pile-up suppression applied (remove cPFOs whose associated track does not come from the PV).
- getHadRecoilPFOVoronoi : uses the sum of all particle flow objects (at LCW scale) to calculate the Voronoi areas. Then, Voronoi pile-up suppression is applied to the neutral PFOs (nPFOs). Eventually, charged-suppressed cPFOs and nPFOs are added.
- getHadRecoilPFOVoronoiSpreading : same as above but using Voronoi spreading instead of simple Voronoi subtraction for nPFOs.
- getHadRecoilPFOSKVoronoi : same as getHadRecoilPFOVoronoi but using softkiller + Voronoi subtraction for nPFOs.

NB : a refined treatment of the leptons has to be done for the case of multijet background evaluation. This is not done at the moment ; at 7 TeV it used to replace the cone hole with the lepton isolation value

- bias_wplusenu.pdf: mean value of bias (|u_par| - pTW) vs pTW, Wplusenu

- bias_zee.pdf: mean value of bias (|u_par| - pTZ) vs pTZ, Zee

- rms_uperpVsSumet_wplusenu.pdf: RMS of u_perp vs CST SumEt, Wplusenu

- rms_uperpVsSumet_zee.pdf: RMS of u_perp vs CST SumEt, Zee

- rmsUperpOverUAvg_wplusenu_ratio.pdf: sigma(u_perp) / |u| vs pTW, Wplusenu, ratio w.r.t. nominal

- rmsUperpOverUAvg_zee_ratio.pdf: sigma(u_perp) / |u| vs pTZ, Zee, ratio w.r.t. nominal

- rmsUperpOverUMU_wplusenu.pdf: sigma(u_perp) / |u| vs NPV, Wplusenu

- rmsUperpOverUMU_zee.pdf: sigma(u_perp) / |u| vs NPV, Zee

- rmsUperpOverUNPV_wplusenu.pdf: sigma(u_perp) / |u| vs mu, Wplusenu

- rmsUperpOverUNPV_zee.pdf: sigma(u_perp) / |u| vs mu, Zee

-- FabriceBalli - 2016-04-27

I | Attachment | History | Action | Size | Date | Who | Comment |
---|---|---|---|---|---|---|---|

bias_wplusenu.pdf | r1 | manage | 22.0 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

bias_zee.pdf | r1 | manage | 24.4 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rmsUperpOverUAvg_wplusenu_ratio.pdf | r1 | manage | 15.0 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rmsUperpOverUAvg_zee_ratio.pdf | r1 | manage | 14.9 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rmsUperpOverUMU_wplusenu.pdf | r1 | manage | 18.9 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rmsUperpOverUMU_zee.pdf | r1 | manage | 20.7 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rmsUperpOverUNPV_wplusenu.pdf | r1 | manage | 21.3 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rmsUperpOverUNPV_zee.pdf | r1 | manage | 22.6 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rms_uperpVsSumet_wplusenu.pdf | r1 | manage | 23.7 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu | |

rms_uperpVsSumet_zee.pdf | r1 | manage | 23.7 K | 2016-04-27 - 12:31 | FabriceBalli | performance plots for hadronic recoil package , Zee and Wplusenu |

Edit | Attach | ~~Watch~~ | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions

Topic revision: r1 - 2016-04-27 - FabriceBalli

Copyright &© 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.

or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback

or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback