|
[Rivet-svn] r3965 - branches/2012-06-aidarivet/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgMon Oct 15 10:28:22 BST 2012
Author: hoeth Date: Mon Oct 15 10:28:22 2012 New Revision: 3965 Log: NaN bugfix in CMS_2012_I1107658 Modified: branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1107658.cc Modified: branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1107658.cc ============================================================================== --- branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1107658.cc Mon Oct 15 09:59:45 2012 (r3964) +++ branches/2012-06-aidarivet/src/Analyses/CMS_2012_I1107658.cc Mon Oct 15 10:28:22 2012 (r3965) @@ -101,9 +101,9 @@ _h_pTsum_towards_pTmumu-> fill(Zpt, 1./area * ptSumTowards, weight); _h_pTsum_transverse_pTmumu-> fill(Zpt, 1./area * ptSumTransverse, weight); _h_pTsum_away_pTmumu-> fill(Zpt, 1./area * ptSumAway, weight); - _h_avgpT_towards_pTmumu-> fill(Zpt, 1.*nTowards/ptSumTowards, weight); - _h_avgpT_transverse_pTmumu-> fill(Zpt, 1.*nTransverse/ptSumTransverse, weight); - _h_avgpT_away_pTmumu-> fill(Zpt, 1.*nAway/ptSumAway, weight); + if (nTowards > 0) _h_avgpT_towards_pTmumu-> fill(Zpt, ptSumTowards/nTowards, weight); + if (nTransverse > 0) _h_avgpT_transverse_pTmumu-> fill(Zpt, ptSumTransverse/nTransverse, weight); + if (nAway > 0) _h_avgpT_away_pTmumu-> fill(Zpt, ptSumAway/nAway, weight); _h_Nchg_towards_zmass_81_101-> fill(nTowards, weight); _h_Nchg_transverse_zmass_81_101->fill(nTransverse, weight); _h_Nchg_away_zmass_81_101-> fill(nAway, weight); @@ -112,7 +112,7 @@ if (Zpt < 5.) { _h_Nchg_towards_plus_transverse_Mmumu->fill(Zmass, (nTowards + nTransverse)/(2.*area), weight); _h_pTsum_towards_plus_transverse_Mmumu->fill(Zmass, (ptSumTowards + ptSumTransverse)/(2.*area), weight); - _h_avgpT_towards_plus_transverse_Mmumu->fill(Zmass, 1.*(nTowards + nTransverse)/(ptSumTowards + ptSumTransverse), weight); + if ((nTowards + nTransverse) > 0) _h_avgpT_towards_plus_transverse_Mmumu->fill(Zmass, (ptSumTowards + ptSumTransverse)/(nTowards + nTransverse), weight); _h_Nchg_transverse_zpt_5->fill(nTransverse, weight); }
More information about the Rivet-svn mailing list |