message: PWGHF/jetsHF/AliAnalysisTaskHFJetIPQA.cxx: change from hard cut on R=0.4 to using radius given by jet container
diff: @@ -559,20 +559,35 @@ Bool_t AliAnalysisTaskHFJetIPQA::Run(){
  if (jetprob >0.9)FillHist(Form("fh1dJetRecPt_0_9JP_%sAccepted","all"),jetpt,1); //this->fXsectionWeightingFactor );
  if (jetprob >0.95)FillHist(Form("fh1dJetRecPt_0_95JP_%sAccepted","all"),jetpt,1); //this->fXsectionWeightingFactor );
  }
  + AliVParticle *vp=0x0;
  - for(Int_t itrack = 0; itrack < InputEvent()->GetNumberOfTracks(); ++itrack)
  - {
  + for(UInt_t i = 0; i < jetrec->GetNumberOfTracks(); i++) {
  TrackWeight=1;
  Double_t xyzatcda[3];
  - AliAODTrack * trackV = (AliAODTrack *) InputEvent()->GetTrack(itrack);
  +
  + vp = static_cast<AliVParticle*>(jetrec->TrackAt(i, jetconrec->GetParticleContainer()->GetArray()));
  + if (!vp){
  + Printf("AliVParticle associated to constituent not found");
  + continue;
  + }
  +
  + AliVTrack *vtrack = dynamic_cast<AliVTrack*>(vp);
  + if (!vtrack) {
  + printf("ERROR: Could not receive track%d\n", i);
  + continue;
  + }
  +
  + AliAODTrack *trackV = dynamic_cast<AliAODTrack*>(vtrack);
  +
  if (!trackV || !jetrec) continue;
  - if (jetrec->DeltaR(trackV) > 0.4) continue;
  +
  if (!IsTrackAccepted((AliAODTrack*)trackV,3)) continue;
  if(GetImpactParameterWrtToJet((AliAODTrack*)trackV,(AliAODEvent*)InputEvent(),jetrec,dca,cov,xyzatcda,sign)){
  if(fEventVertex) {
  delete fEventVertex;
  fEventVertex =nullptr;
  }
  +
  Int_t corridx=-1;double ppt;
  fIsPythia ? TrackWeight = GetMonteCarloCorrectionFactor(trackV,corridx,ppt) : TrackWeight =1;
  dca[0]=fabs(dca[0]);
  @@ -769,7 +784,7 @@ Bool_t AliAnalysisTaskHFJetIPQA::Run(){
  SJetIpPati d(cursImParXYZSig, TrackWeight,kFALSE,kFALSE,corridx);sImpParXYZSig.push_back(d);
  }
  - }
  + }//end trackloop
  std::sort(sImpParXY.begin(),sImpParXY.end(), AliAnalysisTaskHFJetIPQA::mysort);
  std::sort(sImpParXYSig.begin(),sImpParXYSig.end(), AliAnalysisTaskHFJetIPQA::mysort);
  std::sort(sImpParXYZ.begin(),sImpParXYZ.end(), AliAnalysisTaskHFJetIPQA::mysort);
  