changes from fzhou

diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskPi0V2.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskPi0V2.cxx
index 0a3053679cfcb97c610a45bc2af2a66e7d62f573..76bd5560a516cff247643b0133673fa619cecda0 100644 (file)
--- a/PWGGA/EMCALTasks/AliAnalysisTaskPi0V2.cxx
+++ b/PWGGA/EMCALTasks/AliAnalysisTaskPi0V2.cxx
@@ -405,7 +405,7 @@ void AliAnalysisTaskPi0V2::UserCreateOutputObjects()
     fOutput->SetOwner();  // IMPORTANT!
 
     hEvtCount = new TH1F("hEvtCount", " Event Plane", 10, 0.5, 10.5);
-    hEvtCount->GetXaxis()->SetBinLabel(1,"All");
+    hEvtCount->GetXaxis()->SetBinLabel(1,"SemiMB");
     hEvtCount->GetXaxis()->SetBinLabel(2,"vert");
     hEvtCount->GetXaxis()->SetBinLabel(3,"cent");
     hEvtCount->GetXaxis()->SetBinLabel(4,"EPtask");
@@ -505,11 +505,15 @@ void AliAnalysisTaskPi0V2::UserExec(Option_t *)
     // Main loop
     // Called for each event
         
-        
     // Create pointer to reconstructed event
    AliVEvent *event = InputEvent();
    if (!event) { Printf("ERROR: Could not retrieve event"); return; }
 
+  Bool_t isSelected =0;      
+  isSelected = (((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & (AliVEvent::kMB | AliVEvent::kSemiCentral));
+    if(!isSelected )
+        return; 
+
     // create pointer to event
     fESD = dynamic_cast<AliESDEvent*>(event);
     if (!fESD) {

diff --git a/PWGGA/EMCALTasks/macros/AddTaskEMCALPi0V2.C b/PWGGA/EMCALTasks/macros/AddTaskEMCALPi0V2.C
index 6c3c2a5210613d6a16492433b09790f636e02a3b..c8d08da5e922ee5847a08ef9314a14654637da2f 100644 (file)
--- a/PWGGA/EMCALTasks/macros/AddTaskEMCALPi0V2.C
+++ b/PWGGA/EMCALTasks/macros/AddTaskEMCALPi0V2.C
@@ -1,6 +1,6 @@
 // $Id: AddTaskEMCALPi0V2.C 56081 2012-05-01 08:57:08Z loizides $
 
-AliAnalysisTaskPi0V2 *AddTaskEMCALPi0V2()
+AliAnalysisTask *AddTaskEMCALPi0V2()
 {
   // Get the pointer to the existing analysis manager via the static access method.
   //==============================================================================
@@ -18,7 +18,6 @@ AliAnalysisTaskPi0V2 *AddTaskEMCALPi0V2()
   // Create the task and configure it.
   //===========================================================================
   AliAnalysisTaskPi0V2* ana = new  AliAnalysisTaskPi0V2("Pi0v2Task");
-  ana->SelectCollisionCandidates( AliVEvent::kAnyINT | AliVEvent::kCentral | AliVEvent::kSemiCentral);
   
   mgr->AddTask(ana);
   
@@ -27,7 +26,7 @@ AliAnalysisTaskPi0V2 *AddTaskEMCALPi0V2()
   //==============================================================================
   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("histosEMCALP0v2", 
                                                            TList::Class(),AliAnalysisManager::kOutputContainer,
-                                                           Form("%s", AliAnalysisManager::GetCommonFileName()));
+                                                           Form("%s_v2task", AliAnalysisManager::GetCommonFileName()));
   
   mgr->ConnectInput  (ana, 0, mgr->GetCommonInputContainer());
   mgr->ConnectOutput (ana, 1, coutput1 );

diff --git a/PWGGA/EMCALTasks/macros/AddTaskEMCALPi0V2hardCodeEP.C b/PWGGA/EMCALTasks/macros/AddTaskEMCALPi0V2hardCodeEP.C
new file mode 100644 (file)
index 0000000..1e2d76e
--- /dev/null
+++ b/PWGGA/EMCALTasks/macros/AddTaskEMCALPi0V2hardCodeEP.C
@@ -0,0 +1,53 @@
+AliAnalysisTask *AddTaskEMCALPi0V2hardCodeEP()
+{
+  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+  if (!mgr) {
+    Error("AddTaskEMCALPi0V2hardCodeEP", "No analysis manager found.");
+    return NULL;
+  }
+
+  if (!mgr->GetInputEventHandler()) {
+    ::Error("AddTaskEMCALPi0V2hardCodeEP", "This task requires an input event handler");
+    return NULL;
+  }
+  TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
+  if (type=="AOD"){
+    ::Error("AddTaskEMCALPi0V2hardCodeEP", "The tasks exits because AODs are in input");
+    return NULL;
+  }
+  
+  //Event plane task
+  AliEPSelectionTask *eventplaneTask = new AliEPSelectionTask("EventplaneSelection");
+  eventplaneTask->SetTrackType("TPC");
+  eventplaneTask->SetUsePtWeight();
+  eventplaneTask->SetUsePhiWeight();
+  eventplaneTask->SetSaveTrackContribution();
+  
+  AliESDtrackCuts* epTrackCuts = new AliESDtrackCuts("AliESDtrackCuts", "Standard");
+  epTrackCuts->SetPtRange(0.1, 4);
+  eventplaneTask->SetPersonalESDtrackCuts(epTrackCuts);
+
+  mgr->AddTask(eventplaneTask);
+
+  TString containerName3 = mgr->GetCommonFileName();
+  containerName3 += ":PWGGA_pi0v2CalEventPlane";
+  
+  AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
+  AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("EPStatTPC",TList::Class(), AliAnalysisManager::kOutputContainer,containerName3.Data());
+  mgr->ConnectInput(eventplaneTask, 0, mgr->GetCommonInputContainer());
+  mgr->ConnectOutput(eventplaneTask,1,coutput1);
+
+  //analysis task 
+
+ AliAnalysisTaskSE* taskMB = new  AliAnalysisTaskPi0V2("Pi0v2Task"); 
+  
+  TString containerName = mgr->GetCommonFileName();
+  containerName += ":PWGGA_pi0v2CalSemiCentral";
+  
+  AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
+  AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histv2task", TList::Class(),AliAnalysisManager::kOutputContainer, containerName.Data());
+  mgr->ConnectInput(taskMB, 0, cinput);
+  mgr->ConnectOutput(taskMB, 1, coutput2);
+  
+  return NULL;
+}