fixing psi in MC from the header
[u/mrichter/AliRoot.git] / PWGCF / EBYE / BalanceFunctions / AliAnalysisTaskBFPsi.cxx
index 030e735..63aa6a4 100755 (executable)
@@ -235,7 +235,7 @@ void AliAnalysisTaskBFPsi::UserCreateOutputObjects() {
   fList->Add(fHistVz);\r
 \r
   //Event plane\r
-  fHistEventPlane = new TH1F("fHistEventPlane",";#Psi_{2} [rad];Counts",100,0,2.*TMath::Pi());\r
+  fHistEventPlane = new TH1F("fHistEventPlane",";#Psi_{2} [deg.];Counts",100,0,360.);\r
   fList->Add(fHistEventPlane);\r
 \r
   // QA histograms\r
@@ -444,7 +444,7 @@ void AliAnalysisTaskBFPsi::UserExec(Option_t *) {
                    gVZEROEventPlane = ep->CalculateVZEROEventPlane(gESD,10,2,qxTot,qyTot);\r
                  if(gVZEROEventPlane < 0.) gVZEROEventPlane += TMath::Pi();\r
                  fHistEventPlane->Fill(gVZEROEventPlane);\r
-                 gReactionPlane = gVZEROEventPlane;\r
+                 gReactionPlane = gVZEROEventPlane*TMath::RadToDeg();\r
                  //========Get the VZERO event plane========//\r
 \r
                  //Printf("There are %d tracks in this event", gESD->GetNumberOfTracks());\r
@@ -711,7 +711,7 @@ void AliAnalysisTaskBFPsi::UserExec(Option_t *) {
                    gVZEROEventPlane = ep->CalculateVZEROEventPlane(gAOD,10,2,qxTot,qyTot);\r
                  if(gVZEROEventPlane < 0.) gVZEROEventPlane += TMath::Pi();\r
                  fHistEventPlane->Fill(gVZEROEventPlane);\r
-                 gReactionPlane = gVZEROEventPlane;\r
+                 gReactionPlane = gVZEROEventPlane*TMath::RadToDeg();\r
                  //========Get the VZERO event plane========//\r
 \r
                  //Printf("There are %d tracks in this event", gAOD->GetNumberOfTracks());\r
@@ -865,7 +865,7 @@ void AliAnalysisTaskBFPsi::UserExec(Option_t *) {
                    gVZEROEventPlane = ep->CalculateVZEROEventPlane(gESD,10,2,qxTot,qyTot);\r
                  if(gVZEROEventPlane < 0.) gVZEROEventPlane += TMath::Pi();\r
                  fHistEventPlane->Fill(gVZEROEventPlane);\r
-                 gReactionPlane = gVZEROEventPlane;\r
+                 gReactionPlane = gVZEROEventPlane*TMath::RadToDeg();\r
                  //========Get the VZERO event plane========//\r
 \r
                  //Printf("There are %d tracks in this event", gESD->GetNumberOfTracks());\r
@@ -984,7 +984,7 @@ void AliAnalysisTaskBFPsi::UserExec(Option_t *) {
        fCentrality = gImpactParameter;\r
       }\r
       fCentrality = gImpactParameter;\r
-      fHistEventPlane->Fill(gReactionPlane);\r
+      fHistEventPlane->Fill(gReactionPlane*TMath::RadToDeg());\r
 \r
       // take only events inside centrality class (DIDN'T CHANGE THIS UP TO NOW)\r
       if((fImpactParameterMin > gImpactParameter) || (fImpactParameterMax < gImpactParameter))\r
@@ -1147,9 +1147,9 @@ void AliAnalysisTaskBFPsi::UserExec(Option_t *) {
   \r
   // calculate balance function\r
   if(fUseMultiplicity) \r
-    fBalance->CalculateBalance(gNumberOfAcceptedTracks,gReactionPlane,chargeVector);\r
+    fBalance->CalculateBalance(gNumberOfAcceptedTracks,gReactionPlane*TMath::RadToDeg(),chargeVector);\r
   else                 \r
-    fBalance->CalculateBalance(fCentrality,gReactionPlane,chargeVector);\r
+    fBalance->CalculateBalance(fCentrality,gReactionPlane*TMath::RadToDeg(),chargeVector);\r
 \r
   if(fRunShuffling) {\r
     // shuffle charges\r