fixing psi in MC from the header
authorpchrist <pchrist@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 4 Apr 2012 09:42:12 +0000 (09:42 +0000)
committerpchrist <pchrist@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 4 Apr 2012 09:42:12 +0000 (09:42 +0000)
PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBFPsi.cxx
PWGCF/EBYE/macros/configBalanceFunctionPsiAnalysis.C

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
index ed2b8d3..34d76e4 100644 (file)
@@ -13,9 +13,9 @@ AliBalancePsi *GetBalanceFunctionObject(const char* analysisLevel = "ESD",
 \r
   //Set all analyses separately\r
   //Rapidity\r
-  gBalance->SetInterval(AliBalance::kRapidity,-0.8,0.8,16,-0.8,0.8,15.);  \r
+  gBalance->SetInterval(AliBalance::kRapidity,-0.8,0.8,32,-1.6,1.6,15.);  \r
   //Eta\r
-  gBalance->SetInterval(AliBalance::kEta,-0.8,0.8,16,-0.8,0.8,15);\r
+  gBalance->SetInterval(AliBalance::kEta,-0.8,0.8,32,-1.6,1.6,15);\r
   //Qlong\r
   gBalance->SetInterval(AliBalance::kQlong,-1,1,200,0.0,4.0,15);\r
   //Qout\r