@@ -459,15 +459,16 @@ struct TauThreeProngEventTableProducer {
459459 return -1 ;
460460 int nITSbits = 0 ;
461461 int firstThreeLayers = 0 ;
462+ const int threeLayers = 3 ;
462463 uint32_t clusterSizes = track.itsClusterSizes ();
463464 for (int layer = 0 ; layer < 7 ; layer++) {
464465 if ((clusterSizes >> (layer * 4 )) & 0xf ) {
465466 nITSbits++;
466- if (layer < 3 )
467+ if (layer < threeLayers) // 3
467468 firstThreeLayers++;
468469 }
469470 } // end of loop over ITS bits
470- if (firstThreeLayers == 3 )
471+ if (firstThreeLayers == threeLayers) // 3
471472 nITSbits += 10 ;
472473
473474 return nITSbits;
@@ -1012,6 +1013,10 @@ struct TauThreeProngEventTableProducer {
10121013 FullMCUDTracks const & tracks)
10131014 {
10141015 // registrySkim.get<TH1>(HIST("skim/efficiencyMC"))->Fill(0., 1.);
1016+
1017+ const int fourTracks = 4 ;
1018+ const int sixTracks = 6 ;
1019+
10151020 if (verbose)
10161021 LOGF (info, " 0. <MC> UDMcCollision size %d, Collisions size %d, UDtracks %d, UDMcParticles %d" , mcCollisions.size (), collisions.size (), tracks.size (), mcParticles.size ());
10171022
@@ -1139,15 +1144,15 @@ struct TauThreeProngEventTableProducer {
11391144
11401145 registrySkim.get <TH1>(HIST (" skim/nChPartMC" ))->Fill (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ]); // N charged particles from taus
11411146 // check number of charged particles in MC event
1142- if ((nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] != 4 ) && (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] != 6 )) {
1147+ if ((nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] != fourTracks ) && (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] != sixTracks )) {
11431148 if (verbose)
11441149 LOGF (info, " Different from 4/6 charged particles (%d) from both taus. Jump to the next MC event." , nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ]);
11451150 continue ;
11461151 }
1147- registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (3 ., 1 .); // 1+3 (3+3) topology
1148- if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == 4 ) {
1152+ registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (3 ., 1 .); // 1+3 (3+3) topology
1153+ if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == fourTracks ) { // 4
11491154 registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (4 ., 1 .);
1150- } else if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == 6 ) {
1155+ } else if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == sixTracks ) { // 6
11511156 registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (5 ., 1 .);
11521157 }
11531158
@@ -1156,10 +1161,10 @@ struct TauThreeProngEventTableProducer {
11561161 LOGF (info, " At least one daughter particle from taus out of pseudo-rapidity (|eta|<0.9). Jump to the next MC event." );
11571162 continue ;
11581163 }
1159- registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (6 ., 1 .); // particles from tau in |eta|<0.9
1160- if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == 4 ) {
1164+ registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (6 ., 1 .); // particles from tau in |eta|<0.9
1165+ if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == fourTracks ) { // 4
11611166 registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (7 ., 1 .);
1162- } else if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == 6 ) {
1167+ } else if (nChargedDaughtersTau[0 ] + nChargedDaughtersTau[1 ] == sixTracks ) { // 6
11631168 registrySkim.get <TH1>(HIST (" skim/efficiencyMC" ))->Fill (8 ., 1 .);
11641169 }
11651170
@@ -1350,7 +1355,7 @@ struct TauThreeProngEventTableProducer {
13501355 countPi0++;
13511356
13521357 // check whether 1+3 or 3+3 topology is present
1353- if (countDaughters > 6 ) {
1358+ if (countDaughters > sixTracks ) { // 6
13541359 if (verbose)
13551360 LOGF (info, " Truth collision has more than 6 charged daughters from 2 taus. Breaking the daughter loop." );
13561361 // histos.get<TH1>(HIST("Truth/hTroubles"))->Fill(3);
@@ -1470,7 +1475,7 @@ struct TauThreeProngEventTableProducer {
14701475 countPi0++;
14711476
14721477 // check whether 1+3 or 3+3 topology is present
1473- if (countDaughters > 6 ) {
1478+ if (countDaughters > sixTracks ) { // 6
14741479 if (verbose)
14751480 LOGF (info, " Truth collision has more than 6 charged daughters from taus. Breaking the daughter loop." );
14761481 // histos.get<TH1>(HIST("Truth/hTroubles"))->Fill(13);
0 commit comments