@@ -100,7 +100,7 @@ struct cascademcfinder {
100100 histos.add (" hPtXiPlusReconstructed" , " hPtXiPlusReconstructed" , kTH1F , {axisPt});
101101 histos.add (" hPtOmegaMinusReconstructed" , " hPtOmegaMinusReconstructed" , kTH1F , {axisPt});
102102 histos.add (" hPtOmegaPlusReconstructed" , " hPtOmegaPlusReconstructed" , kTH1F , {axisPt});
103-
103+
104104 histos.add (" hPtXiMinusGlobal" , " hPtXiMinusGlobal" , kTH1F , {axisPt});
105105 histos.add (" hPtXiPlusGlobal" , " hPtXiPlusGlobal" , kTH1F , {axisPt});
106106 histos.add (" hPtOmegaMinusGlobal" , " hPtOmegaMinusGlobal" , kTH1F , {axisPt});
@@ -207,69 +207,69 @@ struct cascademcfinder {
207207 for (auto & mcParticle : mcParticles) {
208208 if (mcParticle.pdgCode () == 3312 && findXiMinus) {
209209 reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS, positiveTPC, negativeTPC, bachelorTPC);
210- bool reconstructed = false ;
211- for (auto & mcParticle : mcParticles) {
212- if (mcParticle.pdgCode () == 3312 && findXiMinus) {
213- reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS);
214- if (fabs (mcParticle.y ()) < 0.5 ) {
215- histos.fill (HIST (" hPtXiMinusGenerated" ), mcParticle.pt ());
216- if (reconstructed)
217- histos.fill (HIST (" hPtXiMinusReconstructed" ), mcParticle.pt ());
218- if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
219- histos.fill (HIST (" hPtXiMinusGlobal" ), mcParticle.pt ());
220- if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
221- histos.fill (HIST (" hPtXiMinusGlobalWithPV" ), mcParticle.pt ());
222- }
223- }
224- if (mcParticle.pdgCode () == -3312 && findXiPlus) {
225- reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS, positiveTPC, negativeTPC, bachelorTPC);
226- if (fabs (mcParticle.y ()) < 0.5 ) {
227- histos.fill (HIST (" hPtXiPlusGenerated" ), mcParticle.pt ());
228- if (reconstructed)
229- histos.fill (HIST (" hPtXiPlusReconstructed" ), mcParticle.pt ());
230- if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
231- histos.fill (HIST (" hPtXiPlusGlobal" ), mcParticle.pt ());
232- if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
233- histos.fill (HIST (" hPtXiPlusGlobalWithPV" ), mcParticle.pt ());
234- }
235- }
236- if (mcParticle.pdgCode () == 3334 && findOmegaMinus) {
237- reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS, positiveTPC, negativeTPC, bachelorTPC);
238- if (fabs (mcParticle.y ()) < 0.5 ) {
239- histos.fill (HIST (" hPtOmegaMinusGenerated" ), mcParticle.pt ());
240- if (reconstructed)
241- histos.fill (HIST (" hPtOmegaMinusReconstructed" ), mcParticle.pt ());
242- if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
243- histos.fill (HIST (" hPtOmegaMinusGlobal" ), mcParticle.pt ());
244- if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
245- histos.fill (HIST (" hPtOmegaMinusGlobalWithPV" ), mcParticle.pt ());
246- }
247- }
248- if (mcParticle.pdgCode () == -3334 && findOmegaPlus) {
249- reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS, positiveTPC, negativeTPC, bachelorTPC);
250- if (fabs (mcParticle.y ()) < 0.5 ) {
251- histos.fill (HIST (" hPtOmegaPlusGenerated" ), mcParticle.pt ());
252- if (reconstructed)
253- histos.fill (HIST (" hPtOmegaPlusReconstructed" ), mcParticle.pt ());
254- if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
255- histos.fill (HIST (" hPtOmegaPlusGlobal" ), mcParticle.pt ());
256- if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
257- histos.fill (HIST (" hPtOmegaPlusGlobalWithPV" ), mcParticle.pt ());
210+ bool reconstructed = false ;
211+ for (auto & mcParticle : mcParticles) {
212+ if (mcParticle.pdgCode () == 3312 && findXiMinus) {
213+ reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS);
214+ if (fabs (mcParticle.y ()) < 0.5 ) {
215+ histos.fill (HIST (" hPtXiMinusGenerated" ), mcParticle.pt ());
216+ if (reconstructed)
217+ histos.fill (HIST (" hPtXiMinusReconstructed" ), mcParticle.pt ());
218+ if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
219+ histos.fill (HIST (" hPtXiMinusGlobal" ), mcParticle.pt ());
220+ if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
221+ histos.fill (HIST (" hPtXiMinusGlobalWithPV" ), mcParticle.pt ());
222+ }
223+ }
224+ if (mcParticle.pdgCode () == -3312 && findXiPlus) {
225+ reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS, positiveTPC, negativeTPC, bachelorTPC);
226+ if (fabs (mcParticle.y ()) < 0.5 ) {
227+ histos.fill (HIST (" hPtXiPlusGenerated" ), mcParticle.pt ());
228+ if (reconstructed)
229+ histos.fill (HIST (" hPtXiPlusReconstructed" ), mcParticle.pt ());
230+ if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
231+ histos.fill (HIST (" hPtXiPlusGlobal" ), mcParticle.pt ());
232+ if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
233+ histos.fill (HIST (" hPtXiPlusGlobalWithPV" ), mcParticle.pt ());
234+ }
235+ }
236+ if (mcParticle.pdgCode () == 3334 && findOmegaMinus) {
237+ reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS, positiveTPC, negativeTPC, bachelorTPC);
238+ if (fabs (mcParticle.y ()) < 0.5 ) {
239+ histos.fill (HIST (" hPtOmegaMinusGenerated" ), mcParticle.pt ());
240+ if (reconstructed)
241+ histos.fill (HIST (" hPtOmegaMinusReconstructed" ), mcParticle.pt ());
242+ if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
243+ histos.fill (HIST (" hPtOmegaMinusGlobal" ), mcParticle.pt ());
244+ if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
245+ histos.fill (HIST (" hPtOmegaMinusGlobalWithPV" ), mcParticle.pt ());
246+ }
247+ }
248+ if (mcParticle.pdgCode () == -3334 && findOmegaPlus) {
249+ reconstructed = ProcessCascade (mcParticle, tracks, v0s, bestCollisionIndex, positiveITS, negativeITS, bachelorITS, positiveTPC, negativeTPC, bachelorTPC);
250+ if (fabs (mcParticle.y ()) < 0.5 ) {
251+ histos.fill (HIST (" hPtOmegaPlusGenerated" ), mcParticle.pt ());
252+ if (reconstructed)
253+ histos.fill (HIST (" hPtOmegaPlusReconstructed" ), mcParticle.pt ());
254+ if (reconstructed && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
255+ histos.fill (HIST (" hPtOmegaPlusGlobal" ), mcParticle.pt ());
256+ if (reconstructed && bestCollisionIndex >= 0 && positiveITS && negativeITS && bachelorITS && positiveTPC && negativeTPC && bachelorTPC)
257+ histos.fill (HIST (" hPtOmegaPlusGlobalWithPV" ), mcParticle.pt ());
258+ }
259+ }
258260 }
259261 }
260- }
261- }
262262
263- // sort according to collision ID
264- auto sortedIndices = sort_indices (casccollisionId);
263+ // sort according to collision ID
264+ auto sortedIndices = sort_indices (casccollisionId);
265265
266- // V0 list established, populate
267- for (auto ic : sortedIndices) {
268- if (casccollisionId[ic] >= 0 ) {
269- cascades (casccollisionId[ic], cascv0Index[ic], cascbachelorIndex[ic]);
266+ // V0 list established, populate
267+ for (auto ic : sortedIndices) {
268+ if (casccollisionId[ic] >= 0 ) {
269+ cascades (casccollisionId[ic], cascv0Index[ic], cascbachelorIndex[ic]);
270+ }
271+ }
270272 }
271- }
272- }
273273};
274274
275275WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
0 commit comments