add DYNAMIC_COLUMN for T0AC

Author: AllaMaevskaya

Common/DataModel/T0infoTable.h

@@ -1,3 +1,5 @@
+
+
 // Copyright 2019-2020 CERN and copyright holders of ALICE O2.
 // See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
 // All rights not expressly granted are reserved.
@@ -12,19 +14,26 @@
 #define O2_T0_INFO_TABLE_H_
 
 #include "Framework/AnalysisDataModel.h"
+//#include "CommonConstants/PhysicsConstants.h"
 
 namespace o2::aod
 {
 namespace infoFT0
 {
-DECLARE_SOA_COLUMN(T0A, t0A, float);           //! Collision time A-side, corrected with primary vertex
-DECLARE_SOA_COLUMN(T0C, t0C, float);           //! Collision time C-side, corrected with primary vertex
-DECLARE_SOA_COLUMN(T0AC, t0AC, float);         //! Collision time (A+C)/2
-DECLARE_SOA_COLUMN(T0vertex, t0vertex, float); //! FT0 vertex
+DECLARE_SOA_COLUMN(T0A, t0A, float);
+//! Collision time A-side, corrected with primary vertex
+DECLARE_SOA_COLUMN(T0C, t0C, float);
+//! Collision time C-side, corrected with primary vertex
+DECLARE_SOA_COLUMN(T0vertex, t0vertex, float);
+//! FT0 vertex
+DECLARE_SOA_DYNAMIC_COLUMN(T0AC, t0AC,
+                           [](float t0C, float t0A) -> float { return 0.5 * (t0A + t0C); });
+//! Collision time (A+C)/2
 } // namespace infoFT0
 DECLARE_SOA_TABLE(T0info, "AOD", "T0info",
                   infoFT0::T0A, infoFT0::T0C,
-                  infoFT0::T0AC, infoFT0::T0vertex);
+                  infoFT0::T0AC<infoFT0::T0A, infoFT0::T0C>,
+                  infoFT0::T0vertex);
 using infosFT0 = T0info::iterator;
 } // namespace o2::aod
 

Common/TableProducer/FT0Table.cxx

@@ -40,8 +40,8 @@ struct FT0Table {
       int64_t foundFT0 = collision.foundFT0();
       float vertexPV = collision.posZ();
       float vertex_corr = vertexPV / o2::constants::physics::LightSpeedCm2NS;
-      float t0A, t0C, t0AC, vertexT0;
-      t0A = t0C = t0AC = vertexT0 = 32000;
+      float t0A, t0C, vertexT0;
+      t0A = t0C = vertexT0 = 32000;
       if (foundFT0 != -1) {
         auto ft0 = ft0s.iteratorAt(foundFT0);
         t0A = ft0.timeA();
@@ -54,7 +54,6 @@ struct FT0Table {
         LOGF(debug, " T0A = %f, T0C %f, vertex_corr %f, triggersignals %i", t0A, t0C, vertex_corr, triggersignals);
         if (vertexTrigger) {
           vertexT0 = o2::constants::physics::LightSpeedCm2NS * (t0C - t0A) / 2;
-          t0AC = (t0A + t0C) / 2;
         }
         if (ora) {
           t0A += vertex_corr;
@@ -63,8 +62,8 @@ struct FT0Table {
           t0C -= vertex_corr;
         }
       }
-      LOGF(debug, " T0 collision time = %f, T0A = %f, T0C = %f, T0 vertex = %f PV = %f", t0A, t0C, t0AC, vertexT0, vertexPV);
-      table(t0A, t0C, t0AC, vertexT0);
+      LOGF(debug, " T0 collision time T0A = %f, T0C = %f, T0 vertex = %f PV = %f", t0A, t0C, vertexT0, vertexPV);
+      table(t0A, t0C, vertexT0);
     }
   }
 };

Common/Tasks/ft0infoQaTask.cxx

@@ -44,7 +44,7 @@ struct FT0InfoQaTask {
   {
     int64_t foundFT0 = col.foundFT0();
     if (foundFT0 != -1) {
-      LOGF(info, "multV0A=%5.0f;  multT0A=%5.0f; PV=%f; T0vertex=%f;T0A=%f; T0C=%f; T0AC=%f ", col.multV0A(), col.multT0A(), col.posZ(), col.t0A(), col.t0C(), col.t0AC(), col.t0vertex());
+      LOGF(debug, "multV0A=%5.0f;  multT0A=%5.0f; PV=%f; T0vertex=%f;T0A=%f; T0C=%f; T0AC=%f ", col.multV0A(), col.multT0A(), col.posZ(), col.t0A(), col.t0C(), col.t0AC(), col.t0vertex());
       hT0A->Fill(col.t0A());
       hT0C->Fill(col.t0C());
       hT0AC->Fill(col.t0AC());