using pT() function of pythia8/event

cms-sw · Sep 5, 2016 · 2e88f79 · 2e88f79
1 parent c694e9a
commit 2e88f79
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/GeneratorInterface/Pythia8Interface/src/PTFilterHook.cc b/GeneratorInterface/Pythia8Interface/src/PTFilterHook.cc
@@ -25,8 +25,7 @@ bool PTFilterHook::checkVetoPT( int iPos, const Pythia8::Event& event) {
   //look for quark
   for (int i = 0; i < event.size(); ++i) {
     if ( (abs(event[i].id()) == quark_ ) && ( abs(event[i].y()) <= quarkY_ ) ) {
-       double pT = sqrt(pow(event[i].px(),2) + pow(event[i].py(),2));
-       if (pT >= quarkPt_) { 
+       if (event[i].pT() >= quarkPt_) { 
           foundQuark = true;
           break;
        }