|
[Rivet-svn] rivet: ATLAS done up to analyses/pluginATLAS/ATLAS_2015_I1394679.ccRivet Mercurial rivet at projects.hepforge.orgWed Sep 27 13:30:02 BST 2017
details: https://rivet.hepforge.org/hg/rivet/rev/2901986bb2e0 branches: multiweight changeset: 6051:2901986bb2e0 user: Chris Pollard <cspollard at gmail.com> date: Wed Sep 27 14:28:56 2017 +0200 description: ATLAS done up to analyses/pluginATLAS/ATLAS_2015_I1394679.cc diffs (truncated from 1326 to 50 lines): --- a/analyses/pluginATLAS/ATLAS_2014_I1298023.cc Wed Sep 27 14:07:01 2017 +0200 +++ b/analyses/pluginATLAS/ATLAS_2014_I1298023.cc Wed Sep 27 14:28:56 2017 +0200 @@ -55,8 +55,6 @@ /// Perform the per-event analysis void analyze(const Event& event) { - const double weight = 1.0; - const vector<DressedLepton>& leptons = apply<DressedLeptons>(event, "leptons").dressedLeptons(); if ( leptons.size() < 2 ) vetoEvent; @@ -98,10 +96,10 @@ if ( dijet.mass() <= 500*GeV ) vetoEvent; // inclusive region - _hist->fill(0.5, weight); + _hist->fill(0.5); // VBS region - if ( deltaRap(jets[0], jets[1]) > 2.4 ) _hist->fill(1.5, weight); + if ( deltaRap(jets[0], jets[1]) > 2.4 ) _hist->fill(1.5); } /// Normalise histograms etc., after the run --- a/analyses/pluginATLAS/ATLAS_2014_I1298811.cc Wed Sep 27 14:07:01 2017 +0200 +++ b/analyses/pluginATLAS/ATLAS_2014_I1298811.cc Wed Sep 27 14:28:56 2017 +0200 @@ -57,9 +57,6 @@ // Require at least one jet in the event if (jets.empty()) vetoEvent; - // Get the event weight since we will be filling some histos - const double weight = 1.0; - // Identify the leading jet and its phi and pT const FourMomentum plead = jets[0].momentum(); const double philead = plead.phi(); @@ -131,27 +128,27 @@ } // Plot profiles and distributions which have no max/min region definition - _p_etsum25_vs_ptlead_trans[itopo]->fill(ptlead/GeV, etsum25[0]/5.0/dphi[0]/GeV, weight); - _p_etsum48_vs_ptlead_trans[itopo]->fill(ptlead/GeV, etsum48[0]/9.6/dphi[0]/GeV, weight); + _p_etsum25_vs_ptlead_trans[itopo]->fill(ptlead/GeV, etsum25[0]/5.0/dphi[0]/GeV); + _p_etsum48_vs_ptlead_trans[itopo]->fill(ptlead/GeV, etsum48[0]/9.6/dphi[0]/GeV); if (etsum25[0] > 0) { - _p_chratio_vs_ptlead_trans[itopo]->fill(ptlead/GeV, ptsum[0]/etsum25[0], weight); + _p_chratio_vs_ptlead_trans[itopo]->fill(ptlead/GeV, ptsum[0]/etsum25[0]); } const double ptmean = safediv(ptsum[0], nch[0], -1); ///< Return -1 if div by zero if (ptmean >= 0) {
More information about the Rivet-svn mailing list |