@@ -266,7 +266,7 @@ struct kstarInOO {
266266 histos.add (" hGen_pT_GoodEv" , " Gen_pT_GoodTrig (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
267267 histos.add (" hGen_pT_GoodTrig" , " Gen_pT_GoodTrig (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
268268 histos.add (" hGen_pT_GoodEvTrig" , " Gen_pT_GoodEvTrig (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
269-
269+
270270 histos.add (" hEffRec_pT" , " EffRec_pT (GeV/c)" , kTH1F , {{1600 , 0 ., 80 .}});
271271 histos.add (" hEffRecTest1_pT" , " EffRecTest1_pT (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
272272 histos.add (" hEffRecTest2_pT" , " EffRecTest2_pT (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
@@ -277,7 +277,7 @@ struct kstarInOO {
277277 histos.add (" hEffRecTest7_pT" , " EffRecTest7_pT (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
278278 histos.add (" hEffRecTest8_pT" , " EffRecTest8_pT (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
279279 histos.add (" hEffGen_pT" , " EffGen_pT (GeV/c)" , kTH1F , {{800 , 0 ., 40 .}});
280-
280+
281281 histos.add (" hMotherPdg1" , " hMotherPdg1" , kTH1F , {{5000 , 0 ., 5000 .}});
282282 histos.add (" hMotherPdg2" , " hMotherPdg2" , kTH1F , {{5000 , 0 ., 5000 .}});
283283 }
@@ -1510,58 +1510,58 @@ struct kstarInOO {
15101510 std::cout << " Processed MC (GEN) Events: " << nprocessEffEvents << std::endl;
15111511 }
15121512 }
1513- if (cfgJetMCHistos){
1513+ if (cfgJetMCHistos) {
15141514 histos.fill (HIST (" nEvents_Gen" ), 0.5 );
15151515 }
1516-
1516+
15171517 for (auto & particle : mcParticles) {
15181518 if (particle.pdgCode () != 313 )
15191519 continue ;
15201520 if (std::fabs (particle.eta ()) > cfgTrackMaxEta)
15211521 continue ;
15221522 if (fabs (collision.posZ ()) > cfgEventVtxCut)
1523- break ;
1523+ break ;
15241524
1525- if (cfgJetMCHistos){
1526- histos.fill (HIST (" hGen_pT_Raw" ), particle.pt ());
1525+ if (cfgJetMCHistos) {
1526+ histos.fill (HIST (" hGen_pT_Raw" ), particle.pt ());
15271527 }
1528- }// Unreconstructed collisions(=Raw coll) for correction
1528+ } // Unreconstructed collisions(=Raw coll) for correction
15291529
15301530 if (recocolls.size () <= 0 ) { // not reconstructed
15311531 return ;
15321532 }
15331533
15341534 for (auto & recocoll : recocolls) { // poorly reconstructed
1535- if (recocoll.posZ () > cfgEventVtxCut)
1536- continue ;
1535+ if (recocoll.posZ () > cfgEventVtxCut)
1536+ continue ;
15371537 auto goodEv = jetderiveddatautilities::selectCollision (recocoll, eventSelectionBits);
15381538 auto goodTrig = jetderiveddatautilities::selectTrigger (recocoll, RealTriggerMaskBits);
15391539 for (auto & particle : mcParticles) {
1540- if (particle.pdgCode () != 313 )
1541- continue ;
1542- if (std::fabs (particle.eta ()) > cfgTrackMaxEta)
1543- continue ;
1544- if (cfgJetMCHistos) {
1545- // check K* PID
1546- if (goodEv) {
1547- histos.fill (HIST (" hGen_pT_GoodEv" ), particle.pt ());
1548- }
1549- if (goodTrig){
1550- histos.fill (HIST (" hGen_pT_GoodTrig" ), particle.pt ());
1551- }
1552- if (goodEv && goodTrig) {
1553- histos.fill (HIST (" hGen_pT_GoodEvTrig" ), particle.pt ());
1554- }
1555- } // cfgJetMCHistos
1556- }// mcParticles
1557- }// recocolls (=reconstructed collisions)
1540+ if (particle.pdgCode () != 313 )
1541+ continue ;
1542+ if (std::fabs (particle.eta ()) > cfgTrackMaxEta)
1543+ continue ;
1544+ if (cfgJetMCHistos) {
1545+ // check K* PID
1546+ if (goodEv) {
1547+ histos.fill (HIST (" hGen_pT_GoodEv" ), particle.pt ());
1548+ }
1549+ if (goodTrig) {
1550+ histos.fill (HIST (" hGen_pT_GoodTrig" ), particle.pt ());
1551+ }
1552+ if (goodEv && goodTrig) {
1553+ histos.fill (HIST (" hGen_pT_GoodEvTrig" ), particle.pt ());
1554+ }
1555+ } // cfgJetMCHistos
1556+ } // mcParticles
1557+ } // recocolls (=reconstructed collisions)
15581558
15591559 // =================
15601560 // || Efficiency
15611561 // =================
15621562 if (fabs (collision.posZ ()) > cfgEventVtxCut)
15631563 return ;
1564-
1564+
15651565 for (auto & recocoll : recocolls) { // poorly reconstructed
15661566 auto goodEv = jetderiveddatautilities::selectCollision (recocoll, eventSelectionBits);
15671567 if (goodEv) {
@@ -1580,24 +1580,24 @@ struct kstarInOO {
15801580 continue ;
15811581
15821582 /* // Not Yet
1583- if (cfg_Force_BR) {
1584- bool baddecay = false;
1585- for (auto& phidaughter : particle.daughters_as<aod::McParticles>()) {
1586- if (std::fabs(phidaughter.pdgCode()) != 321) {
1587- baddecay = true;
1588- break;
1589- }
1590- if (cfg_Force_Kaon_Acceptence) {
1591- if (std::fabs(phidaughter.eta()) > cfg_Track_MaxEta) {
1592- baddecay = true;
1593- break;
1594- }
1595- }
1596- } // loop over daughters
1597-
1598- if (baddecay)
1599- continue;
1600- } // enforce BR restriction
1583+ if (cfg_Force_BR) {
1584+ bool baddecay = false;
1585+ for (auto& phidaughter : particle.daughters_as<aod::McParticles>()) {
1586+ if (std::fabs(phidaughter.pdgCode()) != 321) {
1587+ baddecay = true;
1588+ break;
1589+ }
1590+ if (cfg_Force_Kaon_Acceptence) {
1591+ if (std::fabs(phidaughter.eta()) > cfg_Track_MaxEta) {
1592+ baddecay = true;
1593+ break;
1594+ }
1595+ }
1596+ } // loop over daughters
1597+
1598+ if (baddecay)
1599+ continue;
1600+ } // enforce BR restriction
16011601 */
16021602
16031603 if (cfgJetMCHistos) {
0 commit comments