Change the task name from AliAnalysisTaskParticleCorrelation to AliAnalysisTaskCaloTr...

[u/mrichter/AliRoot.git] / ANALYSIS / AliEPSelectionTask.cxx

diff --git a/ANALYSIS/AliEPSelectionTask.cxx b/ANALYSIS/AliEPSelectionTask.cxx
index 905b746ef20c9ec7cd993c04dd93db113e80b1fe..2d6007b4fc398b7a52e98a8db0e2a51e0c87155f 100644 (file)
--- a/ANALYSIS/AliEPSelectionTask.cxx
+++ b/ANALYSIS/AliEPSelectionTask.cxx
@@ -316,48 +316,46 @@ void AliEPSelectionTask::UserExec(Option_t */*option*/)
         fRunNumber = aod->GetRunNumber();
         SetPhiDist();      
       }
-    }
 
-    if (fUseMCRP) {
-      AliAODMCHeader *headerH = dynamic_cast<AliAODMCHeader*>(aod->GetList()->FindObject(AliAODMCHeader::StdBranchName()));
-      if (headerH) fRP = headerH->GetReactionPlaneAngle();
-    }
+      if (fUseMCRP) {
+       AliAODMCHeader *headerH = dynamic_cast<AliAODMCHeader*>(aod->GetList()->FindObject(AliAODMCHeader::StdBranchName()));
+       if (headerH) fRP = headerH->GetReactionPlaneAngle();
+      }
   
-    if (aod){
       esdEP = aod->GetHeader()->GetEventplaneP();
       esdEP->Reset(); 
      
       Int_t maxID = 0;
       TObjArray* tracklist = GetAODTracksAndMaxID(aod,maxID);
        
-    if (fSaveTrackContribution) {
-      esdEP->GetQContributionXArray()->Set(maxID+1);
-      esdEP->GetQContributionYArray()->Set(maxID+1);
-      esdEP->GetQContributionXArraysub1()->Set(maxID+1);
-      esdEP->GetQContributionYArraysub1()->Set(maxID+1);
-      esdEP->GetQContributionXArraysub2()->Set(maxID+1);
-      esdEP->GetQContributionYArraysub2()->Set(maxID+1);
-    }
+      if (fSaveTrackContribution) {
+       esdEP->GetQContributionXArray()->Set(maxID+1);
+       esdEP->GetQContributionYArray()->Set(maxID+1);
+       esdEP->GetQContributionXArraysub1()->Set(maxID+1);
+       esdEP->GetQContributionYArraysub1()->Set(maxID+1);
+       esdEP->GetQContributionXArraysub2()->Set(maxID+1);
+       esdEP->GetQContributionYArraysub2()->Set(maxID+1);
+      }
        
-    const int NT = tracklist->GetEntries();
+      const int NT = tracklist->GetEntries();
       
-  if (NT>4){
-    fQVector = new TVector2(GetQ(esdEP,tracklist));
-    fEventplaneQ = fQVector->Phi()/2.; 
-    GetQsub(qq1, qq2, tracklist, esdEP);
-    fQsub1 = new TVector2(qq1);
-    fQsub2 = new TVector2(qq2);
-    fQsubRes = (fQsub1->Phi()/2. - fQsub2->Phi()/2.);
+      if (NT>4){
+       fQVector = new TVector2(GetQ(esdEP,tracklist));
+       fEventplaneQ = fQVector->Phi()/2.; 
+       GetQsub(qq1, qq2, tracklist, esdEP);
+       fQsub1 = new TVector2(qq1);
+       fQsub2 = new TVector2(qq2);
+       fQsubRes = (fQsub1->Phi()/2. - fQsub2->Phi()/2.);
        
-    esdEP->SetQVector(fQVector);
-    esdEP->SetEventplaneQ(fEventplaneQ);
-    esdEP->SetQsub(fQsub1,fQsub2);
-    esdEP->SetQsubRes(fQsubRes);
+       esdEP->SetQVector(fQVector);
+       esdEP->SetEventplaneQ(fEventplaneQ);
+       esdEP->SetQsub(fQsub1,fQsub2);
+       esdEP->SetQsubRes(fQsubRes);
+       
+       fHOutEventplaneQ->Fill(fEventplaneQ);
+       fHOutsub1sub2->Fill(fQsub1->Phi()/2.,fQsub2->Phi()/2.);
+       fHOutNTEPRes->Fill(NT,fQsubRes);
        
-    fHOutEventplaneQ->Fill(fEventplaneQ);
-    fHOutsub1sub2->Fill(fQsub1->Phi()/2.,fQsub2->Phi()/2.);
-    fHOutNTEPRes->Fill(NT,fQsubRes);
-
        if (fUseMCRP) fHOutDiff->Fill(fEventplaneQ, fRP);
        
        for (int iter = 0; iter<NT;iter++){