K0s Analysis update (Matt Steinpries)
authordgangadh <dgangadh@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 31 Oct 2013 09:10:26 +0000 (09:10 +0000)
committerdgangadh <dgangadh@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 31 Oct 2013 09:10:26 +0000 (09:10 +0000)
PWGCF/FEMTOSCOPY/K0Analysis/AliFemtoK0Analysis.cxx
PWGCF/FEMTOSCOPY/K0Analysis/AliFemtoK0Analysis.h

index bdd94b5..f377649 100644 (file)
@@ -45,6 +45,7 @@
 //     - added merit cut choice, pass as argument (10/16/13)\r
 //             - 1-mass, 2-v0dca, 3-dddca, 4-combination (used to be v0dca)\r
 //     - added passable argument for two-track minimum separation (10/16/13)\r
+//     - added boolean to turn off field-sign dependence for train (10/30/13)\r
 ////////////////////////////////////////////////////////////////////////////////\r
 \r
 \r
@@ -89,6 +90,7 @@ ClassImp(AliFemtoK0Analysis)
 //________________________________________________________________________\r
 AliFemtoK0Analysis::AliFemtoK0Analysis():\r
 AliAnalysisTaskSE(),\r
+  fSignDep(kFALSE),\r
   fFieldPos(kTRUE),\r
   fOnlineCase(kTRUE),\r
   fMeritCase(kTRUE),\r
@@ -106,8 +108,9 @@ AliAnalysisTaskSE(),
 {\r
 }\r
 //________________________________________________________________________\r
-AliFemtoK0Analysis::AliFemtoK0Analysis(const char *name, bool FieldPositive, bool OnlineCase, bool MeritCase, float MinDL, int MeritCutChoice, float MinSep) \r
-: AliAnalysisTaskSE(name), \r
+AliFemtoK0Analysis::AliFemtoK0Analysis(const char *name, bool SignDep, bool FieldPositive, bool OnlineCase, bool MeritCase, float MinDL, int MeritCutChoice, float MinSep) \r
+: AliAnalysisTaskSE(name),\r
+  fSignDep(SignDep),\r
   fFieldPos(FieldPositive),\r
   fOnlineCase(OnlineCase),\r
   fMeritCase(MeritCase),\r
@@ -124,6 +127,7 @@ AliFemtoK0Analysis::AliFemtoK0Analysis(const char *name, bool FieldPositive, boo
   fPidAOD(0x0)\r
 {\r
   //main constructor\r
+  fSignDep             = SignDep;\r
   fFieldPos    = FieldPositive;\r
   fOnlineCase  = OnlineCase;\r
   fMeritCase   = MeritCase;\r
@@ -139,6 +143,7 @@ AliFemtoK0Analysis::AliFemtoK0Analysis(const char *name, bool FieldPositive, boo
 //________________________________________________________________________\r
 AliFemtoK0Analysis::AliFemtoK0Analysis(const AliFemtoK0Analysis &obj)\r
 : AliAnalysisTaskSE(obj.fName),\r
+  fSignDep(obj.fSignDep),\r
   fFieldPos(obj.fFieldPos),\r
   fOnlineCase(obj.fOnlineCase),\r
   fMeritCase(obj.fMeritCase),\r
@@ -161,6 +166,7 @@ AliFemtoK0Analysis &AliFemtoK0Analysis::operator=(const AliFemtoK0Analysis &obj)
  //Assignment operator\r
  if (this == &obj) return *this;\r
  \r
+ fSignDep   = obj.fSignDep;\r
  fFieldPos     = obj.fFieldPos;\r
  fOnlineCase   = obj.fOnlineCase;\r
  fMeritCase    = obj.fMeritCase;\r
@@ -452,9 +458,12 @@ void AliFemtoK0Analysis::Exec(Option_t *)
 \r
   float bField=0;\r
   bField = fAOD->GetMagneticField();\r
-  if(fFieldPos && bField < 0) return;\r
-  if(!fFieldPos && bField > 0) return;\r
   if(bField == 0) return;\r
+  if(fSignDep){\r
+   if(fFieldPos && bField < 0) return;\r
+   if(!fFieldPos && bField > 0) return;\r
+  }\r
+\r
   \r
   int zBin=0;\r
   double zStep=2*10/double(kZVertexBins), zstart=-10.;\r
@@ -810,7 +819,7 @@ void AliFemtoK0Analysis::Exec(Option_t *)
   float pairPt, pairMt, pairKt;                        //pair momentum values\r
   float pairMInv, pairPDotQ;\r
   float qinv, q0, qx, qy, qz;                          //pair q values\r
-  float qLength, thetaSH, thetaSHCos, phiSH;            //Spherical Harmonics values\r
+  //float qLength, thetaSH, thetaSHCos, phiSH;            //Spherical Harmonics values\r
   float am12, epm, h1, p12, p112, ppx, ppy, ppz, ks;   //PRF\r
   //float qOutLCMS;\r
   float qOutPRF, qSide, qLong;                         //relative momentum in LCMS/PRF frame\r
@@ -938,10 +947,10 @@ void AliFemtoK0Analysis::Exec(Option_t *)
                *histname3D += ktBin;\r
 \r
         //Spherical harmonics\r
-        qLength = sqrt(qLong*qLong + qSide*qSide + qOutPRF*qOutPRF);\r
-        thetaSHCos = qLong/qLength;\r
-        thetaSH = acos(thetaSHCos);\r
-        phiSH = acos(qOutPRF/(qLength*sin(thetaSH)));\r
+        //qLength = sqrt(qLong*qLong + qSide*qSide + qOutPRF*qOutPRF);\r
+        //thetaSHCos = qLong/qLength;\r
+        //thetaSH = acos(thetaSHCos);\r
+        //phiSH = acos(qOutPRF/(qLength*sin(thetaSH)));\r
 \r
         //Finding average separation of daughters throughout TPC - two-track cut\r
         float posPositions1[9][3] = {{0}};\r
index 6f86550..1e127fb 100644 (file)
@@ -30,7 +30,7 @@ class AliESDpid;
 class AliFemtoK0Analysis : public AliAnalysisTaskSE {
  public:
   AliFemtoK0Analysis();
-  AliFemtoK0Analysis(const char *name, bool FieldPositive = kTRUE, bool OnlineCase = kTRUE, bool MeritCase = kTRUE, float MinDL = 0.0, int MeritCutChoice = 2, float MinSep = 5.0);
+  AliFemtoK0Analysis(const char *name, bool SignDep = kFALSE, bool FieldPositive = kTRUE, bool OnlineCase = kTRUE, bool MeritCase = kTRUE, float MinDL = 0.0, int MeritCutChoice = 2, float MinSep = 5.0);
   virtual ~AliFemtoK0Analysis();
   AliFemtoK0Analysis(const AliFemtoK0Analysis&);
   AliFemtoK0Analysis& operator=(const AliFemtoK0Analysis&);
@@ -56,6 +56,7 @@ class AliFemtoK0Analysis : public AliAnalysisTaskSE {
     nphibins    = 72
   };
 
+  bool fSignDep;
   bool fFieldPos;
   bool fOnlineCase;
   bool fMeritCase;