diff --git a/ALICE3/Tools/handleParamTOFResoALICE3.cxx b/ALICE3/Tools/handleParamTOFResoALICE3.cxx index 8e309a0d0f8..83f4cbe9425 100644 --- a/ALICE3/Tools/handleParamTOFResoALICE3.cxx +++ b/ALICE3/Tools/handleParamTOFResoALICE3.cxx @@ -70,7 +70,7 @@ int main(int argc, char* argv[]) const unsigned int mode = vm["mode"].as(); const std::string path = vm["ccdb-path"].as(); std::map metadata; - std::map* headers; + std::map* headers = nullptr; o2::ccdb::CcdbApi api; const std::string url = vm["url"].as(); api.init(url); diff --git a/Common/Core/EventMixing.h b/Common/Core/EventMixing.h index 8d227ae2afd..9317a1213af 100644 --- a/Common/Core/EventMixing.h +++ b/Common/Core/EventMixing.h @@ -33,9 +33,9 @@ static int getMixingBin(const T1& vtxBins, const T1& multBins, const T2& vtx, co return -1; } - for (int i = 1; i < vtxBins.size(); i++) { + for (unsigned int i = 1; i < vtxBins.size(); i++) { if (vtx < vtxBins.at(i)) { - for (int j = 1; j < multBins.size(); j++) { + for (unsigned int j = 1; j < multBins.size(); j++) { if (mult < multBins.at(j)) { return i + j * (vtxBins.size() + 1); } diff --git a/Common/TableProducer/eventSelection.cxx b/Common/TableProducer/eventSelection.cxx index b9a7d8e887f..7f9c6e7a610 100644 --- a/Common/TableProducer/eventSelection.cxx +++ b/Common/TableProducer/eventSelection.cxx @@ -431,15 +431,15 @@ struct EventSelectionTask { int64_t foundFT0 = bc.foundFT0(); if (foundFT0 < 0) { // search in +/-4 sigma around meanBC - uint64_t apprBC = bc.globalBC(); - uint64_t meanBC = apprBC - std::lround(col.collisionTime() / o2::constants::lhc::LHCBunchSpacingNS); + int64_t apprBC = bc.globalBC(); + int64_t meanBC = apprBC - std::lround(col.collisionTime() / o2::constants::lhc::LHCBunchSpacingNS); int64_t deltaBC = std::ceil(col.collisionTimeRes() / o2::constants::lhc::LHCBunchSpacingNS * 4); // search forward int forwardMoveCount = 0; int64_t forwardBcDist = deltaBC + 1; - for (; bc != bcs.end() && bc.globalBC() - meanBC <= deltaBC; ++bc, ++forwardMoveCount) { + for (; bc != bcs.end() && apprBC - meanBC <= deltaBC; ++bc, ++forwardMoveCount) { if (bc.foundFT0() >= 0) { - forwardBcDist = bc.globalBC() - meanBC; + forwardBcDist = apprBC - meanBC; break; } } @@ -447,9 +447,9 @@ struct EventSelectionTask { // search backward int backwardMoveCount = 0; int64_t backwardBcDist = deltaBC + 1; - for (; bc != bcs.begin() && bc.globalBC() - meanBC >= -deltaBC; --bc, --backwardMoveCount) { + for (; bc != bcs.begin() && apprBC - meanBC >= -deltaBC; --bc, --backwardMoveCount) { if (bc.foundFT0() >= 0) { - backwardBcDist = meanBC - bc.globalBC(); + backwardBcDist = meanBC - apprBC; break; } } diff --git a/Common/Tools/handleParamTOFReso.cxx b/Common/Tools/handleParamTOFReso.cxx index 68521848e40..86df6b9a532 100644 --- a/Common/Tools/handleParamTOFReso.cxx +++ b/Common/Tools/handleParamTOFReso.cxx @@ -74,7 +74,7 @@ int main(int argc, char* argv[]) const unsigned int mode = vm["mode"].as(); const std::string path = vm["ccdb-path"].as(); std::map metadata; - std::map* headers; + std::map* headers = nullptr; o2::ccdb::CcdbApi api; const std::string url = vm["url"].as(); api.init(url); diff --git a/PWGLF/TableProducer/cascadebuilder.cxx b/PWGLF/TableProducer/cascadebuilder.cxx index 8c1fd0c3d28..ea3fc57bc98 100644 --- a/PWGLF/TableProducer/cascadebuilder.cxx +++ b/PWGLF/TableProducer/cascadebuilder.cxx @@ -75,7 +75,7 @@ struct cascadebuilder { OutputObj hCascFiltered{TH1F("hCascFiltered", "", 15, 0, 15)}; OutputObj hCascCandidate{TH1F("hCascCandidate", "", 10, 0, 10)}; - //Configurables + // Configurables Configurable d_bz{"d_bz", -5.0, "bz field"}; Configurable d_UseAbsDCA{"d_UseAbsDCA", true, "Use Abs DCAs"}; @@ -92,7 +92,7 @@ struct cascadebuilder { void process(aod::MatchedV0Cascades const& MatchedV0Cascades, aod::V0Datas const&, aod::Cascades const&, aod::Collisions const&, soa::Join const&) { - //Define o2 fitter, 2-prong + // Define o2 fitter, 2-prong o2::vertexing::DCAFitterN<2> fitterV0, fitterCasc; fitterV0.setBz(d_bz); fitterV0.setPropagateToPCA(true); @@ -120,19 +120,19 @@ struct cascadebuilder { std::array pVtx = {v0.collision().posX(), v0.collision().posY(), v0.collision().posZ()}; - auto b = casc.bachelor_as(); + // auto b = casc.bachelor_as(); if (tpcrefit) { if (!(v0.posTrack_as().trackType() & o2::aod::track::TPCrefit)) { - continue; //TPC refit + continue; // TPC refit } hCascFiltered->Fill(1.5); if (!(v0.negTrack_as().trackType() & o2::aod::track::TPCrefit)) { - continue; //TPC refit + continue; // TPC refit } hCascFiltered->Fill(2.5); if (!(casc.bachelor_as().trackType() & o2::aod::track::TPCrefit)) { - continue; //TPC refit + continue; // TPC refit } hCascFiltered->Fill(3.5); } @@ -161,7 +161,7 @@ struct cascadebuilder { } hCascFiltered->Fill(9.5); - //V0 selections + // V0 selections if (fabs(v0.mLambda() - 1.116) > lambdamasswindow && fabs(v0.mAntiLambda() - 1.116) > lambdamasswindow) { continue; } @@ -192,7 +192,7 @@ struct cascadebuilder { hCascCandidate->Fill(0.5); - //Acquire basic tracks + // Acquire basic tracks auto pTrack = getTrackParCov(v0.posTrack_as()); auto nTrack = getTrackParCov(v0.negTrack_as()); auto bTrack = getTrackParCov(casc.bachelor_as()); @@ -213,7 +213,7 @@ struct cascadebuilder { std::array cov1 = {0}; std::array covV0 = {0}; - //Covariance matrix calculation + // Covariance matrix calculation const int momInd[6] = {9, 13, 14, 18, 19, 20}; // cov matrix elements for momentum component fitterV0.getTrack(0).getPxPyPzGlo(pvecpos); fitterV0.getTrack(1).getPxPyPzGlo(pvecneg); @@ -235,7 +235,7 @@ struct cascadebuilder { const std::array momentum = {pvecpos[0] + pvecneg[0], pvecpos[1] + pvecneg[1], pvecpos[2] + pvecneg[2]}; auto tV0 = o2::track::TrackParCov(vertex, momentum, covV0, 0); - tV0.setQ2Pt(0); //No bending, please + tV0.setQ2Pt(0); // No bending, please int nCand2 = fitterCasc.process(tV0, bTrack); if (nCand2 != 0) { fitterCasc.propagateTracksToVertex(); @@ -245,9 +245,9 @@ struct cascadebuilder { posXi[i] = cascvtx[i]; } fitterCasc.getTrack(1).getPxPyPzGlo(pvecbach); - } //end if cascade recoed - } //end if v0 recoed - //Fill table, please + } // end if cascade recoed + } // end if v0 recoed + // Fill table, please cascdata( v0.globalIndex(), casc.bachelor_as().globalIndex(), diff --git a/PWGLF/TableProducer/cascadefinder.cxx b/PWGLF/TableProducer/cascadefinder.cxx index 264780e4888..4bd3c0e538b 100644 --- a/PWGLF/TableProducer/cascadefinder.cxx +++ b/PWGLF/TableProducer/cascadefinder.cxx @@ -122,7 +122,7 @@ struct cascadeprefilter { { for (auto& t0 : goodPosTracks) { if (!(t0.trackType() & o2::aod::track::TPCrefit)) { - continue; //TPC refit + continue; // TPC refit } if (t0.tpcNClsCrossedRows() < mincrossedrows) { continue; @@ -131,7 +131,7 @@ struct cascadeprefilter { } for (auto& t0 : goodNegTracks) { if (!(t0.trackType() & o2::aod::track::TPCrefit)) { - continue; //TPC refit + continue; // TPC refit } if (t0.tpcNClsCrossedRows() < mincrossedrows) { continue; @@ -167,16 +167,16 @@ struct cascadefinder { OutputObj hCandPerEvent{TH1F("hCandPerEvent", "", 100, 0, 100)}; - //Configurables + // Configurables Configurable d_bz{"d_bz", +5.0, "bz field"}; Configurable d_UseAbsDCA{"d_UseAbsDCA", kTRUE, "Use Abs DCAs"}; - //Selection criteria - Configurable v0cospa{"casccospa", 0.998, "Casc CosPA"}; //double -> N.B. dcos(x)/dx = 0 at x=0) + // Selection criteria + Configurable v0cospa{"casccospa", 0.998, "Casc CosPA"}; // double -> N.B. dcos(x)/dx = 0 at x=0) Configurable dcav0dau{"dcacascdau", 1.0, "DCA Casc Daughters"}; Configurable v0radius{"cascradius", 1.0, "cascradius"}; - //Process: subscribes to a lot of things! + // Process: subscribes to a lot of things! void process(aod::Collision const& collision, soa::Join const& tracks, aod::V0Datas const& V0s, @@ -185,7 +185,7 @@ struct cascadefinder { aod::CascGoodPosTracks const& pBachtracks, aod::CascGoodNegTracks const& nBachtracks) { - //Define o2 fitter, 2-prong + // Define o2 fitter, 2-prong o2::vertexing::DCAFitterN<2> fitterV0, fitterCasc; fitterV0.setBz(d_bz); fitterV0.setPropagateToPCA(true); @@ -207,20 +207,20 @@ struct cascadefinder { Long_t lNCand = 0; - std::array pVtx = {collision.posX(), collision.posY(), collision.posZ()}; + // std::array pVtx = {collision.posX(), collision.posY(), collision.posZ()}; std::array pos = {0.}; std::array posXi = {0.}; std::array pvecpos = {0.}; std::array pvecneg = {0.}; std::array pvecbach = {0.}; - //Cascades first + // Cascades first for (auto& v0id : lambdas) { - //required: de-reference the tracks for cascade building + // required: de-reference the tracks for cascade building auto v0 = v0id.goodLambda(); auto pTrack = getTrackParCov(v0.posTrack_as>()); auto nTrack = getTrackParCov(v0.negTrack_as>()); - //Let's do the slow part first: the V0 recalculation from scratch + // Let's do the slow part first: the V0 recalculation from scratch int nCand = fitterV0.process(pTrack, nTrack); if (nCand != 0) { fitterV0.propagateTracksToVertex(); @@ -233,7 +233,7 @@ struct cascadefinder { std::array cov1 = {0}; std::array covV0 = {0}; - //Covariance matrix calculation + // Covariance matrix calculation const int momInd[6] = {9, 13, 14, 18, 19, 20}; // cov matrix elements for momentum component fitterV0.getTrack(0).getPxPyPzGlo(pvecpos); fitterV0.getTrack(1).getPxPyPzGlo(pvecneg); @@ -255,7 +255,7 @@ struct cascadefinder { const std::array momentum = {pvecpos[0] + pvecneg[0], pvecpos[1] + pvecneg[1], pvecpos[2] + pvecneg[2]}; auto tV0 = o2::track::TrackParCov(vertex, momentum, covV0, 0); - tV0.setQ2Pt(0); //No bending, please + tV0.setQ2Pt(0); // No bending, please for (auto& t0id : nBachtracks) { auto t0 = t0id.goodNegTrack_as>(); @@ -271,7 +271,7 @@ struct cascadefinder { fitterCasc.getTrack(1).getPxPyPzGlo(pvecbach); lNCand++; - //If we got here, it means this is a good candidate! + // If we got here, it means this is a good candidate! cascdata(v0.globalIndex(), v0.posTrack().globalIndex(), v0.negTrack().collisionId(), -1, posXi[0], posXi[1], posXi[2], pos[0], pos[1], pos[2], pvecpos[0], pvecpos[1], pvecpos[2], @@ -281,18 +281,18 @@ struct cascadefinder { v0.dcapostopv(), v0.dcanegtopv(), t0id.dcaXY()); - } //end if cascade recoed - } //end loop over bachelor - } //end if v0 recoed - } //end loop over cascades + } // end if cascade recoed + } // end loop over bachelor + } // end if v0 recoed + } // end loop over cascades - //Anticascades + // Anticascades for (auto& v0id : antiLambdas) { - //required: de-reference the tracks for cascade building + // required: de-reference the tracks for cascade building auto v0 = v0id.goodAntiLambda(); auto pTrack = getTrackParCov(v0.posTrack_as>()); auto nTrack = getTrackParCov(v0.negTrack_as>()); - //Let's do the slow part first: the V0 recalculation from scratch + // Let's do the slow part first: the V0 recalculation from scratch int nCand = fitterV0.process(pTrack, nTrack); if (nCand != 0) { fitterV0.propagateTracksToVertex(); @@ -305,7 +305,7 @@ struct cascadefinder { std::array cov1 = {0}; std::array covV0 = {0}; - //Covariance matrix calculation + // Covariance matrix calculation const int momInd[6] = {9, 13, 14, 18, 19, 20}; // cov matrix elements for momentum component fitterV0.getTrack(0).getPxPyPzGlo(pvecpos); fitterV0.getTrack(1).getPxPyPzGlo(pvecneg); @@ -327,7 +327,7 @@ struct cascadefinder { const std::array momentum = {pvecpos[0] + pvecneg[0], pvecpos[1] + pvecneg[1], pvecpos[2] + pvecneg[2]}; auto tV0 = o2::track::TrackParCov(vertex, momentum, covV0, 0); - tV0.setQ2Pt(0); //No bending, please + tV0.setQ2Pt(0); // No bending, please for (auto& t0id : pBachtracks) { auto t0 = t0id.goodPosTrack_as>(); @@ -343,7 +343,7 @@ struct cascadefinder { fitterCasc.getTrack(1).getPxPyPzGlo(pvecbach); lNCand++; - //If we got here, it means this is a good candidate! + // If we got here, it means this is a good candidate! cascdata(v0.globalIndex(), v0.posTrack().globalIndex(), v0.negTrack().collisionId(), +1, posXi[0], posXi[1], posXi[2], pos[0], pos[1], pos[2], pvecpos[0], pvecpos[1], pvecpos[2], @@ -353,25 +353,25 @@ struct cascadefinder { v0.dcapostopv(), v0.dcanegtopv(), t0id.dcaXY()); - } //end if cascade recoed - } //end loop over bachelor - } //end if v0 recoed - } //end loop over anticascades + } // end if cascade recoed + } // end loop over bachelor + } // end if v0 recoed + } // end loop over anticascades hCandPerEvent->Fill(lNCand); } }; struct cascadefinderQA { - //Basic checks + // Basic checks OutputObj h3dMassXiMinus{TH3F("h3dMassXiMinus", "", 20, 0, 100, 200, 0, 10, 200, 1.322 - 0.100, 1.322 + 0.100)}; OutputObj h3dMassXiPlus{TH3F("h3dMassXiPlus", "", 20, 0, 100, 200, 0, 10, 200, 1.322 - 0.100, 1.322 + 0.100)}; OutputObj h3dMassOmegaMinus{TH3F("h3dMassOmegaMinus", "", 20, 0, 100, 200, 0, 10, 200, 1.672 - 0.100, 1.672 + 0.100)}; OutputObj h3dMassOmegaPlus{TH3F("h3dMassOmegaPlus", "", 20, 0, 100, 200, 0, 10, 200, 1.672 - 0.100, 1.672 + 0.100)}; - //Selection criteria - Configurable v0cospa{"v0cospa", 0.999, "V0 CosPA"}; //double -> N.B. dcos(x)/dx = 0 at x=0) - Configurable casccospa{"casccospa", 0.999, "Casc CosPA"}; //double -> N.B. dcos(x)/dx = 0 at x=0) + // Selection criteria + Configurable v0cospa{"v0cospa", 0.999, "V0 CosPA"}; // double -> N.B. dcos(x)/dx = 0 at x=0) + Configurable casccospa{"casccospa", 0.999, "Casc CosPA"}; // double -> N.B. dcos(x)/dx = 0 at x=0) Configurable dcav0dau{"dcav0dau", 1.0, "DCA V0 Daughters"}; Configurable dcacascdau{"dcacascdau", .3, "DCA Casc Daughters"}; Configurable dcanegtopv{"dcanegtopv", .1, "DCA Neg To PV"}; @@ -387,7 +387,7 @@ struct cascadefinderQA { aod::cascdata::dcabachtopv > dcabachtopv&& aod::cascdata::dcaV0daughters < dcav0dau&& aod::cascdata::dcacascdaughters < dcacascdau; - ///Connect to CascFinderData: newly indexed, note: CascDataExt table incompatible with standard V0 table! + /// Connect to CascFinderData: newly indexed, note: CascDataExt table incompatible with standard V0 table! void process(soa::Join::iterator const& collision, soa::Filtered const& Cascades) { if (!collision.alias()[kINT7]) { @@ -397,13 +397,13 @@ struct cascadefinderQA { return; } for (auto& casc : Cascades) { - //FIXME: dynamic columns cannot be filtered on? + // FIXME: dynamic columns cannot be filtered on? if (casc.v0radius() > v0radius && casc.cascradius() > cascradius && casc.v0cosPA(collision.posX(), collision.posY(), collision.posZ()) > v0cospa && casc.casccosPA(collision.posX(), collision.posY(), collision.posZ()) > casccospa && casc.dcav0topv(collision.posX(), collision.posY(), collision.posZ()) > dcav0topv) { - if (casc.sign() < 0) { //FIXME: could be done better... + if (casc.sign() < 0) { // FIXME: could be done better... if (TMath::Abs(casc.yXi()) < 0.5) { h3dMassXiMinus->Fill(collision.centV0M(), casc.pt(), casc.mXi()); } diff --git a/PWGLF/TableProducer/lambdakzerobuilder.cxx b/PWGLF/TableProducer/lambdakzerobuilder.cxx index 74690bb5208..a5314128827 100644 --- a/PWGLF/TableProducer/lambdakzerobuilder.cxx +++ b/PWGLF/TableProducer/lambdakzerobuilder.cxx @@ -61,7 +61,7 @@ using namespace o2::framework; using namespace o2::framework::expressions; using std::array; -//This table stores a filtered list of valid V0 indices +// This table stores a filtered list of valid V0 indices namespace o2::aod { namespace v0goodindices @@ -89,9 +89,9 @@ using MyTracks = FullTracksExt; #define MY_DEBUG_MSG(condition, cmd) #endif -//This prefilter creates a skimmed list of good V0s to re-reconstruct so that -//CPU is saved in case there are specific selections that are to be done -//Note: more configurables, more options to be added as needed +// This prefilter creates a skimmed list of good V0s to re-reconstruct so that +// CPU is saved in case there are specific selections that are to be done +// Note: more configurables, more options to be added as needed struct lambdakzeroprefilterpairs { Configurable dcanegtopv{"dcanegtopv", .1, "DCA Neg To PV"}; Configurable dcapostopv{"dcapostopv", .1, "DCA Pos To PV"}; @@ -134,11 +134,11 @@ struct lambdakzeroprefilterpairs { if (isRun2) { if (!(V0.posTrack_as().trackType() & o2::aod::track::TPCrefit)) { MY_DEBUG_MSG(isK0SfromLc, LOG(info) << "posTrack " << labelPos << " has no TPC refit"); - continue; //TPC refit + continue; // TPC refit } if (!(V0.negTrack_as().trackType() & o2::aod::track::TPCrefit)) { MY_DEBUG_MSG(isK0SfromLc, LOG(info) << "negTrack " << labelNeg << " has no TPC refit"); - continue; //TPC refit + continue; // TPC refit } } registry.fill(HIST("hGoodIndices"), 1.5); @@ -181,12 +181,12 @@ struct lambdakzerobuilder { }, }; - //Configurables + // Configurables Configurable d_bz{"d_bz", -5.0, "bz field"}; // Configurable d_UseAbsDCA{"d_UseAbsDCA", 1, "Use Abs DCAs"}; uncomment this once we want to use the weighted DCA - //Selection criteria - Configurable v0cospa{"v0cospa", 0.995, "V0 CosPA"}; //double -> N.B. dcos(x)/dx = 0 at x=0) + // Selection criteria + Configurable v0cospa{"v0cospa", 0.995, "V0 CosPA"}; // double -> N.B. dcos(x)/dx = 0 at x=0) Configurable dcav0dau{"dcav0dau", 1.0, "DCA V0 Daughters"}; Configurable v0radius{"v0radius", 5.0, "v0radius"}; @@ -207,7 +207,7 @@ struct lambdakzerobuilder { #endif ) { - //Define o2 fitter, 2-prong + // Define o2 fitter, 2-prong o2::vertexing::DCAFitterN<2> fitter; fitter.setBz(d_bz); fitter.setPropagateToPCA(true); @@ -219,7 +219,7 @@ struct lambdakzerobuilder { fitter.setUseAbsDCA(true); // use d_UseAbsDCA once we want to use the weighted DCA registry.fill(HIST("hEventCounter"), 0.5); - std::array pVtx = {collision.posX(), collision.posY(), collision.posZ()}; + // std::array pVtx = {collision.posX(), collision.posY(), collision.posZ()}; for (auto& V0 : V0s) { std::array pos = {0.}; @@ -254,7 +254,7 @@ struct lambdakzerobuilder { MY_DEBUG_MSG(isK0SfromLc, LOG(info) << "in builder 0: posTrack --> " << labelPos << ", negTrack --> " << labelNeg); - //Apply selections so a skimmed table is created only + // Apply selections so a skimmed table is created only if (fitter.getChi2AtPCACandidate() > dcav0dau) { MY_DEBUG_MSG(isK0SfromLc, LOG(info) << "posTrack --> " << labelPos << ", negTrack --> " << labelNeg << " will be skipped due to dca cut"); continue; @@ -266,7 +266,7 @@ struct lambdakzerobuilder { continue; } - auto V0radius = RecoDecay::sqrtSumOfSquares(pos[0], pos[1]); //probably find better name to differentiate the cut from the variable + auto V0radius = RecoDecay::sqrtSumOfSquares(pos[0], pos[1]); // probably find better name to differentiate the cut from the variable if (V0radius < v0radius) { MY_DEBUG_MSG(isK0SfromLc, LOG(info) << "posTrack --> " << labelPos << ", negTrack --> " << labelNeg << " will be skipped due to radius cut"); continue; diff --git a/PWGLF/Tasks/NucleiSpectraTask.cxx b/PWGLF/Tasks/NucleiSpectraTask.cxx index 1977ef10f1d..7aa2e982e4b 100644 --- a/PWGLF/Tasks/NucleiSpectraTask.cxx +++ b/PWGLF/Tasks/NucleiSpectraTask.cxx @@ -142,9 +142,9 @@ struct NucleiSpectraTask { // spectra.fill(HIST("histPairCount"), negTracks.size(), posTracks.size()); // - for (Int_t iPos = 0; iPos < posTracks.size(); iPos++) { + for (unsigned int iPos = 0; iPos < posTracks.size(); iPos++) { TLorentzVector& vecPos = posTracks[iPos]; - for (Int_t jNeg = 0; jNeg < negTracks.size(); jNeg++) { + for (unsigned int jNeg = 0; jNeg < negTracks.size(); jNeg++) { TLorentzVector& vecNeg = negTracks[jNeg]; TLorentzVector vecMother = vecPos + vecNeg; spectra.fill(HIST("histInvMassData"), vecMother.Pt(), vecMother.M()); diff --git a/PWGUD/Tasks/upcForward.cxx b/PWGUD/Tasks/upcForward.cxx index 2de8d5647f2..99459e66bb9 100644 --- a/PWGUD/Tasks/upcForward.cxx +++ b/PWGUD/Tasks/upcForward.cxx @@ -75,7 +75,7 @@ struct UPCForward { // V0 and FD information bool isBeamBeamV0A = bc.bbV0A(); bool isBeamGasV0A = bc.bgV0A(); - bool isBeamBeamV0C = bc.bbV0C(); + // bool isBeamBeamV0C = bc.bbV0C(); bool isBeamGasV0C = bc.bgV0C(); bool isBeamBeamFDA = bc.bbFDA();