@@ -164,17 +164,17 @@ struct strangenesstofpid {
164164 ConfigurableAxis axisNSigma{" axisNSigma" , {200 , -10 .0f , +10 .0f }, " N(#sigma)" };
165165 ConfigurableAxis axisRatioMethods{" axisRatioMethods" , {400 , 0 .9f , 1 .9f }, " T_{method 1}/T_{method 0}" };
166166 ConfigurableAxis axisSnp{" axisSnp" , {220 , -1 .1f , 1 .1f }, " snp" };
167-
167+
168168 // master p axis
169169 ConfigurableAxis axisP{" axisP" , {VARIABLE_WIDTH, 0 .0f , 0 .1f , 0 .2f , 0 .3f , 0 .4f , 0 .5f , 0 .6f , 0 .7f , 0 .8f , 0 .9f , 1 .0f , 1 .1f , 1 .2f , 1 .3f , 1 .4f , 1 .5f , 1 .6f , 1 .7f , 1 .8f , 1 .9f , 2 .0f , 2 .2f , 2 .4f , 2 .6f , 2 .8f , 3 .0f , 3 .2f , 3 .4f , 3 .6f , 3 .8f , 4 .0f , 4 .4f , 4 .8f , 5 .2f , 5 .6f , 6 .0f , 6 .5f , 7 .0f , 7 .5f , 8 .0f , 9 .0f , 10 .0f , 11 .0f , 12 .0f , 13 .0f , 14 .0f , 15 .0f , 17 .0f , 19 .0f , 21 .0f , 23 .0f , 25 .0f , 30 .0f , 35 .0f , 40 .0f , 50 .0f }, " p_{T} (GeV/c)" };
170170
171171 // for zooming in at low values only (e-loss studies and effective correction)
172172 ConfigurableAxis axisSmallP{" axisSmallP" , {250 , 0 .0f , 2 .5f }, " p_{T} (GeV/c)" };
173173
174- // for BC shift QA plots
174+ // for BC shift QA plots
175175 // binning to actually match BC shifts but show in picoseconds
176- const double bcShiftValuePS = o2::constants::lhc::LHCBunchSpacingNS* 1000 .0f ;
177- ConfigurableAxis axisBCshift{" axisBCshift" , {130 , -120 .5f * bcShiftValuePS, 9 .5f * bcShiftValuePS}, " time shift (ps)" };
176+ const double bcShiftValuePS = o2::constants::lhc::LHCBunchSpacingNS * 1000 .0f ;
177+ ConfigurableAxis axisBCshift{" axisBCshift" , {130 , -120 .5f * bcShiftValuePS, 9 .5f * bcShiftValuePS}, " time shift (ps)" };
178178
179179 // very broad time axis
180180 ConfigurableAxis axisTimeLong{" axisTimeLong" , {3000 , -1500000 .0f , 1500000 .0f }, " time (ps)" };
@@ -422,11 +422,11 @@ struct strangenesstofpid {
422422 // per event
423423 histos.add (" hCandidateCounter" , " hCandidateCounter" , kTH1F , {{500 , -0 .5f , 499 .5f }});
424424
425- histos.add (" hV0PositiveBCShift" , " hV0PositiveBCShift" , kTH1F , {axes.axisBCshift });
426- histos.add (" hV0NegativeBCShift" , " hV0NegativeBCShift" , kTH1F , {axes.axisBCshift });
427- histos.add (" hCascadePositiveBCShift" , " hCascadePositiveBCShift" , kTH1F , {axes.axisBCshift });
428- histos.add (" hCascadeNegativeBCShift" , " hCascadeNegativeBCShift" , kTH1F , {axes.axisBCshift });
429- histos.add (" hCascadeBachelorBCShift" , " hCascadeBachelorBCShift" , kTH1F , {axes.axisBCshift });
425+ histos.add (" hV0PositiveBCShift" , " hV0PositiveBCShift" , kTH1F , {axes.axisBCshift });
426+ histos.add (" hV0NegativeBCShift" , " hV0NegativeBCShift" , kTH1F , {axes.axisBCshift });
427+ histos.add (" hCascadePositiveBCShift" , " hCascadePositiveBCShift" , kTH1F , {axes.axisBCshift });
428+ histos.add (" hCascadeNegativeBCShift" , " hCascadeNegativeBCShift" , kTH1F , {axes.axisBCshift });
429+ histos.add (" hCascadeBachelorBCShift" , " hCascadeBachelorBCShift" , kTH1F , {axes.axisBCshift });
430430
431431 histos.add (" hTOFSignalPositive" , " hTOFSignalPositive" , kTH1F , {axes.axisTimeLong });
432432 histos.add (" hTOFSignalNegative" , " hTOFSignalNegative" , kTH1F , {axes.axisTimeLong });
@@ -951,14 +951,14 @@ struct strangenesstofpid {
951951 // int incompatibilityReason = 0;
952952 // if (std::abs(v0tof.nSigmaPositiveK0ShortPi - o2::aod::v0data::kNoTOFValue) > o2::aod::v0data::kEpsilon && std::abs(v0tof.nSigmaPositiveK0ShortPi) > 4) {
953953 // compatibleK0Short = false; // reject only if info present and incompatible
954- // incompatibilityReason += 1;
954+ // incompatibilityReason += 1;
955955 // }
956956 // if (std::abs(v0tof.nSigmaNegativeK0ShortPi - o2::aod::v0data::kNoTOFValue) > o2::aod::v0data::kEpsilon && std::abs(v0tof.nSigmaNegativeK0ShortPi) > 4) {
957957 // compatibleK0Short = false; // reject only if info present and incompatible
958- // incompatibilityReason += 2;
958+ // incompatibilityReason += 2;
959959 // }
960960
961- // if(!compatibleK0Short && passesQAcuts && k0ShortCandidate){
961+ // if(!compatibleK0Short && passesQAcuts && k0ShortCandidate){
962962 // histos.fill(HIST("hIncompatibilityReason"), incompatibilityReason);
963963 // // LOGF(info, "Incompatible K0, sigmas = (%.2f %.2f), lengths = (%.2f %.2f) tofSignals = (%.2f %.2f) evtimes = (%.2f %.2f)", v0tof.nSigmaPositiveK0ShortPi, v0tof.nSigmaNegativeK0ShortPi, lengthPositive, lengthNegative, pTof.tofSignal, nTof.tofSignal, pTof.tofEvTime, nTof.tofEvTime);
964964 // }
@@ -1367,34 +1367,34 @@ struct strangenesstofpid {
13671367
13681368 auto pTra = V0.posTrack_as <TracksWithAllExtras>();
13691369 auto nTra = V0.negTrack_as <TracksWithAllExtras>();
1370- double deltaTimePos = 0 .0f ;
1371- double deltaTimeNeg = 0 .0f ;
1370+ double deltaTimePos = 0 .0f ;
1371+ double deltaTimeNeg = 0 .0f ;
13721372
13731373 auto collisionV0 = collisions.rawIteratorAt (V0.collisionId ());
13741374 auto bcV0 = bcs.rawIteratorAt (collisionV0.bcId ());
13751375
1376- if (pTra.collisionId ()>= 0 ){
1376+ if (pTra.collisionId () >= 0 ) {
13771377 auto collisionPos = collisions.rawIteratorAt (pTra.collisionId ());
13781378 auto bcPos = bcs.rawIteratorAt (collisionPos.bcId ());
13791379 const int64_t deltaBcPos = bcPos.globalBC () - bcV0.globalBC ();
1380- deltaTimePos = o2::constants::lhc::LHCBunchSpacingNS* deltaBcPos* 1000 .0f ;
1380+ deltaTimePos = o2::constants::lhc::LHCBunchSpacingNS * deltaBcPos * 1000 .0f ;
13811381 histos.fill (HIST (" hV0PositiveBCShift" ), deltaTimePos);
13821382 }
13831383
1384- if (nTra.collisionId ()>= 0 ){
1384+ if (nTra.collisionId () >= 0 ) {
13851385 auto collisionNeg = collisions.rawIteratorAt (nTra.collisionId ());
13861386 auto bcNeg = bcs.rawIteratorAt (collisionNeg.bcId ());
13871387 const int64_t deltaBcNeg = bcNeg.globalBC () - bcV0.globalBC ();
1388- deltaTimeNeg = o2::constants::lhc::LHCBunchSpacingNS* deltaBcNeg* 1000 .0f ;
1388+ deltaTimeNeg = o2::constants::lhc::LHCBunchSpacingNS * deltaBcNeg * 1000 .0f ;
13891389 histos.fill (HIST (" hV0NegativeBCShift" ), deltaTimeNeg);
13901390 }
13911391
1392- if (pTof.hasTOF ){
1392+ if (pTof.hasTOF ) {
13931393 histos.fill (HIST (" hTOFSignalPositive" ), pTof.tofSignal );
13941394 histos.fill (HIST (" h2dTOFSignalPositive" ), pTof.tofSignal , deltaTimePos);
13951395 }
1396-
1397- if (nTof.hasTOF ){
1396+
1397+ if (nTof.hasTOF ) {
13981398 histos.fill (HIST (" hTOFSignalNegative" ), nTof.tofSignal );
13991399 histos.fill (HIST (" h2dTOFSignalNegative" ), nTof.tofSignal , deltaTimeNeg);
14001400 }
@@ -1454,34 +1454,34 @@ struct strangenesstofpid {
14541454 auto nTra = cascade.negTrack_as <TracksWithAllExtras>();
14551455 auto bTra = cascade.bachelor_as <TracksWithAllExtras>();
14561456
1457- double deltaTimePos = 0 .0f ;
1458- double deltaTimeNeg = 0 .0f ;
1459- double deltaTimeBach = 0 .0f ;
1457+ double deltaTimePos = 0 .0f ;
1458+ double deltaTimeNeg = 0 .0f ;
1459+ double deltaTimeBach = 0 .0f ;
14601460
14611461 auto collisionCascade = collisions.rawIteratorAt (cascade.collisionId ());
14621462 auto bcV0 = bcs.rawIteratorAt (collisionCascade.bcId ());
14631463
1464- if (pTra.collisionId ()>= 0 ){
1464+ if (pTra.collisionId () >= 0 ) {
14651465 auto collisionPos = collisions.rawIteratorAt (pTra.collisionId ());
14661466 auto bcPos = bcs.rawIteratorAt (collisionPos.bcId ());
14671467 const int64_t deltaBcPos = bcPos.globalBC () - bcV0.globalBC ();
1468- deltaTimePos = o2::constants::lhc::LHCBunchSpacingNS* deltaBcPos* 1000 .0f ;
1468+ deltaTimePos = o2::constants::lhc::LHCBunchSpacingNS * deltaBcPos * 1000 .0f ;
14691469 histos.fill (HIST (" hCascadePositiveBCShift" ), deltaTimePos);
14701470 }
14711471
1472- if (nTra.collisionId ()>= 0 ){
1472+ if (nTra.collisionId () >= 0 ) {
14731473 auto collisionNeg = collisions.rawIteratorAt (nTra.collisionId ());
14741474 auto bcNeg = bcs.rawIteratorAt (collisionNeg.bcId ());
14751475 const int64_t deltaBcNeg = bcNeg.globalBC () - bcV0.globalBC ();
1476- deltaTimeNeg = o2::constants::lhc::LHCBunchSpacingNS* deltaBcNeg* 1000 .0f ;
1476+ deltaTimeNeg = o2::constants::lhc::LHCBunchSpacingNS * deltaBcNeg * 1000 .0f ;
14771477 histos.fill (HIST (" hCascadeNegativeBCShift" ), deltaTimeNeg);
14781478 }
14791479
1480- if (bTra.collisionId ()>= 0 ){
1480+ if (bTra.collisionId () >= 0 ) {
14811481 auto collisionBach = collisions.rawIteratorAt (bTra.collisionId ());
14821482 auto bcBach = bcs.rawIteratorAt (collisionBach.bcId ());
14831483 const int64_t deltaBcBach = bcBach.globalBC () - bcV0.globalBC ();
1484- deltaTimeBach = o2::constants::lhc::LHCBunchSpacingNS* deltaBcBach* 1000 .0f ;
1484+ deltaTimeBach = o2::constants::lhc::LHCBunchSpacingNS * deltaBcBach * 1000 .0f ;
14851485 histos.fill (HIST (" hCascadeNegativeBCShift" ), deltaTimeBach);
14861486 }
14871487
@@ -1594,11 +1594,11 @@ struct strangenesstofpid {
15941594 if (tofIndices[V0.posTrackExtraId ()] >= 0 && collision.eventTime () > -1e+5 ) {
15951595 auto pTofExt = dauTrackTOFPIDs.rawIteratorAt (tofIndices[V0.posTrackExtraId ()]);
15961596
1597- if (pTofExt.straCollisionId ()>= 0 ) {
1598- // extract BC for BC time shift
1597+ if (pTofExt.straCollisionId () >= 0 ) {
1598+ // extract BC for BC time shift
15991599 auto collisionTrack = collisions.rawIteratorAt (pTofExt.straCollisionId ());
16001600 const int64_t deltaBc = collisionTrack.globalBC () - collision.globalBC ();
1601- const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS* deltaBc* 1000 .0f ;
1601+ const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS * deltaBc * 1000 .0f ;
16021602 histos.fill (HIST (" hV0PositiveBCShift" ), deltaTimeBc);
16031603 deltaTimeBcPos = deltaTimeBc;
16041604
@@ -1619,11 +1619,11 @@ struct strangenesstofpid {
16191619 if (tofIndices[V0.negTrackExtraId ()] >= 0 && collision.eventTime () > -1e+5 ) {
16201620 auto nTofExt = dauTrackTOFPIDs.rawIteratorAt (tofIndices[V0.negTrackExtraId ()]);
16211621
1622- if (nTofExt.straCollisionId ()>= 0 ) {
1623- // extract BC for BC time shift
1622+ if (nTofExt.straCollisionId () >= 0 ) {
1623+ // extract BC for BC time shift
16241624 auto collisionTrack = collisions.rawIteratorAt (nTofExt.straCollisionId ());
16251625 const int64_t deltaBc = collisionTrack.globalBC () - collision.globalBC ();
1626- const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS* deltaBc* 1000 .0f ;
1626+ const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS * deltaBc * 1000 .0f ;
16271627 histos.fill (HIST (" hV0NegativeBCShift" ), deltaTimeBc);
16281628 deltaTimeBcNeg = deltaTimeBc;
16291629
@@ -1635,12 +1635,12 @@ struct strangenesstofpid {
16351635 nTof.length = nTofExt.length ();
16361636 }
16371637 }
1638- if (pTof.hasTOF ){
1638+ if (pTof.hasTOF ) {
16391639 histos.fill (HIST (" hTOFSignalPositive" ), pTof.tofSignal );
16401640 histos.fill (HIST (" h2dTOFSignalPositive" ), pTof.tofSignal , deltaTimeBcPos);
16411641 }
1642-
1643- if (nTof.hasTOF ){
1642+
1643+ if (nTof.hasTOF ) {
16441644 histos.fill (HIST (" hTOFSignalNegative" ), nTof.tofSignal );
16451645 histos.fill (HIST (" h2dTOFSignalNegative" ), nTof.tofSignal , deltaTimeBcNeg);
16461646 }
@@ -1687,11 +1687,11 @@ struct strangenesstofpid {
16871687 if (tofIndices[cascade.posTrackExtraId ()] >= 0 && collision.eventTime () > -1e+5 ) {
16881688 auto pTofExt = dauTrackTOFPIDs.rawIteratorAt (tofIndices[cascade.posTrackExtraId ()]);
16891689
1690- if (pTofExt.straCollisionId ()>= 0 ) {
1691- // extract BC for BC time shift
1690+ if (pTofExt.straCollisionId () >= 0 ) {
1691+ // extract BC for BC time shift
16921692 auto collisionTrack = collisions.rawIteratorAt (pTofExt.straCollisionId ());
16931693 const int64_t deltaBc = collisionTrack.globalBC () - collision.globalBC ();
1694- const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS* deltaBc* 1000 .0f ;
1694+ const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS * deltaBc * 1000 .0f ;
16951695 histos.fill (HIST (" hCascadePositiveBCShift" ), deltaTimeBc);
16961696
16971697 pTof.collisionId = pTofExt.straCollisionId ();
@@ -1710,11 +1710,11 @@ struct strangenesstofpid {
17101710 if (tofIndices[cascade.negTrackExtraId ()] >= 0 && collision.eventTime () > -1e+5 ) {
17111711 auto nTofExt = dauTrackTOFPIDs.rawIteratorAt (tofIndices[cascade.negTrackExtraId ()]);
17121712
1713- if (nTofExt.straCollisionId ()>= 0 ) {
1714- // extract BC for BC time shift
1713+ if (nTofExt.straCollisionId () >= 0 ) {
1714+ // extract BC for BC time shift
17151715 auto collisionTrack = collisions.rawIteratorAt (nTofExt.straCollisionId ());
17161716 const int64_t deltaBc = collisionTrack.globalBC () - collision.globalBC ();
1717- const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS* deltaBc* 1000 .0f ;
1717+ const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS * deltaBc * 1000 .0f ;
17181718 histos.fill (HIST (" hCascadeNegativeBCShift" ), deltaTimeBc);
17191719
17201720 nTof.collisionId = nTofExt.straCollisionId ();
@@ -1733,11 +1733,11 @@ struct strangenesstofpid {
17331733 if (tofIndices[cascade.bachTrackExtraId ()] >= 0 && collision.eventTime () > -1e+5 ) {
17341734 auto bTofExt = dauTrackTOFPIDs.rawIteratorAt (tofIndices[cascade.bachTrackExtraId ()]);
17351735
1736- if (bTofExt.straCollisionId ()>= 0 ) {
1737- // extract BC for BC time shift
1736+ if (bTofExt.straCollisionId () >= 0 ) {
1737+ // extract BC for BC time shift
17381738 auto collisionTrack = collisions.rawIteratorAt (bTofExt.straCollisionId ());
17391739 const int64_t deltaBc = collisionTrack.globalBC () - collision.globalBC ();
1740- const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS* deltaBc* 1000 .0f ;
1740+ const double deltaTimeBc = o2::constants::lhc::LHCBunchSpacingNS * deltaBc * 1000 .0f ;
17411741 histos.fill (HIST (" hCascadeBachelorBCShift" ), deltaTimeBc);
17421742
17431743 bTof.collisionId = bTofExt.straCollisionId ();
0 commit comments