Skip to content

Commit 28bb7de

Browse files
Mingze129Mingze129
authored
Merge pull request #23 from alibuild/alibot-cleanup-12436
Please consider the following formatting changes to #12436
2 parents 9878b6a + a04eb65 commit 28bb7de

File tree

1 file changed

+61
-64
lines changed

1 file changed

+61
-64
lines changed

‎PWGHF/HFJ/Tasks/taskDstarPolarisationInJet.cxx‎

Lines changed: 61 additions & 64 deletions
Original file line numberDiff line numberDiff line change
@@ -83,7 +83,7 @@ struct HfTaskDstarPolarisationInJet {
8383
float invMassMax{1000.f};
8484
float bkgRotationAngleStep{0.f};
8585
uint8_t nMassHypos{0u};
86-
struct HfHistoInput{
86+
struct HfHistoInput {
8787
float invMassCharmHad;
8888
float ptCharmHad;
8989
float rapCharmHad;
@@ -300,8 +300,8 @@ struct HfTaskDstarPolarisationInJet {
300300
{
301301

302302
if (cosThetaStarType == charm_polarisation::CosThetaStarType::Helicity) { // Helicity
303-
if (!doMc) { // data
304-
if (withMl) { // with ML
303+
if (!doMc) { // data
304+
if (withMl) { // with ML
305305
if (activateTrackingSys) {
306306
if (nBkgRotations > 0) {
307307
registry.fill(HIST("hHelicity"), recoHistoInput.invMassCharmHad, recoHistoInput.ptCharmHad, recoHistoInput.rapCharmHad, recoHistoInput.invMassD0, recoHistoInput.cosThetaStar, recoHistoInput.absEtaMin, recoHistoInput.numItsClsMin, recoHistoInput.numTpcClsMin, recoHistoInput.isRotatedCandidate, recoHistoInput.outputMl[0], /*recoHistoInput.outputMl[1],*/ recoHistoInput.outputMl[2], recoHistoInput.zParallel, recoHistoInput.jetPt);
@@ -330,8 +330,8 @@ struct HfTaskDstarPolarisationInJet {
330330
}
331331
}
332332
}
333-
} else { // MC --> no distinction among channels, since rotational bkg not supported
334-
if (withMl) { // with ML
333+
} else { // MC --> no distinction among channels, since rotational bkg not supported
334+
if (withMl) { // with ML
335335
if (recoHistoInput.origin == RecoDecay::OriginType::Prompt) { // prompt
336336
if (activateTrackingSys) {
337337
if (!isPartRecoDstar) {
@@ -361,7 +361,7 @@ struct HfTaskDstarPolarisationInJet {
361361
}
362362
}
363363
}
364-
} else { // without ML
364+
} else { // without ML
365365
if (recoHistoInput.origin == RecoDecay::OriginType::Prompt) { // prompt
366366
if (activateTrackingSys) {
367367
if (!isPartRecoDstar) {
@@ -395,8 +395,8 @@ struct HfTaskDstarPolarisationInJet {
395395
}
396396
}
397397
} else if (cosThetaStarType == charm_polarisation::CosThetaStarType::Production) { // Production
398-
if (!doMc) { // data
399-
if (withMl) { // with ML
398+
if (!doMc) { // data
399+
if (withMl) { // with ML
400400
if (activateTrackingSys) {
401401
if (nBkgRotations > 0) {
402402
registry.fill(HIST("hProduction"), recoHistoInput.invMassCharmHad, recoHistoInput.ptCharmHad, recoHistoInput.rapCharmHad, recoHistoInput.invMassD0, recoHistoInput.cosThetaStar, recoHistoInput.absEtaMin, recoHistoInput.numItsClsMin, recoHistoInput.numTpcClsMin, recoHistoInput.isRotatedCandidate, recoHistoInput.outputMl[0], /*recoHistoInput.outputMl[1],*/ recoHistoInput.outputMl[2], recoHistoInput.zParallel, recoHistoInput.jetPt);
@@ -425,8 +425,8 @@ struct HfTaskDstarPolarisationInJet {
425425
}
426426
}
427427
}
428-
} else { // MC --> no distinction among channels, since rotational bkg not supported
429-
if (withMl) { // with ML
428+
} else { // MC --> no distinction among channels, since rotational bkg not supported
429+
if (withMl) { // with ML
430430
if (recoHistoInput.origin == RecoDecay::OriginType::Prompt) { // prompt
431431
if (activateTrackingSys) {
432432
if (!isPartRecoDstar) {
@@ -456,7 +456,7 @@ struct HfTaskDstarPolarisationInJet {
456456
}
457457
}
458458
}
459-
} else { // without ML
459+
} else { // without ML
460460
if (recoHistoInput.origin == RecoDecay::OriginType::Prompt) { // prompt
461461
if (activateTrackingSys) {
462462
if (!isPartRecoDstar) {
@@ -490,8 +490,8 @@ struct HfTaskDstarPolarisationInJet {
490490
}
491491
}
492492
} else if (cosThetaStarType == charm_polarisation::CosThetaStarType::JetAxis) { // JetAxis
493-
if (!doMc) { // data
494-
if (withMl) { // with ML
493+
if (!doMc) { // data
494+
if (withMl) { // with ML
495495
if (activateTrackingSys) {
496496
if (nBkgRotations > 0) {
497497
registry.fill(HIST("hJetAxis"), recoHistoInput.invMassCharmHad, recoHistoInput.ptCharmHad, recoHistoInput.rapCharmHad, recoHistoInput.invMassD0, recoHistoInput.cosThetaStar, recoHistoInput.absEtaMin, recoHistoInput.numItsClsMin, recoHistoInput.numTpcClsMin, recoHistoInput.isRotatedCandidate, recoHistoInput.outputMl[0], /*recoHistoInput.outputMl[1],*/ recoHistoInput.outputMl[2], recoHistoInput.zParallel, recoHistoInput.jetPt);
@@ -520,8 +520,8 @@ struct HfTaskDstarPolarisationInJet {
520520
}
521521
}
522522
}
523-
} else { // MC --> no distinction among channels, since rotational bkg not supported
524-
if (withMl) { // with ML
523+
} else { // MC --> no distinction among channels, since rotational bkg not supported
524+
if (withMl) { // with ML
525525
if (recoHistoInput.origin == RecoDecay::OriginType::Prompt) { // prompt
526526
if (activateTrackingSys) {
527527
if (!isPartRecoDstar) {
@@ -551,7 +551,7 @@ struct HfTaskDstarPolarisationInJet {
551551
}
552552
}
553553
}
554-
} else { // without ML
554+
} else { // without ML
555555
if (recoHistoInput.origin == RecoDecay::OriginType::Prompt) { // prompt
556556
if (activateTrackingSys) {
557557
if (!isPartRecoDstar) {
@@ -810,66 +810,63 @@ struct HfTaskDstarPolarisationInJet {
810810
// helicity
811811
cosThetaStarHelicity = helicityVec.Dot(threeVecDauCM) / std::sqrt(threeVecDauCM.Mag2()) / std::sqrt(helicityVec.Mag2());
812812
HfHistoInput helicityInput{
813-
.invMassCharmHad = invMassCharmHadForSparse,
814-
.ptCharmHad = ptCharmHad,
815-
.rapCharmHad = rapidity,
816-
.invMassD0 = invMassD0,
817-
.cosThetaStar = cosThetaStarHelicity,
818-
.outputMl = outputMl,
819-
.isRotatedCandidate = isRotatedCandidate,
820-
.origin = origin,
821-
.ptBhadMother = ptBhadMother,
822-
.absEtaMin = absEtaTrackMin,
823-
.numItsClsMin = numItsClsMin,
824-
.numTpcClsMin = numTpcClsMin,
825-
.nMuons = nMuons,
826-
.zParallel = zParallel,
827-
.jetPt = jetPt
828-
};
813+
.invMassCharmHad = invMassCharmHadForSparse,
814+
.ptCharmHad = ptCharmHad,
815+
.rapCharmHad = rapidity,
816+
.invMassD0 = invMassD0,
817+
.cosThetaStar = cosThetaStarHelicity,
818+
.outputMl = outputMl,
819+
.isRotatedCandidate = isRotatedCandidate,
820+
.origin = origin,
821+
.ptBhadMother = ptBhadMother,
822+
.absEtaMin = absEtaTrackMin,
823+
.numItsClsMin = numItsClsMin,
824+
.numTpcClsMin = numTpcClsMin,
825+
.nMuons = nMuons,
826+
.zParallel = zParallel,
827+
.jetPt = jetPt};
829828
fillRecoHistos(charm_polarisation::CosThetaStarType::Helicity, withMl, doMc, partRecoDstar, helicityInput);
830829
}
831830
if (activateTHnSparseCosThStarProduction) {
832831
// production
833832
cosThetaStarProduction = normalVec.Dot(threeVecDauCM) / std::sqrt(threeVecDauCM.Mag2()) / std::sqrt(normalVec.Mag2());
834833
HfHistoInput productionInput{
835-
.invMassCharmHad = invMassCharmHadForSparse,
836-
.ptCharmHad = ptCharmHad,
837-
.rapCharmHad = rapidity,
838-
.invMassD0 = invMassD0,
839-
.cosThetaStar = cosThetaStarProduction,
840-
.outputMl = outputMl,
841-
.isRotatedCandidate = isRotatedCandidate,
842-
.origin = origin,
843-
.ptBhadMother = ptBhadMother,
844-
.absEtaMin = absEtaTrackMin,
845-
.numItsClsMin = numItsClsMin,
846-
.numTpcClsMin = numTpcClsMin,
847-
.nMuons = nMuons,
848-
.zParallel = zParallel,
849-
.jetPt = jetPt
850-
};
834+
.invMassCharmHad = invMassCharmHadForSparse,
835+
.ptCharmHad = ptCharmHad,
836+
.rapCharmHad = rapidity,
837+
.invMassD0 = invMassD0,
838+
.cosThetaStar = cosThetaStarProduction,
839+
.outputMl = outputMl,
840+
.isRotatedCandidate = isRotatedCandidate,
841+
.origin = origin,
842+
.ptBhadMother = ptBhadMother,
843+
.absEtaMin = absEtaTrackMin,
844+
.numItsClsMin = numItsClsMin,
845+
.numTpcClsMin = numTpcClsMin,
846+
.nMuons = nMuons,
847+
.zParallel = zParallel,
848+
.jetPt = jetPt};
851849
fillRecoHistos(charm_polarisation::CosThetaStarType::Production, withMl, doMc, partRecoDstar, productionInput);
852850
}
853851
if (activateTHnSparseCosThStarJetAxis) {
854852
// jet axis
855853
cosThetaStarJet = jetaxisVec.Dot(threeVecDauCM) / std::sqrt(threeVecDauCM.Mag2()) / std::sqrt(jetaxisVec.Mag2());
856854
HfHistoInput jetAxisInput{
857-
.invMassCharmHad = invMassCharmHadForSparse,
858-
.ptCharmHad = ptCharmHad,
859-
.rapCharmHad = rapidity,
860-
.invMassD0 = invMassD0,
861-
.cosThetaStar = cosThetaStarJet,
862-
.outputMl = outputMl,
863-
.isRotatedCandidate = isRotatedCandidate,
864-
.origin = origin,
865-
.ptBhadMother = ptBhadMother,
866-
.absEtaMin = absEtaTrackMin,
867-
.numItsClsMin = numItsClsMin,
868-
.numTpcClsMin = numTpcClsMin,
869-
.nMuons = nMuons,
870-
.zParallel = zParallel,
871-
.jetPt = jetPt
872-
};
855+
.invMassCharmHad = invMassCharmHadForSparse,
856+
.ptCharmHad = ptCharmHad,
857+
.rapCharmHad = rapidity,
858+
.invMassD0 = invMassD0,
859+
.cosThetaStar = cosThetaStarJet,
860+
.outputMl = outputMl,
861+
.isRotatedCandidate = isRotatedCandidate,
862+
.origin = origin,
863+
.ptBhadMother = ptBhadMother,
864+
.absEtaMin = absEtaTrackMin,
865+
.numItsClsMin = numItsClsMin,
866+
.numTpcClsMin = numTpcClsMin,
867+
.nMuons = nMuons,
868+
.zParallel = zParallel,
869+
.jetPt = jetPt};
873870
fillRecoHistos(charm_polarisation::CosThetaStarType::JetAxis, withMl, doMc, partRecoDstar, jetAxisInput);
874871
}
875872
} /// end loop over mass hypotheses

0 commit comments

Comments
 (0)