123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115 |
- #ifndef MPD_FLOW_NAMESPACE_H
- #define MPD_FLOW_NAMESPACE_H
- #include <iostream>
- #include <vector>
- namespace mfdca
- {
- const float ptBins[] = {0., 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.5, 3.};
- const int NptBins = 12;
- const float etaBins[] = {-1.5, -1.2, -1., -0.8, -0.6, -0.4, -0.2, 0., 0.2, 0.4, 0.6, 0.8, 1., 1.2, 1.5};
- const int NetaBins = 14;
- const int Ndim = 3;
- const int NdcaFitIter = 3;
- const float NsigmaDcaCut = 2.;
- int GetPtBin(float pt);
- int GetEtaBin(float eta);
- }; // namespace mfdca
- namespace mfcent
- {
- const float ptRefMultCut = 0.15;
- const float etaRefMultCut = 0.5;
- float CentB(float bimp);
- // bool isRefMult(float pt, float eta);
- }; // namespace mfcent
- namespace mfEp
- {
- const int NetaGaps = 3;
- enum EpNames
- {
- kTpc_R_EtaGap1,
- kTpc_R_EtaGap2,
- kTpc_R_EtaGap3,
- kTpc_L_EtaGap1,
- kTpc_L_EtaGap2,
- kTpc_L_EtaGap3,
- kSpTpc_R_EtaGap1,
- kSpTpc_R_EtaGap2,
- kSpTpc_R_EtaGap3,
- kSpTpc_L_EtaGap1,
- kSpTpc_L_EtaGap2,
- kSpTpc_L_EtaGap3,
- kFHCal_R,
- kFHCal_L,
- kFHCal,
- };
- const int NsubEvents = kFHCal + 1;
- const int NsubEventsTpc = kSpTpc_L_EtaGap3 + 1;
- const float TpcEtaGap[NetaGaps] = {0.05, 0.1, 0.25};
- const int TpcREp[NetaGaps] = {kTpc_R_EtaGap1, kTpc_R_EtaGap2, kTpc_R_EtaGap3};
- const int TpcLEp[NetaGaps] = {kTpc_L_EtaGap1, kTpc_L_EtaGap2, kTpc_L_EtaGap3};
- const int SpTpcREp[NetaGaps] = {kSpTpc_R_EtaGap1, kSpTpc_R_EtaGap2, kSpTpc_R_EtaGap3};
- const int SpTpcLEp[NetaGaps] = {kSpTpc_L_EtaGap1, kSpTpc_L_EtaGap2, kSpTpc_L_EtaGap3};
- const int Nmodules = 45;
- double GetFHCalPhi(int iModule);
- std::pair<double, double> GetFHCalCoord(int iModule);
- float RangeAngle_02Pi(float phi, float harm = 1.); // make it from 0 to 2*Pi
- float RangeAngle(float phi, float harm = 1.); // make it from -Pi to Pi
- const int MinTpcEpMult = 4;
- const float MinFHCalEpEnergy = 0.;
- const std::vector<int> BadFHCalModules = {22, 67};
- const int Nflattening = 12;
- }; // namespace mfEp
- namespace mfRes
- {
- const int NresBin = 10;
- const int NresChiIter = 50;
- double ResEventPlane(double chi, double harmonic);
- double ResChi(double res, double harmonic, int Niter);
- }; // namespace mfRes
- namespace mfTrack
- {
- const float MaxTrackPt = 3.;
- const int MinTrackNhits = 32;
- const float MaxTrackEta = 1.5;
- const int motherId = -1;
- }; // namespace mfTrack
- namespace mfPid
- {
- const float MinPidProbability = 0.9;
- const int Npid = 8;
- const float Mass2[] = {0.01948, 0.2437, 0.88035};
- const int pdg[] = {211, 321, 2212, -211, -321, -2212};
- const int BadTofFlags[] = {0, 4};
- }; // namespace mfPid
- namespace mfFlow
- {
- const float MinV1PtEta = 0.2;
- }; // namespace mfFlow
- namespace mfAcceptance
- {
- // Bad TPC sector (in 0 to 2Pi range)
- const float AcceptanceTpcPhiMin = 0.26179939; // 15 deg
- const float AcceptanceTpcPhiMax = 0.78539816; // 45 deg
- // Bad FHCal modules (not 22, 67)
- const std::vector<int> AcceptanceFHCalModules = {15, 73};
- }
- #endif
|