Merge pull request #12 from JETSCAPE/pt_broadening_fix

applies Ritoban's pT-broadening fix from JETSCAPE

Author: latessa

src/jet/LBT.cc

@@ -191,6 +191,11 @@ void LBT::DoEnergyLoss(double deltaT, double time, double Q2,
       return;
     }
 
+    if (pIn[0].pstat()==101) {
+      // pOut.push_back(pIn[0]);
+       // JSINFO << BOLDYELLOW << " broadenend parton rejected by LBT ";
+        return;
+    }
     // pass particle infomation to LBT array (only pass one particle each time)
 
     jetClean();

src/jet/Matter.cc

@@ -1681,13 +1681,17 @@ void Matter::DoEnergyLoss(double deltaT, double time, double Q2,
               energy -= drag;
               pOut[iout].reset_momentum(px, py, pz, energy);
             }
-
+            pOut[iout].set_stat(101);
             VERBOSE(8) << BOLDYELLOW << " p after b & d, E = " << energy
                        << " pz = " << pz << " px = " << px << " py = " << py;
           }
-
+          else{
+            pOut.push_back(pIn[i]);
+          }
         } // end if(broadening_on)
+        else{
           pOut.push_back(pIn[i]);
+        }
       }
     } else { // virtuality too low lets broaden it
 
@@ -1837,13 +1841,18 @@ void Matter::DoEnergyLoss(double deltaT, double time, double Q2,
             energy -= drag;
             pOut[iout].reset_momentum(px, py, pz, energy);
           }
-
+          pOut[iout].set_stat(101);
           VERBOSE(8) << BOLDYELLOW << " p after b & d, E = " << energy
                      << " pz = " << pz << " px = " << px << " py = " << py;
         }
-
+        else{
+            pOut.push_back(pIn[i]);
+        }
         //pOut.push_back(pIn[i]);
       }
+      else{
+        pOut.push_back(pIn[i]);
+      }
     }
 
   } // particle loop