message: adding pt1 to the pp and fixing the EP
diff: @@ -207,6 +207,7 @@ void AliAnalysisTaskNewJetSubstructure::UserCreateOutputObjects() {
  {
  fShapesVarNames[18] = "subjet1";
  fShapesVarNames[19] = "subjet2";
  +
  }
  if (fJetShapeType == kPythiaDef)
  {
  @@ -214,6 +215,8 @@ void AliAnalysisTaskNewJetSubstructure::UserCreateOutputObjects() {
  fShapesVarNames[13] = "subjet2";
  }
  }
  +
  + if (fJetShapeType == kData && fCentSelectOn == false)fShapesVarNames[14] = "ptsub1";
  if (fDoFlow) fShapesVarNames[20] = "EP";
  for (Int_t ivar = 0; ivar < nVar; ivar++) {
  @@ -566,7 +569,14 @@ Bool_t AliAnalysisTaskNewJetSubstructure::FillHistograms() {
  {
  Double_t fRPAngle = ((AliVAODHeader*)InputEvent()->GetHeader())->GetEventplane();
  Double_t phiBinT = RelativePhi(jet1->Phi(),fRPAngle);
  - fShapesVar[20] = phiBinT;
  + AliEventplane *aliEP = InputEvent()->GetEventplane();
  + Double_t fEPV0 = 0;
  + Double_t phiBinT2 = 0;
  + if (aliEP) {
  + fEPV0 = aliEP->GetEventplane("V0" ,InputEvent());
  + phiBinT2 = RelativePhi(jet1->Phi(),fEPV0);
  + }
  + fShapesVar[20] = phiBinT2;
  }
  fTreeSubstructure->Fill();
  @@ -963,6 +973,7 @@ void AliAnalysisTaskNewJetSubstructure::IterativeParentsPP(
  double zg = 0;
  double xktg = 0;
  double cumtf = 0;
  + double ptone = 0;
  while (jj.has_parents(j1, j2)) {
  nall = nall + 1;
  if (j1.perp() < j2.perp())
  @@ -989,6 +1000,7 @@ void AliAnalysisTaskNewJetSubstructure::IterativeParentsPP(
  j2first =j2;
  *sub2 = j2first;
  flagSubjet = 1;
  + ptone = j1.perp();
  }
  if (lnpt_rel > 0) {
  cumtf = cumtf + form;
  @@ -1012,6 +1024,7 @@ void AliAnalysisTaskNewJetSubstructure::IterativeParentsPP(
  fShapesVar[2] = nsd;
  fShapesVar[3] = zg;
  fShapesVar[4] = Rg;
  + fShapesVar[14] = ptone;
  } catch (fastjet::Error) {
  AliError(" [w] FJ Exception caught.");
  