message: Add matching distance to the response THnSparse
diff: @@ -82,10 +82,10 @@ void AliAnalysisTaskEmcalJetEnergyScale::UserCreateOutputObjects(){
  AliAnalysisTaskEmcal::UserCreateOutputObjects();
  TLinearBinning jetPtBinningDet(300, 0., 300.), jetPtBinningPart(500, 0., 500), nefbinning(100, 0., 1.), ptdiffbinning(200, -1., 1.), jetEtaBinning(100, -0.9, 0.9), jetPhiBinning(100, 0., TMath::TwoPi()),
  - subsampleBinning(2, -0.5, 1.5);
  + subsampleBinning(2, -0.5, 1.5), deltaRbinning(20, 0., 1.);
  const TBinning *diffbinning[3] = {&jetPtBinningPart, &nefbinning, &ptdiffbinning},
  - *corrbinning[4] = {&jetPtBinningPart, &jetPtBinningDet, &nefbinning, &subsampleBinning},
  + *corrbinning[5] = {&jetPtBinningPart, &jetPtBinningDet, &nefbinning, &deltaRbinning,&subsampleBinning},
  *effbinning[3] = {&jetPtBinningPart, &jetEtaBinning, &jetPhiBinning};
  TCustomBinning jetPtBinningCoarseDet, jetPtBinningCoarsePart;
  @@ -107,7 +107,7 @@ void AliAnalysisTaskEmcalJetEnergyScale::UserCreateOutputObjects(){
  fHistos->CreateTH1("hSpectrumTrueTruncated", "jet pt spectrum particle level, truncated", jetPtBinningCoarsePart);
  fHistos->CreateTH2("hJetResponseCoarse", "Response matrix, coarse binning", jetPtBinningCoarseDet, jetPtBinningCoarsePart);
  fHistos->CreateTHnSparse("hPtDiff", "pt diff det/part", 3, diffbinning, "s");
  - fHistos->CreateTHnSparse("hPtCorr", "Correlation det pt / part pt", 4, corrbinning, "s");
  + fHistos->CreateTHnSparse("hPtCorr", "Correlation det pt / part pt", 5, corrbinning, "s");
  fHistos->CreateTHnSparse("hPartJetsAccepted", "Accepted particle level jets", 3, effbinning, "s");
  fHistos->CreateTHnSparse("hPartJetsReconstructed", "Accepted and reconstructed particle level jets", 3, effbinning, "s");
  for(auto h : *(fHistos->GetListOfHistograms())) fOutput->Add(h);
  @@ -148,8 +148,11 @@ Bool_t AliAnalysisTaskEmcalJetEnergyScale::Run(){
  AliDebugStream(2) << "No tagged jet" << std::endl;
  continue;
  }
  + TVector3 basevec, tagvec;
  + basevec.SetPtEtaPhi(detjet->Pt(), detjet->Eta(), detjet->Phi());
  + tagvec.SetPtEtaPhi(partjet->Pt(), partjet->Eta(), partjet->Phi());
  taggedjets.emplace_back(partjet);
  - double pointCorr[4] = {partjet->Pt(), detjet->Pt(), detjet->NEF(), fSampleSplitter->Uniform() < fFractionResponseClosure ? 0. : 1.},
  + double pointCorr[5] = {partjet->Pt(), detjet->Pt(), detjet->NEF(), basevec.DeltaR(tagvec), fSampleSplitter->Uniform() < fFractionResponseClosure ? 0. : 1.},
  pointDiff[3] = {partjet->Pt(), detjet->NEF(), (detjet->Pt()-partjet->Pt())/partjet->Pt()};
  fHistos->FillTHnSparse("hPtDiff", pointDiff);
  fHistos->FillTHnSparse("hPtCorr", pointCorr);
  