o add Delta vars for flow (Julian)

diff --git a/PWGDQ/dielectron/AliDielectronVarManager.cxx b/PWGDQ/dielectron/AliDielectronVarManager.cxx
index 43bbb4904d2698b9ed41527098c1d2ec93a8e0a2..9024d44a9aa116e3fa855a54603a4354e71a2e5d 100644 (file)
--- a/PWGDQ/dielectron/AliDielectronVarManager.cxx
+++ b/PWGDQ/dielectron/AliDielectronVarManager.cxx
@@ -169,6 +169,9 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   "v0ArpH2",
   "v0CrpH2",
   "v0ACrpH2",
+  "DeltaPhiv0ArpH2",        
+  "DeltaPhiv0CrpH2",        
+  "DeltaPhiv0ACrpH2",       
   "v0ATPCDiffH2",
   "v0CTPCDiffH2",
   "v0Av0CDiffH2",

diff --git a/PWGDQ/dielectron/AliDielectronVarManager.h b/PWGDQ/dielectron/AliDielectronVarManager.h
index 204b0ce88a24102a653ae146d5f02892b98ed9d8..d0a6cf8e7ae8911d4579a33699b439d71af0374e 100644 (file)
--- a/PWGDQ/dielectron/AliDielectronVarManager.h
+++ b/PWGDQ/dielectron/AliDielectronVarManager.h
@@ -221,6 +221,9 @@ public:
     kv0ArpH2,                  // VZERO-A reaction plane of the Q vector for 2nd harmonic
     kv0CrpH2,                  //         reaction plane
     kv0ACrpH2,                 // VZERO-AC reaction plane of the Q vector for 2nd harmonic
+    kDeltaPhiv0ArpH2,          // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-A
+    kDeltaPhiv0CrpH2,          // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-C
+    kDeltaPhiv0ACrpH2,         // Delta phi of the pair with respect to the 2nd order harmonic reaction plane from V0-AC
     kv0ATPCDiffH2,             // V0A-TPC reaction plane difference for 2nd harmonic
     kv0CTPCDiffH2,             // V0C-TPC reaction plane difference for 2nd harmonic
     kv0Av0CDiffH2,             // V0A-V0C reaction plane difference for 2nd harmonic
@@ -1093,6 +1096,26 @@ inline void AliDielectronVarManager::FillVarDielectronPair(const AliDielectronPa
   values[AliDielectronVarManager::kV0ACrpH2FlowV2] = TMath::Cos(2.0*delta);  // 2nd harmonic flow coefficient
   values[AliDielectronVarManager::kDeltaPhiV0ACrpH2] = delta;
 
+
+  // quantities using the values of  AliEPSelectionTask
+  values[AliDielectronVarManager::kDeltaPhiv0ArpH2] = values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0ArpH2];
+  values[AliDielectronVarManager::kDeltaPhiv0CrpH2] = values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0CrpH2];
+  values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] = values[AliDielectronVarManager::kPhi] - values[AliDielectronVarManager::kv0ACrpH2];
+  // keep the interval [-pi,+pi]                                                                                 
+  if ( values[AliDielectronVarManager::kDeltaPhiv0ArpH2] > TMath::Pi() ) 
+    values[AliDielectronVarManager::kDeltaPhiv0ArpH2] -= TMath::TwoPi(); 
+  if ( values[AliDielectronVarManager::kDeltaPhiv0CrpH2] > TMath::Pi() ) 
+    values[AliDielectronVarManager::kDeltaPhiv0CrpH2] -= TMath::TwoPi(); 
+  if ( values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] > TMath::Pi() ) 
+    values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] -= TMath::TwoPi(); 
+
+  if ( values[AliDielectronVarManager::kDeltaPhiv0ArpH2] < -1.*TMath::Pi() ) 
+    values[AliDielectronVarManager::kDeltaPhiv0ArpH2] += TMath::TwoPi(); 
+  if ( values[AliDielectronVarManager::kDeltaPhiv0CrpH2] < -1.*TMath::Pi() ) 
+    values[AliDielectronVarManager::kDeltaPhiv0CrpH2] += TMath::TwoPi(); 
+  if ( values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] < -1.*TMath::Pi() )
+    values[AliDielectronVarManager::kDeltaPhiv0ACrpH2] += TMath::TwoPi(); 
+
              
   AliDielectronMC *mc=AliDielectronMC::Instance();