<html> <head> <meta http-equiv="content-type" content="text/html; charset=utf-8"> </head> <body bgcolor="#FFFFFF" text="#000000"> Hi Michael, during testing we found that your analysis plugin crashes for certain events when trying to calculate the median pTDensity. Could I please ask you to have a look at these code snippets. <tt> // Get the jet pT densities</tt><tt> </tt><tt> vector< vector<double> > ptDensities(_eta_bins_areaoffset.size()-1);</tt><tt> </tt><tt> FastJets fastjets = apply<FastJets>(event, "KtJetsD05");</tt><tt> </tt><tt> const shared_ptr<fastjet::ClusterSequenceArea> clust_seq_area = fastjets.clusterSeqArea();</tt><tt> </tt><tt> for (const Jet& jet : fastjets.jets()) {</tt><tt> </tt><tt> const double area = clust_seq_area->area(jet); //< Implicit call to pseudojet()</tt><tt> </tt><tt> /// @todo Should be 1e-4?</tt><tt> </tt><tt> if (area > 10e-4 && jet.abseta() < _eta_bins_areaoffset.back()) {</tt><tt> </tt><tt> ptDensities.at(getEtaBin(jet.abseta(), true)) += jet.pT()/area;</tt><tt> </tt><tt> }</tt><tt> </tt><tt> }</tt><tt> </tt><tt> </tt><tt> // Now compute the median energy densities</tt><tt> </tt><tt> vector<double> ptDensity;</tt><tt> </tt><tt> for (size_t b = 0; b < _eta_bins_areaoffset.size()-1; ++b) {</tt><tt> </tt><tt> if (ptDensities[b].size() >0 ) ptDensity += median(ptDensities[b]);</tt><tt> </tt><tt> </tt><tt> }</tt><tt> </tt> The red bit was added by us to prevent calculating the median of an empty vector --- do you remember what the logic in your analysis code was for cases like that? If you find the time, could you please comment on the blue bit as well? You can find the version of the code in the current release here: <blockquote> <a class="moz-txt-link-freetext" href="https://rivet.hepforge.org/hg/rivet/file/0175188f89df/src/Analyses/ATLAS_2010_S8914702.cc">https://rivet.hepforge.org/hg/rivet/file/0175188f89df/src/Analyses/ATLAS_2010_S8914702.cc</a> </blockquote> Cheers, Holger </body> </html>