- fShapesVar[2] = GetJetpTD(jet1);
- fShapesVar[3] = GetJetMass(jet1);
- fShapesVar[4] = 1.*GetJetNumberOfConstituents(jet1);
- fShapesVar[5] = GetJetAngularity(jet1);
- fShapesVar[6] = GetJetCircularity(jet1);
- fShapesVar[7] = GetJetLeSub(jet1);
- fShapesVar[8] = GetSigma2(jet1);
+ fShapesVar[2] = GetJetpTD(jet1,0);
+ fShapesVar[3] = GetJetMass(jet1,0);
+ fShapesVar[4] = 1.*GetJetNumberOfConstituents(jet1,0);
+ fShapesVar[5] = GetJetAngularity(jet1,0);
+ fShapesVar[6] = GetJetCircularity(jet1,0);
+ fShapesVar[7] = GetJetLeSub(jet1,0);
+ fShapesVar[8] = GetSigma2(jet1,0);
+
+ Float_t ptMatch=0., ptDMatch=0., massMatch=0., constMatch=0.,angulMatch=0.,circMatch=0., lesubMatch=0., sigma2Match=0.;
+ Int_t kMatched = 0;
+ if (fJetShapeType == kTrueDet || fJetShapeType == kDetEmb) {
+ kMatched = 1;
+ ptMatch=jet2->Pt();
+ ptDMatch=GetJetpTD(jet2, kMatched);
+ massMatch=GetJetMass(jet2,kMatched);
+ constMatch=1.*GetJetNumberOfConstituents(jet2,kMatched);
+ angulMatch=GetJetAngularity(jet2, kMatched);
+ circMatch=GetJetCircularity(jet2, kMatched);
+ lesubMatch=GetJetLeSub(jet2, kMatched);
+ sigma2Match = GetSigma2(jet2, kMatched);
+
+ }
+
+ if (fJetShapeType == kTrue || fJetShapeType == kData) {
+ kMatched = 0;
+ ptMatch=0.;
+ ptDMatch=0.;
+ massMatch=0.;
+ constMatch=0.;
+ angulMatch=0.;
+ circMatch=0.;
+ lesubMatch=0.;
+ sigma2Match =0.;
+
+ }
+
+ fShapesVar[9] = ptMatch;
+ fShapesVar[10] = ptDMatch;
+ fShapesVar[11] = massMatch;
+ fShapesVar[12] = constMatch;
+ fShapesVar[13] = angulMatch;
+ fShapesVar[14] = circMatch;
+ fShapesVar[15] = lesubMatch;
+ fShapesVar[16] = sigma2Match;
+ fShapesVar[17] = kWeight;