6161#include " CCDB/BasicCCDBManager.h"
6262
6363// from phi
64- #include " PWGLF/DataModel/EPCalibrationTables.h"
6564#include " Common/DataModel/PIDResponseITS.h"
6665
6766using namespace o2 ;
@@ -85,8 +84,8 @@ struct F0980pbpbanalysis {
8584 Configurable<float > cfgCutVertex{" cfgCutVertex" , 10.0 , " PV selection" };
8685 Configurable<bool > cfgQvecSel{" cfgQvecSel" , true , " Reject events when no QVector" };
8786 Configurable<bool > cfgOccupancySel{" cfgOccupancySel" , false , " Occupancy selection" };
88- Configurable<int > cfgMaxOccupancy{ " cfgMaxOccupancy " , 999999 , " maximum occupancy of tracks in neighbouring collisions in a given time range" };
89- Configurable<int > cfgMinOccupancy{ " cfgMinOccupancy " , -100 , " maximum occupancy of tracks in neighbouring collisions in a given time range" };
87+ Configurable<int > cfgOccupancyMax{ " cfgOccupancyMax " , 999999 , " maximum occupancy of tracks in neighbouring collisions in a given time range" };
88+ Configurable<int > cfgOccupancyMin{ " cfgOccupancyMin " , -100 , " minimum occupancy of tracks in neighbouring collisions in a given time range" };
9089 Configurable<bool > cfgNCollinTR{" cfgNCollinTR" , false , " Additional selection for the number of coll in time range" };
9190 Configurable<bool > cfgPVSel{" cfgPVSel" , false , " Additional PV selection flag for syst" };
9291 Configurable<float > cfgPV{" cfgPV" , 8.0 , " Additional PV selection range for syst" };
@@ -101,8 +100,8 @@ struct F0980pbpbanalysis {
101100 Configurable<int > cfgTPCcluster{" cfgTPCcluster" , 70 , " Number of TPC cluster" };
102101 Configurable<float > cfgRatioTPCRowsOverFindableCls{" cfgRatioTPCRowsOverFindableCls" , 0.8 , " TPC Crossed Rows to Findable Clusters" };
103102
104- Configurable<float > cfgMinRap{ " cfgMinRap " , -0.5 , " Minimum rapidity for pair" };
105- Configurable<float > cfgMaxRap{ " cfgMaxRap " , 0.5 , " Maximum rapidity for pair" };
103+ Configurable<float > cfgRapMin{ " cfgRapMin " , -0.5 , " Minimum rapidity for pair" };
104+ Configurable<float > cfgRapMax{ " cfgRapMax " , 0.5 , " Maximum rapidity for pair" };
106105
107106 Configurable<bool > cfgPrimaryTrack{" cfgPrimaryTrack" , true , " Primary track selection" };
108107 Configurable<bool > cfgGlobalWoDCATrack{" cfgGlobalWoDCATrack" , true , " Global track selection without DCA" };
@@ -115,24 +114,23 @@ struct F0980pbpbanalysis {
115114 Configurable<int > cfgSelectPID{" cfgSelectPID" , 0 , " PID selection type" };
116115 Configurable<int > cfgSelectPtl{" cfgSelectPtl" , 0 , " Particle selection type" };
117116
118- Configurable<int > cfgnMods{ " cfgnMods " , 1 , " The number of modulations of interest starting from 2" };
117+ Configurable<int > cfgNMods{ " cfgNMods " , 1 , " The number of modulations of interest starting from 2" };
119118 Configurable<int > cfgNQvec{" cfgNQvec" , 7 , " The number of total Qvectors for looping over the task" };
120119
121120 Configurable<std::string> cfgQvecDetName{" cfgQvecDetName" , " FT0C" , " The name of detector to be analyzed" };
122121 Configurable<std::string> cfgQvecRefAName{" cfgQvecRefAName" , " TPCpos" , " The name of detector for reference A" };
123122 Configurable<std::string> cfgQvecRefBName{" cfgQvecRefBName" , " TPCneg" , " The name of detector for reference B" };
124123
125- Configurable<bool > cfgRotBkg{ " cfgRotBkg " , true , " flag to construct rotational backgrounds" };
126- Configurable<int > cfgNRotBkg{ " cfgNRotBkg " , 10 , " the number of rotational backgrounds" };
124+ Configurable<bool > cfgRotBkgSel{ " cfgRotBkgSel " , true , " flag to construct rotational backgrounds" };
125+ Configurable<int > cfgRotBkgNum{ " cfgRotBkgNum " , 10 , " the number of rotational backgrounds" };
127126
128127 // for phi test
129128 Configurable<bool > cfgTPCFinableClsSel{" cfgTPCFinableClsSel" , true , " TPC Crossed Rows to Findable Clusters selection flag" };
130129 Configurable<bool > cfgITSClsSel{" cfgITSClsSel" , false , " ITS cluster selection flag" };
131130 Configurable<int > cfgITScluster{" cfgITScluster" , 0 , " Number of ITS cluster" };
132- Configurable<bool > cfgpTDepPID{" cfgpTDepPID" , false , " pT dependent PID" };
133- Configurable<bool > cfgBetaCutSel{" cfgBetaCutSel" , false , " TOF beta cut selection flag" };
134- Configurable<float > cfgCutTOFBeta{" cfgCutTOFBeta" , 0.0 , " cut TOF beta" };
135- Configurable<bool > isDeepAngle{" isDeepAngle" , true , " Deep Angle cut" };
131+ Configurable<bool > cfgTOFBetaSel{" cfgTOFBetaSel" , false , " TOF beta cut selection flag" };
132+ Configurable<float > cfgTOFBetaCut{" cfgTOFBetaCut" , 0.0 , " cut TOF beta" };
133+ Configurable<bool > cfgDeepAngleSel{" cfgDeepAngleSel" , true , " Deep Angle cut" };
136134 Configurable<double > cfgDeepAngle{" cfgDeepAngle" , 0.04 , " Deep Angle cut value" };
137135
138136 ConfigurableAxis massAxis{" massAxis" , {400 , 0.2 , 2.2 }, " Invariant mass axis" };
@@ -186,8 +184,7 @@ struct F0980pbpbanalysis {
186184 // Filter PIDcutFilter = nabs(aod::pidtpc::tpcNSigmaKa) < cMaxTPCnSigmaPion;
187185 // Filter PIDcutFilter = nabs(aod::pidTPCFullKa::tpcNSigmaKa) < cMaxTPCnSigmaPion;
188186
189- using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::TPCMults, aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0Cs, aod::CentFT0As, aod::Mults, aod::Qvectors, aod::EPCalibrationTables>>;
190- // aod::EPCalibrationTables 추가됨
187+ using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::TPCMults, aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0Cs, aod::CentFT0As, aod::Mults, aod::Qvectors>>;
191188 using TrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTPCFullKa, aod::pidTOFbeta>>;
192189 // aod::pidTOFbeta 추가됨
193190
@@ -239,7 +236,7 @@ struct F0980pbpbanalysis {
239236 if (cfgQvecSel && (collision.qvecAmp ()[detId] < 1e-4 || collision.qvecAmp ()[refAId] < 1e-4 || collision.qvecAmp ()[refBId] < 1e-4 )) {
240237 return 0 ;
241238 }
242- if (cfgOccupancySel && (collision.trackOccupancyInTimeRange () > cfgMaxOccupancy || collision.trackOccupancyInTimeRange () < cfgMinOccupancy )) {
239+ if (cfgOccupancySel && (collision.trackOccupancyInTimeRange () > cfgOccupancyMax || collision.trackOccupancyInTimeRange () < cfgOccupancyMin )) {
243240 return 0 ;
244241 }
245242 if (cfgNCollinTR && !collision.selection_bit (o2::aod::evsel::kNoCollInTimeRangeStandard )) {
@@ -349,7 +346,7 @@ struct F0980pbpbanalysis {
349346 p1 = track1.p ();
350347 p2 = track2.p ();
351348 angle = std::acos ((pt1 * pt2 + pz1 * pz2) / (p1 * p2));
352- if (isDeepAngle && angle < cfgDeepAngle) {
349+ if (cfgDeepAngleSel && angle < cfgDeepAngle) {
353350 return 0 ;
354351 }
355352 return 1 ;
@@ -434,7 +431,7 @@ struct F0980pbpbanalysis {
434431 pion2.SetXYZM (trk2.px (), trk2.py (), trk2.pz (), massPtl);
435432 reco = pion1 + pion2;
436433
437- if (reco.Rapidity () > cfgMaxRap || reco.Rapidity () < cfgMinRap ) {
434+ if (reco.Rapidity () > cfgRapMax || reco.Rapidity () < cfgRapMin ) {
438435 continue ;
439436 }
440437
@@ -448,8 +445,8 @@ struct F0980pbpbanalysis {
448445 histos.fill (HIST (" hInvMass_f0980_LSmm_EPA" ), reco.M (), reco.Pt (), centrality, relPhi);
449446 }
450447
451- if (cfgRotBkg && trk1.sign () * trk2.sign () < 0 ) {
452- for (int nr = 0 ; nr < cfgNRotBkg ; nr++) {
448+ if (cfgRotBkgSel && trk1.sign () * trk2.sign () < 0 ) {
449+ for (int nr = 0 ; nr < cfgRotBkgNum ; nr++) {
453450 auto randomPhi = rn->Uniform (o2::constants::math::PI * 5.0 / 6.0 , o2::constants::math::PI * 7.0 / 6.0 );
454451 randomPhi += pion2.Phi ();
455452 pion2Rot.SetXYZM (pion2.Pt () * std::cos (randomPhi), pion2.Pt () * std::sin (randomPhi), trk2.pz (), massPtl);
0 commit comments