@@ -80,6 +80,7 @@ struct qVectorsTable {
8080 Configurable<float > cfgMaxPtOnTPC{" cfgMaxPtOnTPC" , 5 ., " maximum transverse momentum selection for TPC tracks participating in Q-vector reconstruction" };
8181 Configurable<int > cfgCorrLevel{" cfgCorrLevel" , 4 , " calibration step: 0 = no corr, 1 = gain corr, 2 = rectr, 3 = twist, 4 = full" };
8282 Configurable<std::vector<int >> cfgnMods{" cfgnMods" , {2 , 3 }, " Modulation of interest" };
83+ Configurable<float > cfgMaxCentrality{" cfgMaxCentrality" , 100 .f , " max. centrality for Q vector calibration" };
8384
8485 Configurable<std::string> cfgGainEqPath{" cfgGainEqPath" , " Users/j/junlee/Qvector/GainEq" , " CCDB path for gain equalization constants" };
8586 Configurable<std::string> cfgQvecCalibPath{" cfgQvecCalibPath" , " Analysis/EventPlane/QVecCorrections" , " CCDB pasth for Q-vecteor calibration constants" };
@@ -522,14 +523,14 @@ struct qVectorsTable {
522523 coll.centFV0A ()};
523524 cent = centAllEstim[cfgCentEsti];
524525 bool IsCalibrated = true ;
525- if (cent < 0 . || cent > 90 . ) {
526+ if (cent < 0 . || cent > cfgMaxCentrality ) {
526527 cent = 110 .;
527528 IsCalibrated = false ;
528529 }
529530 for (std::size_t id = 0 ; id < cfgnMods->size (); id++) {
530531 int ind = cfgnMods->at (id);
531532 CalQvec (ind, coll, tracks, qvecRe, qvecIm, qvecAmp, TrkBPosLabel, TrkBNegLabel, TrkBTotLabel);
532- if (cent < 90 . ) {
533+ if (cent < cfgMaxCentrality ) {
533534 for (auto i{0u }; i < kBTot + 1 ; i++) {
534535 helperEP.DoRecenter (qvecRe[(kBTot + 1 ) * 4 * id + i * 4 + 1 ], qvecIm[(kBTot + 1 ) * 4 * id + i * 4 + 1 ],
535536 objQvec.at (id)->GetBinContent (static_cast <int >(cent) + 1 , 1 , i + 1 ), objQvec.at (id)->GetBinContent (static_cast <int >(cent) + 1 , 2 , i + 1 ));
0 commit comments