change library dependencies and paths to the new places of the analysis under PWGGA...
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 24 Jan 2012 10:57:30 +0000 (10:57 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 24 Jan 2012 10:57:30 +0000 (10:57 +0000)
26 files changed:
PWGGA/CaloTasks/macros/anaCaloFilter.C [moved from PWGGA/EMCALTasks/macros/anaCaloFilter.C with 95% similarity]
PWGGA/CaloTrackCorrelations/macros/ana.C
PWGGA/CaloTrackCorrelations/macros/anaM.C
PWGGA/EMCALTasks/macros/AddTaskEMCALClusterize.C
PWGGA/EMCALTasks/macros/anaEMCALCalib.C
PWGGA/EMCALTasks/macros/emcalReclusterize.C
PWGGA/GammaConv/macros/AddTaskCaloConv.C
PWGGA/GammaConv/macros/AddTaskGammaConvDalitz.C
PWGGA/GammaConv/macros/ConfigGammaConversion.C
PWGGA/GammaConv/macros/Cuts_Events_new.C
PWGGA/GammaConv/macros/ExtractSignal.C
PWGGA/GammaConv/macros/Extract_IntegratedPi0Yield.C
PWGGA/GammaConv/macros/Extract_Pi0_Characteristics.C
PWGGA/GammaConv/macros/MakeCutLog.C
PWGGA/GammaConv/macros/Photon_Characteristics_Events.C
PWGGA/GammaConv/macros/Photon_Resolution.C
PWGGA/GammaConv/macros/Plot_Mapping_Histos_Events.C
PWGGA/GammaConv/macros/PlottingGammaConversionAdditional.h
PWGGA/GammaConv/macros/PlottingGammaConversionHistos.h
PWGGA/GammaConv/macros/start_GammaConversionAnalysis.sh
PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C
PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskPHOSQA.C
PWGGA/PHOSTasks/CaloCellQA/macros/ana_example.C
PWGGA/PHOSTasks/PHOS_PbPbQA/macros/PHOSPbPb.C
PWGGA/PHOSTasks/omega3pi/AddTaskomega3pi.C
PWGGA/PHOSTasks/omega3pi/grid/AnaTaskOmega3pi.C

similarity index 95%
rename from PWGGA/EMCALTasks/macros/anaCaloFilter.C
rename to PWGGA/CaloTasks/macros/anaCaloFilter.C
index d76d550..9892f5b 100644 (file)
@@ -170,14 +170,8 @@ void  LoadLibraries(const anaModes mode) {
     gSystem->Load("libCORRFW.so");
     gSystem->Load("libPHOSUtils.so"); 
     gSystem->Load("libEMCALUtils.so");
-    gSystem->Load("libPWG4CaloCalib.so");
+    gSystem->Load("libPWGGACaloTasks.so");
        
-    /*
-      //     gSystem->Load("libPWG4omega3pi.so");
-      //     gSystem->Load("libCORRFW.so");
-      //     gSystem->Load("libPWGCFbase.so");
-      //     gSystem->Load("libPWGmuon.so");
- */
     //--------------------------------------------------------
     //If you want to use root and par files from aliroot
     //--------------------------------------------------------  
@@ -191,7 +185,7 @@ void  LoadLibraries(const anaModes mode) {
           SetupPar("EMCALUtils");
           //Create Geometry
           TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
-          SetupPar("PWG4CaloCalib");
+          SetupPar("PWGGACaloTasks");
 */
   }
 
@@ -209,9 +203,7 @@ void  LoadLibraries(const anaModes mode) {
     //    gProof->ClearPackage("ESD");
     //    gProof->ClearPackage("AOD");
     //    gProof->ClearPackage("ANALYSIS");   
-    //    gProof->ClearPackage("PWG4PartCorrBase");
-    //    gProof->ClearPackage("PWG4PartCorrDep");
-    
+     
     // Enable the STEERBase Package
     gProof->UploadPackage("STEERBase.par");
     gProof->EnablePackage("STEERBase");
@@ -227,11 +219,6 @@ void  LoadLibraries(const anaModes mode) {
     // Enable the PHOS geometry Package
     //gProof->UploadPackage("PHOSUtils.par");
     //gProof->EnablePackage("PHOSUtils");
-    // Enable PartCorr analysis
-    gProof->UploadPackage("PWG4PartCorrBase.par");
-    gProof->EnablePackage("PWG4PartCorrBase");
-    gProof->UploadPackage("PWG4PartCorrDep.par");
-    gProof->EnablePackage("PWG4PartCorrDep");    
     gProof->ShowEnabledPackages();
   }  
   
index db54c49..3832cf8 100644 (file)
@@ -214,7 +214,7 @@ void ana(Int_t mode=mGRID)
   if(kInputData=="ESD"){
     printf("* Configure photon conversion analysis in macro \n");
     TString arguments = "-run-on-train -use-own-xyz  -force-aod -mc-off ";
-    gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/ConfigGammaConversion.C");
+    gROOT->LoadMacro("$ALICE_ROOT/PWGGA/GammaConv/macros/ConfigGammaConversion.C");
     AliAnalysisTaskGammaConversion * taskGammaConversion = 
     ConfigGammaConversion(arguments,mgr->GetCommonInputContainer());
     taskGammaConversion->SelectCollisionCandidates();
@@ -230,8 +230,8 @@ void ana(Int_t mode=mGRID)
   
   Bool_t kPrint   = kFALSE;
   Bool_t deltaAOD = kFALSE;
-  gROOT->LoadMacro("AddTaskCaloTrackCorr.C");   // $ALICE_ROOT/PWG4/macros
-  gROOT->LoadMacro("AddTaskEMCALClusterize.C"); // $ALICE_ROOT/PWG4/CaloCalib/macros  
+  gROOT->LoadMacro("AddTaskCaloTrackCorr.C");   // $ALICE_ROOT/PWGGA/CaloTrackCorrelations/macros
+  gROOT->LoadMacro("AddTaskEMCALClusterize.C"); // $ALICE_ROOT/PWGGA/EMCALTasks/macros  
   
   
   // ------
@@ -351,7 +351,7 @@ void  LoadLibraries(Int_t mode)
   if (mode == mPROOF) {
     //TProof::Mgr("ccalpmaster")->SetROOTVersion("ALICE_v5-27-06b");
     gROOT->LoadMacro("/afs/in2p3.fr/group/alice/laf/EnableAliRootForLAF.C");
-    TProof* proof = EnableAliRootForLAF("ccaplmaster",nPROOFWorkers.Data(),ccin2p3UserName.Data(),alienUserName.Data(),"",kFALSE,kTRUE,kTRUE,"OADB:ANALYSIS:ANALYSISalice:AOD:ESD:CORRFW:STEERBase:EMCALUtils:PHOSUtils:PWG4PartCorrBase:PWG4PartCorrDep:PWG4CaloCalib");
+    TProof* proof = EnableAliRootForLAF("ccaplmaster",nPROOFWorkers.Data(),ccin2p3UserName.Data(),alienUserName.Data(),"",kFALSE,kTRUE,kTRUE,"OADB:ANALYSIS:ANALYSISalice:AOD:ESD:CORRFW:STEERBase:EMCALUtils:PHOSUtils:PWGCaloTrackCorrBase:PWGGACaloTrackCorrelations:PWGGAEMCALTasks");
     
     //  TProof* proof = TProof::Open("ccaplmaster",Form("workers=%s",nPROOFWorkers.Data()));
     
@@ -366,9 +366,9 @@ void  LoadLibraries(Int_t mode)
     //     //proof->UploadPackage("JETAN");
     //     proof->UploadPackage("PHOSUtils");
     //     proof->UploadPackage("EMCALUtils");
-    //     proof->UploadPackage("PWG4PartCorrBase");
-    //     proof->UploadPackage("PWG4PartCorrDep");
-    //     proof->UploadPackage("PWG4CaloCalib");
+    //     proof->UploadPackage("PWGCaloTrackBase");
+    //     proof->UploadPackage("PWGGACaloTrackCorrelations");
+    //     proof->UploadPackage("PWGGAEMCALTasks");
     
     //     proof->EnablePackage("STEERBase");
     //     proof->EnablePackage("ESD");
@@ -380,9 +380,9 @@ void  LoadLibraries(Int_t mode)
     //     //proof->EnablePackage("JETAN");
     //     proof->EnablePackage("PHOSUtils");
     //     proof->EnablePackage("EMCALUtils");
-    //     proof->EnablePackage("PWG4PartCorrBase");
-    //     proof->EnablePackage("PWG4PartCorrDep");
-    //     proof->EnablePackage("PWG4CaloCalib");
+    //     proof->EnablePackage("PWGCaloTrackBase");
+    //     proof->EnablePackage("PWGGACaloTrackCorrelations");
+    //     proof->EnablePackage("PWGGAEMCALTasks");
     return;
   }  
   
@@ -429,20 +429,20 @@ void  LoadLibraries(Int_t mode)
   
   gSystem->Load("libPHOSUtils");
   gSystem->Load("libEMCALUtils");
-  gSystem->Load("libPWG4PartCorrBase");
-  gSystem->Load("libPWG4PartCorrDep");
-  gSystem->Load("libPWG4CaloCalib");
-  //SetupPar("PWG4PartCorrBase");
-  //SetupPar("PWG4PartCorrDep");
-  //SetupPar("PWG4CaloCalib");
+  gSystem->Load("libPWGCaloTrackBase");
+  gSystem->Load("libPWGGACaloCorrelations");
+  gSystem->Load("libPWGGAEMCALTasks");
+  //SetupPar("PWGCaloTrackBase");
+  //SetupPar("PWGGACaloCorrelations");
+  //SetupPar("PWGGAEMCALTasks");
   
   //gSystem->Load("libJETAN");
   //gSystem->Load("FASTJETAN");
-  //gSystem->Load("PWG4JetTasks");
+  //gSystem->Load("PWGJE");
 
   gSystem->Load("libCORRFW.so");
-  gSystem->Load("libPWG4GammaConv.so"); 
-  //SetupPar("PWG4GammaConv"); 
+  gSystem->Load("libPWGGAGammaConv.so"); 
+  //SetupPar("PWGGAGammaConv"); 
   
   // needed for plugin?
   gSystem->AddIncludePath("-I$ALICE_ROOT");
index 2d4f685..f7dc8a9 100644 (file)
@@ -110,16 +110,16 @@ void anaM()
     //correlation analysis
     gROOT->LoadMacro("AddTaskCaloTrackCorrM.C");
     
-    AliAnalysisTaskParticleCorrelationM *taskEMCAL = AddTaskCaloTrackCorrM(kInputData,"EMCAL",kFALSE);
+    AliAnalysisTaskCaloTrackCorrelationM *taskEMCAL = AddTaskCaloTrackCorrM(kInputData,"EMCAL",kFALSE);
 
     mgr->AddTask(taskEMCAL);
     
-    AliAnalysisTaskParticleCorrelationM *taskPHOS  = AddTaskCaloTrackCorrM(kInputData,"PHOS", kFALSE);
+    AliAnalysisTaskCaloTrackCorrelationM *taskPHOS  = AddTaskCaloTrackCorrM(kInputData,"PHOS", kFALSE);
 
     mgr->AddTask(taskPHOS);
     
     //gROOT->LoadMacro("AddTaskChargeCorr.C");
-    AliAnalysisTaskParticleCorrelationM *taskCharge  = AddTaskCaloTrackCorrM(kInputData, "CTS",kFALSE);
+    AliAnalysisTaskCaloTrackCorrelationM *taskCharge  = AddTaskCaloTrackCorrM(kInputData, "CTS",kFALSE);
 //    if(!kMC)
 //      taskCharge->SelectCollisionCandidates();
     mgr->AddTask(taskCharge);
@@ -158,8 +158,8 @@ void  LoadLibraries() {
     SetupPar("PHOSUtils");
     SetupPar("EMCALUtils");
     
-    SetupPar("PWG4CaloTrackCorrBase");
-    SetupPar("PWG4CaloTrackCorrDep");
+    SetupPar("PWGCaloTrackCorrBase");
+    SetupPar("PWGGACaloTrackCorrelations");
   }
   else{
     //--------------------------------------------------------
@@ -173,8 +173,8 @@ void  LoadLibraries() {
     gSystem->Load("libANALYSISalice");
          gSystem->Load("libPHOSUtils");
          gSystem->Load("libEMCALUtils");
-    gSystem->Load("libPWG4CaloTrackCorrBase");
-    gSystem->Load("libPWG4CaloTrackCorrDep");
+    gSystem->Load("libPWGCaloTrackCorrBase");
+    gSystem->Load("libPWGGACaloTrackCorrelations");
   }
   
 }
index fd20af8..74f529c 100644 (file)
@@ -116,7 +116,7 @@ AliAnalysisTaskEMCALClusterize* AddTaskEMCALClusterize(
   TGeoHMatrix* matrix[10];
   AliEMCALRecoUtils * reco = clusterize->GetRecoUtils();
   
-  gROOT->LoadMacro("ConfigureEMCALRecoUtils.C"); // $ALICE_ROOT/PWG4/CaloCalib/macros
+  gROOT->LoadMacro("ConfigureEMCALRecoUtils.C"); // $ALICE_ROOT/PWGGA/EMCALTasks/macros
   ConfigureEMCALRecoUtils(
                           reco,
                           bMC, 
index 0796e53..33bc1ab 100644 (file)
@@ -1,7 +1,6 @@
 /* $Id:  $ */
 //--------------------------------------------------
-// Example macro to do analysis with the 
-// analysis classes in PWG4PartCorr
+// Example macro to do EMCAL calibration analysis
 // Can be executed with Root and AliRoot
 //
 // Pay attention to the options and definitions
@@ -42,11 +41,11 @@ void anaEMCALCalib(Int_t mode=mLocal)
   // change whatever you need for your analysis case
   // ------------------------------------------------------------------
   LoadLibraries(mode) ;
-  //gSystem->Unload("libPWG4CaloCalib.so");
+  //gSystem->Unload("libPWGGAEMCALTasks.so");
   //Try to set the new library
-  //gSystem->Load("./PWG4CaloCalib/libPWG4CaloCalib.so");
+  //gSystem->Load("./PWGGAEMCALTasks/libPWGGAEMCALTasks.so");
   //gSystem->ListLibraries();
-
+  
   //-------------------------------------------------------------------------------------------------
   //Create chain from ESD and from cross sections files, look below for options.
   //-------------------------------------------------------------------------------------------------
@@ -56,10 +55,10 @@ void anaEMCALCalib(Int_t mode=mLocal)
     cout<<"Wrong  data type "<<kInputData<<endl;
     break;
   }
-
+  
   TChain *chain       = new TChain(kTreeName) ;
   CreateChain(mode, chain);  
-
+  
   if(chain){
     AliLog::SetGlobalLogLevel(AliLog::kError);//Minimum prints on screen
     
@@ -67,24 +66,24 @@ void anaEMCALCalib(Int_t mode=mLocal)
     // Make the analysis manager
     //-------------------------------------
     AliAnalysisManager *mgr  = new AliAnalysisManager("Manager", "Manager");
-
+    
     //input
     if(kInputData == "ESD")
-      {
-       // ESD handler
-       AliESDInputHandler *esdHandler = new AliESDInputHandler();
-       mgr->SetInputEventHandler(esdHandler);
-       esdHandler->SetReadFriends(kFALSE);
-       cout<<"ESD handler "<<mgr->GetInputEventHandler()<<endl;
-      }
+    {
+      // ESD handler
+      AliESDInputHandler *esdHandler = new AliESDInputHandler();
+      mgr->SetInputEventHandler(esdHandler);
+      esdHandler->SetReadFriends(kFALSE);
+      cout<<"ESD handler "<<mgr->GetInputEventHandler()<<endl;
+    }
     if(kInputData == "AOD")
-      {
-       // AOD handler
-       AliAODInputHandler *aodHandler = new AliAODInputHandler();
-       mgr->SetInputEventHandler(aodHandler);
-       cout<<"AOD handler "<<mgr->GetInputEventHandler()<<endl;
-       
-      }
+    {
+      // AOD handler
+      AliAODInputHandler *aodHandler = new AliAODInputHandler();
+      mgr->SetInputEventHandler(aodHandler);
+      cout<<"AOD handler "<<mgr->GetInputEventHandler()<<endl;
+      
+    }
     
     // mgr->SetDebugLevel(1);
     
@@ -92,19 +91,20 @@ void anaEMCALCalib(Int_t mode=mLocal)
     //Define task, put here any other task that you want to use.
     //-------------------------------------------------------------------------
     
-    // ESD filter task
+/*    
     if(kInputData == "ESD"){
-
-      //gROOT->LoadMacro("AddTaskPhysicsSelection.C");
-      //AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
-
+      
+      gROOT->LoadMacro("AddTaskPhysicsSelection.C");
+      AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
+      
     }
+*/    
     
     // Create containers for input/output
     AliAnalysisDataContainer *cinput1  = mgr->GetCommonInputContainer();
-
+    
     AliAnalysisTaskEMCALPi0CalibSelection * pi0calib = new AliAnalysisTaskEMCALPi0CalibSelection ("EMCALPi0Calibration");
-     if(kInputData == "ESD") pi0calib->SelectCollisionCandidates(); 
+    if(kInputData == "ESD") pi0calib->SelectCollisionCandidates(); 
     //pi0calib->SetDebugLevel(10); 
     //pi0calib->UseFilteredEventAsInput();
     pi0calib->SetClusterMinEnergy(0.5);
@@ -114,74 +114,96 @@ void anaEMCALCalib(Int_t mode=mLocal)
     pi0calib->SetNCellsGroup(0);
     pi0calib->SwitchOnSameSM();
     //pi0calib->SwitchOnOldAODs();
-
+    
+/*
     TGeoHMatrix *matrix[4];
     
     double rotationMatrix[4][9] = {-0.014585, -0.999892, -0.002031, 0.999892, -0.014589,  0.001950, -0.001979, -0.002003,  0.999996,
-                                  -0.014585,  0.999892,  0.002031, 0.999892,  0.014589, -0.001950, -0.001979,  0.002003, -0.999996,
-                                  -0.345861, -0.938280, -0.003412, 0.938281, -0.345869,  0.001950, -0.003010, -0.002527,  0.999992,
-                                  -0.345861,  0.938280,  0.003412, 0.938281,  0.345869, -0.001950, -0.003010,  0.002527, -0.999992};
+      -0.014585,  0.999892,  0.002031, 0.999892,  0.014589, -0.001950, -0.001979,  0.002003, -0.999996,
+      -0.345861, -0.938280, -0.003412, 0.938281, -0.345869,  0.001950, -0.003010, -0.002527,  0.999992,
+      -0.345861,  0.938280,  0.003412, 0.938281,  0.345869, -0.001950, -0.003010,  0.002527, -0.999992};
     
     double translationMatrix[4][3] = {0.367264,    446.508738,  175.97185+0.3,
-                                     1.078181,    445.826258, -174.026758+0.3,
-                                     -153.843916, 418.304256,  175.956905+0.8,
-                                     -152.649580, 417.621779, -174.040392+0.8};
+      1.078181,    445.826258, -174.026758+0.3,
+      -153.843916, 418.304256,  175.956905+0.8,
+      -152.649580, 417.621779, -174.040392+0.8};
     for(int j=0; j<4; j++)
-      {
-       matrix[j] = new TGeoHMatrix();
-       matrix[j]->SetRotation(rotationMatrix[j]);
-       matrix[j]->SetTranslation(translationMatrix[j]);
-       matrix[j]->Print();
-       pi0calib->SetGeometryMatrixInSM(matrix[j],j);
-      }
+    {
+      matrix[j] = new TGeoHMatrix();
+      matrix[j]->SetRotation(rotationMatrix[j]);
+      matrix[j]->SetTranslation(translationMatrix[j]);
+      matrix[j]->Print();
+      pi0calib->SetGeometryMatrixInSM(matrix[j],j);
+    }
+    
     
-
     pi0calib->SwitchOnLoadOwnGeometryMatrices();
-
+*/
+    
     pi0calib->SwitchOnClusterCorrection();
     AliEMCALRecoUtils * reco = pi0calib->GetEMCALRecoUtils();
-        reco->SetParticleType(AliEMCALRecoUtils::kPhoton);
+    reco->SetParticleType(AliEMCALRecoUtils::kPhoton);
     reco->SetW0(4.5);
+    
+    //---------------------
+    // Geometry alignment
+    //---------------------
 
     //reco->SetPositionAlgorithm(AliEMCALRecoUtils::kUnchanged);
-
     reco->SetPositionAlgorithm(AliEMCALRecoUtils::kPosTowerGlobal);
+    
     //reco->SetMisalTransShift(0,1.134);   reco->SetMisalTransShift(1,8.2); reco->SetMisalTransShift(2,1.197);
     //reco->SetMisalTransShift(3,-3.093);  reco->SetMisalTransShift(4,6.82);reco->SetMisalTransShift(5,1.635);
-
+    
     //reco->SetPositionAlgorithm(AliEMCALRecoUtils::kPosTowerIndex);
     //reco->SetMisalTransShift(0,1.08);   reco->SetMisalTransShift(1,8.35); reco->SetMisalTransShift(2,1.12);
     //reco->SetMisalRotShift(3,-8.05);    reco->SetMisalRotShift(4,8.05);  
     //reco->SetMisalTransShift(3,-0.42);  reco->SetMisalTransShift(5,1.55);
-
+    
+    //---------------------
+    // Non linearity
+    //---------------------
+    
     reco->SetNonLinearityFunction(AliEMCALRecoUtils::kNoCorrection);
     //reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0GammaGamma);
     //reco->SetNonLinearityParam(0,1.04);     reco->SetNonLinearityParam(1,-0.1445);
     //reco->SetNonLinearityParam(2,1.046);    
+    
+    //     reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0GammaConversion);
+    //     reco->SetNonLinearityParam(0,1.033);     reco->SetNonLinearityParam(1,0.0566186);
+    //     reco->SetNonLinearityParam(2,0.982133);    
+    
+    
+    //      reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0MC);
+    //      reco->SetNonLinearityParam(0,1.001);     reco->SetNonLinearityParam(1,-0.01264);
+    //      reco->SetNonLinearityParam(2,-0.03632);    
+    //      reco->SetNonLinearityParam(3,0.1798);     reco->SetNonLinearityParam(4,-0.522);
+    
+    
+    //---------------------
+    // Recalibration
+    //---------------------
+
+    // 
+    //     reco->SwitchOnRecalibration();
+    //      TFile * f = new TFile("RecalibrationFactors.root","read");
+    //      TH2F * h0 = (TH2F*)f->Get("EMCALRecalFactors_SM0")->Clone();
+    //      TH2F * h1 = (TH2F*)f->Get("EMCALRecalFactors_SM1")->Clone();
+    //      TH2F * h2 = (TH2F*)f->Get("EMCALRecalFactors_SM2")->Clone();
+    //      TH2F * h3 = (TH2F*)f->Get("EMCALRecalFactors_SM3")->Clone();
+    
+    //      reco->SetEMCALChannelRecalibrationFactors(0,h0);
+    //      reco->SetEMCALChannelRecalibrationFactors(1,h1);
+    //      reco->SetEMCALChannelRecalibrationFactors(2,h2);
+    //      reco->SetEMCALChannelRecalibrationFactors(3,h3);
+    
+    //---------------------
+    // Bad channels removal
+    //---------------------
 
-//     reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0GammaConversion);
-//     reco->SetNonLinearityParam(0,1.033);     reco->SetNonLinearityParam(1,0.0566186);
-//     reco->SetNonLinearityParam(2,0.982133);    
-
-
-//      reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0MC);
-//      reco->SetNonLinearityParam(0,1.001);     reco->SetNonLinearityParam(1,-0.01264);
-//      reco->SetNonLinearityParam(2,-0.03632);    
-//      reco->SetNonLinearityParam(3,0.1798);     reco->SetNonLinearityParam(4,-0.522);
-
-//     reco->SwitchOnRecalibration();
-//      TFile * f = new TFile("RecalibrationFactors.root","read");
-//      TH2F * h0 = (TH2F*)f->Get("EMCALRecalFactors_SM0")->Clone();
-//      TH2F * h1 = (TH2F*)f->Get("EMCALRecalFactors_SM1")->Clone();
-//      TH2F * h2 = (TH2F*)f->Get("EMCALRecalFactors_SM2")->Clone();
-//      TH2F * h3 = (TH2F*)f->Get("EMCALRecalFactors_SM3")->Clone();
-
-//      reco->SetEMCALChannelRecalibrationFactors(0,h0);
-//      reco->SetEMCALChannelRecalibrationFactors(1,h1);
-//      reco->SetEMCALChannelRecalibrationFactors(2,h2);
-//      reco->SetEMCALChannelRecalibrationFactors(3,h3);
-
+/*    
     reco->SwitchOnBadChannelsRemoval();
+    
     // SM0
     reco->SetEMCALChannelStatus(0,3,13);  reco->SetEMCALChannelStatus(0,44,1); reco->SetEMCALChannelStatus(0,3,13); 
     reco->SetEMCALChannelStatus(0,20,7);  reco->SetEMCALChannelStatus(0,38,2);   
@@ -199,24 +221,25 @@ void anaEMCALCalib(Int_t mode=mLocal)
     reco->SetEMCALChannelStatus(2,19,22);
     //SM3
     reco->SetEMCALChannelStatus(3,4,7);
-
+*/
+    
     reco->SetNumberOfCellsFromEMCALBorder(1);
-
+    
     //reco->Print("");
-
+    
     pi0calib->PrintInfo();
     mgr->AddTask(pi0calib);
-
+    
     AliAnalysisDataContainer *coutput2 = 
     mgr->CreateContainer("pi0calib", TList::Class(), AliAnalysisManager::kOutputContainer, "pi0calib.root");
     
     AliAnalysisDataContainer *cout_cuts = mgr->CreateContainer("Cuts", TList::Class(), 
-                                                      AliAnalysisManager::kOutputContainer, "pi0calib.root");
+                                                               AliAnalysisManager::kOutputContainer, "pi0calib.root");
     
     mgr->ConnectInput  (pi0calib,     0, cinput1);
     mgr->ConnectOutput (pi0calib, 1, coutput2 );
     mgr->ConnectOutput (pi0calib, 2, cout_cuts);
-
+    
     //-----------------------
     // Run the analysis
     //-----------------------    
@@ -231,9 +254,9 @@ void anaEMCALCalib(Int_t mode=mLocal)
     mgr->InitAnalysis();
     mgr->PrintStatus();
     mgr->StartAnalysis(smode.Data(),chain);
-
-cout <<" Analysis ended sucessfully "<< endl ;
-
+    
+    cout <<" Analysis ended sucessfully "<< endl ;
+    
   }
   else cout << "Chain was not produced ! "<<endl;
   
@@ -257,7 +280,7 @@ void  LoadLibraries(const anaModes mode) {
     // If you want to use already compiled libraries 
     // in the aliroot distribution
     //--------------------------------------------------------
-
+    
     gSystem->Load("libSTEERBase.so");
     gSystem->Load("libESD.so");
     gSystem->Load("libAOD.so");
@@ -265,14 +288,9 @@ void  LoadLibraries(const anaModes mode) {
     gSystem->Load("libANALYSISalice.so");
     gSystem->Load("libANALYSISalice.so");
     //TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
-    gSystem->Load("libPWG4CaloCalib.so");
-    //SetupPar("PWG4CaloCalib");
-    /*
-      //     gSystem->Load("libPWG4omega3pi.so");
-      //     gSystem->Load("libCORRFW.so");
-      //     gSystem->Load("libPWGHFbase.so");
-      //     gSystem->Load("libPWGmuon.so");
- */
+    gSystem->Load("libPWGGAEMCALTasks.so");
+    //SetupPar("PWGGAEMCALTasks");
+    
     //--------------------------------------------------------
     //If you want to use root and par files from aliroot
     //--------------------------------------------------------  
@@ -286,10 +304,10 @@ void  LoadLibraries(const anaModes mode) {
           SetupPar("EMCALUtils");
           //Create Geometry
           TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
-          SetupPar("PWG4CaloCalib");
-*/
+          SetupPar("PWGGAEMCALTasks");
+     */
   }
-
+  
   //---------------------------------------------------------
   // <<<<<<<<<< PROOF mode >>>>>>>>>>>>
   //---------------------------------------------------------
@@ -304,8 +322,6 @@ void  LoadLibraries(const anaModes mode) {
     //    gProof->ClearPackage("ESD");
     //    gProof->ClearPackage("AOD");
     //    gProof->ClearPackage("ANALYSIS");   
-    //    gProof->ClearPackage("PWG4PartCorrBase");
-    //    gProof->ClearPackage("PWG4PartCorrDep");
     
     // Enable the STEERBase Package
     gProof->UploadPackage("STEERBase.par");
@@ -322,11 +338,6 @@ void  LoadLibraries(const anaModes mode) {
     // Enable the PHOS geometry Package
     //gProof->UploadPackage("PHOSUtils.par");
     //gProof->EnablePackage("PHOSUtils");
-    // Enable PartCorr analysis
-    gProof->UploadPackage("PWG4PartCorrBase.par");
-    gProof->EnablePackage("PWG4PartCorrBase");
-    gProof->UploadPackage("PWG4PartCorrDep.par");
-    gProof->EnablePackage("PWG4PartCorrDep");    
     gProof->ShowEnabledPackages();
   }  
   
@@ -337,7 +348,7 @@ void SetupPar(char* pararchivename)
   //Load par files, create analysis libraries
   //For testing, if par file already decompressed and modified
   //classes then do not decompress.
+  
   TString cdir(Form("%s", gSystem->WorkingDirectory() )) ; 
   TString parpar(Form("%s.par", pararchivename)) ; 
   if ( gSystem->AccessPathName(parpar.Data()) ) {
@@ -404,7 +415,7 @@ void CreateChain(const anaModes mode, TChain * chain){
     //If you want to add several ESD files sitting in a common directory INDIR
     //Specify as environmental variables the directory (INDIR), the number of files 
     //to analyze (NFILES) and the pattern name of the directories with files (PATTERN)
-
+    
     if(gSystem->Getenv("INDIR"))  
       kInDir = gSystem->Getenv("INDIR") ; 
     else cout<<"INDIR not set, use default: "<<kInDir<<endl;   
@@ -421,15 +432,15 @@ void CreateChain(const anaModes mode, TChain * chain){
     if ( kInDir  && kFile) {
       printf("Get %d files from directory %s\n",kFile,kInDir);
       if ( ! gSystem->cd(kInDir) ) {//check if ESDs directory exist
-       printf("%s does not exist\n", kInDir) ;
-       return ;
+        printf("%s does not exist\n", kInDir) ;
+        return ;
       }
-
-  
+      
+      
       cout<<"INDIR   : "<<kInDir<<endl;
       cout<<"NFILES  : "<<kFile<<endl;
       cout<<"PATTERN : " <<kPattern<<endl;
-
+      
       TString datafile="";
       if(kInputData == "ESD") datafile = "AliESDs.root" ;
       else if(kInputData == "AOD") datafile = "aod.root" ;
@@ -441,19 +452,19 @@ void CreateChain(const anaModes mode, TChain * chain){
       char file[120] ;
       
       for (event = 0 ; event < kFile ; event++) {
-       sprintf(file, "%s/%s%d/%s", kInDir,kPattern,event,datafile.Data()) ; 
-       TFile * fESD = 0 ; 
-       //Check if file exists and add it, if not skip it
-       if ( fESD = TFile::Open(file)) {
-         if ( fESD->Get(kTreeName) ) { 
-           printf("++++ Adding %s\n", file) ;
-           chain->AddFile(file);
-         }
-       }
-       else { 
-         printf("---- Skipping %s\n", file) ;
-         skipped++ ;
-       }
+        sprintf(file, "%s/%s%d/%s", kInDir,kPattern,event,datafile.Data()) ; 
+        TFile * fESD = 0 ; 
+        //Check if file exists and add it, if not skip it
+        if ( fESD = TFile::Open(file)) {
+          if ( fESD->Get(kTreeName) ) { 
+            printf("++++ Adding %s\n", file) ;
+            chain->AddFile(file);
+          }
+        }
+        else { 
+          printf("---- Skipping %s\n", file) ;
+          skipped++ ;
+        }
       }
       printf("number of entries # %lld, skipped %d\n", chain->GetEntries(), skipped*100) ;     
     }
@@ -471,7 +482,7 @@ void CreateChain(const anaModes mode, TChain * chain){
   else if(mode == mGRID){
     //Get colection file. It is specified by the environmental
     //variable XML
-
+    
     if(gSystem->Getenv("XML") )
       kXML = gSystem->Getenv("XML");
     else
@@ -482,12 +493,12 @@ void CreateChain(const anaModes mode, TChain * chain){
       return ;
     }
     else cout<<"XML file "<<kXML<<endl;
-
+    
     //Load necessary libraries and connect to the GRID
     gSystem->Load("libNetx.so") ; 
     gSystem->Load("libRAliEn.so"); 
     TGrid::Connect("alien://") ;
-
+    
     //Feed Grid with collection file
     //TGridCollection * collection =  (TGridCollection*)gROOT->ProcessLine(Form("TAlienCollection::Open(\"%s\", 0)", kXML));
     TGridCollection * collection = (TGridCollection*) TAlienCollection::Open(kXML);
@@ -496,7 +507,7 @@ void CreateChain(const anaModes mode, TChain * chain){
       return kFALSE ; 
     }
     TGridResult* result = collection->GetGridResult("",0 ,0);
-   
+    
     // Makes the ESD chain 
     printf("*** Getting the Chain       ***\n");
     for (Int_t index = 0; index < result->GetEntries(); index++) {
index c6dea75..b0cbcac 100644 (file)
@@ -180,13 +180,12 @@ void  LoadLibraries(const anaModes mode) {
     gSystem->Load("libEMCALUtils.so");  
 
     //SetupPar("EMCALUtils"); 
-    //SetupPar("PWG4CaloCalib"); 
+    //SetupPar("PWGGAEMCALTasks"); 
     
     //TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
-    gSystem->Load("libPWG4CaloCalib.so");
+    gSystem->Load("libPWGGAEMCALTasks.so");
    
     /*
-      //     gSystem->Load("libPWG4omega3pi.so");
       //     gSystem->Load("libCORRFW.so");
       //     gSystem->Load("libPWGHFbase.so");
       //     gSystem->Load("libPWGmuon.so");
@@ -204,7 +203,7 @@ void  LoadLibraries(const anaModes mode) {
           SetupPar("EMCALUtils");
           //Create Geometry
           TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
-          SetupPar("PWG4CaloCalib");
+          SetupPar("PWGGAEMCALTasks");
 */
   }
 
@@ -222,8 +221,6 @@ void  LoadLibraries(const anaModes mode) {
     //    gProof->ClearPackage("ESD");
     //    gProof->ClearPackage("AOD");
     //    gProof->ClearPackage("ANALYSIS");   
-    //    gProof->ClearPackage("PWG4PartCorrBase");
-    //    gProof->ClearPackage("PWG4PartCorrDep");
     
     // Enable the STEERBase Package
     gProof->UploadPackage("STEERBase.par");
@@ -240,11 +237,6 @@ void  LoadLibraries(const anaModes mode) {
     // Enable the PHOS geometry Package
     //gProof->UploadPackage("PHOSUtils.par");
     //gProof->EnablePackage("PHOSUtils");
-    // Enable PartCorr analysis
-    gProof->UploadPackage("PWG4PartCorrBase.par");
-    gProof->EnablePackage("PWG4PartCorrBase");
-    gProof->UploadPackage("PWG4PartCorrDep.par");
-    gProof->EnablePackage("PWG4PartCorrDep");    
     gProof->ShowEnabledPackages();
   }  
   
index f2ff41c..5826819 100644 (file)
@@ -55,9 +55,9 @@ AliAnalysisTaskCaloConv * AddTaskCaloConv(){
   mgr->ConnectInput(task, 0, cinput);
 
   TString outputfile = AliAnalysisManager::GetCommonFileName();
-  AliAnalysisDataContainer *coutput = mgr->CreateContainer("CaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWG4_CaloConv",outputfile.Data()));
+  AliAnalysisDataContainer *coutput = mgr->CreateContainer("CaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWGGA_CaloConv",outputfile.Data()));
   mgr->ConnectOutput(task, 1, coutput);
-  AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("CFCaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWG4_CFCaloConv",outputfile.Data()));
+  AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("CFCaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWGGA_CFCaloConv",outputfile.Data()));
   mgr->ConnectOutput(task, 2, coutput2);
 
   return task ;
index d996515..24b9122 100644 (file)
@@ -106,7 +106,7 @@ AliAnalysisTaskGammaConvDalitz* AddTaskGammaConvDalitz(
 
     // Connect I/O to the task
     TString outputfile = AliAnalysisManager::GetCommonFileName();
-    outputfile += ":PWG4_GammaConversion";
+    outputfile += ":PWGGA_GammaConversion";
     outputfile += Form("_%s",kGCAnalysisCutSelectionId.Data());
     AliAnalysisDataContainer *coutput3 = mgr->CreateContainer("histogramsAliGammaConversionDalitz", TList::Class(),AliAnalysisManager::kOutputContainer, outputfile);
     mgr->ConnectInput ( dalitzTask, 0, mgr->GetCommonInputContainer() );
index 436acd3..9336074 100644 (file)
@@ -915,8 +915,8 @@ TString kGCoutputFileName = "histogramsGammaConversion";
 TString kGCoutputFileAppendix = "";
 TString kGCdataList = "";
 Bool_t kGCwriteNtuple = kFALSE;
-// WE DOO NOT NEED TO CHANGE THIS (kGCusePWG4PartCorr) ANYMORE SINCE IT IS TAKEN CARE OF AUTOMATICALLY NOW
-Bool_t kGCusePWG4PartCorr = kTRUE;
+// WE DOO NOT NEED TO CHANGE THIS (kGCusePWGGACaloTrackCorrelations) ANYMORE SINCE IT IS TAKEN CARE OF AUTOMATICALLY NOW
+Bool_t kGCusePWGGACaloTrackCorrelations = kTRUE;
 
 /** Flag to enable running on train  */
 Bool_t kGCrunOnTrain = kFALSE;
@@ -1183,21 +1183,21 @@ Bool_t scanArguments(TString arguments){
 }
 
 void SetVersionLibrary(){
-  // Check if the file $ALICE_ROOT/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx exists.
-  // If yes, we set kGCusePWG4PartCorr to false since we have a newer version
-  // If no, kGCusePWG4PartCorr is true.
+  // Check if the file $ALICE_ROOT/PWGGA/GammaConv/AliAnalysisTaskGammaConversion.cxx exists.
+  // If yes, we set kGCusePWGGACaloTrackCorrelations to false since we have a newer version
+  // If no, kGCusePWGGACaloTrackCorrelations is true.
        
   TString file = gSystem->Getenv("ALICE_ROOT");
-  file+="/PWG4/PartCorr/AliAnalysisTaskGammaConversion.cxx";
+  file+="/PWGGA/PartCorr/AliAnalysisTaskGammaConversion.cxx";
        
   ifstream stream;
   stream.open(file.Data());
        
   if(!stream){
-    kGCusePWG4PartCorr=kFALSE;
+    kGCusePWGGACaloTrackCorrelations=kFALSE;
   }
   else{
-    kGCusePWG4PartCorr=kTRUE;
+    kGCusePWGGACaloTrackCorrelations=kTRUE;
   }
   stream.close();
 }
@@ -1212,7 +1212,7 @@ AliAnalysisTaskGammaConversion* ConfigGammaConversion(TString arguments, AliAnal
     break;
   }
        
-  SetVersionLibrary(); // checks if PWG4GammaConv or PWG4PartCorr is used
+  SetVersionLibrary(); // checks if PWGGAGammaConv or PWGGACaloTrackCorrelations is used
        
   if(cin_esd == NULL && kGCrunOnTrain == kTRUE){
     cout<<"Error: kGCrunOnTrain flag is set to true but the input AliAnalysisDataContainer is NULL"<<endl;
@@ -1379,7 +1379,7 @@ AliAnalysisTaskGammaConversion* ConfigGammaConversion(TString arguments, AliAnal
   // Declare Common Input Tchain
   AliAnalysisDataContainer *cinput1 = NULL;
   
-  if(kGCusePWG4PartCorr){
+  if(kGCusePWGGACaloTrackCorrelations){
   
     if(kGCrunOnTrain == kFALSE){
       cinput1 = mgr->CreateContainer("GammaConvChain",TChain::Class(),AliAnalysisManager::kInputContainer);
@@ -1394,7 +1394,7 @@ AliAnalysisTaskGammaConversion* ConfigGammaConversion(TString arguments, AliAnal
   }
        
   // Common Output Tree in common ??????default?????? output file
-  // CKB kGCusePWG4PartCorr and writestandard are not mutually exclusive?
+  // CKB kGCusePWGGACaloTrackCorrelations and writestandard are not mutually exclusive?
   AliAnalysisDataContainer *coutput1 = NULL;
 
   if(kGCWriteAOD) {
@@ -1417,8 +1417,8 @@ AliAnalysisTaskGammaConversion* ConfigGammaConversion(TString arguments, AliAnal
   
   TString outputfile = AliAnalysisManager::GetCommonFileName();
   cout<<"Analyis cut selection ID is: "<<kGCAnalysisCutSelectionId.Data()<<endl;
-  //  outputfile += Form(":PWG4_GammaConversion_%s",kGCAnalysisCutSelectionId.Data());
-  outputfile += Form(":PWG4_GammaConversion");
+  //  outputfile += Form(":PWGGA_GammaConversion_%s",kGCAnalysisCutSelectionId.Data());
+  outputfile += Form(":PWGGA_GammaConversion");
 
 //   if(kGCrunNeutralMeson==kTRUE) outputfile +="1";  else outputfile +="0";
 
@@ -1693,7 +1693,7 @@ AliAnalysisTaskGammaConversion* ConfigGammaConversion(TString arguments, AliAnal
 
   if( kGCrunDalitz ){
    
-   gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskGammaConvDalitz.C");
+   gROOT->LoadMacro("$ALICE_ROOT/PWGGA/macros/AddTaskGammaConvDalitz.C");
    AddTaskGammaConvDalitz( v0Reader, kGCcalculateBackground, kGCRunStandalone );
 
   }
@@ -1742,7 +1742,7 @@ void LoadLibraries() {
   gSystem->Load("libANALYSIS.so");
   gSystem->Load("libANALYSISalice.so");
   gSystem->Load("libCORRFW.so");
-  gSystem->Load("libPWG4GammaConv.so");
+  gSystem->Load("libPWGGAGammaConv.so");
   
   //  gSystem->ChangeDirectory(pwd.Data());
   
@@ -1802,11 +1802,11 @@ void build() {
   gSystem->Load("libCORRFW.so");
        
   ////
-  //Setting up PWG4GammaConv.par//
+  //Setting up PWGGAGammaConv.par//
   ////
-  cout<<"compiling PWG4GammaConv"<<endl;
-  setupPar("PWG4GammaConv");
-  gSystem->Load("libPWG4GammaConv.so");
+  cout<<"compiling PWGGAGammaConv"<<endl;
+  setupPar("PWGGAGammaConv");
+  gSystem->Load("libPWGGAGammaConv.so");
   
   gSystem->ChangeDirectory(pwd.Data());
 }
index 20e68b5..474917c 100644 (file)
@@ -1,5 +1,5 @@
 /****************************************************************************************************************************
-******         provided by Gamma Conversion Group, PWG4,                                                                                                               *****
+******         provided by Gamma Conversion Group, PWGGA,                                                                                                              *****
 ******         Friederike Bock, friederike.bock@cern.ch                                                                                                        *****
 *****************************************************************************************************************************
 *** This macro can be used to observe the effects of the different cuts used in the GammaConversion analysis, it can   *****
@@ -55,12 +55,12 @@ void  Cuts_Events_new(const char *data = "myOutput", const char *MCfile = "",con
        char *StandardYAxis = "#gamma/ event scaled by multiplicity";
 
        if(cutsel != ""){
-               char *GammaDirectory = Form("PWG4_GammaConversion_%s",  cutsel);
+               char *GammaDirectory = Form("PWGGA_GammaConversion_%s",  cutsel);
                cout << GammaDirectory << endl;
                char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
                cout << GammaList << endl;
        }else{
-               char *GammaDirectory = "PWG4_GammaConversion";
+               char *GammaDirectory = "PWGGA_GammaConversion";
                cout << GammaDirectory << endl;
                char *GammaList = "histogramsAliGammaConversion";
                cout << GammaList << endl;
@@ -76,13 +76,13 @@ void  Cuts_Events_new(const char *data = "myOutput", const char *MCfile = "",con
        Char_t filename_data[200] = (Form("%s%s",path,data));   
        TFile f(filename_data);  
 
-        TDirectory *fPWG4GammaConversion_data = new TDirectory(); // definition of first folder / list 
+        TDirectory *fPWGGAGammaConversion_data = new TDirectory(); // definition of first folder / list        
         TList *fHistosGammaConversion_data = new TList(); // definition of first folder / list
         TList *fESDContainer_data = new TList();  // definition of following folder / list
         TList *fMappingContainer_data = new TList();    
 
-        if(!(fPWG4GammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl; 
-        if(!(fHistosGammaConversion_data = (TList*)fPWG4GammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
+        if(!(fPWGGAGammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl; 
+        if(!(fHistosGammaConversion_data = (TList*)fPWGGAGammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
         if(!(fESDContainer_data = (TList*)fHistosGammaConversion_data->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl; 
        fMappingContainer_data = (TList*)fHistosGammaConversion_data->FindObject("Mapping histograms"); 
        if(MC1){
@@ -313,13 +313,13 @@ void  Cuts_Events_new(const char *data = "myOutput", const char *MCfile = "",con
                
                
                // for new versions
-                TDirectory *fPWG4GammaConversion_montecarlo = new TDirectory(); // definition of first folder / list   
+                TDirectory *fPWGGAGammaConversion_montecarlo = new TDirectory(); // definition of first folder / list  
                 TList *fHistosGammaConversion_montecarlo = new TList(); // definition of first folder / list
                 TList *fESDContainer_montecarlo = new TList();  // definition of following folder / list
                 TList *fMappingContainer_montecarlo = new TList();  // definition of following folder / list
 
-                if(!(fPWG4GammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl; 
-                if(!(fHistosGammaConversion_montecarlo = (TList*)fPWG4GammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
+                if(!(fPWGGAGammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl; 
+                if(!(fHistosGammaConversion_montecarlo = (TList*)fPWGGAGammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
                 if(!(fESDContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl; 
                fMappingContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("Mapping histograms");
 
index 52d8e23..6d12c43 100644 (file)
@@ -1,4 +1,4 @@
-// provided by Gamma Conversion Group, $ALICE_ROOT/PWG4/GammaConv ;https://twiki.cern.ch/twiki/bin/view/ALICE/PWG4GammaConversion
+// provided by Gamma Conversion Group, $ALICE_ROOT/PWGGA/GammaConv ;https://twiki.cern.ch/twiki/bin/view/ALICE/PWG4GammaConversion
 
 #include <stdlib.h>
 #include <iostream>
@@ -390,9 +390,9 @@ void ExtractSignal(TString meson, TString file, TString cutSelection, TString su
        TFile f(file.Data());
        fCutSelection = cutSelection;
 
-       TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWG4_GammaConversion_%s",fCutSelection.Data()));
+       TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWGGA_GammaConversion_%s",fCutSelection.Data()));
        if(pwg4dir == NULL){
-               cout<<Form("ERROR: PWG4_GammaConversion_%s",fCutSelection.Data())<<" is not found in the file"<<endl;
+               cout<<Form("ERROR: PWGGA_GammaConversion_%s",fCutSelection.Data())<<" is not found in the file"<<endl;
                return;
        }
        TList *fHistosGammaConversion = (TList*)pwg4dir->Get(Form("histogramsAliGammaConversion_%s",fCutSelection.Data()));
index a83ed41..d4ba887 100644 (file)
@@ -1,4 +1,4 @@
-// provided by Gamma Conversion Group, PWG4, Kathrin Koch, kkoch@physi.uni-heidelberg.de
+// provided by Gamma Conversion Group, PWGGA, Kathrin Koch, kkoch@physi.uni-heidelberg.de
 
 #include <fstream>
 #include <Riostream.h>
@@ -44,7 +44,7 @@ void Extract_IntegratedPi0Yield(const char* cutSelection="", const char *inputRo
   c_mass->SetFillColor(0);
 
       
-  TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWG4_GammaConversion_%s",cutSelection));
+  TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWGGA_GammaConversion_%s",cutSelection));
   TList *fHistosGammaConversion = (TList*)pwg4dir->Get(Form("histogramsAliGammaConversion_%s",cutSelection));
   TList *fESDContainer = (TList*)fHistosGammaConversion->FindObject("ESD histograms");
   TList *fMCContainer = (TList*)fHistosGammaConversion->FindObject("MC histograms");
index f1a0d6d..463892a 100644 (file)
@@ -43,7 +43,7 @@ void Extract_Pi0_Characteristics(const char *cutSelection="", const char *inputR
   if(makeMappingPlots==kTRUE){
     outputFile = new TFile(Form("%sPi0Characteristics.root",outputDir),"UPDATE");
   }
-  TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWG4_GammaConversion_%s",cutSelection));
+  TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWGGA_GammaConversion_%s",cutSelection));
   TList *fHistosGammaConversion = (TList*)pwg4dir->Get(Form("histogramsAliGammaConversion_%s",cutSelection));
   TList *fESDContainer = (TList*)fHistosGammaConversion->FindObject("ESD histograms");
   TList *fMCContainer = (TList*)fHistosGammaConversion->FindObject("MC histograms");
index 2849292..cf01921 100644 (file)
@@ -1,4 +1,4 @@
-// provided by Gamma Conversion Group, PWG4, Kathrin Koch, kkoch@physi.uni-heidelberg.de
+// provided by Gamma Conversion Group, PWGGA, Kathrin Koch, kkoch@physi.uni-heidelberg.de
 
 #include <Riostream>
 #include <fstream>
@@ -22,13 +22,13 @@ void MakeCutLog(const char *inputRootFile = "AnalysisResults",const char *path =
 
     TObject * o = f.Get(directories->At(entFile)->GetName()); // get the object in the base directory
 
-    if(TString(o->IsA()->GetName())=="TDirectoryFile"){ // means that this is a directory (PWG4......)
+    if(TString(o->IsA()->GetName())=="TDirectoryFile"){ // means that this is a directory (PWGGA......)
       
       TDirectory *pwg4dir =(TDirectory*)o;
  
       TString baseDirName = pwg4dir->GetName();
       
-      TString reconstructionFlagString = ""; // this is for new scheme where also the flags are coded in numbers in the PWG4.... name
+      TString reconstructionFlagString = ""; // this is for new scheme where also the flags are coded in numbers in the PWGGA.... name
 
       if(baseDirName.Length()>31){
        reconstructionFlagString = baseDirName(baseDirName.Index("GammaConversion_")+16,8);
index ac8808e..9f74ee4 100644 (file)
@@ -1,5 +1,5 @@
 /****************************************************************************************************************************
-******         provided by Gamma Conversion Group, PWG4,                                                                                                               *****
+******         provided by Gamma Conversion Group, PWGGA,                                                                                                              *****
 ******         Kathrin Koch, kkoch@physi.uni-heidelberg.de                                                                                                     *****
 ******         Friederike Bock, friederike.bock@cern.ch                                                                                                        *****
 *****************************************************************************************************************************
@@ -59,12 +59,12 @@ void  Photon_Characteristics_Events(const char *data = "myOutput", const char *M
                    Scaling->SetLineWidth(2);   
 
        if(cutsel != ""){
-               char *GammaDirectory = Form("PWG4_GammaConversion_%s",  cutsel);
+               char *GammaDirectory = Form("PWGGA_GammaConversion_%s",  cutsel);
                cout << GammaDirectory << endl;
                char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
                cout << GammaList << endl;
        }else{
-               char *GammaDirectory = "PWG4_GammaConversion";
+               char *GammaDirectory = "PWGGA_GammaConversion";
                cout << GammaDirectory << endl;
                char *GammaList = "histogramsAliGammaConversion";
                cout << GammaList << endl;
@@ -127,13 +127,13 @@ void  Photon_Characteristics_Events(const char *data = "myOutput", const char *M
        TFile f(filename_data);  
        
        //for new versions
-        TDirectory *fPWG4GammaConversion_data = new TDirectory(); // definition of first folder / list 
+        TDirectory *fPWGGAGammaConversion_data = new TDirectory(); // definition of first folder / list        
         TList *fHistosGammaConversion_data = new TList(); // definition of first folder / list
         TList *fESDContainer_data = new TList();  // definition of following folder / list
        TList *fMappingContainer_data = new TList();
         
-        if(!(fPWG4GammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl; 
-        if(!(fHistosGammaConversion_data = (TList*)fPWG4GammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
+        if(!(fPWGGAGammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl; 
+        if(!(fHistosGammaConversion_data = (TList*)fPWGGAGammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
         if(!(fESDContainer_data = (TList*)fHistosGammaConversion_data->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl; 
        fMappingContainer_data = (TList*)fHistosGammaConversion_data->FindObject("Mapping histograms"); 
 
@@ -222,13 +222,13 @@ void  Photon_Characteristics_Events(const char *data = "myOutput", const char *M
                
                
                // for new versions
-                TDirectory *fPWG4GammaConversion_montecarlo = new TDirectory(); // definition of first folder / list   
+                TDirectory *fPWGGAGammaConversion_montecarlo = new TDirectory(); // definition of first folder / list  
                 TList *fHistosGammaConversion_montecarlo = new TList(); // definition of first folder / list
                 TList *fESDContainer_montecarlo = new TList();  // definition of following folder / list
                 TList *fMappingContainer_montecarlo = new TList();  // definition of following folder / list
                 TList *fMCtruth_montecarlo = new TList();
-                if(!(fPWG4GammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl; 
-                if(!(fHistosGammaConversion_montecarlo = (TList*)fPWG4GammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
+                if(!(fPWGGAGammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl; 
+                if(!(fHistosGammaConversion_montecarlo = (TList*)fPWGGAGammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
                 if(!(fESDContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl; 
                fMappingContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("Mapping histograms");              
                fMCtruth_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("MC histograms");            
index 02febda..4e02798 100644 (file)
@@ -1,5 +1,5 @@
 /***********************************************************************************************
-*** provided by Gamma Conversion Group, PWG4,                                                                  ******
+*** provided by Gamma Conversion Group, PWGGA,                                                                         ******
 ***        Friederike Bock, fbock@physi.uni-heidelberg.de ***                                                  ******
 ************************************************************************************************
 ************************************************************************************************
@@ -74,12 +74,12 @@ void  Photon_Resolution(const char *MCfile = "", const char *cutsel = "",const c
        Char_t filename_MC[100] = (Form("%s%s",path,MCfile));
 
        if(cutsel != ""){
-               char *GammaDirectory = Form("PWG4_GammaConversion_%s",  cutsel);
+               char *GammaDirectory = Form("PWGGA_GammaConversion_%s",  cutsel);
                cout << GammaDirectory << endl;
                char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
                cout << GammaList << endl;
        }else{
-               char *GammaDirectory = "PWG4_GammaConversion";
+               char *GammaDirectory = "PWGGA_GammaConversion";
                cout << GammaDirectory << endl;
                char *GammaList = "histogramsAliGammaConversion";
                cout << GammaList << endl;
@@ -92,12 +92,12 @@ void  Photon_Resolution(const char *MCfile = "", const char *cutsel = "",const c
        TFile f(filename_MC);  
        
        //************************** Container Loading ********************************************************************
-       TDirectory *fPWG4GammaConversion_montecarlo = new TDirectory(); 
+       TDirectory *fPWGGAGammaConversion_montecarlo = new TDirectory(); 
        TList *fHistosGammaConversion_montecarlo = new TList(); 
        TList *fResolutionContainer_montecarlo = new TList();  
      TList *fESDContainer_montecarlo = new TList();
-       if(!(fPWG4GammaConversion_montecarlo = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl; 
-       if(!(fHistosGammaConversion_montecarlo = (TList*)fPWG4GammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
+       if(!(fPWGGAGammaConversion_montecarlo = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl; 
+       if(!(fHistosGammaConversion_montecarlo = (TList*)fPWGGAGammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl; 
        if(fResolutionContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("Resolution histograms")) cout<<"list loaded correctly!"<<endl; 
        if(!(fESDContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl; 
        TH1F * ESD_ConvGamma_Pt=fESDContainer_montecarlo->FindObject("ESD_ConvGamma_Pt");
index 6178f97..35f1046 100644 (file)
@@ -1,4 +1,4 @@
-// provided by Gamma Conversion Group, PWG4, Kathrin Koch, kkoch@physi.uni-heidelberg.de and Friederike Bock, fbock@physi.uni-heidelberg.de
+// provided by Gamma Conversion Group, PWGGA, Kathrin Koch, kkoch@physi.uni-heidelberg.de and Friederike Bock, fbock@physi.uni-heidelberg.de
 
 #include <Riostream.h>
 #include <fstream>
@@ -70,12 +70,12 @@ void  Plot_Mapping_Histos_Events(const char *input1 = "myOutput", const char *se
 
        // choice of dateset
        if(cutsel != ""){
-               char *GammaDirectory = Form("PWG4_GammaConversion_%s",  cutsel);
+               char *GammaDirectory = Form("PWGGA_GammaConversion_%s",  cutsel);
                cout << GammaDirectory << endl;
                char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
                cout << GammaList << endl;
        }else{
-               char *GammaDirectory = "PWG4_GammaConversion";
+               char *GammaDirectory = "PWGGA_GammaConversion";
                cout << GammaDirectory << endl;
                char *GammaList = "histogramsAliGammaConversion";
                cout << GammaList << endl;
@@ -91,13 +91,13 @@ void  Plot_Mapping_Histos_Events(const char *input1 = "myOutput", const char *se
                    EtaRange->SetLineWidth(6);  
 
 //------------------------------- Reading FILES        ----------------------------------------------------------------------
-        TDirectory *fPWG4GammaConversion_input1 = new TDirectory(); // definition of first folder / list       
+        TDirectory *fPWGGAGammaConversion_input1 = new TDirectory(); // definition of first folder / list      
         TList *fHistosGammaConversion_input1 = new TList(); // definition of first folder / list
         TList *fESDContainer_input1 = new TList();  // definition of following folder / list
         TList *fMappingContainer_input1 = new TList();
 
-       fPWG4GammaConversion_input1 = (TDirectory*)f1.Get(GammaDirectory); 
-       fHistosGammaConversion_input1 = (TList*)fPWG4GammaConversion_input1->Get(GammaList); 
+       fPWGGAGammaConversion_input1 = (TDirectory*)f1.Get(GammaDirectory); 
+       fHistosGammaConversion_input1 = (TList*)fPWGGAGammaConversion_input1->Get(GammaList); 
        fMappingContainer_input1 = (TList*)fHistosGammaConversion_input1->FindObject("Mapping histograms"); 
        fESDContainer_input1 = (TList*)fHistosGammaConversion_input1->FindObject("ESD histograms"); 
        
@@ -232,13 +232,13 @@ void  Plot_Mapping_Histos_Events(const char *input1 = "myOutput", const char *se
                
                input2 = new TFile(Form("%s%s",path, secondinput));     
 
-                TDirectory *fPWG4GammaConversion_input2 = new TDirectory(); // definition of first folder / list       
+                TDirectory *fPWGGAGammaConversion_input2 = new TDirectory(); // definition of first folder / list      
                 TList *fHistosGammaConversion_input2 = new TList(); // definition of first folder / list
                 TList *fMappingContainer_input2 = new TList();  // definition of following folder / list
                 TList *fESDContainer_input2 = new TList();  // definition of following folder / list
                 
-               fPWG4GammaConversion_input2 = (TDirectory*)input2->Get(GammaDirectory); 
-               fHistosGammaConversion_input2 = (TList*)fPWG4GammaConversion_input2->Get(GammaList); 
+               fPWGGAGammaConversion_input2 = (TDirectory*)input2->Get(GammaDirectory); 
+               fHistosGammaConversion_input2 = (TList*)fPWGGAGammaConversion_input2->Get(GammaList); 
                fMappingContainer_input2 = (TList*)fHistosGammaConversion_input2->FindObject("Mapping histograms");
                fESDContainer_input2 = (TList*)fHistosGammaConversion_input2->FindObject("ESD histograms"); 
 
index a09c8b4..6a44df5 100644 (file)
@@ -1,5 +1,5 @@
 /***********************************************************************************************
-*** provided by Gamma Conversion Group, PWG4,                                                                  ******
+*** provided by Gamma Conversion Group, PWGGA,                                                                         ******
 ***        Friederike Bock, fbock@physi.uni-heidelberg.de ***                                                  ******
 ************************************************************************************************/
 
index 596df6d..5693548 100644 (file)
@@ -1,5 +1,5 @@
 /***********************************************************************************************
-*** provided by Gamma Conversion Group, PWG4,                                                                  ******
+*** provided by Gamma Conversion Group, PWGGA,                                                                         ******
 ***        Friederike Bock, fbock@physi.uni-heidelberg.de ***                                                  ******
 ************************************************************************************************/
 
index c44984e..e1ed58e 100755 (executable)
@@ -49,7 +49,7 @@ if [ -f $OutputDirectory/$DatFilenameBase.dat ]; then
 fi
 
 
-root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/MakeCutLog.C\(\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\)
+root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/MakeCutLog.C\(\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\)
 
 # Read the different cuts form the Cut selection log file
 exec<"$OutputDirectory/CutSelection.log"
@@ -57,12 +57,12 @@ exec<"$OutputDirectory/CutSelection.log"
 while read cutSelection
 do
     echo CutSelection is $cutSelection;
-    root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Extract_IntegratedPi0Yield.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\"\);
+    root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Extract_IntegratedPi0Yield.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\"\);
     
-    root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Extract_Pi0_Characteristics.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\);
+    root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Extract_Pi0_Characteristics.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\);
 
 done
 
-    root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Plot_IntegratedPi0Yield.C\(\"$DatFilenameBase\"\,\"$OutputDirectory\"\,\"$Suffix\"\)
+    root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Plot_IntegratedPi0Yield.C\(\"$DatFilenameBase\"\,\"$OutputDirectory\"\,\"$Suffix\"\)
 
-    root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Plot_Pi0_Characteristics.C
\ No newline at end of file
+    root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Plot_Pi0_Characteristics.C
\ No newline at end of file
index 5796f38..0eefc67 100644 (file)
@@ -5,8 +5,8 @@ AliAnalysisTaskCaloCellsQA* AddTaskCaloCellsQA(Int_t nmods = 10, Int_t det = 0,
   //
   // Usage example for EMCAL:
   //
-  //   gSystem->Load("libPWG4UserTasks.so");
-  //   gROOT->LoadMacro("$ALICE_ROOT/PWG4/UserTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
+  //   gSystem->Load("libPWGGAPHOSTasks.so");
+  //   gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
   //   AliAnalysisTaskCaloCellsQA *taskQA = AddTaskCaloCellsQA(10); // 10 supermodules
   //   taskQA->SelectCollisionCandidates(AliVEvent::kMB); // if necessary
   //   // taskQA->SetAvoidPileup(kFALSE); // some customization
@@ -16,8 +16,8 @@ AliAnalysisTaskCaloCellsQA* AddTaskCaloCellsQA(Int_t nmods = 10, Int_t det = 0,
   //
   // Usage example for PHOS:
   //
-  //   gSystem->Load("libPWG4UserTasks.so");
-  //   gROOT->LoadMacro("$ALICE_ROOT/PWG4/UserTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
+  //   gSystem->Load("libPWGAPHOSTasks.so");
+  //   gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
   //   AliAnalysisTaskCaloCellsQA *taskQA = AddTaskCaloCellsQA(4, 1);
   //   taskQA->SelectCollisionCandidates(AliVEvent::kMB); // if necessary
   //   taskQA->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.1); // increase statistics
index 7ebfdd2..01b2919 100644 (file)
@@ -1,6 +1,6 @@
 AliAnalysisTaskSE * AddTaskPHOSQA()
 {
-  gROOT->LoadMacro("$ALICE_ROOT/PWG4/UserTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C"); 
+  gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C"); 
 
   AliAnalysisTaskCaloCellsQA *taskPHOSCellQA = AddTaskCaloCellsQA(4, 1, "CaloCellsQA.root","PHOSCellsQA");
   taskPHOSCellQA->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.3,1.0); 
index f8ce5ea..b646665 100644 (file)
@@ -3,7 +3,7 @@ void ana_example()
   // Example macro to run QA
 
   // load relevant library
-  gSystem->Load("libPWG4UserTasks.so");
+  gSystem->Load("libPWGGAPHOSTasks.so");
 
   // change next line to a working code
   TChain* chain = NULL; //CreateChain("wn.xml");
index c6b964c..ad2cc3a 100644 (file)
@@ -9,7 +9,7 @@ void PHOSPbPbQA(const char* dataset="collection.xml")
   gSystem->Load("libANALYSIS");
   gSystem->Load("libANALYSISalice"); //AliAnalysisTaskSE
 
-  gSystem->Load("libPWG4UserTasks.so");
+  gSystem->Load("libPWGGAPHOSTasks.so");
   
   gSystem->AddIncludePath("-I$ALICE_ROOT/include -I$ALICE_ROOT/PHOS");
 
index a0987c8..f5df112 100644 (file)
@@ -26,7 +26,7 @@ AliAnalysisTaskOmegaPi0PiPi* AddTaskomega3pi()
 
 
   TString outputfile = AliAnalysisManager::GetCommonFileName();
-  AliAnalysisDataContainer *coutput = mgr->CreateContainer("omega3pi", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWG4_omega3pi",outputfile.Data()));
+  AliAnalysisDataContainer *coutput = mgr->CreateContainer("omega3pi", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWGGA_omega3pi",outputfile.Data()));
   
   // Connect input/output
   mgr->ConnectInput(omegaTask, 0, cinput);
index 53ce78c..768ea34 100644 (file)
@@ -9,7 +9,7 @@ void AnaTask(const char* dataset="minbias_LHC09a4_81040_81050.xml")
     //load analysis framework
     gSystem->Load("libANALYSIS");
     gSystem->Load("libANALYSISalice"); //AliAnalysisTaskSE
-    gSystem->Load("libPWG4Omega3pi");
+    gSystem->Load("libPWGGAPHOSTasks");
     gSystem->AddIncludePath("-I$ALICE_ROOT/include -I$ALICE_ROOT/PHOS");
 
     // A task can be compiled dynamically with AClic