@@ -67,6 +67,7 @@ struct PhosPi0 {
6767 Configurable<float > nonlinB{" nonlinB" , 0 ., " nonlinsrity param B (a+b*exp(-e/c))" };
6868 Configurable<float > nonlinC{" nonlinC" , 1 ., " nonlinsrity param C (a+b*exp(-e/c))" };
6969 Configurable<int > tofEffParam{" tofEffParam" , 0 , " parameterization of TOF cut efficiency" };
70+ Configurable<float > timeOffset{" timeOffset" , 0 ., " time offset to compensate imperfection of time calibration" };
7071
7172 using SelCollisions = soa::Join<aod::Collisions, aod::EvSels>;
7273 using SelCollisionsMC = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels>;
@@ -261,16 +262,16 @@ struct PhosPi0 {
261262 hMiBoth = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiBoth" , " inv mass for centrality" ,
262263 HistType::kTH2F , {mggAxis, ptAxis}))
263264 .get ();
264- hMiOneAll = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiAll " , " inv mass for centrality" ,
265+ hMiOneAll = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiOneAll " , " inv mass for centrality" ,
265266 HistType::kTH2F , {mggAxis, ptAxis}))
266267 .get ();
267- hMiOneCPV = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiCPV " , " inv mass for centrality" ,
268+ hMiOneCPV = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiOneCPV " , " inv mass for centrality" ,
268269 HistType::kTH2F , {mggAxis, ptAxis}))
269270 .get ();
270- hMiOneDisp = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiDisp " , " inv mass for centrality" ,
271+ hMiOneDisp = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiOneDisp " , " inv mass for centrality" ,
271272 HistType::kTH2F , {mggAxis, ptAxis}))
272273 .get ();
273- hMiOneBoth = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiBoth " , " inv mass for centrality" ,
274+ hMiOneBoth = std::get<std::shared_ptr<TH2>>(mHistManager .add (" mggMiOneBoth " , " inv mass for centrality" ,
274275 HistType::kTH2F , {mggAxis, ptAxis}))
275276 .get ();
276277 if (isMC) {
@@ -549,28 +550,28 @@ struct PhosPi0 {
549550 }
550551 }
551552 // Test time eff
552- if (std::abs (ph1.time ) < 12.5e-9 ) { // strict cut on first photon
553- if (std::abs (ph2.time ) < 100 .e -9 ) {
553+ if (std::abs (ph1.time - timeOffset ) < 12.5e-9 ) { // strict cut on first photon
554+ if (std::abs (ph2.time - timeOffset ) < 100 .e -9 ) {
554555 hReTime100->Fill (m, ph2.pt ());
555- if (std::abs (ph2.time ) < 50 .e -9 ) {
556+ if (std::abs (ph2.time - timeOffset ) < 50 .e -9 ) {
556557 hReTime50->Fill (m, ph2.pt ());
557- if (std::abs (ph2.time ) < 30 .e -9 ) {
558+ if (std::abs (ph2.time - timeOffset ) < 30 .e -9 ) {
558559 hReTime30->Fill (m, ph2.pt ());
559- if (std::abs (ph2.time ) < 12.5e-9 ) {
560+ if (std::abs (ph2.time - timeOffset ) < 12.5e-9 ) {
560561 hReTime12->Fill (m, ph2.pt ());
561562 }
562563 }
563564 }
564565 }
565566 }
566- if (std::abs (ph2.time ) < 12.5e-9 ) { // strict cut on first photon
567- if (std::abs (ph1.time ) < 100 .e -9 ) {
567+ if (std::abs (ph2.time - timeOffset ) < 12.5e-9 ) { // strict cut on first photon
568+ if (std::abs (ph1.time - timeOffset ) < 100 .e -9 ) {
568569 hReTime100->Fill (m, ph1.pt ());
569- if (std::abs (ph1.time ) < 50 .e -9 ) {
570+ if (std::abs (ph1.time - timeOffset ) < 50 .e -9 ) {
570571 hReTime50->Fill (m, ph1.pt ());
571- if (std::abs (ph1.time ) < 30 .e -9 ) {
572+ if (std::abs (ph1.time - timeOffset ) < 30 .e -9 ) {
572573 hReTime30->Fill (m, ph1.pt ());
573- if (std::abs (ph1.time ) < 12.5e-9 ) {
574+ if (std::abs (ph1.time - timeOffset ) < 12.5e-9 ) {
574575 hReTime12->Fill (m, ph1.pt ());
575576 }
576577 }
0 commit comments