2 more Double32_t stored in the AOD header. These are the Qx and Qy components of

the TPC event plane vector. With this we are able to do all the necessary correction on the TPC event
plane calculation in nanoAODs (where the huge event plane object is not stored).
Julian Book <jbook@ikf.uni-frankfurt.de>

diff --git a/STEER/AOD/AliAODHeader.cxx b/STEER/AOD/AliAODHeader.cxx
index 653e4a3b92c85123682aa09489d0159125b6088f..1c26d0a6cdd09b2379bef68aafbe76f4b4b003ab 100644 (file)
--- a/STEER/AOD/AliAODHeader.cxx
+++ b/STEER/AOD/AliAODHeader.cxx
@@ -36,6 +36,8 @@ AliAODHeader::AliAODHeader() :
   fCentrality(-999.),
   fEventplane(-999.),
   fEventplaneMag(-999.),
+  fEventplaneQx(-999.),
+  fEventplaneQy(-999.),
   fZDCN1Energy(-999.),
   fZDCP1Energy(-999.),
   fZDCN2Energy(-999.),
@@ -97,6 +99,8 @@ AliAODHeader::AliAODHeader(Int_t nRun,
   fCentrality(-999.),
   fEventplane(-999.),
   fEventplaneMag(-999.),
+  fEventplaneQx(-999.),
+  fEventplaneQy(-999.),
   fZDCN1Energy(-999.),
   fZDCP1Energy(-999.),
   fZDCN2Energy(-999.),
@@ -178,6 +182,8 @@ AliAODHeader::AliAODHeader(Int_t nRun,
   fCentrality(cent),
   fEventplane(eventplane),
   fEventplaneMag(0),
+  fEventplaneQx(0),
+  fEventplaneQy(0),
   fZDCN1Energy(n1Energy),
   fZDCP1Energy(p1Energy),
   fZDCN2Energy(n2Energy),
@@ -245,6 +251,8 @@ AliAODHeader::AliAODHeader(const AliAODHeader& hdr) :
   fCentrality(hdr.fCentrality),
   fEventplane(hdr.fEventplane),
   fEventplaneMag(hdr.fEventplaneMag),
+  fEventplaneQx(hdr.fEventplaneQx),
+  fEventplaneQy(hdr.fEventplaneQy),
   fZDCN1Energy(hdr.fZDCN1Energy),
   fZDCP1Energy(hdr.fZDCP1Energy),
   fZDCN2Energy(hdr.fZDCN2Energy),
@@ -322,6 +330,8 @@ AliAODHeader& AliAODHeader::operator=(const AliAODHeader& hdr)
     fCentrality       = hdr.fCentrality;
     fEventplane       = hdr.fEventplane;
     fEventplaneMag    = hdr.fEventplaneMag;
+    fEventplaneQx     = hdr.fEventplaneQx;
+    fEventplaneQy     = hdr.fEventplaneQy;
     fZDCN1Energy      = hdr.fZDCN1Energy;
     fZDCP1Energy      = hdr.fZDCP1Energy;
     fZDCN2Energy      = hdr.fZDCN2Energy;
@@ -453,6 +463,8 @@ void AliAODHeader::Clear(Option_t* /*opt*/)
     fEventplaneP = 0;
     fEventplane = -999;
     fEventplaneMag = -999.;
+    fEventplaneQx = -999.;
+    fEventplaneQy = -999.;
   }
   return;
 }
@@ -473,8 +485,10 @@ void AliAODHeader::Print(Option_t* /*option*/) const
   printf("Muon mag. field scale   : %f\n", fMuonMagFieldScale);
   
   printf("Centrality              : %f\n", fCentrality);
-  printf("Event plane             : %f\n", fEventplane);
-  printf("Event plane             : %f\n", fEventplaneMag);
+  printf("Event plane Ang         : %f\n", fEventplane);
+  printf("Event plane Mag         : %f\n", fEventplaneMag);
+  printf("Event plane Qx          : %f\n", fEventplaneQx);
+  printf("Event plane Qy          : %f\n", fEventplaneQy);
   printf("ZDC N1 Energy           : %f\n", fZDCN1Energy);
   printf("ZDC P1 Energy           : %f\n", fZDCP1Energy);
   printf("ZDC N2 Energy           : %f\n", fZDCN2Energy);

diff --git a/STEER/AOD/AliAODHeader.h b/STEER/AOD/AliAODHeader.h
index 1e7623ec0dd1dbd20d3632b5ed8588f0edc0af6a..016b121b03749c89d7d4c4a07ccf482e345dc469 100644 (file)
--- a/STEER/AOD/AliAODHeader.h
+++ b/STEER/AOD/AliAODHeader.h
@@ -74,6 +74,8 @@ class AliAODHeader : public AliVHeader {
   Double_t  GetCentrality()         const { return fCentrality; }
   Double_t  GetEventplane()         const { return fEventplane; }
   Double_t  GetEventplaneMag()      const { return fEventplaneMag; }
+  Double_t  GetEventplaneQx()       const { return fEventplaneQx; }
+  Double_t  GetEventplaneQy()       const { return fEventplaneQy; }
   Double_t  GetZDCN1Energy()        const { return fZDCN1Energy; }
   Double_t  GetZDCP1Energy()        const { return fZDCP1Energy; }
   Double_t  GetZDCN2Energy()        const { return fZDCN2Energy; }
@@ -204,6 +206,8 @@ class AliAODHeader : public AliVHeader {
   Double32_t  fCentrality;          // Centrality
   Double32_t  fEventplane;          // Event plane angle
   Double32_t  fEventplaneMag;       // Length of Q vector from TPC event plance
+  Double32_t  fEventplaneQx;        // Q vector component x from TPC event plance
+  Double32_t  fEventplaneQy;        // Q vector component y from TPC event plance
   Double32_t  fZDCN1Energy;         // reconstructed energy in the neutron1 ZDC
   Double32_t  fZDCP1Energy;         // reconstructed energy in the proton1 ZDC
   Double32_t  fZDCN2Energy;         // reconstructed energy in the neutron2 ZDC
@@ -267,11 +271,19 @@ void AliAODHeader::SetEventplane(AliEventplane* eventplane)      {
        fEventplane = eventplane->GetEventplane("Q");
         const TVector2* qvect=eventplane->GetQVector();
         fEventplaneMag = -999;
-        if (qvect) fEventplaneMag=qvect->Mod();
+       fEventplaneQx = -999;
+       fEventplaneQy = -999;
+        if (qvect) {
+         fEventplaneMag=qvect->Mod();
+         fEventplaneQx=qvect->X();
+         fEventplaneQy=qvect->Y();
+       }
     }
     else{
        fEventplane = -999;
         fEventplaneMag = -999;
+       fEventplaneQx = -999;
+       fEventplaneQy = -999;
     }
 }
 inline