Obsolete
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 29 Mar 2010 07:45:06 +0000 (07:45 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 29 Mar 2010 07:45:06 +0000 (07:45 +0000)
41 files changed:
prod/acrcaf/qa_pp.sh [deleted file]
prod/acrcaf/qa_pp/AddTaskCalorimeterQA.C [deleted file]
prod/acrcaf/qa_pp/AddTaskQAsym.C [deleted file]
prod/acrcaf/qa_pp/AddTaskVZEROQA.C [deleted file]
prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx [deleted file]
prod/acrcaf/qa_pp/AliAnalysisTaskQASym.h [deleted file]
prod/acrcaf/qa_pp/plot_macros/plot.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotD.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotDca.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotDcaVz.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotDeltaPhi.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotEta-over-pt.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotEta.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotEtaPtVz.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotEtaVz.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotLadderDca.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotLadderEta.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotLadderPhi.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotLadderPt.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotLadderSigma.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotPhi.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotPhiVz.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotPhi_Eta.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotPt.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotPtVz.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotPt_Eta.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotQPtVz.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotQPt_Eta.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorDca.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorDcaRatio.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorEta.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPt.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPtRatio.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorEtaRatio.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorPt.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorPtRatio.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorQPt.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/plotSectorQPtRatio.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/pos-NegInv.C [deleted file]
prod/acrcaf/qa_pp/plot_macros/sigma.C [deleted file]
prod/acrcaf/qa_pp/qa_pp.C [deleted file]

diff --git a/prod/acrcaf/qa_pp.sh b/prod/acrcaf/qa_pp.sh
deleted file mode 100755 (executable)
index c8bb119..0000000
+++ /dev/null
@@ -1,14 +0,0 @@
-#!/bin/bash
-
-source /afs/cern.ch/alice/caf/caf-lxplus.sh -alien v4-18-12-AN
-
-  [ -d qa_pp ] || mkdir qa_pp
-
-  [ -z $1 ] && { echo "Usage: qa_pp.sh <run_number>"; exit 1; }
-
-  cd qa_pp
- root.exe qa_pp.C\($1\)
-
-  rm plot_macros/QAsym.proof.root
-  ln -s /home/alishift/acrcaf/qa_pp/run$1.root plot_macros/QAsym.proof.root
-
diff --git a/prod/acrcaf/qa_pp/AddTaskCalorimeterQA.C b/prod/acrcaf/qa_pp/AddTaskCalorimeterQA.C
deleted file mode 100644 (file)
index ac70c8c..0000000
+++ /dev/null
@@ -1,113 +0,0 @@
-AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kUseKinematics = kFALSE, Bool_t kPrintSettings = kFALSE)\r
-{\r
-  // Creates a PartCorr task for calorimeters performance studies, configures it and adds it to the analysis manager.\r
-  \r
-  // Get the pointer to the existing analysis manager via the static access method.\r
-  //==============================================================================\r
-  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
-  if (!mgr) {\r
-    ::Error("AddTaskPartCorr", "No analysis manager to connect to.");\r
-    return NULL;\r
-  }  \r
-  \r
-  // Check the analysis type using the event handlers connected to the analysis manager.\r
-  //==============================================================================\r
-  if (!mgr->GetInputEventHandler()) {\r
-    ::Error("AddTaskPartCorr", "This task requires an input event handler");\r
-    return NULL;\r
-  }\r
-   TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
-  \r
-   // Configure analysis\r
-   //===========================================================================\r
-  \r
-  //Reader\r
-  //For this particular analysis few things done by the reader.\r
-  //Nothing else needs to be set.\r
-  AliCaloTrackReader * reader = 0x0;\r
-  if(data=="AOD")      reader = new AliCaloTrackAODReader();\r
-  else if(data=="ESD") reader = new AliCaloTrackESDReader();\r
-  //reader->SetDebug(10);//10 for lots of messages\r
-   if(kPrintSettings) reader->Print("");\r
-  \r
-  if(kUseKinematics){\r
-               if(inputDataType == "ESD"){\r
-                       reader->SwitchOnStack();          \r
-                       reader->SwitchOffAODMCParticles(); \r
-               }\r
-               else if(inputDataType == "AOD"){\r
-                       reader->SwitchOffStack();          \r
-                       reader->SwitchOnAODMCParticles(); \r
-               }\r
-   }\r
-       \r
-   reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches.\r
-\r
-       \r
-  // ##### Analysis algorithm settings ####\r
-   //Only needed now for MC data\r
-   //AliFiducialCut * fidCut = new AliFiducialCut();\r
-   //fidCut->DoCTSFiducialCut(kFALSE) ;\r
-   //fidCut->DoEMCALFiducialCut(kTRUE) ;\r
-   //fidCut->DoPHOSFiducialCut(kTRUE) ;\r
-               \r
-  AliAnaCalorimeterQA *emcalQA = new AliAnaCalorimeterQA();\r
-  //emcalQA->SetDebug(2); //10 for lots of messages\r
-  emcalQA->SetCalorimeter("EMCAL");\r
-  if(kUseKinematics) emcalQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.\r
-  else  emcalQA->SwitchOffDataMC() ;\r
-  emcalQA->AddToHistogramsName("EMCAL_"); //Begining of histograms name\r
-  //emcalQA->SetFiducialCut(fidCut);\r
-  emcalQA->SwitchOffFidutialCut();\r
-  if(kPrintSettings) emcalQA->Print("");       \r
-  //emcalQA->SwithOnPlotsMaking();\r
-  //emcalQA->GetMCAnalysisUtils()->SetDebug(10);\r
-       \r
-   AliAnaCalorimeterQA *phosQA = new AliAnaCalorimeterQA();\r
-   //phosQA->SetDebug(2); //10 for lots of messages\r
-   phosQA->SetCalorimeter("PHOS");\r
-   if(kUseKinematics) phosQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.\r
-   else  phosQA->SwitchOffDataMC() ;  \r
-   phosQA->AddToHistogramsName("PHOS_");//Begining of histograms name\r
-   //phosQA->SetFiducialCut(fidCut);\r
-   phosQA->SwitchOffFidutialCut();\r
-   if(kPrintSettings)phosQA->Print("");        \r
-   //phosQA->GetMCAnalysisUtils()->SetDebug(10);\r
-   //phosQA->SwithOnPlotsMaking();\r
-  \r
-  // #### Configure Maker ####\r
-  AliAnaPartCorrMaker * maker = new AliAnaPartCorrMaker();\r
-  maker->SetReader(reader);//pointer to reader\r
-  maker->AddAnalysis(emcalQA,0);\r
-  maker->AddAnalysis(phosQA,1);\r
-  maker->SetAnaDebug(0)  ;\r
-  maker->SwitchOnHistogramsMaker()  ;\r
-  maker->SwitchOffAODsMaker()  ;\r
-  if(kPrintSettings) maker->Print("");\r
\r
-\r
-  printf("======================== \n");\r
-  printf(" End Configuration of Calorimeter QA \n");\r
-  printf("======================== \n");\r
-  \r
-   // Create task\r
-   //===========================================================================\r
-  AliAnalysisTaskParticleCorrelation * task = new AliAnalysisTaskParticleCorrelation ("CalorimeterPerformance");\r
-  task->SetConfigFileName(""); //Don't configure the analysis via configuration file.\r
-  task->SetDebugLevel(0);\r
-  task->SetAnalysisMaker(maker);                               \r
-  mgr->AddTask(task);\r
-  \r
-  AliAnalysisDataContainer *cout_pc = mgr->CreateContainer("Calo.Performance",TList::Class(),\r
-                                                          AliAnalysisManager::kOutputContainer, "Calo.Performance.root");\r
-  \r
-  // Create ONLY the output containers for the data produced by the task.\r
-  // Get and connect other common input/output containers via the manager as below\r
-  //==============================================================================\r
-  mgr->ConnectInput  (task, 0, mgr->GetCommonInputContainer());\r
-  mgr->ConnectOutput (task, 1, cout_pc);\r
-  \r
-  return task;\r
-}\r
-\r
-\r
diff --git a/prod/acrcaf/qa_pp/AddTaskQAsym.C b/prod/acrcaf/qa_pp/AddTaskQAsym.C
deleted file mode 100644 (file)
index b16cf4a..0000000
+++ /dev/null
@@ -1,120 +0,0 @@
-AliAnalysisTaskQASym * AddTaskQAsym(Int_t runNumber)\r
-\r
-{\r
-  // Creates a QA task exploiting simple symmetries phi, eta +/-, charge ...\r
-  \r
-  // Get the pointer to the existing analysis manager via the static access method.\r
-  //==============================================================================\r
-  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
-  if (!mgr) {\r
-    ::Error("AddTaskQAsym", "No analysis manager to connect to.");\r
-    return NULL;\r
-  }  \r
-  \r
-  // Check the analysis type using the event handlers connected to the analysis manager.\r
-  //==============================================================================\r
-  if (!mgr->GetInputEventHandler()) {\r
-    ::Error("AddTasQAsym", "This task requires an input event handler");\r
-    return NULL;\r
-  }\r
-   TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
-  \r
-   // Configure analysis\r
-   //===========================================================================\r
-   \r
\r
-   //Task for global tracks\r
-   AliAnalysisTaskQASym *task0 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_Global");\r
-   task0->SetTrackType(0);\r
-   task0->SelectCollisionCandidates();\r
-   //Task for ITS tracks \r
-   AliAnalysisTaskQASym *task1 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_ITS");\r
-   task1->SetTrackType(1);\r
-   task1->SelectCollisionCandidates();\r
-   //Task for TPC tracks \r
-   AliAnalysisTaskQASym *task2 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_TPC");\r
-   task2->SetTrackType(2);\r
-   task2->SelectCollisionCandidates();\r
-\r
-   //cuts for global tracks\r
-   AliESDtrackCuts* esdTrackCutsL0 = new AliESDtrackCuts("AliESDtrackCuts0","Global");\r
-   esdTrackCutsL0->SetMinNClustersTPC(70);\r
-   esdTrackCutsL0->SetRequireTPCRefit(kTRUE);\r
-   esdTrackCutsL0->SetMaxDCAToVertexXY(3.);\r
-   esdTrackCutsL0->SetMaxDCAToVertexZ(3.);\r
-   esdTrackCutsL0->SetAcceptKinkDaughters(kFALSE);\r
-\r
-   AliESDtrackCuts* esdTrackCutsL1 = new AliESDtrackCuts("AliESDtrackCuts1","ITS");\r
-   esdTrackCutsL1->SetMaxDCAToVertexXY(3.);\r
-   esdTrackCutsL1->SetMaxDCAToVertexZ(3.);\r
-   esdTrackCutsL1->SetAcceptKinkDaughters(kFALSE);\r
-   esdTrackCutsL1->SetRequireITSRefit(kTRUE);\r
-   esdTrackCutsL1->SetRequireITSStandAlone(kTRUE);\r
-   \r
-   //cuts for TPC tracks\r
-   AliESDtrackCuts* esdTrackCutsL2 = new AliESDtrackCuts("AliESDtrackCuts2","TPC");\r
-   esdTrackCutsL2->SetRequireTPCRefit(kFALSE);\r
-   esdTrackCutsL2->SetAcceptKinkDaughters(kFALSE);\r
-   //jacek's cuts:\r
-   esdTrackCutsL2->SetMinNClustersTPC(70);\r
-   // cut on max ncl=160 in Task\r
-   esdTrackCutsL2->SetMaxDCAToVertexXY(3.);\r
-   esdTrackCutsL2->SetMaxDCAToVertexZ(3.);\r
-   esdTrackCutsL2->SetMaxChi2PerClusterTPC(3.999);\r
-   //cut minChi=0 in task\r
-   //esdTrackCutsL2->SetPRange(0.15,16); // not needed for QA\r
-   //esdTrackCutsL2->SetEtaRange(-0.8, 0.7999); // not needed for QA\r
-  \r
-\r
-   task0->SetCuts(esdTrackCutsL0);\r
-   task1->SetCuts(esdTrackCutsL1);\r
-   task2->SetCuts(esdTrackCutsL2);\r
-\r
-   mgr->AddTask(task0);\r
-   mgr->AddTask(task1);\r
-   mgr->AddTask(task2);\r
-  \r
-   AliAnalysisDataContainer *cout0  = 0;\r
-   AliAnalysisDataContainer *cout1  = 0;\r
-   AliAnalysisDataContainer *cout2  = 0;\r
-   \r
-   if(runNumber>0){ \r
-    cout0 =  mgr->CreateContainer("QAsymHists_Global",TList::Class(),\r
-                                 AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber));\r
-    cout1 =  mgr->CreateContainer("QAsymHists_ITS",TList::Class(),\r
-                                 AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber));\r
-    cout2 =  mgr->CreateContainer("QAsymHists_TPC",TList::Class(),\r
-                                 AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber));\r
-   }\r
-   \r
-   else{\r
-      cout0 = mgr->CreateContainer("QAsymHists_Global",TList::Class(),\r
-                                AliAnalysisManager::kOutputContainer, \r
-                                Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName()));\r
-      cout1 = mgr->CreateContainer("QAsymHists_ITS",TList::Class(),\r
-                                  AliAnalysisManager::kOutputContainer, \r
-                                Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName()));\r
-      cout2 = mgr->CreateContainer("QAsymHists_TPC",TList::Class(),\r
-                                  AliAnalysisManager::kOutputContainer, \r
-                                Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName()));\r
-   }\r
-\r
-\r
-   mgr->ConnectInput  (task0, 0, mgr->GetCommonInputContainer());\r
-   mgr->ConnectInput  (task1, 0, mgr->GetCommonInputContainer());\r
-   mgr->ConnectInput  (task2, 0, mgr->GetCommonInputContainer());\r
-   \r
-   mgr->ConnectOutput (task0, 0, mgr->GetCommonOutputContainer());\r
-   mgr->ConnectOutput (task1, 0, mgr->GetCommonOutputContainer());\r
-   mgr->ConnectOutput (task2, 0, mgr->GetCommonOutputContainer());\r
-\r
-   mgr->ConnectOutput (task0, 1, cout0);\r
-   mgr->ConnectOutput (task1, 1, cout1);\r
-   mgr->ConnectOutput (task2, 1, cout2);\r
-\r
-  \r
-   return task0;\r
-\r
-}\r
-\r
-\r
diff --git a/prod/acrcaf/qa_pp/AddTaskVZEROQA.C b/prod/acrcaf/qa_pp/AddTaskVZEROQA.C
deleted file mode 100644 (file)
index a8ed634..0000000
+++ /dev/null
@@ -1,40 +0,0 @@
-AliAnalysisTaskSE* AddTaskVZEROQA(Int_t runNumber)\r
-{\r
-  // Creates a QA task exploiting simple symmetries phi, eta +/-, charge ...\r
-  \r
-  // Get the pointer to the existing analysis manager via the static access method.\r
-  //==============================================================================\r
-  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
-  if (!mgr) {\r
-    ::Error("AddTaskQAsym", "No analysis manager to connect to.");\r
-    return NULL;\r
-  }  \r
-  \r
-  // Check the analysis type using the event handlers connected to the analysis manager.\r
-  //==============================================================================\r
-  if (!mgr->GetInputEventHandler()) {\r
-    ::Error("AddTasQAsym", "This task requires an input event handler");\r
-    return NULL;\r
-  }\r
-   TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
-  \r
-   // Configure analysis\r
-   //===========================================================================\r
-   \r
\r
\r
-   AliAnaVZEROQA* task = new AliAnaVZEROQA("AliAnaVZEROQA");\r
-   mgr->AddTask(task);\r
-  \r
-   AliAnalysisDataContainer *cout  = mgr->CreateContainer("QAVZEROHists",TList::Class(),\r
-                                                       AliAnalysisManager::kOutputContainer, "VZERO.Performance.root");\r
-\r
-   mgr->ConnectInput  (task, 0, mgr->GetCommonInputContainer());\r
-   mgr->ConnectOutput (task, 1, cout);\r
-\r
-   return task;\r
-   \r
-  \r
-}\r
-\r
-\r
diff --git a/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx b/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx
deleted file mode 100644 (file)
index 7d64234..0000000
+++ /dev/null
@@ -1,1298 +0,0 @@
-#include "TChain.h"
-#include "TTree.h"
-#include "TH1F.h"
-#include "TH2F.h"
-#include "TH3F.h"
-#include "TCanvas.h"
-#include "TList.h"
-#include "TParticle.h"
-#include "TParticlePDG.h"
-#include "TProfile.h"
-#include "TNtuple.h"
-#include "TFile.h"
-
-#include "AliAnalysisTask.h"
-#include "AliAnalysisManager.h"
-
-#include "AliESDEvent.h"
-#include "AliLog.h"
-#include "AliESDVertex.h"
-#include "AliESDInputHandler.h"
-#include "AliESDtrackCuts.h"
-#include "AliMultiplicity.h"
-
-#include "AliAnalysisTaskQASym.h"
-#include "AliExternalTrackParam.h"
-#include "AliTrackReference.h"
-#include "AliHeader.h"
-#include "AliGenEventHeader.h"
-#include "AliGenDPMjetEventHeader.h"
-
-// Analysis Task for basic QA on the ESD
-
-// Authors: Jan Fiete Grosse-Oetringhaus, Christian Klein-Boesing,
-//          Andreas Morsch, Eva Sicking
-
-ClassImp(AliAnalysisTaskQASym)
-
-  //________________________________________________________________________
-  AliAnalysisTaskQASym::AliAnalysisTaskQASym(const char *name) 
-    : AliAnalysisTaskSE(name) 
-    ,fTrackType(0)
-    ,fFieldOn(kTRUE)
-
-    ,fHists(0)
-
-    ,fHistRECpt(0)
-    ,fEta(0)
-    ,fEtaPhi(0)
-    ,fEtaPt(0)
-    ,fQPt(0)
-    ,fDca(0)
-    ,fqRec(0)
-    ,fsigmaPt(0)
-  
-    ,fRecPtPos(0)
-    ,fRecPtNeg(0)
-    ,fRecPhiPos(0)
-    ,fRecPhiNeg(0)
-    ,fRecEtaPos(0)
-    ,fRecEtaNeg(0)
-    ,fRecEtaPtPos(0)
-    ,fRecEtaPtNeg(0)
-    ,fRecDcaPos(0)
-    ,fRecDcaNeg(0)
-    ,fRecDcaNegInv(0)
-    ,fRecDPos(0)
-    ,fRecDNeg(0)
-
-
-    ,fRecQPtPosEta(0)
-    ,fRecQPtNegEta(0)
-    ,fRecPtPosEta(0)
-    ,fRecPtNegEta(0)
-    ,fRecPhiPosEta(0)
-    ,fRecPhiNegEta(0)
-    ,fRecDcaPosEta(0)
-    ,fRecDcaNegEta(0)
-    ,fRecDPosEta(0)
-    ,fRecDNegEta(0)
-  
-    ,fRecPtPosVz(0)
-    ,fRecPtNegVz(0)
-    ,fRecEtaPosVz(0)
-    ,fRecEtaNegVz(0)
-    ,fRecPhiPosVz(0)
-    ,fRecPhiNegVz(0)
-    ,fSignedDcaPosVz(0)
-    ,fSignedDcaNegVz(0)
-    ,fRecQPtPosEtaVz(0)
-    ,fRecQPtNegEtaVz(0)
-    ,fRecEtaPtPosVz(0)
-    ,fRecEtaPtNegVz(0)
-  
-  
-    ,fDeltaPhiAll(0)
-    ,fDeltaPhiLeading(0) 
-    ,fDiffDcaD(0)
-
-    ,fPhiRec(0)
-    ,fThetaRec(0)
-    ,fNumber(0)
-    ,fVx(0)
-    ,fVy(0)
-    ,fVz(0)
-    ,test(0)
-  
-    ,fRecDcaPosPhi(0)
-    ,fRecDcaNegPhi(0)
-    ,fRecPtPosPhi(0)
-    ,fRecPtNegPhi(0)
-    ,fRecEtaPosPhi(0)
-    ,fRecEtaNegPhi(0)
-    ,fRecQPtPhi(0)
-    ,fRecEtaPtPosPhi(0)
-    ,fRecEtaPtNegPhi(0)
-
-    ,fRecPtPosEtaPos(0)
-    ,fRecPtNegEtaPos(0)
-    ,fRecPtPosEtaNeg(0)
-    ,fRecPtNegEtaNeg(0)
-
-    ,fRec1PtPosEtaPos(0)
-    ,fRec1PtNegEtaPos(0)
-    ,fRec1PtPosEtaNeg(0)
-    ,fRec1PtNegEtaNeg(0)
-
-    ,fRecPhiPosEtaPos(0)
-    ,fRecPhiNegEtaPos(0)
-    ,fRecPhiPosEtaNeg(0)
-    ,fRecPhiNegEtaNeg(0)
-
-    ,fRecDcaPosPhiEtaPos(0)
-    ,fRecDcaNegPhiEtaPos(0) 
-    ,fRecDcaPosPhiEtaNeg(0)  
-    ,fRecDcaNegPhiEtaNeg(0)  
-    ,fRecPtPosPhiEtaPos(0)  
-    ,fRecPtNegPhiEtaPos(0)  
-    ,fRecPtPosPhiEtaNeg(0) 
-    ,fRecPtNegPhiEtaNeg(0) 
-
-
-    ,fRecDcaPhiPtPosEtaPos(0)
-    ,fRecDcaPhiPtNegEtaPos(0)
-    ,fRecDcaPhiPtPosEtaNeg(0)  
-    ,fRecDcaPhiPtNegEtaNeg(0)  
-    ,fEtavPt(0)  
-
-
-    
-    ,sdca(0)
-    ,xy(0)
-    ,z(0)
-    ,xvertexcor(0)
-    ,yvertexcor(0)  
-    ,fCuts(0)
-
-{
-  // Constructor
-  for(Int_t i = 0;i<18;++i){
-    fRecPtTpcSector[i] = 0;
-    fRecEtaTpcSector[i] = 0;
-    fSignedDcaTpcSector[i] = 0;
-    fRecQPtTpcSector[i] = 0;
-    fRecEtaPtTpcSector[i] = 0;
-  }
-
-  for(Int_t i = 0;i< 7;++i){
-    fRecPtPosLadder[i] = 0;
-    fRecPtNegLadder[i] = 0;
-    fRecPhiPosLadder[i] = 0;
-    fRecPhiNegLadder[i] = 0;
-    fRecEtaPosLadder[i] = 0;
-    fRecEtaNegLadder[i] = 0;
-    fSignDcaPos[i] = 0;
-    fSignDcaNeg[i] = 0;
-    fSignDcaNegInv[i] = 0;
-    fPtSigmaPos[i] =0;
-    fPtSigmaNeg[i] =0;
-    fqPtRec[i] =0;
-    fDcaSigmaPos[i] =0;
-    fDcaSigmaNeg[i] =0;
-  }
-
-  DefineOutput(1,  TList::Class()); 
-
-  
-  
-}
-
-
-//________________________________________________________________________
-void AliAnalysisTaskQASym::UserCreateOutputObjects()
-{
-  // Create histograms
-  // Called once
-
-  Bool_t oldStatus = TH1::AddDirectoryStatus();
-  TH1::AddDirectory(kFALSE);
-
-  Double_t range = 0.2;
-  Double_t pt = 20.;
-
-  fHists = new TList();
-  //  test   = new TNtuple("test","test",  
-  //                     "pt:phi:theta:x:y:z:charge");
-  fHistRECpt   = new TH1F("fHistRECpt", 
-                         " p_{T}",
-                         100, 0., 0.6);
-  fEta   = new TH1F("fEta", 
-                   " #eta",
-                   200, -2., 2.);
- fEtavPt   = new TH2F("fEtavPt", 
-                   " #eta -p_{T}",
-                    200, -2., 2.,
-                    100, 0, 1.5);
-  fEtaPhi   = new TH2F("fEtaPhi", 
-                      " #eta - #phi",
-                      200, -2., 2., 128, 0., 2. * TMath::Pi());
-  
-  fThetaRec   = new TH1F("fThetaRec", 
-                        " #theta",
-                        180, 0., TMath::Pi());
-  fPhiRec   = new TH1F("fPhiRec", 
-                      " #phi",
-                      180, 0., 2*TMath::Pi());
-  fNumber   = new TH1F("fNumber", 
-                      "number of tracks per event",
-                      50, 0.5, 49.5);
-  fVx   = new TH1F("fVx", 
-                  "X of vertex",
-                  100, -1., 1.);
-  fVy   = new TH1F("fVy", 
-                  "Y of vertex",
-                  100, -1., 1.);
-  fVz   = new TH1F("fVz", 
-                  "Z of vertex",
-                  200, -50., 50.);
-
-  fEtaPt   = new TH1F("fEtaPt", 
-                     " #eta/p_{T} ",
-                     100, -1., 1.);
-
-  fQPt   = new TH1F("fQPt", 
-                   " charge/p_{T} ",
-                   100, -1., 1.);
-
-  fDca   = new TH1F("fDca", 
-                   " dca ",
-                   200,  -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-
-
-  fqRec    = new TH1F("fqRec",   
-                     " charge all reconstructed particle",
-                     21, -9.5, 10.5);
-  
-  fsigmaPt    = new TH1F("fsigmaPt",   
-                        "Log_{10}(#sigma_{p_{T}})",
-                        200, -4., 8.);
-
-
-
-
-  //------------
-  for(Int_t ITSlayer_case=0;ITSlayer_case<7;ITSlayer_case++){
-
-    fSignDcaPos[ITSlayer_case]   = new TH1F(Form("fSignDcaPos%d", ITSlayer_case),  
-                                           " Signed dca", 
-                                           200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-    fSignDcaPos[ITSlayer_case]->GetXaxis()->SetTitle("dca");
-    fSignDcaPos[ITSlayer_case]->GetYaxis()->SetTitle("");
-   
-    fSignDcaNeg[ITSlayer_case]   = new TH1F(Form("fSignDcaNeg%d", ITSlayer_case),  
-                                           " Signed dcas",
-                                           200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-    fSignDcaNeg[ITSlayer_case]->GetXaxis()->SetTitle("dca");
-    fSignDcaNeg[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-    fSignDcaNegInv[ITSlayer_case]   = new TH1F(Form("fSignDcaNegInv%d", ITSlayer_case),  
-                                              " inverse Signed dca ",
-                                              200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-    fSignDcaNegInv[ITSlayer_case]->GetXaxis()->SetTitle("-dca");
-    fSignDcaNegInv[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-
-
-
-    fPtSigmaPos[ITSlayer_case]   = new TH1F(Form("fPtSigmaPos%d", ITSlayer_case),  
-                                           " #sigma_{pT} ",
-                                           208, -4., 8.);
-    fPtSigmaPos[ITSlayer_case]->GetXaxis()->SetTitle("Log_{10}(#sigma_{pT})");
-    fPtSigmaPos[ITSlayer_case]->GetYaxis()->SetTitle("");
-    
-    
-    fPtSigmaNeg[ITSlayer_case]   = new TH1F(Form("fPtSigmaNeg%d",ITSlayer_case),  
-                                           " #sigma_{pT}",
-                                           208, -4., 8.);
-    fPtSigmaNeg[ITSlayer_case]->GetXaxis()->SetTitle("Log_{10}(#sigma_{pT})");
-    fPtSigmaNeg[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-
-
-
-
-    fqPtRec[ITSlayer_case]   = new TH1F(Form("fqPtRec%d",ITSlayer_case),  
-                                       "q/ p_{T}",
-                                       200, -100., 100.);
-    fqPtRec[ITSlayer_case]->GetXaxis()->SetTitle("q_{tr}/p_{T, tr} (GeV/c)");
-    fqPtRec[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-  
-   
-
-
-    fDcaSigmaPos[ITSlayer_case]   = new TH2F(Form("fDcaSigmaPos%d", ITSlayer_case),  
-                                            " p_{T} shift vs #sigma_{pT} ",
-                                            200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200, -4., 4. );
-    fDcaSigmaPos[ITSlayer_case]->GetXaxis()->SetTitle("signed DCA)");
-    fDcaSigmaPos[ITSlayer_case]->GetYaxis()->SetTitle("log_{10}(#sigma_{pT})");
-    
-    
-    fDcaSigmaNeg[ITSlayer_case]   = new TH2F(Form("fDcaSigmaNeg%d", ITSlayer_case),  
-                                            " p_{T} shift vs #sigma_{pT} ",
-                                            200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200, -4., 4. );
-    fDcaSigmaNeg[ITSlayer_case]->GetXaxis()->SetTitle("signed DCA");
-    fDcaSigmaNeg[ITSlayer_case]->GetYaxis()->SetTitle("log_{10}(#sigma_{pT})");
-
-
-  }
-    
-    
-  // YIELDs---------- positive and negative particles
-    
-  fRecPtPos   = new TH1F("fRecPtPos", 
-                        " p_{T}",
-                        100, 0.,pt);
-  fRecPtPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-  fRecPtNeg   = new TH1F("fRecPtNeg", 
-                        " p_{T} ",
-                        100, 0., pt);
-  fRecPtNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-    
-  fRecPhiPos   = new TH1F("fRecPhiPos", 
-                         "#phi",
-                         361, 0., 360.);
-  fRecPhiPos->GetXaxis()->SetTitle("#phi (deg)");
-  
-  fRecPhiNeg   = new TH1F("fRecPhiNeg", 
-                         "#phi ",
-                         361, 0., 360.);
-  fRecPhiNeg->GetXaxis()->SetTitle("#phi (deg)");
-    
-  fRecEtaPos   = new TH1F("fRecEtaPos", 
-                         "#eta",
-                         200, -2., 2.);
-  fRecEtaPos->GetXaxis()->SetTitle("#eta");
-
-  fRecEtaNeg   = new TH1F("fRecEtaNeg", 
-                         "#eta",
-                         200, -2., 2.);
-  fRecEtaNeg->GetXaxis()->SetTitle("#eta");
-    
-  fRecEtaPtPos   = new TH1F("fRecEtaPtPos", 
-                           "#eta/p_{T}",
-                           200, -0.1, .1);
-  fRecEtaPtPos->GetXaxis()->SetTitle("#eta/p_{T}");
-
-  fRecEtaPtNeg   = new TH1F("fRecEtaPtNeg", 
-                           "#eta/p_{T}",
-                           200, -.1, .1);
-  fRecEtaPtNeg->GetXaxis()->SetTitle("#eta/p_{T}");
-
-  fRecDcaPos   = new TH1F("fRecDcaPos", 
-                         " dca",
-                         100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaPos->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaNeg   = new TH1F("fRecDcaNeg", 
-                         " dca",
-                         100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaNeg->GetXaxis()->SetTitle("dca (cm)");
-
-  fRecDcaNegInv   = new TH1F("fRecDcaNegInv", 
-                            " dca",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaNegInv->GetXaxis()->SetTitle("dca (cm)");
-
-
-  fRecDPos   = new TH1F("fRecDPos", 
-                       " d",
-                       100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDPos->GetXaxis()->SetTitle("d (cm)");
-  fRecDNeg   = new TH1F("fRecDNeg", 
-                       "d",
-                       100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDNeg->GetXaxis()->SetTitle("d (cm)");
-
-
-  //  YIELDs ---------------- positive and negative eta
-    
-    
-  fRecQPtPosEta   = new TH1F("fRecQPtPosEta", 
-                            "q/p_{T}",
-                            200, -0.5, 0.5);
-  fRecQPtPosEta->GetXaxis()->SetTitle("q/p_{T} ");
-
-  fRecQPtNegEta   = new TH1F("fRecQPtNegEta", 
-                            "q/p_{T}",
-                            200, -0.5, 0.5);
-  fRecQPtNegEta->GetXaxis()->SetTitle("q/p_{T}");
-    
-  fRecPtPosEta   = new TH1F("fRecPtPosEta", 
-                           " p_{T} ",
-                           100, 0., pt);
-  fRecPtPosEta->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-  fRecPtNegEta   = new TH1F("fRecPtNegEta", 
-                           " p_{T}",
-                           100, 0., pt);
-  fRecPtNegEta->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-    
-  fRecPhiPosEta   = new TH1F("fRecPhiPosEta", 
-                            "#phi",
-                            361, 0., 360);
-  fRecPhiPosEta->GetXaxis()->SetTitle("#phi (deg)");
-
-  fRecPhiNegEta   = new TH1F("fRecPhiNegEta", 
-                            "#phi ",
-                            361, 0, 360);
-  fRecPhiNegEta->GetXaxis()->SetTitle("#phi (deg)");
-
-  fRecDcaPosEta   = new TH1F("fRecDcaPosEta", 
-                            " dca ",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaPosEta->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaNegEta   = new TH1F("fRecDcaNegEta", 
-                            " dca",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaNegEta->GetXaxis()->SetTitle("dca (cm)");
-
-  fRecDPosEta   = new TH1F("fRecDPosEta", 
-                          " d",
-                          100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDPosEta->GetXaxis()->SetTitle("d (cm)");
-  fRecDNegEta   = new TH1F("fRecDNegEta", 
-                          "d",
-                          100, -5., 5.);
-  fRecDNegEta->GetXaxis()->SetTitle("d (cm)");
-
-  fRecDcaPosPhi   = new TH2F("fRecDcaPosPhi", 
-                            " dca vs. phi",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
-  fRecDcaPosPhi->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecDcaNegPhi   = new TH2F("fRecDcaNegPhi", 
-                            " dca vs. phi",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
-  fRecDcaNegPhi->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
-  fRecPtPosPhi   = new TH2F("fRecPtPosPhi", 
-                            " log(p_T) vs. phi",
-                            100, -2.5, 2., 180, 0, TMath::Pi()*2);
-  fRecPtPosPhi->GetXaxis()->SetTitle("log_{10}(p_{T})");
-  fRecPtPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecPtNegPhi   = new TH2F("fRecPtNegPhi", 
-                            " log(p_T) vs. phi",
-                            100,-2.5 , 2., 180, 0, TMath::Pi()*2);
-  fRecPtNegPhi->GetXaxis()->SetTitle("log_{10}(p_{T})");
-  fRecPtNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
-  fRecEtaPosPhi   = new TH2F("fRecEtaPosPhi", 
-                            " eta vs. phi",
-                            100, -1.5, 1.5, 180, 0, TMath::Pi()*2);
-  fRecEtaPosPhi->GetXaxis()->SetTitle("#eta");
-  fRecEtaPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecEtaNegPhi   = new TH2F("fRecEtaNegPhi", 
-                            " eta vs. phi",
-                            100, -1.5, 1.5, 180, 0, TMath::Pi()*2);
-  fRecEtaNegPhi->GetXaxis()->SetTitle("#eta");
-  fRecEtaNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
-  fRecQPtPhi   = new TH2F("fRecQPtPhi", 
-                            " charge/p_T vs. phi",
-                            100,-1. , 1., 180, 0, TMath::Pi()*2);
-  fRecQPtPhi->GetXaxis()->SetTitle("charge/p_{T}");
-  fRecQPtPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
-  fRecEtaPtPosPhi   = new TH2F("fRecEtaPtPosPhi", 
-                            " eta/p_T vs. phi",
-                            100, -5, 5., 180, 0, TMath::Pi()*2);
-  fRecEtaPtPosPhi->GetXaxis()->SetTitle("#eta/p_{T}");
-  fRecEtaPtPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecEtaPtNegPhi   = new TH2F("fRecEtaPtNegPhi", 
-                            " eta/p_T vs. phi",
-                            100,-5 , 5., 180, 0, TMath::Pi()*2);
-  fRecEtaPtNegPhi->GetXaxis()->SetTitle("#eta/p_{T}");
-  fRecEtaPtNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
-
-
-
-
-  fRecDcaPosPhiEtaPos   = new TH2F("fRecDcaPosPhiEtaPos", 
-                            " dca vs. phi",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
-  fRecDcaPosPhiEtaPos->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaPosPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecDcaNegPhiEtaPos   = new TH2F("fRecDcaNegPhiEtaPos", 
-                            " dca vs. phi",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
-  fRecDcaNegPhiEtaPos->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaNegPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
-
-  fRecPtPosPhiEtaPos   = new TH2F("fRecPtPosPhiEtaPos", 
-                            " log(p_T) vs. phi",
-                            100, -2.5, 2., 180, 0, TMath::Pi()*2);
-  fRecPtPosPhiEtaPos->GetXaxis()->SetTitle("log_{10}(p_{T})");
-  fRecPtPosPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecPtNegPhiEtaPos   = new TH2F("fRecPtNegPhiEtaPos", 
-                            " log(p_T) vs. phi",
-                            100,-2.5 , 2., 180, 0, TMath::Pi()*2);
-  fRecPtNegPhiEtaPos->GetXaxis()->SetTitle("log_{10}(p_{T})");
-  fRecPtNegPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
-
-
-  fRecDcaPosPhiEtaNeg   = new TH2F("fRecDcaPosPhiEtaNeg", 
-                            " dca vs. phi",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
-  fRecDcaPosPhiEtaNeg->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaPosPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecDcaNegPhiEtaNeg   = new TH2F("fRecDcaNegPhiEtaNeg", 
-                            " dca vs. phi",
-                            100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
-  fRecDcaNegPhiEtaNeg->GetXaxis()->SetTitle("dca (cm)");
-  fRecDcaNegPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
-
-  fRecPtPosPhiEtaNeg   = new TH2F("fRecPtPosPhiEtaNeg", 
-                            " log(p_T) vs. phi",
-                            100, -2.5, 2., 180, 0, TMath::Pi()*2);
-  fRecPtPosPhiEtaNeg->GetXaxis()->SetTitle("log_{10}(p_{T})");
-  fRecPtPosPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecPtNegPhiEtaNeg   = new TH2F("fRecPtNegPhiEtaNeg", 
-                            " log(p_T) vs. phi",
-                            100,-2.5 , 2., 180, 0, TMath::Pi()*2);
-  fRecPtNegPhiEtaNeg->GetXaxis()->SetTitle("log_{10}(p_{T})");
-  fRecPtNegPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
-
-
-
-    
-  //  YIELDs ---------------- for TPC sectors
-  for(Int_t sector=0; sector<18;sector++){
-      
-
-    fRecPtTpcSector[sector]   = new TH1F(Form("fRecPtTpcSector%02d",sector), 
-                                        Form("p_{T} distribution: TPC sector %d",
-                                             sector),100, 0., pt);
-    fRecPtTpcSector[sector]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-    fRecEtaTpcSector[sector]   = new TH1F(Form("fRecEtaTpcSector%02d",sector), 
-                                         Form("#eta distribution: TPC sector %d",
-                                              sector),200, -2., 2.);
-    fRecEtaTpcSector[sector]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-     
-
-    fSignedDcaTpcSector[sector]   = new TH1F(Form("fSignedDcaTpcSector%02d",sector), 
-                                            Form("dca distribution: TPC sector %d",
-                                                 sector),200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9) );
-    fSignedDcaTpcSector[sector]->GetXaxis()->SetTitle("dca");
-
-    fRecQPtTpcSector[sector]   = new TH1F(Form("fRecQPtTpcSector%02d",sector), 
-                                         Form("Q/ p_{T} distribution: TPC sector %d",
-                                              sector),100, -1., 1.);
-    fRecQPtTpcSector[sector]->GetXaxis()->SetTitle("Q/p_{T} (GeV/c)");
-
-    fRecEtaPtTpcSector[sector]   = new TH1F(Form("fRecEtaPtTpcSector%02d",sector), 
-                                           Form("#eta/ p_{T} distribution: TPC sector %d",
-                                                sector),100, -1., 1.);
-    fRecEtaPtTpcSector[sector]->GetXaxis()->SetTitle("#eta/p_{T} (GeV/c)");
-  }
-  // YIELDS ITS ladder
-  for(Int_t i=0;i<7;i++){
-    fRecPtPosLadder[i]   = new TH1F(Form("fRecPtPosLadder%d", i), 
-                                   " p_{T} distribution",
-                                   100, 0., pt);
-    fRecPtPosLadder[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-    fRecPtNegLadder[i]   = new TH1F(Form("fRecPtNegLadder%d",i), 
-                                   " p_{T} distribution ",
-                                   100, 0., pt);
-    fRecPtNegLadder[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
-    fRecPhiPosLadder[i]   = new TH1F(Form("fRecPhiPosLadder%d",i), 
-                                    "#phi distribution: all pos eta",
-                                    361, 0., 360);
-    fRecPhiPosLadder[i]->GetXaxis()->SetTitle("#phi (deg)");
-      
-    fRecPhiNegLadder[i]   = new TH1F(Form("fRecPhiNegLadder%d", i),
-                                    "#phi distribution: all neg eta",
-                                    361, 0, 360);
-    fRecPhiNegLadder[i]->GetXaxis()->SetTitle("#phi (deg)");
-
-
-
-    fRecEtaPosLadder[i]   = new TH1F(Form("fRecEtaPosLadder%d",i), 
-                                    "#eta distribution",
-                                    200, -2., 2.);
-    fRecEtaPosLadder[i]->GetXaxis()->SetTitle("#eta)");
-      
-    fRecEtaNegLadder[i]   = new TH1F(Form("fRecEtaNegLadder%d", i),
-                                    "#eta distribution",
-                                    200, -2., 2.);
-    fRecEtaNegLadder[i]->GetXaxis()->SetTitle("#eta");
-  }
-
-  Double_t vzmax = 15;
-
-  fRecPtPosVz = new TH2F("fRecPtPosVz", 
-                        "p_{T} distribution vs Vz()",
-                        100, -1., 2., 200,-vzmax,vzmax);
-  fRecPtPosVz->GetXaxis()->SetTitle("log_{10}(p_{T})");
-    
-  fRecPtNegVz = new TH2F("fRecPtNegVz",
-                        "p_{T} distribution vs Vz()",
-                        100, -1., 2.,200,-vzmax,vzmax);
-  fRecPtNegVz->GetXaxis()->SetTitle("Log_{10}(p_{T})");
-    
-   
-  fRecEtaPosVz= new TH2F("fRecEtaPosVz", 
-                        "#eta distribution vs Vz()",
-                        100, -2., 2., 200,-vzmax,vzmax);
-  fRecEtaPosVz->GetXaxis()->SetTitle("#eta");
-  fRecEtaNegVz = new TH2F("fRecEtaNegVz",
-                         "#eta distribution vs Vz()",
-                         100, -2., 2.,200,-vzmax,vzmax);
-  fRecEtaNegVz->GetXaxis()->SetTitle("#eta");
-
-  fRecPhiPosVz= new TH2F("fRecPhiPosVz", 
-                        "#eta distribution vs Vz()",
-                        361, 0., 360., 200,-vzmax,vzmax);
-  fRecPhiPosVz->GetXaxis()->SetTitle("#phi (deg)");
-  fRecPhiNegVz = new TH2F("fRecPhiNegVz",
-                         "dca vs Vz()",
-                         361, 0., 360.,200,-vzmax,vzmax);
-  fRecPhiNegVz->GetXaxis()->SetTitle("#phi (deg)");
-
-  fSignedDcaPosVz= new TH2F("fSignedDcaPosVz", 
-                           "#eta distribution vs Vz()",
-                           200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 200,-vzmax,vzmax);
-  fSignedDcaPosVz->GetXaxis()->SetTitle("dca (cm)");
-  fSignedDcaNegVz = new TH2F("fSignedDcaNegVz",
-                            "dca vs Vz()",
-                            200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200,-vzmax,vzmax);
-  fSignedDcaNegVz->GetXaxis()->SetTitle("dca (cm)");
-
-  fRecQPtPosEtaVz= new TH2F("fRecQPtPosEtaVz",
-                           " Q/p_{T} distribution vs Vz()",
-                           100, -1., 1., 200,-vzmax,vzmax);
-  fRecQPtPosEtaVz->GetXaxis()->SetTitle("Q/p_{T}");
-  fRecQPtNegEtaVz = new TH2F("fRecQPtNegEtaVz",
-                            " Q/p_{T} distribution vs Vz()",
-                            100, -1., 1.,200,-vzmax,vzmax);
-  fRecQPtNegEtaVz->GetXaxis()->SetTitle("Q/p_{T}");
-
-  fRecEtaPtPosVz= new TH2F("fRecEtaPtPosVz",
-                          " #eta/p_{T} distribution vs Vz()",
-                          100, -1., 1., 200,-vzmax,vzmax);
-  fRecEtaPtPosVz->GetXaxis()->SetTitle("#eta/p_{T");
-  fRecEtaPtNegVz = new TH2F("fRecEtaPtNegVz",
-                           " #eta/p_{T} distribution vs Vz()",
-                           100, -1., 1.,200,-vzmax,vzmax);
-  fRecEtaPtNegVz->GetXaxis()->SetTitle("#eta/p_{T}");
-
-  //new
-  fDeltaPhiAll = new TH1F("fDeltaPhiAll",
-                         " #Delta #phi",200,-360,360);
-  fDeltaPhiAll->GetXaxis()->SetTitle("#Delta #phi");
-
-
-  fDeltaPhiLeading = new TH2F("fDeltaPhiLeading",
-                             " #Delta #phi",361,-360,360, 361,0, 360);
-  fDeltaPhiLeading->GetXaxis()->SetTitle("#Delta #phi (deg.)");
-  fDeltaPhiLeading->GetYaxis()->SetTitle("#phi_{leading particle} (deg.)");
-
-  fDiffDcaD    = new TH1F("fDiffDcaD",   
-                         "dca-d",
-                         200, -1., 1.);
-
-  
-  fRecPtPosEtaPos = new TH1F("fRecPtPosEtaPos",
-                            "p_{T} distribution",100,0,pt);
-  fRecPtPosEtaPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-  fRecPtNegEtaPos = new TH1F("fRecPtNegEtaPos",
-                            "p_{T} distribution",100,0,pt);
-  fRecPtNegEtaPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-  fRecPtPosEtaNeg = new TH1F("fRecPtPosEtaNeg",
-                            "p_{T} distribution",100,0,pt);
-  fRecPtPosEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-  fRecPtNegEtaNeg = new TH1F("fRecPtNegEtaNeg",
-                            "p_{T} distribution",100,0,pt);
-  fRecPtNegEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
-
-  fRec1PtPosEtaPos = new TH1F("fRec1PtPosEtaPos",
-                            "1/p_{T} distribution",100,0,0.5);
-  fRec1PtPosEtaPos->GetXaxis()->SetTitle("p_{T} (c/GeV)");
-
-  fRec1PtNegEtaPos = new TH1F("fRec1PtNegEtaPos",
-                            "1/p_{T} distribution",100,0,0.5);
-  fRec1PtNegEtaPos->GetXaxis()->SetTitle("p_{T} (c/GeV)");
-
-  fRec1PtPosEtaNeg = new TH1F("fRec1PtPosEtaNeg",
-                            "1/p_{T} distribution",100,0,0.5);
-  fRec1PtPosEtaNeg->GetXaxis()->SetTitle("p_{T} (c/GeV)");
-
-  fRec1PtNegEtaNeg = new TH1F("fRec1PtNegEtaNeg",
-                            "1/p_{T} distribution",100,0,0.5);
-  fRec1PtNegEtaNeg->GetXaxis()->SetTitle("1/p_{T} (c/GeV)");
-
-
-  fRecPhiPosEtaPos = new TH1F("fRecPhiPosEtaPos",
-                            "#phi",180,0,2*TMath::Pi());
-  fRecPhiPosEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
-
-  fRecPhiNegEtaPos = new TH1F("fRecPhiNegEtaPos",
-                            "#phi",180,0,2*TMath::Pi());
-  fRecPhiNegEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
-
-  fRecPhiPosEtaNeg = new TH1F("fRecPhiPosEtaNeg",
-                            "#phi",180,0,2*TMath::Pi());
-  fRecPhiPosEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
-
-  fRecPhiNegEtaNeg = new TH1F("fRecPhiNegEtaNeg",
-                            "#phi",180,0,2*TMath::Pi());
-  fRecPhiNegEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
-
-
-
-  fRecDcaPhiPtPosEtaPos = new TH3F("fRecDcaPhiPtPosEtaPos",
-                                  "#phi- p_{T} - DCA",
-                                  180,0,2*TMath::Pi(),
-                                  100,0,pt,
-                                  100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaPhiPtPosEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
-  fRecDcaPhiPtPosEtaPos->GetYaxis()->SetTitle("p_{T} (GeV/c)");
-  fRecDcaPhiPtPosEtaPos->GetZaxis()->SetTitle("dca (cm)");
-
-  fRecDcaPhiPtPosEtaNeg = new TH3F("fRecDcaPhiPtPosEtaNeg",
-                                  "#phi- p_{T} - DCA",
-                                  180,0,2*TMath::Pi(),
-                                  100,0,pt,
-                                  100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaPhiPtPosEtaNeg->GetZaxis()->SetTitle("dca (cm)");
-  fRecDcaPhiPtPosEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
-  fRecDcaPhiPtPosEtaNeg->GetYaxis()->SetTitle("p_{T} (GeV/c)");
-
-  fRecDcaPhiPtNegEtaPos = new TH3F("fRecDcaPhiPtNegEtaPos",
-                                  "#phi- p_{T} - DCA",
-                                  180,0,2*TMath::Pi(),
-                                  100,0,pt,
-                                  100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaPhiPtNegEtaPos->GetZaxis()->SetTitle("dca (cm)");
-  fRecDcaPhiPtNegEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
-  fRecDcaPhiPtNegEtaPos->GetYaxis()->SetTitle("p_{T} (GeV/c)");
-
-  fRecDcaPhiPtNegEtaNeg = new TH3F("fRecDcaPhiPtNegEtaNeg",
-                                  "#phi- p_{T} - DCA",
-                                  180,0,2*TMath::Pi(),
-                                  100,0,pt,
-                                  100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-  fRecDcaPhiPtNegEtaNeg->GetZaxis()->SetTitle("dca (cm)");
-  fRecDcaPhiPtNegEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
-  fRecDcaPhiPtNegEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
-
-
-  fHists->SetOwner();
-
-  fHists->Add(fHistRECpt);
-  fHists->Add(fEta);
-  fHists->Add(fEtavPt);
-  fHists->Add(fEtaPhi);
-  fHists->Add(fThetaRec);
-  fHists->Add(fPhiRec);
-  fHists->Add(fNumber);
-  fHists->Add(fVx);
-  fHists->Add(fVy);
-  fHists->Add(fVz);
-
-  fHists->Add(fEtaPt);
-  fHists->Add(fQPt);
-  fHists->Add(fDca);
-
-  fHists->Add(fDeltaPhiAll);
-  fHists->Add(fDeltaPhiLeading);
-  fHists->Add(fDiffDcaD);
-
-  fHists->Add(fqRec);
-  fHists->Add(fsigmaPt);
-
-  fHists->Add(fRecPtPos);
-  fHists->Add(fRecPtNeg);
-  fHists->Add(fRecPhiPos);
-  fHists->Add(fRecPhiNeg);
-  fHists->Add(fRecEtaPos);
-  fHists->Add(fRecEtaNeg);
-  fHists->Add(fRecEtaPtPos);
-  fHists->Add(fRecEtaPtNeg);
-  fHists->Add(fRecDcaPos);
-  fHists->Add(fRecDcaNeg);
-  fHists->Add(fRecDcaNegInv);
-  fHists->Add(fRecDPos);
-  fHists->Add(fRecDNeg);
-
-
-  fHists->Add(fRecQPtPosEta);
-  fHists->Add(fRecQPtNegEta);
-  fHists->Add(fRecPtPosEta);
-  fHists->Add(fRecPtNegEta);
-  fHists->Add(fRecPhiPosEta);
-  fHists->Add(fRecPhiNegEta);
-  fHists->Add(fRecDcaPosEta);
-  fHists->Add(fRecDcaNegEta);
-  fHists->Add(fRecDPosEta);
-  fHists->Add(fRecDNegEta);
-
-
-  for(Int_t i=0;i<18;i++){
-    fHists->Add(fRecPtTpcSector[i]);
-    fHists->Add(fRecEtaTpcSector[i]);
-    fHists->Add(fSignedDcaTpcSector[i]);
-    fHists->Add(fRecQPtTpcSector[i]);
-    fHists->Add(fRecEtaPtTpcSector[i]);
-  }
-
-  for(Int_t i=0;i<7;i++){
-    fHists->Add(fRecPtPosLadder[i]);
-    fHists->Add(fRecPtNegLadder[i]);
-    fHists->Add(fRecPhiPosLadder[i]);
-    fHists->Add(fRecPhiNegLadder[i]);
-    fHists->Add(fRecEtaPosLadder[i]);
-    fHists->Add(fRecEtaNegLadder[i]);
-  }
-
-  fHists->Add(fRecPtPosVz);
-  fHists->Add(fRecPtNegVz);
-  fHists->Add(fRecEtaPosVz);
-  fHists->Add(fRecEtaNegVz);
-  fHists->Add(fRecPhiPosVz);
-  fHists->Add(fRecPhiNegVz);
-  fHists->Add(fSignedDcaPosVz);
-  fHists->Add(fSignedDcaNegVz);
-  fHists->Add(fRecQPtPosEtaVz);
-  fHists->Add(fRecQPtNegEtaVz);
-  fHists->Add(fRecEtaPtPosVz);
-  fHists->Add(fRecEtaPtNegVz);
-
-
-  for(Int_t i=0;i<7;i++){
-    fHists->Add(fSignDcaPos[i]);
-    fHists->Add(fSignDcaNeg[i]);
-    fHists->Add(fSignDcaNegInv[i]);
-    fHists->Add(fPtSigmaPos[i]);
-    fHists->Add(fPtSigmaNeg[i]);
-    fHists->Add(fqPtRec[i]);
-  
-    fHists->Add(fDcaSigmaPos[i]);
-    fHists->Add(fDcaSigmaNeg[i]);
-
-  } 
-  
-  fHists->Add(fRecDcaPosPhi);
-  fHists->Add(fRecDcaNegPhi);   
-  fHists->Add(fRecPtPosPhi);
-  fHists->Add(fRecPtNegPhi);   
-  fHists->Add(fRecEtaPosPhi);
-  fHists->Add(fRecEtaNegPhi);  
-  fHists->Add(fRecQPtPhi);   
-  fHists->Add(fRecEtaPtPosPhi);
-  fHists->Add(fRecEtaPtNegPhi);   
-
-  fHists->Add(fRecPtPosEtaPos);   
-  fHists->Add(fRecPtNegEtaPos);   
-  fHists->Add(fRecPtPosEtaNeg);   
-  fHists->Add(fRecPtNegEtaNeg); 
-
-  fHists->Add(fRec1PtPosEtaPos);   
-  fHists->Add(fRec1PtNegEtaPos);   
-  fHists->Add(fRec1PtPosEtaNeg);   
-  fHists->Add(fRec1PtNegEtaNeg);   
-
-  fHists->Add(fRecPhiPosEtaPos);   
-  fHists->Add(fRecPhiNegEtaPos);   
-  fHists->Add(fRecPhiPosEtaNeg);   
-  fHists->Add(fRecPhiNegEtaNeg);   
-
-  fHists->Add(fRecDcaPosPhiEtaPos);
-  fHists->Add(fRecDcaNegPhiEtaPos);   
-  fHists->Add(fRecPtPosPhiEtaPos);
-  fHists->Add(fRecPtNegPhiEtaPos);  
-  fHists->Add(fRecDcaPosPhiEtaNeg);
-  fHists->Add(fRecDcaNegPhiEtaNeg);   
-  fHists->Add(fRecPtPosPhiEtaNeg);
-  fHists->Add(fRecPtNegPhiEtaNeg); 
-
-  fHists->Add(fRecDcaPhiPtPosEtaPos); 
-  fHists->Add(fRecDcaPhiPtPosEtaNeg); 
-  fHists->Add(fRecDcaPhiPtNegEtaPos); 
-  fHists->Add(fRecDcaPhiPtNegEtaNeg); 
-
-
-
-
-    
-//   for (Int_t i=0; i<fHists->GetEntries(); ++i) {
-//     TH1 *h1 = dynamic_cast<TH1*>(fHists->At(i));
-//     if (h1){
-//     //  Printf("%s ",h1->GetName());
-//       h1->Sumw2();
-//     }
-//   }
-  // BKC
-
-  TH1::AddDirectory(oldStatus);
-}
-
-//__________________________________________________________
-
-void AliAnalysisTaskQASym::UserExec(Option_t *) 
-{
-  AliVEvent *event = InputEvent();
-  if (!event) {
-    Printf("ERROR: Could not retrieve event");
-    return;
-  }
-
-
-  if(Entry()==0){
-    AliESDEvent* esd = dynamic_cast<AliESDEvent*>(event);
-    if(esd){
-      Printf("We are reading from ESD");
-    }
-   
-  }
-
-
-
-  if(fDebug>1)Printf("There are %d tracks in this event", event->GetNumberOfTracks());
-
-  
-  Int_t   leadingTrack  =   0;
-  Float_t leadingEnergy = -20.;
-  Float_t leadingPhi    =   0;//TMath::Pi();
-
-
-  if(event->GetNumberOfTracks()!=0) fNumber->Fill(event->GetNumberOfTracks());
-
-  const AliVVertex* vertex = event->GetPrimaryVertex();
-  Float_t vz = vertex->GetZ();
-  if (TMath::Abs(vz) > 10.) return;
-
-  for (Int_t iTrack = 0; iTrack < event->GetNumberOfTracks(); iTrack++) {
-    
-
-    AliVParticle *track = event->GetTrack(iTrack);
-    AliESDtrack *esdtrack =  dynamic_cast<AliESDtrack*>(track);
-    esdtrack->PropagateToDCA(event->GetPrimaryVertex(),
-                            event->GetMagneticField(), 10000.);
-
-    if (!track) {
-      Printf("ERROR: Could not receive track %d", iTrack);
-      continue;
-    }
-    //__________
-    // run Task for global tracks or ITS tracks or TPC tracks
-    AliESDtrack *tpcP = 0x0;
-    const AliExternalTrackParam *tpcPin = 0x0;
-    Double_t phiIn=0;
-    if(fTrackType==0){
-      //Fill all histograms with global tracks
-      tpcP = esdtrack;
-      phiIn = tpcP->Phi();
-      if (!tpcP) continue;
-      if (!fCuts->AcceptTrack(tpcP)) continue;
-    }
-    else if(fTrackType==1){
-      //Fill all histograms with ITS tracks
-      tpcP = esdtrack;
-      phiIn = tpcP->Phi();
-      if (!tpcP) continue;
-      if (!fCuts->AcceptTrack(tpcP)) continue;
-    }
-    else if(fTrackType==2){     
-      //Fill all histograms with TPC track information
-      tpcPin = esdtrack->GetInnerParam();
-      if (!tpcPin) continue;
-      phiIn=tpcPin->Phi();
-
-      tpcP = AliESDtrackCuts::GetTPCOnlyTrack(dynamic_cast<AliESDEvent*>(event),esdtrack->GetID());
-      if (!tpcP) continue;
-      if (!fCuts->AcceptTrack(tpcP)) continue;
-      if(tpcP->GetNcls(1)>160)continue;//jacek's track cut
-      if(tpcP->GetConstrainedChi2TPC()<0)continue; // jacek's track cut
-    }
-    else{
-      Printf("ERROR: wrong track type \n");
-      continue;
-    }
-    //___________
-    //
-  
-    if(tpcP->E()>leadingEnergy){
-      leadingTrack=iTrack;
-      leadingEnergy=tpcP->E();
-      leadingPhi=phiIn;
-    }
-   
-
-    fqRec->Fill(tpcP->Charge());
-  
-
-    Double_t sigmapt = tpcP->GetSigma1Pt2();
-    sigmapt= sqrt(sigmapt);
-    sigmapt= sigmapt *(tpcP->Pt()*tpcP->Pt()); 
-
-    if(sigmapt == 0.)continue;
-    fsigmaPt->Fill(TMath::Log10(sigmapt));
-
-    // hits in ITS layer
-    Int_t cas=-1;
-    if(tpcP->HasPointOnITSLayer(0)) 
-      cas=0;
-    else if(!tpcP->HasPointOnITSLayer(0)
-           &&  tpcP->HasPointOnITSLayer(1)) 
-      cas=1;
-    else if(!tpcP->HasPointOnITSLayer(0)
-           && !tpcP->HasPointOnITSLayer(1) 
-           &&  tpcP->HasPointOnITSLayer(2)) 
-      cas=2;
-    else if(!tpcP->HasPointOnITSLayer(0)
-           && !tpcP->HasPointOnITSLayer(1) 
-           && !tpcP->HasPointOnITSLayer(2)
-           &&  tpcP->HasPointOnITSLayer(3)) 
-      cas=3;
-    else if(!tpcP->HasPointOnITSLayer(0)
-           && !tpcP->HasPointOnITSLayer(1) 
-           && !tpcP->HasPointOnITSLayer(2)
-           && !tpcP->HasPointOnITSLayer(3)
-           &&  tpcP->HasPointOnITSLayer(4)) 
-      cas=4;
-    else if(   !tpcP->HasPointOnITSLayer(0)
-              && !tpcP->HasPointOnITSLayer(1)
-              && !tpcP->HasPointOnITSLayer(2)
-              && !tpcP->HasPointOnITSLayer(3)
-              && !tpcP->HasPointOnITSLayer(4) 
-              &&  tpcP->HasPointOnITSLayer(5)) 
-      cas=5;
-    else 
-      cas=6;
-  
-   
-   
-    //------------------- 
-
-    xvertexcor = tpcP->Xv() - vertex->GetX(); // coordinate corrected for vertex position
-    yvertexcor = tpcP->Yv() - vertex->GetY(); // "
-    sdca = (tpcP->Py()*xvertexcor - tpcP->Px()*yvertexcor)/tpcP->Pt();
-
-
-    fqPtRec[cas]->Fill(tpcP->Charge()/tpcP->Pt());
-    
-    
-
-    fHistRECpt->Fill(tpcP->Pt());
-    fEta->Fill(tpcP->Eta());
-    fEtavPt->Fill(tpcP->Eta(), tpcP->Pt());
-    fEtaPhi->Fill(tpcP->Eta(), phiIn);
-    fThetaRec->Fill(tpcP->Theta());
-    fPhiRec->Fill(phiIn);
-    fVx->Fill(tpcP->Xv());
-    fVy->Fill(tpcP->Yv());
-    fVz->Fill(tpcP->Zv());
-  
-
-    fEtaPt->Fill(tpcP->Eta()/tpcP->Pt());
-    fQPt->Fill(tpcP->Charge()/tpcP->Pt());
-    fDca->Fill(sdca);
-    fRecQPtPhi->Fill(tpcP->Charge()/tpcP->Pt(), phiIn);
-
-
-    tpcP->GetImpactParameters(xy,z);
-    fDiffDcaD->Fill(sdca+xy);
-
-   
-    
-    //for positive particles
-
-    if(tpcP->Charge()>0){
-      fRecPtPos->Fill(tpcP->Pt());
-      fRecPtPosLadder[cas]->Fill(tpcP->Pt());
-      fRecPtPosVz->Fill(TMath::Log10(tpcP->Pt()),tpcP->Zv());
-      fRecPhiPos->Fill(TMath::RadToDeg()*phiIn);
-    
-     
-      fRecPhiPosLadder[cas]->Fill(TMath::RadToDeg()*phiIn);
-      fRecPhiPosVz->Fill(TMath::RadToDeg()*phiIn,tpcP->Zv());
-      fSignedDcaPosVz->Fill(sdca,tpcP->Zv());
-
-      fRecEtaPos->Fill(tpcP->Eta());
-      fRecEtaPosLadder[cas]->Fill(tpcP->Eta());
-      fRecEtaPtPos->Fill(tpcP->Eta()/tpcP->Pt());
-      fRecEtaPosVz->Fill(tpcP->Eta(),tpcP->Zv());
-      fRecEtaPtPosVz->Fill(tpcP->Eta()/tpcP->Pt(),tpcP->Zv());
-     
-      fRecDcaPos->Fill(sdca);
-      fRecDcaPosPhi->Fill(sdca, phiIn);
-      fRecPtPosPhi->Fill(TMath::Log10(tpcP->Pt()), phiIn);
-      fRecEtaPtPosPhi->Fill(tpcP->Eta()/tpcP->Pt(), phiIn);
-      fRecEtaPosPhi->Fill(tpcP->Eta(), phiIn);
-      fRecDPos->Fill(xy);
-      fSignDcaPos[cas]->Fill(sdca);
-    
-     
-      fDcaSigmaPos[cas]->Fill(sdca, TMath::Log10(sigmapt));
-    
-      fPtSigmaPos[cas]->Fill(TMath::Log10(sigmapt));
-      //pos eta
-      if(tpcP->Eta()>0){
-       fRecPtPosEtaPos->Fill(tpcP->Pt());
-       fRec1PtPosEtaPos->Fill(1/tpcP->Pt());
-       fRecPhiPosEtaPos->Fill(phiIn);
-       fRecDcaPosPhiEtaPos->Fill(sdca, phiIn);
-       fRecPtPosPhiEtaPos->Fill(TMath::Log10(tpcP->Pt()), phiIn);
-       fRecDcaPhiPtPosEtaPos->Fill(phiIn, tpcP->Pt(), sdca);
-      }
-      //neg eta
-      else{
-       fRecPtPosEtaNeg->Fill(tpcP->Pt());
-       fRec1PtPosEtaNeg->Fill(1/tpcP->Pt());
-       fRecPhiPosEtaNeg->Fill(phiIn);
-       fRecDcaPosPhiEtaNeg->Fill(sdca, phiIn);
-       fRecPtPosPhiEtaNeg->Fill(TMath::Log10(tpcP->Pt()), phiIn);
-       fRecDcaPhiPtPosEtaNeg->Fill(phiIn, tpcP->Pt(), sdca);
-      }
-      
-    }
-    //and negative particles
-    else {
-      fRecPtNeg->Fill(tpcP->Pt());
-      fRecPtNegLadder[cas]->Fill(tpcP->Pt());
-      fRecPtNegVz->Fill(TMath::Log10(tpcP->Pt()),tpcP->Zv());
-           
-      fRecPhiNeg->Fill(TMath::RadToDeg()*phiIn);
-      fRecPhiNegLadder[cas]->Fill(TMath::RadToDeg()*phiIn);
-      fRecPhiNegVz->Fill(TMath::RadToDeg()*phiIn,tpcP->Zv());
-      fSignedDcaNegVz->Fill(sdca,tpcP->Zv());
-      fRecEtaPtNegVz->Fill(tpcP->Eta()/tpcP->Pt(),tpcP->Zv());
-
-      fRecEtaNeg->Fill(tpcP->Eta());
-      fRecEtaNegLadder[cas]->Fill(tpcP->Eta());
-      fRecEtaPtNeg->Fill(tpcP->Eta()/tpcP->Pt());
-      fRecEtaNegVz->Fill(tpcP->Eta(),tpcP->Zv());
-     
-      fRecDcaNeg->Fill(sdca);
-      fRecDcaNegInv->Fill(-sdca);
-      fRecDcaNegPhi->Fill(sdca, phiIn);
-      fRecPtNegPhi->Fill(TMath::Log10(tpcP->Pt()), phiIn);
-      fRecEtaNegPhi->Fill(tpcP->Eta(), phiIn);
-      fRecEtaPtNegPhi->Fill(tpcP->Eta()/tpcP->Pt(), phiIn);
-      fRecDNeg->Fill(xy);
-      fSignDcaNeg[cas]->Fill(sdca);
-      fSignDcaNegInv[cas]->Fill(-sdca);
-     
-     
-      fDcaSigmaNeg[cas]->Fill(sdca,TMath::Log10(sigmapt));
-   
-      fPtSigmaNeg[cas]->Fill(TMath::Log10(sigmapt));
-      
-      //pos eta
-      if(tpcP->Eta()>0){
-       fRecPtNegEtaPos->Fill(tpcP->Pt());
-       fRec1PtNegEtaPos->Fill(1/tpcP->Pt());
-       fRecPhiNegEtaPos->Fill(phiIn);
-       fRecDcaNegPhiEtaPos->Fill(sdca, phiIn);
-       fRecPtNegPhiEtaPos->Fill(TMath::Log10(tpcP->Pt()), phiIn);
-       fRecDcaPhiPtNegEtaPos->Fill(phiIn, tpcP->Pt(), sdca);
-      }
-      //neg eta
-      else{
-       fRecPtNegEtaNeg->Fill(tpcP->Pt());
-       fRec1PtNegEtaNeg->Fill(1/tpcP->Pt());
-       fRecPhiNegEtaNeg->Fill(phiIn);
-       fRecDcaNegPhiEtaNeg->Fill(sdca, phiIn);
-       fRecPtNegPhiEtaNeg->Fill(TMath::Log10(tpcP->Pt()), phiIn);
-       fRecDcaPhiPtNegEtaNeg->Fill(phiIn, tpcP->Pt(), sdca);
-      }
-
-    }
-    
-
-
-    //all particles with positive eta
-    if(tpcP->Eta()>0){
-      fRecQPtPosEta->Fill(tpcP->Charge()/tpcP->Pt());
-      fRecPtPosEta->Fill(tpcP->Pt());
-      fRecPhiPosEta->Fill(TMath::RadToDeg()*phiIn);
-      fRecQPtPosEtaVz->Fill(tpcP->Charge()/tpcP->Pt(),tpcP->Zv());
-      fRecDcaPosEta->Fill(sdca);
-      fRecDPosEta->Fill(xy);
-    }
-    //all particles with negative eta (and eta==0)
-    else{
-      fRecQPtNegEta->Fill(tpcP->Charge()/tpcP->Pt());
-      fRecPtNegEta->Fill(tpcP->Pt());
-      fRecPhiNegEta->Fill(TMath::RadToDeg()*phiIn);
-      fRecQPtNegEtaVz->Fill(tpcP->Charge()/tpcP->Pt(),tpcP->Zv());
-      fRecDcaNegEta->Fill(sdca);
-      fRecDNegEta->Fill(xy);
-
-    }
-
-
-    fRecPtTpcSector[Int_t(phiIn*
-                         TMath::RadToDeg()/20)]->Fill(tpcP->Pt());
-    fRecEtaTpcSector[Int_t(phiIn*
-                          TMath::RadToDeg()/20)]->Fill(tpcP->Eta());
-    fSignedDcaTpcSector[Int_t(phiIn*
-                             TMath::RadToDeg()/20)]->Fill(sdca); 
-    fRecQPtTpcSector[Int_t(phiIn*
-                          TMath::RadToDeg()/20)]->Fill(tpcP->Charge()/tpcP->Pt());
-    fRecEtaPtTpcSector[Int_t(phiIn*
-                            TMath::RadToDeg()/20)]->Fill(tpcP->Eta()/tpcP->Pt());
-     
-
-
-//     // another track loop
-//     for (Int_t iTrack2 = 0; iTrack2 < event->GetNumberOfTracks(); iTrack2++) {
-      
-//       if(LeadingTrack==iTrack2) continue;
-
-//       AliVParticle *track2 = event->GetTrack(iTrack2);
-//       AliESDtrack* esdtrack2 =  dynamic_cast<AliESDtrack*>(track2);
-//       if (!track2) {
-//     Printf("ERROR: Could not receive track %d", iTrack);
-//     continue;
-//       }
-//       if (!fCuts->AcceptTrack(esdtrack2)) continue;
-//       //propagate to dca
-//       esdtrack2->PropagateToDCA(event->GetPrimaryVertex(),
-//                             event->GetMagneticField(), 10000.);
-//       fDeltaPhiLeading->Fill((LeadingPhi-esdtrack2->Phi())*TMath::RadToDeg(),
-//                          LeadingPhi*TMath::RadToDeg() );
-
-     
-
-//     }//second track loop
-  }//first track loop
-
-  
-
-  // Post output data.
-  // PostData(1, fHistPt);
-  PostData(1, fHists);
-}      
-
-
-
-
-
-//________________________________________________________________________
-void AliAnalysisTaskQASym::Terminate(Option_t *) 
-{
-
-
-}  
-
-
-
-
-
diff --git a/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.h b/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.h
deleted file mode 100644 (file)
index d20b5b0..0000000
+++ /dev/null
@@ -1,192 +0,0 @@
-#ifndef AliAnalysisTaskQASym_cxx
-#define AliAnalysisTaskQASym_cxx
-
-class TH1F;
-class TH2F;
-class TH3F;
-class TList;
-class TNtuple;
-
-class AliESDEvent;
-class AliESDtrack;
-class AliESDtrackCuts;
-
-
-#include "AliAnalysisTaskSE.h"
-#include "TFile.h"
-#include "TNtuple.h"
-
-class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
- public:
-  AliAnalysisTaskQASym(const char *name = "AliAnalysisTaskQASym");
-  virtual ~AliAnalysisTaskQASym() {}
-  
-  virtual void   UserCreateOutputObjects();
-  virtual void   UserExec(Option_t *option);
-  virtual void   Terminate(Option_t *);
-  virtual void   SetTrackType(Int_t type) {fTrackType = type;}  
-  
-  
-  virtual void   SetCuts(AliESDtrackCuts* cuts)
-     {fCuts = cuts;}
-
-  virtual void   SetFieldOn(Bool_t b = kTRUE){fFieldOn = b;} 
-
-  
- private:
-
-  Int_t       fTrackType;
-  Bool_t      fFieldOn;
-
-  TList       *fHists;          // List of histos
-
-  //old
-  TH1F        *fHistRECpt;      // pt 
-  TH1F        *fEta;            // eta
-  TH2F        *fEtaPhi;         // eta-phi
-  TH1F        *fEtaPt;          // eta over pt 
-  TH1F        *fQPt;            // charge over pt 
-  TH1F        *fDca;            // distance of closest approach
-  TH1F        *fqRec;           // reconstrcuted charge
-  TH1F        *fsigmaPt;        // sigma_pT
-
-   //positive und negative tracks
-  TH1F        *fRecPtPos;      // pt of pos tracks
-  TH1F        *fRecPtNeg;      // pt of neg tracks
-  TH1F        *fRecPhiPos;     // phi of pos. tracks
-  TH1F        *fRecPhiNeg;     // phi of neg. tracks
-  TH1F        *fRecEtaPos;     // eta of neg. tracks
-  TH1F        *fRecEtaNeg;     // eta of neg. tracks
-  TH1F        *fRecEtaPtPos;   // eta over pt of neg. tracks
-  TH1F        *fRecEtaPtNeg;   // eta over pt of neg. tracks
-  TH1F        *fRecDcaPos;     // distance of closest approach of neg. tracks
-  TH1F        *fRecDcaNeg;     // distance of closest of neg. tracks
-  TH1F        *fRecDcaNegInv;  // invers dca of neg. tracks
-  TH1F        *fRecDPos;       // impact parameter of neg. tracks
-  TH1F        *fRecDNeg;       // impact parameter of neg. tracks
-
-  // two sides of TPC -> Eta/Theta
-  TH1F        *fRecQPtPosEta;   //charge/pT for pos. eta
-  TH1F        *fRecQPtNegEta;   //charge/pT for neg. eta
-  TH1F        *fRecPtPosEta;    //pT        for pos. eta
-  TH1F        *fRecPtNegEta;    //pT        for neg. eta
-  TH1F        *fRecPhiPosEta;   //phi       for pos. eta
-  TH1F        *fRecPhiNegEta;   //phi       for neg. eta
-  TH1F        *fRecDcaPosEta;   //dca       for pos. eta
-  TH1F        *fRecDcaNegEta;   //dca       for neg. eta
-  TH1F        *fRecDPosEta;     //d         for pos. eta
-  TH1F        *fRecDNegEta;     //d         for neg. eta
-
-  // 2D: all measures as function of z of first trackpoint 
-  TH2F        *fRecPtPosVz;     //pt-zfirst of pos tracks
-  TH2F        *fRecPtNegVz;      //pt-zfirst of neg tracks
-  TH2F        *fRecEtaPosVz;    //eta-zfirst of pos tracks
-  TH2F        *fRecEtaNegVz;    //eta-zfirst of neg tracks
-  TH2F        *fRecPhiPosVz;    //phi-zfirst of pos tracks
-  TH2F        *fRecPhiNegVz;    //phi-zfirst of neg tracks
-  TH2F        *fSignedDcaPosVz; //dca-zfirst of pos tracks
-  TH2F        *fSignedDcaNegVz; //dca-zfirst of neg tracks
-  TH2F        *fRecQPtPosEtaVz; //charge/pT-zfirst of pos tracks
-  TH2F        *fRecQPtNegEtaVz; //charge/pT-zfirst of neg tracks
-  TH2F        *fRecEtaPtPosVz;  //eta/pT-zfirst of pos tracks
-  TH2F        *fRecEtaPtNegVz;  //eta/pT-zfirst of neg tracks
-
-
-  //high
-  TH1F * fDeltaPhiAll;         // phiLeaingTracks-phiOthers
-  TH2F * fDeltaPhiLeading;     // phiLeaingTracks-phiOthers vs. phiLeading
-  TH1F * fDiffDcaD;            // d-dca
-
-  //sim
-  TH1F * fPhiRec;              //phi
-  TH1F * fThetaRec;            //theta
-  TH1F * fNumber;              //Number of tracks per event
-  TH1F * fVx;                  // x of first track point
-  TH1F * fVy;                  // y of first track point
-  TH1F * fVz;                  // z of first track point
-  TNtuple * test;
-
-  //new
-  TH2F        *fRecDcaPosPhi;     //dca-phi for pos.
-  TH2F        *fRecDcaNegPhi;     //dca-phi for neg.
-  TH2F        *fRecPtPosPhi;      //pt-phi for pos.
-  TH2F        *fRecPtNegPhi;      //pt-phi for neg.
-  TH2F        *fRecEtaPosPhi;     //eta-phi for pos.
-  TH2F        *fRecEtaNegPhi;     //eta-phi for neg.
-  TH2F        *fRecQPtPhi;        //charge/pt-phi
-  TH2F        *fRecEtaPtPosPhi;   //eta/pt-phi for neg.
-  TH2F        *fRecEtaPtNegPhi;   //eta/pt-phi for pos.
-
-  TH1F        *fRecPtPosEtaPos;   //pt for pos tracks and pos eta
-  TH1F        *fRecPtNegEtaPos;   //pt for neg tracks and pos eta
-  TH1F        *fRecPtPosEtaNeg;   //pt for pos tracks and neg eta
-  TH1F        *fRecPtNegEtaNeg;   //pt for neg tracks and neg eta
-
-  TH1F        *fRec1PtPosEtaPos;   //1/pt for pos tracks and pos eta
-  TH1F        *fRec1PtNegEtaPos;   //1/pt for neg tracks and pos eta
-  TH1F        *fRec1PtPosEtaNeg;   //1/pt for pos tracks and neg eta
-  TH1F        *fRec1PtNegEtaNeg;   //1/pt for neg tracks and neg eta
-
-  TH1F        *fRecPhiPosEtaPos;   //phi for pos tracks and pos eta
-  TH1F        *fRecPhiNegEtaPos;   //phi for neg tracks and pos eta
-  TH1F        *fRecPhiPosEtaNeg;   //phi for pos tracks and neg eta
-  TH1F        *fRecPhiNegEtaNeg;   //phi for neg tracks and neg eta
-
-  TH2F        *fRecDcaPosPhiEtaPos;  //dca-phi for pos tracks and pos eta
-  TH2F        *fRecDcaNegPhiEtaPos;  //dca-phi for neg tracks and pos eta
-  TH2F        *fRecDcaPosPhiEtaNeg;  //dca-phi for pos tracks and neg eta
-  TH2F        *fRecDcaNegPhiEtaNeg;  //dca-phi for neg tracks and neg eta
-
-  TH2F        *fRecPtPosPhiEtaPos;  //pt-phi for pos tracks and pos eta
-  TH2F        *fRecPtNegPhiEtaPos;  //pt-phi for neg tracks and pos eta 
-  TH2F        *fRecPtPosPhiEtaNeg;  //pt-phi for pos tracks and neg eta
-  TH2F        *fRecPtNegPhiEtaNeg;  //pt-phi for neg tracks and neg eta
-
-
-  TH3F        *fRecDcaPhiPtPosEtaPos; //dca-pt-phi for pos tracks and pos eta
-  TH3F        *fRecDcaPhiPtNegEtaPos; //dca-pt-phi for neg tracks and pos eta
-  TH3F        *fRecDcaPhiPtPosEtaNeg; //dca-pt-phi for pos tracks and neg eta
-  TH3F        *fRecDcaPhiPtNegEtaNeg; //dca-pt-phi for neg tracks and neg eta
-
-  TH2F        *fEtavPt;
-
-  Double_t  sdca;
-  Float_t xy, z, xvertexcor, yvertexcor;
-  AliESDtrackCuts* fCuts;                      // List of cuts
-
-  // sectors of TPC 
-  TH1F        *fRecPtTpcSector[18];     //pt for TPC sectors
-  TH1F        *fRecEtaTpcSector[18];    //eta for TPC sectors
-  TH1F        *fRecQPtTpcSector[18];    //charge/pt for TPC sectors
-  TH1F        *fRecEtaPtTpcSector[18];  //eta/pt for TPC sectors
-  TH1F        *fSignedDcaTpcSector[18]; //dca for TPC sectors
-
-
-  // 7 different case of hit in ITS ladders
-  TH1F        *fRecPtPosLadder[7];  //pt for pos tracks
-  TH1F        *fRecPtNegLadder[7];  //pt for neg tracks
-  TH1F        *fRecPhiPosLadder[7]; //phi for pos tracks
-  TH1F        *fRecPhiNegLadder[7]; //phi for neg tracks
-  TH1F        *fRecEtaPosLadder[7]; //eta for pos tracks
-  TH1F        *fRecEtaNegLadder[7]; //eta for neg tracks
-  TH1F        *fSignDcaPos[7];      //dca for pos tracks
-  TH1F        *fSignDcaNeg[7];      //dca for neg tracks
-  TH1F        *fSignDcaNegInv[7];   //-dca for neg tracks
-  TH1F        *fPtSigmaPos[7];      //sigma_pT for pos tracks
-  TH1F        *fPtSigmaNeg[7];      //sigma_pT for neg tracks
-  TH1F        *fqPtRec[7];          // charge/pt 
-  TH2F        *fDcaSigmaPos[7];     // dca - sigma_pT for pos tracks
-  TH2F        *fDcaSigmaNeg[7];     // dca - sigma_pT for neg tracks
-
-  
-  
-
-
-  AliAnalysisTaskQASym(const AliAnalysisTaskQASym&); // not implemented
-  AliAnalysisTaskQASym& operator=(const AliAnalysisTaskQASym&); // not implemented
-  
-  ClassDef(AliAnalysisTaskQASym, 1); // Basic QA exploiting symmetries
-};
-
-#endif
diff --git a/prod/acrcaf/qa_pp/plot_macros/plot.C b/prod/acrcaf/qa_pp/plot_macros/plot.C
deleted file mode 100644 (file)
index eaf3caa..0000000
+++ /dev/null
@@ -1,78 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[7];
-TH1F* fHneg[7];
-
-
-fHpos[0]=  (TH1F*) l->FindObject("fRecPtPos");
-fHneg[0]=  (TH1F*) l->FindObject("fRecPtNeg");
-fHpos[1]=  (TH1F*) l->FindObject("fRecPhiPos");
-fHneg[1]=  (TH1F*) l->FindObject("fRecPhiNeg");
-fHpos[2]=  (TH1F*) l->FindObject("fRecEtaPos");
-fHneg[2]=  (TH1F*) l->FindObject("fRecEtaNeg");
-fHpos[3]=  (TH1F*) l->FindObject("fRecEtaPtPos");
-fHneg[3]=  (TH1F*) l->FindObject("fRecEtaPtNeg");
-
-fHpos[4]=  (TH1F*) l->FindObject("fRecPtPosEta");
-fHneg[4]=  (TH1F*) l->FindObject("fRecPtNegEta");
-fHpos[5]=  (TH1F*) l->FindObject("fRecPhiPosEta");
-fHneg[5]=  (TH1F*) l->FindObject("fRecPhiNegEta");
-fHpos[6]=  (TH1F*) l->FindObject("fRecQPtPosEta");
-fHneg[6]=  (TH1F*) l->FindObject("fRecQPtNegEta");
-
-
-TLegend *legp1;
-legp1= new TLegend(0.9,0.65,0.65,0.9);
-legp1->SetFillColor(kWhite);
-
-TLegend *legp2;
-legp2= new TLegend(0.9,0.65,0.65,0.9);
-legp2->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp1->AddEntry(fun0,"pos. charge","l");
-legp1->AddEntry(fun1,"neg. charge","l");
-
-legp2->AddEntry(fun0,"pos. #eta","l");
-legp2->AddEntry(fun1,"neg. #eta","l");
-
-
-
-for(Int_t i=0; i<7;i++){
-  fHpos[i]->SetLineColor(kBlack);
-  fHpos[i]->SetLineWidth(2);
-  fHneg[i]->SetLineColor(kRed);
-  fHneg[i]->SetLineWidth(2);
-}
-
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(3,3);
-for(Int_t i=0; i<7;i++){
-  dca->cd(i+1);
-  fHpos[i]->Draw();
-  fHneg[i]->Draw("same");
-  if(i==0||i==4)gPad->SetLogy();
-  if (i < 4) 
-    legp1->Draw();
-  else
-    legp2->Draw();
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotD.C b/prod/acrcaf/qa_pp/plot_macros/plotD.C
deleted file mode 100644 (file)
index 4a94302..0000000
+++ /dev/null
@@ -1,73 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecDcaPos");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecDcaNeg");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TH1F* fHposD;
-TH1F* fHnegD;
-
-fHposD=  (TH1F*) l->FindObject("fRecDPos");
-fHposD->SetLineWidth(2);
-
-fHnegD=  (TH1F*) l->FindObject("fRecDNeg");
-fHnegD->SetLineColor(kRed);
-fHnegD->SetLineWidth(2);
-
-
-fH=  (TH2F*) l->FindObject("fDiffDcaD");
-
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,2);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-gPad->SetLogy();
-legp->Draw();
-
-dca->cd(2);
-fHposD->DrawClone();
-fHnegD->DrawClone("same");
-legp->Draw();
-gPad->SetLogy();
-
-
-dca->cd(3);
-fH->Draw("colz");
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotDca.C b/prod/acrcaf/qa_pp/plot_macros/plotDca.C
deleted file mode 100644 (file)
index e5fcba2..0000000
+++ /dev/null
@@ -1,61 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHneg;
-TH1F* fHnegInv;
-
-fHpos=  (TH1F*) l->FindObject("fRecDcaPos");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecDcaNeg");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-fHnegInv=  (TH1F*) l->FindObject("fRecDcaNegInv");
-fHnegInv->SetLineColor(kRed);
-fHnegInv->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHnegInv);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive}/Inverse Yield_{negative}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotDcaVz.C b/prod/acrcaf/qa_pp/plot_macros/plotDcaVz.C
deleted file mode 100644 (file)
index e0bd0cf..0000000
+++ /dev/null
@@ -1,42 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH2F* fHpos;
-TH2F* fHpneg;
-
-fHpos=  (TH2F*) l->FindObject("fSignedDcaPosVz");
-fHneg=  (TH2F*) l->FindObject("fSignedDcaNegVz");
-
-
-
-
-
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 520);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->Draw("colz");
-fHpos->SetTitle("positive tracks");
-fHpos->SetYTitle("esdtrack->Zv()");
-
-dca->cd(2);
-fHneg->Draw("colz");
-fHneg->SetMaximum(fHpos->GetMaximum());
-fHneg->SetTitle("negative tracks");
-fHneg->SetYTitle("esdtrack->Zv()");
-
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotDeltaPhi.C b/prod/acrcaf/qa_pp/plot_macros/plotDeltaPhi.C
deleted file mode 100644 (file)
index e3c9bed..0000000
+++ /dev/null
@@ -1,47 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH2F* fHpos;
-fHpos=  (TH2F*) l->FindObject("fDeltaPhiLeading");
-
-TH1D* px;
-TH1D* py;
-
-//pgood= fHpos->ProjectionY();
-//pgood->Draw();
-
-TCanvas * c1 = new TCanvas("c1", "c1", 100, 100, 1020, 480);
-c1->Divide(2,1);
-// c1->cd(1);
-// fHpos->Draw("colz");
-
-
-c1->cd(1);
-px= fHpos->ProjectionX();
-px->Draw();
-px->SetTitle("");
-
-c1->cd(2);
-py= fHpos->ProjectionY();
-py->Draw();
-py->SetMinimum(0);
-py->SetTitle("");
-
-TCanvas * c2 = new TCanvas("c2", "c2", 100, 100, 620, 480);
-c2->cd();
-gPad->SetLeftMargin(0.13);
-fHpos->SetTitleOffset(1.2, "Y");
-fHpos->Draw("colz");
-fHpos->SetTitle("");
-
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotEta-over-pt.C b/prod/acrcaf/qa_pp/plot_macros/plotEta-over-pt.C
deleted file mode 100644 (file)
index 64eabc9..0000000
+++ /dev/null
@@ -1,56 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecEtaPtPos");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecEtaPtNeg");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-//gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHneg);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive}/Yield_{negative}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotEta.C b/prod/acrcaf/qa_pp/plot_macros/plotEta.C
deleted file mode 100644 (file)
index 8a67c8b..0000000
+++ /dev/null
@@ -1,57 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecEtaPos");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecEtaNeg");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("Hallo");
-fHneg->SetTitle("Hallo");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHneg);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive}/Yield_{negative}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotEtaPtVz.C b/prod/acrcaf/qa_pp/plot_macros/plotEtaPtVz.C
deleted file mode 100644 (file)
index 4870656..0000000
+++ /dev/null
@@ -1,47 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH2F* fHpos;
-TH2F* fHpneg;
-
-fHpos=  (TH2F*) l->FindObject("fRecEtaPtPosVz");
-fHneg=  (TH2F*) l->FindObject("fRecEtaPtNegVz");
-
-
-
-
-
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 520);
-dca->Divide(2,1);
-dca->cd(1);
-gPad->SetRightMargin(0.15);
-fHpos->Draw("colz");
-fHpos->SetTitle("positive tracks");
-fHpos->SetYTitle("esdtrack->Zv()");
-fHpos->SetXTitle("#eta/p_{T}");
-
-
-dca->cd(2);
-gPad->SetRightMargin(0.15);
-fHneg->Draw("colz");
-fHneg->SetMaximum(fHpos->GetMaximum());
-fHneg->SetTitle("negative tracks");
-fHneg->SetYTitle("esdtrack->Zv()");
-fHneg->SetXTitle("#eta/p_{T}");
-
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotEtaVz.C b/prod/acrcaf/qa_pp/plot_macros/plotEtaVz.C
deleted file mode 100644 (file)
index 116fd8b..0000000
+++ /dev/null
@@ -1,42 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecEtaPosVz");
-fHneg=  (TH1F*) l->FindObject("fRecEtaNegVz");
-
-
-
-
-
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 520);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->Draw("colz");
-fHpos->SetTitle("positive tracks");
-fHpos->SetYTitle("esdtrack->Zv()");
-
-dca->cd(2);
-fHneg->Draw("colz");
-fHneg->SetMaximum(fHpos->GetMaximum());
-fHneg->SetTitle("negative tracks");
-fHneg->SetYTitle("esdtrack->Zv()");
-
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotLadderDca.C b/prod/acrcaf/qa_pp/plot_macros/plotLadderDca.C
deleted file mode 100644 (file)
index efa6845..0000000
+++ /dev/null
@@ -1,81 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[7];
-TH1F* fHneg[7];
-
-for(Int_t sector=0;sector<7;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fSignDcaPos%d", sector));
- fHneg[sector]=  (TH1F*) l->FindObject(Form("fSignDcaNegInv%d", sector));
-} 
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<7;i++){
-  fHpos[i]->SetLineColor(kBlack);
-  fHpos[i]->SetLineWidth(2);
-  fHneg[i]->SetLineColor(kRed);
-  fHneg[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 920, 720);
-dca->Divide(3,3);
-for(Int_t i=0; i<7;i++){
-  //fHpos[i]->Divide(fHneg[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
- //  if(i!=1){
-//     fHpos[i]->SetMaximum(fHpos[0]->GetMaximum()*2.5);
-//     fHpos[i]->SetMinimum(fHpos[0]->GetMinimum()-1);
-//   }
-  fHpos[i]->SetMaximum(500);
-  fHpos[i]->SetMinimum(1);
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("p_{T} case %2d",i));
-
-//   gPad->SetLeftMargin(0.15);
-//   gPad->SetBottomMargin(0.15);
-  fHneg[i]->Draw("same");
-  legp->Draw();
-  //gPad->SetGridy();
-  //gPad->SetGridx();
-//   fHpos[i]->SetTitleOffset(1., "X");
-//   fHpos[i]->SetTitleOffset(1., "Y");
-//   fHpos[i]->SetTitleSize(0.08, "X");
-//   fHpos[i]->SetTitleSize(0.05, "Y");
-//   fHpos[i]->SetLabelSize(0.1, "X");
-//   fHpos[i]->SetLabelSize(0.1, "Y");
-//   //fHpos[i]->GetXaxis()->SetNdivision(5);
-//   fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotLadderEta.C b/prod/acrcaf/qa_pp/plot_macros/plotLadderEta.C
deleted file mode 100644 (file)
index 3340bfe..0000000
+++ /dev/null
@@ -1,76 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[7];
-TH1F* fHneg[7];
-
-for(Int_t sector=0;sector<7;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecEtaPosLadder%d", sector));
- fHneg[sector]=  (TH1F*) l->FindObject(Form("fRecEtaNegLadder%d", sector));
-} 
-
-TLegend *legp;
-legp= new TLegend(0.9,0.78,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<7;i++){
-  fHpos[i]->SetLineColor(kBlack);
-  fHpos[i]->SetLineWidth(2);
-  fHneg[i]->SetLineColor(kRed);
-  fHneg[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 920, 720);
-dca->Divide(3,3);
-for(Int_t i=0; i<7;i++){
-  //  fHpos[i]->Divide(fHpos[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  fHpos[i]->SetMaximum(100);
-  fHpos[i]->SetMinimum(0.1);
-
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("#eta for case %02d",i));
-
-//   gPad->SetLeftMargin(0.15);
-//   gPad->SetBottomMargin(0.15);
-  fHneg[i]->Draw("same");
-
-//   fHpos[i]->SetTitleOffset(1., "X");
-//   fHpos[i]->SetTitleOffset(1., "Y");
-//   fHpos[i]->SetTitleSize(0.08, "X");
-//   fHpos[i]->SetTitleSize(0.05, "Y");
-//   fHpos[i]->SetLabelSize(0.1, "X");
-//   fHpos[i]->SetLabelSize(0.1, "Y");
-//   //fHpos[i]->GetXaxis()->SetNdivision(5);
-//   fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-  legp->Draw();
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotLadderPhi.C b/prod/acrcaf/qa_pp/plot_macros/plotLadderPhi.C
deleted file mode 100644 (file)
index a8eca3f..0000000
+++ /dev/null
@@ -1,77 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[6];
-TH1F* fHneg[6];
-
-for(Int_t sector=0;sector<6;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecPhiPosLadder%d", sector));
- fHneg[sector]=  (TH1F*) l->FindObject(Form("fRecPhiNegLadder%d", sector));
-} 
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<6;i++){
-  fHpos[i]->SetLineColor(kBlack);
-  fHpos[i]->SetLineWidth(2);
-  fHneg[i]->SetLineColor(kRed);
-  fHneg[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 920, 720);
-dca->Divide(3,2);
-for(Int_t i=0; i<6;i++){
-  //  fHpos[i]->Divide(fHpos[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  fHpos[i]->SetMaximum(10000);
-  fHpos[i]->SetMinimum(0.1);
-
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("#phi of case %d",i));
-
-//   gPad->SetLeftMargin(0.15);
-//   gPad->SetBottomMargin(0.15);
-  fHneg[i]->Draw("same");
-
-//   fHpos[i]->SetTitleOffset(1., "X");
-//   fHpos[i]->SetTitleOffset(1., "Y");
-//   fHpos[i]->SetTitleSize(0.08, "X");
-//   fHpos[i]->SetTitleSize(0.05, "Y");
-//   fHpos[i]->SetLabelSize(0.1, "X");
-//   fHpos[i]->SetLabelSize(0.1, "Y");
-//   //fHpos[i]->GetXaxis()->SetNdivision(5);
-//   fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotLadderPt.C b/prod/acrcaf/qa_pp/plot_macros/plotLadderPt.C
deleted file mode 100644 (file)
index 9a44c3c..0000000
+++ /dev/null
@@ -1,79 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[7];
-TH1F* fHneg[7];
-
-for(Int_t sector=0;sector<7;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecPtPosLadder%d", sector));
- fHneg[sector]=  (TH1F*) l->FindObject(Form("fRecPtNegLadder%d", sector));
-} 
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<7;i++){
-  fHpos[i]->SetLineColor(kBlack);
-  fHpos[i]->SetLineWidth(2);
-  fHneg[i]->SetLineColor(kRed);
-  fHneg[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 920, 720);
-dca->Divide(3,3);
-for(Int_t i=0; i<7;i++){
-  //fHpos[i]->Divide(fHpos[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
- //  if(i!=1){
-//     fHpos[i]->SetMaximum(fHpos[0]->GetMaximum()*2.5);
-//     fHpos[i]->SetMinimum(fHpos[0]->GetMinimum()-1);
-//   }
-  fHpos[i]->SetMaximum(1000000);
-  fHpos[i]->SetMinimum(0.1);
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("p_{T} case %2d",i));
-
-//   gPad->SetLeftMargin(0.15);
-//   gPad->SetBottomMargin(0.15);
-  fHneg[i]->Draw("same");
-  legp->Draw();
-//   fHpos[i]->SetTitleOffset(1., "X");
-//   fHpos[i]->SetTitleOffset(1., "Y");
-//   fHpos[i]->SetTitleSize(0.08, "X");
-//   fHpos[i]->SetTitleSize(0.05, "Y");
-//   fHpos[i]->SetLabelSize(0.1, "X");
-//   fHpos[i]->SetLabelSize(0.1, "Y");
-//   //fHpos[i]->GetXaxis()->SetNdivision(5);
-//   fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotLadderSigma.C b/prod/acrcaf/qa_pp/plot_macros/plotLadderSigma.C
deleted file mode 100644 (file)
index 38cded4..0000000
+++ /dev/null
@@ -1,80 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[7];
-TH1F* fHneg[7];
-
-for(Int_t sector=0;sector<7;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fPtSigmaPos%d", sector));
- fHneg[sector]=  (TH1F*) l->FindObject(Form("fPtSigmaNeg%d", sector));
-} 
-
-TLegend *legp;
-legp= new TLegend(0.9,0.75,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<7;i++){
-  fHpos[i]->SetLineColor(kBlack);
-  fHpos[i]->SetLineWidth(2);
-  fHneg[i]->SetLineColor(kRed);
-  fHneg[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 920, 720);
-dca->Divide(3,3);
-for(Int_t i=0; i<7;i++){
-  //fHpos[i]->Divide(fHneg[i]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
- //  if(i!=1){
-//     fHpos[i]->SetMaximum(fHpos[0]->GetMaximum()*2.5);
-//     fHpos[i]->SetMinimum(fHpos[0]->GetMinimum()-1);
-//   }
-  fHpos[i]->SetMaximum(1000000);
-  fHpos[i]->SetMinimum(0.1);
-  // fHpos[i]->GetXaxis()->SetRangeUser(-2.5,0.5);
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("#sigma_{p_{T}} case %2d",i));
-
-//   gPad->SetLeftMargin(0.15);
-//   gPad->SetBottomMargin(0.15);
-  fHneg[i]->Draw("same");
-  legp->Draw();
-//   fHpos[i]->SetTitleOffset(1., "X");
-//   fHpos[i]->SetTitleOffset(1., "Y");
-//   fHpos[i]->SetTitleSize(0.08, "X");
-//   fHpos[i]->SetTitleSize(0.05, "Y");
-//   fHpos[i]->SetLabelSize(0.1, "X");
-//   fHpos[i]->SetLabelSize(0.1, "Y");
-//   //fHpos[i]->GetXaxis()->SetNdivision(5);
-//   fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotPhi.C b/prod/acrcaf/qa_pp/plot_macros/plotPhi.C
deleted file mode 100644 (file)
index bec8014..0000000
+++ /dev/null
@@ -1,56 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecPhiPos");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecPhiNeg");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-//gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHneg);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive}/Yield_{negative}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotPhiVz.C b/prod/acrcaf/qa_pp/plot_macros/plotPhiVz.C
deleted file mode 100644 (file)
index 6f89453..0000000
+++ /dev/null
@@ -1,42 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH2F* fHpos;
-TH2F* fHpneg;
-
-fHpos=  (TH2F*) l->FindObject("fRecPhiPosVz");
-fHneg=  (TH2F*) l->FindObject("fRecPhiNegVz");
-
-
-
-
-
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 520);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->Draw("colz");
-fHpos->SetTitle("positive tracks");
-fHpos->SetYTitle("esdtrack->Zv()");
-
-dca->cd(2);
-fHneg->Draw("colz");
-fHneg->SetMaximum(fHpos->GetMaximum());
-fHneg->SetTitle("negative tracks");
-fHneg->SetYTitle("esdtrack->Zv()");
-
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotPhi_Eta.C b/prod/acrcaf/qa_pp/plot_macros/plotPhi_Eta.C
deleted file mode 100644 (file)
index d1bb7c5..0000000
+++ /dev/null
@@ -1,56 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecPhiPosEta");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecPhiNegEta");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive #eta","l");
-legp->AddEntry(fun1,"negative #eta","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-//gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHneg);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive #eta}/Yield_{negative #eta}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotPt.C b/prod/acrcaf/qa_pp/plot_macros/plotPt.C
deleted file mode 100644 (file)
index 989091c..0000000
+++ /dev/null
@@ -1,56 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecPtPos");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecPtNeg");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHneg);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive}/Yield_{negative}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotPtVz.C b/prod/acrcaf/qa_pp/plot_macros/plotPtVz.C
deleted file mode 100644 (file)
index 6a5d9b4..0000000
+++ /dev/null
@@ -1,42 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH2F* fHpos;
-TH2F* fHpneg;
-
-fHpos=  (TH2F*) l->FindObject("fRecPtPosVz");
-fHneg=  (TH2F*) l->FindObject("fRecPtNegVz");
-
-
-
-
-
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 520);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->Draw("colz");
-fHpos->SetTitle("positive tracks");
-fHpos->SetYTitle("esdtrack->Zv()");
-
-dca->cd(2);
-fHneg->Draw("colz");
-fHneg->SetMaximum(fHpos->GetMaximum());
-fHneg->SetTitle("negative tracks");
-fHneg->SetYTitle("esdtrack->Zv()");
-
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotPt_Eta.C b/prod/acrcaf/qa_pp/plot_macros/plotPt_Eta.C
deleted file mode 100644 (file)
index a7f3b7f..0000000
+++ /dev/null
@@ -1,56 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecPtPosEta");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecPtNegEta");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive #eta","l");
-legp->AddEntry(fun1,"negative #eta","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHneg);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive #eta}/Yield_{negative #eta}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotQPtVz.C b/prod/acrcaf/qa_pp/plot_macros/plotQPtVz.C
deleted file mode 100644 (file)
index ca7dc63..0000000
+++ /dev/null
@@ -1,52 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH2F* fHpos;
-TH2F* fHpneg;
-
-fHpos=  (TH2F*) l->FindObject("fRecQPtPosEtaVz");
-fHneg=  (TH2F*) l->FindObject("fRecQPtNegEtaVz");
-
-
-
-
-
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 520);
-dca->Divide(2,1);
-dca->cd(1);
-//fHpos->Divide(fHneg);
-gPad->SetRightMargin(0.15);
-gPad->SetLeftMargin(0.15);
-fHpos->SetTitleOffset(1.8,"Y");
-//fHpos->SetXTitle("(Charge/p_{T})_{pos} / (Charge/p_{T})_{neg}");
-fHpos->SetXTitle("Charge/p_{T}");
-fHpos->Draw("colz");
-fHpos->SetTitle("positive #eta");
-fHpos->SetYTitle("esdtrack->Zv()");
-
-dca->cd(2);
-fHneg->Draw("colz");
-gPad->SetRightMargin(0.15);
-gPad->SetLeftMargin(0.15);
-fHneg->SetTitleOffset(1.8,"Y");
-fHneg->SetXTitle("Charge/p_{T}");
-fHneg->SetMaximum(fHpos->GetMaximum());
-fHneg->SetTitle("negative #eta");
-fHneg->SetYTitle("esdtrack->Zv()");
-
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotQPt_Eta.C b/prod/acrcaf/qa_pp/plot_macros/plotQPt_Eta.C
deleted file mode 100644 (file)
index cb9d8e8..0000000
+++ /dev/null
@@ -1,56 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos;
-TH1F* fHpneg;
-
-fHpos=  (TH1F*) l->FindObject("fRecQPtPosEta");
-fHpos->SetLineWidth(2);
-
-fHneg=  (TH1F*) l->FindObject("fRecQPtNegEta");
-fHneg->SetLineColor(kRed);
-fHneg->SetLineWidth(2);
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive #eta","l");
-legp->AddEntry(fun1,"negative #eta","l");
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
-dca->Divide(2,1);
-dca->cd(1);
-fHpos->DrawClone("");
-fHpos->SetTitle("");
-fHpos->SetYTitle("");
-fHneg->DrawClone("same");
-gPad->SetLogy();
-legp->Draw();
-dca->cd(2);
-fHpos->Divide(fHneg);
-fHpos->Draw("");
-fHpos->SetLineColor(kBlue);
-fHpos->SetTitle("Yield_{positive #eta}/Yield_{negative #eta}");
-fHpos->SetYTitle("");
-
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorDca.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorDca.C
deleted file mode 100644 (file)
index 37ca2de..0000000
+++ /dev/null
@@ -1,76 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
-  fHpos[sector]=  (TH1F*) l->FindObject(Form("fSignedDcaTpcSector%02d", sector));
-
-}
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  //  fHpos[i]->Divide(fHpos[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
- //  if(i!=1){
-//   fHpos[i]->SetMaximum(5);
-//   fHpos[i]->SetMinimum(-3);
-  fHpos[i]->SetMaximum(100);
-  fHpos[i]->SetMinimum(1);
-//   }
-  // fHpos[i]->GetXaxis()->SetRangeUser(-.5,.5);
-  fHpos[i]->SetTitle("");
-  //  fHpos[i]->SetXTitle(Form("dca%d / dca0",i));
-  fHpos[i]->SetXTitle(Form("dca%d",i));
-
-  gPad->SetLeftMargin(0.15);
-  gPad->SetBottomMargin(0.15);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(1., "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  //fHpos[i]->GetXaxis()->SetNdivision(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorDcaRatio.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorDcaRatio.C
deleted file mode 100644 (file)
index ff3b955..0000000
+++ /dev/null
@@ -1,77 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
-  fHpos[sector]=  (TH1F*) l->FindObject(Form("fSignedDcaTpcSector%02d", sector));
-}
-TH1F* fHposref =(TH1F*) l->FindObject("fDca");
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->Divide(fHposref);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
- //  if(i!=1){
-  fHpos[i]->SetMaximum(0.3);
-  fHpos[i]->SetMinimum(-0.1);
-//   fHpos[i]->SetMaximum(100000);
-//   fHpos[i]->SetMinimum(1)
-    ;
-//   }
-  fHpos[i]->GetXaxis()->SetRangeUser(-1,1);
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("dca_{%d} / dca_{all}",i));
-  //fHpos[i]->SetXTitle(Form("dca%d",i));
-
-  gPad->SetLeftMargin(0.20);
-  gPad->SetBottomMargin(0.15);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(1., "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  fHpos[i]->GetYaxis()->SetNdivisions(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-  //  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorEta.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorEta.C
deleted file mode 100644 (file)
index d4ab688..0000000
+++ /dev/null
@@ -1,78 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
-  fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecEtaTpcSector%02d", sector));
-
-}
-//TH1F* fHposref =(TH1F*) l->FindObject("")
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  // fHpos[i]->Divide(fHpos[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-
-  fHpos[i]->SetMaximum(5000);
-  fHpos[i]->SetMinimum(0.1);
-//   fHpos[i]->SetMaximum(3);
-//   fHpos[i]->SetMinimum(-1);
-
-  fHpos[i]->GetXaxis()->SetRangeUser(-1.5,1.5);
-  fHpos[i]->SetTitle("");
-  //fHpos[i]->SetXTitle(Form("#eta_{%d} / #eta_{0} ",i));
-  fHpos[i]->SetXTitle(Form("#eta TPC sector %d ",i));
-
-  gPad->SetLeftMargin(0.15);
-  gPad->SetBottomMargin(0.2);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(1., "Y");
-  fHpos[i]->SetTitleSize(0.1, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  //fHpos[i]->GetXaxis()->SetNdivision(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-  fHpos[i]->GetYaxis()->SetNdivisions(7);
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPt.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPt.C
deleted file mode 100644 (file)
index b08fe24..0000000
+++ /dev/null
@@ -1,74 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecEtaPtTpcSector%02d", sector));
-} 
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  //fHpos[i]->Divide(fHpos[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  // fHpos[i]->GetXaxis()->SetRangeUser(-2,2);
-  fHpos[i]->SetMaximum(1400);
-  fHpos[i]->SetMinimum(200);
-//  fHpos[i]->SetMaximum(5);
-//   fHpos[i]->SetMinimum(-2);
-  fHpos[i]->SetTitle("");
-  //  fHpos[i]->SetXTitle(Form("(#eta/p_{T})_{%d} / (#eta/p_{T})_{0} ",i));
-  fHpos[i]->SetXTitle(Form("(#eta/p_{T})_{%d} ",i));
-  gPad->SetLeftMargin(0.21);
-  gPad->SetBottomMargin(0.17);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(.8, "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  //gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPtRatio.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPtRatio.C
deleted file mode 100644 (file)
index 4fd0dd5..0000000
+++ /dev/null
@@ -1,77 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecEtaPtTpcSector%02d", sector));
-} 
-
-TH1F* fHposref =(TH1F*) l->FindObject("fEtaPt");
-
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->Divide(fHposref);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  // fHpos[i]->GetXaxis()->SetRangeUser(-2,2);
-//   fHpos[i]->SetMaximum(2000);
-//   fHpos[i]->SetMinimum(1);
- fHpos[i]->SetMaximum(0.2);
-  fHpos[i]->SetMinimum(-0.05);
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("(#eta/p_{T})_{%d} / (#eta/p_{T})_{all} ",i));
-  //fHpos[i]->SetXTitle(Form("(#eta/p_{T})_{%d} ",i));
-  gPad->SetLeftMargin(0.20);
-  gPad->SetBottomMargin(0.17);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(.8, "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  //fHpos[i]->GetXaxis()->SetNdivision(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  // gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaRatio.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaRatio.C
deleted file mode 100644 (file)
index 604be1c..0000000
+++ /dev/null
@@ -1,78 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
-  fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecEtaTpcSector%02d", sector));
-
-}
-TH1F* fHposref =(TH1F*) l->FindObject("fEta");
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->Divide(fHposref);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-
-  // fHpos[i]->SetMaximum(1000);
-  //fHpos[i]->SetMinimum(0.1);
-  fHpos[i]->SetMaximum(0.3);
-  fHpos[i]->SetMinimum(-0.1);
-
-  // fHpos[i]->GetXaxis()->SetRangeUser(-1.5,1.5);
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("#eta_{%d} / #eta_{all} ",i));
-  //fHpos[i]->SetXTitle(Form("#eta TPC sector %d ",i));
-
-  gPad->SetLeftMargin(0.15);
-  gPad->SetBottomMargin(0.2);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(1., "Y");
-  fHpos[i]->SetTitleSize(0.1, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  //fHpos[i]->GetXaxis()->SetNdivision(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-  fHpos[i]->GetYaxis()->SetNdivisions(7);
-
-  // gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorPt.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorPt.C
deleted file mode 100644 (file)
index c0a8fad..0000000
+++ /dev/null
@@ -1,78 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecPtTpcSector%02d", sector));
-} 
-
-TH1F* fHposref =  (TH1F*) l->FindObject("fHistRECpt");
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  // fHpos[i]->Divide(fHposref);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  fHpos[i]->SetMaximum(100);
-  fHpos[i]->SetMinimum(1);
-  // fHpos[i]->GetXaxis()->SetRangeUser(0,3);
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("p_{T} sector %2d",i));
-  //fHpos[i]->SetXTitle(Form("p_{T,%d}/p_{T,all}",i));
-
-  gPad->SetLeftMargin(0.2);
-  gPad->SetBottomMargin(0.17);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(.8, "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  //fHpos[i]->GetXaxis()->SetNdivision(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-  fHpos[i]->GetYaxis()->SetNdivisions(5);
-
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorPtRatio.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorPtRatio.C
deleted file mode 100644 (file)
index 7c61ad9..0000000
+++ /dev/null
@@ -1,78 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecPtTpcSector%02d", sector));
-} 
-
-TH1F* fHposref =  (TH1F*) l->FindObject("fHistRECpt");
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->Divide(fHposref);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  fHpos[i]->SetMaximum(0.1);
-  fHpos[i]->SetMinimum(0);
-  fHpos[i]->GetXaxis()->SetRangeUser(0,5);
-  fHpos[i]->SetTitle("");
-  //fHpos[i]->SetXTitle(Form("p_{T} sector %2d",i));
-  fHpos[i]->SetXTitle(Form("p_{T,%d}/p_{T,all}",i));
-
-  gPad->SetLeftMargin(0.2);
-  gPad->SetBottomMargin(0.17);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(.8, "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  //fHpos[i]->GetXaxis()->SetNdivision(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-  fHpos[i]->GetYaxis()->SetNdivisions(5);
-
-
-  //gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorQPt.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorQPt.C
deleted file mode 100644 (file)
index e9b231b..0000000
+++ /dev/null
@@ -1,75 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecQPtTpcSector%02d", sector));
-} 
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  // fHpos[i]->Divide(fHpos[0]);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  // fHpos[i]->SetMaximum(2);
-//   fHpos[i]->SetMinimum(0.);
-  fHpos[i]->SetMaximum(1000);
-  fHpos[i]->SetMinimum(0.1);
-  fHpos[i]->SetTitle("");
-  //fHpos[i]->SetXTitle(Form("(Q/p_{T})_{%d}/(Q/p_{T})_{0}",i));
-  fHpos[i]->SetXTitle(Form("Q/p_{T} sector %2d",i));
-
-  gPad->SetLeftMargin(0.15);
-  gPad->SetBottomMargin(0.17);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(.8, "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  //fHpos[i]->GetXaxis()->SetNdivision(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/plotSectorQPtRatio.C b/prod/acrcaf/qa_pp/plot_macros/plotSectorQPtRatio.C
deleted file mode 100644 (file)
index 1c8b823..0000000
+++ /dev/null
@@ -1,78 +0,0 @@
-{
-gROOT->SetStyle("Plain");
-gStyle->SetOptStat(0);
-gStyle->SetPalette(1);
-
-//open files and get lists
-TList *l;
-TFile f0("QAsym.proof.root");
-l = (TList *)f0.Get("QAsymHists");
-
-
-TH1F* fHpos[18];
-
-for(Int_t sector=0;sector<18;sector++){
- fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecQPtTpcSector%02d", sector));
-} 
-
-TH1F* fHposref =(TH1F*) l->FindObject("fQPt");
-
-TLegend *legp;
-legp= new TLegend(0.9,0.65,0.65,0.9);
-legp->SetFillColor(kWhite);
-
-TF1 *fun0, *fun1;
-fun0= new TF1("fun0","gaus",-5.0,5.0);
-fun1= new TF1("fun1","gaus",-5.0,5.0);
-fun0->SetLineColor(kBlack);
-fun1->SetLineColor(kRed);
-
-legp->AddEntry(fun0,"positive","l");
-legp->AddEntry(fun1,"negative","l");
-
-
-
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->SetLineColor(kRed);
-  // fHpos[i]->SetLineWidth(2);
-}
-      
-
-TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
-dca->Divide(6,3);
-for(Int_t i=0; i<18;i++){
-  fHpos[i]->Divide(fHposref);
-  dca->cd(i+1);
-
-  fHpos[i]->Draw();
-  fHpos[i]->SetMaximum(0.3);
-  fHpos[i]->SetMinimum(-0.1);
-//   fHpos[i]->SetMaximum(100);
-//   fHpos[i]->SetMinimum(1)
-    ;
-  fHpos[i]->SetTitle("");
-  fHpos[i]->SetXTitle(Form("(Q/p_{T})_{%d}/(Q/p_{T})_{all}",i));
-  //fHpos[i]->SetXTitle(Form("Q/p_{T} sector %2d",i));
-
-  gPad->SetLeftMargin(0.20);
-  gPad->SetBottomMargin(0.17);
-
-  fHpos[i]->SetTitleOffset(1., "X");
-  fHpos[i]->SetTitleOffset(.8, "Y");
-  fHpos[i]->SetTitleSize(0.08, "X");
-  fHpos[i]->SetTitleSize(0.05, "Y");
-  fHpos[i]->SetLabelSize(0.1, "X");
-  fHpos[i]->SetLabelSize(0.1, "Y");
-  fHpos[i]->GetYaxis()->SetNdivisions(5);
-  fHpos[i]->GetXaxis()->SetNdivisions(5);
-
-
-  //gPad->SetLogy();
-
-}
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/plot_macros/pos-NegInv.C b/prod/acrcaf/qa_pp/plot_macros/pos-NegInv.C
deleted file mode 100644 (file)
index d92089b..0000000
+++ /dev/null
@@ -1,74 +0,0 @@
-//void posNegInv(TString fname)
-{
-
-  gROOT->SetStyle("Plain");
-  gStyle->SetOptStat(0);
-  gStyle->SetPalette(1);
-  
-  //open files and get lists
-  TList *l0;
-  // TFile f0(fname);
-  TFile f0("QAsym.proof.root");
-  l0 = (TList *)f0.Get("QAsymHists");
-
-  // TList *l1;
-  // TFile f1("QAsym.proof.fp.all.ITSrefit.root");
-  // l1 = (TList *)f1.Get("QAsymHists");
-
-  TH1F* fHpos[6];
-  TH1F* fHneg[6];
-
-  for(int i=0;i<6;i++){
-    fHpos[i]=  (TH1F*) l0->FindObject(Form("fSignDcaPos%d", i));
-    fHpos[i]->SetLineColor(kBlack);
-    fHpos[i]->SetLineWidth(2);
-    fHneg[i]=  (TH1F*) l0->FindObject(Form("fSignDcaNegInv%d", i));
-    fHneg[i]->SetLineColor(kRed);
-    fHneg[i]->SetLineWidth(2);
-  }
-
-
-     
-  TLegend *legp[6];
-  for(int i=0;i<6;i++){
-    legp[i]= new TLegend(0.9,0.65,0.55,0.9);
-    legp[i]->SetFillColor(kWhite);
-    TF1 *fun[2];
-    fun[0]= new TF1(Form("fun%d", i),"gaus",-5.0,5.0);
-    fun[0]->SetLineColor(kBlack);
-    fun[1]= new TF1(Form("fun%d", i),"gaus",-5.0,5.0);
-    fun[1]->SetLineColor(kRed);
-    legp[i]->AddEntry(fun[0],"Positive Particles  f( x)","l");  
-    legp[i]->AddEntry(fun[1],"Negative partivcles f(-x)","l");  
-  }
-
-
-
-
-  Double_t range= 4.;
-
-  TCanvas * dca = new TCanvas("dca", "dca", 100, 100, 880, 680);
-  dca->Divide(2,3);
-  for(Int_t i=0;i<6;i++){
-
-    dca->cd(i+1);
-    fHpos[i]->SetTitle(Form("Case %d",i));
-    fHpos[i]->SetMaximum(1000000);
-    fHpos[i]->SetMinimum(0.1);
-    fHpos[i]->Draw();
-    fHneg[i]->Draw("same");
-    gPad->SetLogy();
-    legp[i]->Draw();
-    gPad->SetGridx();
-    gPad->SetGridy();
-    // fHpos[i]->GetXaxis()->SetRangeUser(-range,range);
-  }
-
-
-}
-  
-
-
diff --git a/prod/acrcaf/qa_pp/plot_macros/sigma.C b/prod/acrcaf/qa_pp/plot_macros/sigma.C
deleted file mode 100644 (file)
index 82b4ab8..0000000
+++ /dev/null
@@ -1,95 +0,0 @@
-//void sigma( TString fname)
-{
-  gROOT->SetStyle("Plain");
-  gStyle->SetOptStat(0);
-  gStyle->SetPalette(1);
-  
-  //open files and get lists
-  TList *l;
-  
-  //TFile f0(fname);
-  TFile f0("QAsym.proof.root");
-  l = (TList *)f0.Get("QAsymHists");
-  
-  
-  TH1F* fHpos[7];
-  TH1F* fHneg[7];
-  
-  for(int i=0;i<7;i++){
-    fHpos[i]=  (TH1F*) l->FindObject(Form("fPtSigmaPos%d", i));
-    fHpos[i]->SetLineColor(i+1);
-    fHpos[i]->SetLineWidth(2);
-    fHpos[i]->GetXaxis()->SetLabelSize(0.06);
-    fHpos[i]->GetYaxis()->SetLabelSize(0.06);
-    fHpos[i]->GetXaxis()->SetTitleSize(0.06);
-    fHpos[i]->GetYaxis()->SetTitleSize(0.06);
-  
-    fHneg[i]=  (TH1F*) l->FindObject(Form("fPtSigmaNeg%d", i));
-    fHneg[i]->SetLineColor(i+1);
-    fHneg[i]->SetLineWidth(2);
-    fHneg[i]->GetXaxis()->SetLabelSize(0.06);
-    fHneg[i]->GetYaxis()->SetLabelSize(0.06);
-    fHneg[i]->GetXaxis()->SetTitleSize(0.06);
-    fHneg[i]->GetYaxis()->SetTitleSize(0.06);
-  }
-      
-  TLegend *legp;
-  legp= new TLegend(0.9,0.45,0.83,0.9);
-  legp->SetFillColor(kWhite);
-  TF1 *fun[7];
-  for(int i=0;i<7;i++){
-    fun[i]= new TF1(Form("fun%d", i),"gaus",-5.0,5.0);
-    fun[i]->SetLineColor(i+1);
-    legp->AddEntry(fun[i],Form("case %d",i),"l");  
-  }
-
-
-
-
-
-
-  TCanvas * dca = new TCanvas("sigma", "sigma", 100, 100, 1020, 720);
-  dca->Divide(1,2);
-
-  for(int i=0;i<7;i++){
-
-    dca->cd(1);
-    gPad->SetBottomMargin(0.18);
-    // fHpos[i]->GetXaxis()->SetRangeUser(-3.,3.);
-    fHpos[i]->SetTitle("positive particles");
-    fHpos[i]->SetMaximum(1000);
-    fHpos[i]->SetMinimum(0.1);
-    if(i==0)fHpos[i]->Draw();
-    else fHpos[i]->Draw("same");
-    gPad->SetLogy();
-    //gPad->SetLogx();
-    legp->Draw();
-    gPad->SetGridx();
-    gPad->SetGridy();
-    //fHpos[i]->GetXaxis()->SetRangeUser(-4,0);
-
-    dca->cd(2);
-    gPad->SetBottomMargin(0.18);
-    // fHneg[i]->GetXaxis()->SetRangeUser(-3.,3.);
-    fHneg[i]->SetTitle("negative particles");
-    fHneg[i]->SetMaximum(1000);
-    fHneg[i]->SetMinimum(0.1);
-    if(i==0)fHneg[i]->Draw();
-    else fHneg[i]->Draw("same");
-
-    gPad->SetLogy();
-    //gPad->SetLogx();
-    legp->Draw();
-    gPad->SetGridx();
-    gPad->SetGridy();
-    //fHneg[i]->GetXaxis()->SetRangeUser(-4,0);
-
-
-  }
-  
-
-
-}
diff --git a/prod/acrcaf/qa_pp/qa_pp.C b/prod/acrcaf/qa_pp/qa_pp.C
deleted file mode 100644 (file)
index 488fa64..0000000
+++ /dev/null
@@ -1,69 +0,0 @@
-
-void AddTaskQAsym  (Int_t runNumber);
-void AddTaskVZEROQA(Int_t runNumber);
-
-void qa_pp(Int_t runNumber) {
-  TStopwatch timer;
-  timer.Start();
-
-  gEnv->SetValue("XSec.GSI.DelegProxy","2");
-  /// Select ROOT version
-  TProof::Mgr("proof02@alicecaf:31093")->SetROOTVersion("v5-24-00a");
-  // Login to CAF
-  TProof::Open("proof02@alicecaf:31093");
-
-  // Enable AliRoot
-  gProof->UploadPackage("/afs/cern.ch/alice/caf/sw/ALICE/PARs/v4-18-12-AN/AF-v4-18-12-AN.par");
-  gProof->EnablePackage("AF-v4-18-12-AN.par");
-
-  
-  // Enable analysis libs
-  gSystem->Load("libANALYSIS.so");
-  gSystem->Load("libANALYSISalice.so");
-  gSystem->Load("libPHOSUtils.so");
-  gSystem->Load("libEMCALUtils.so");
-  //  gSystem->Load("libPWG4PartCorrBase.so");
-  //  gSystem->Load("libPWG4PartCorrDep.so");
-
-  gProof->Exec("gSystem->Load(\"libANALYSIS.so\");",        kTRUE);
-  gProof->Exec("gSystem->Load(\"libANALYSISalice.so\");",   kTRUE);
-  gProof->Exec("gSystem->Load(\"libPHOSUtils.so\");",       kTRUE);
-  gProof->Exec("gSystem->Load(\"libEMCALUtils.so\");",      kTRUE);
-  // gProof->Exec("gSystem->Load(\"libPWG4PartCorrBase.so\");",kTRUE);
-  // gProof->Exec("gSystem->Load(\"libPWG4PartCorrDep.so\");", kTRUE);
-
-  gProof->Load("AliAnalysisTaskQASym.cxx++g");
-  gProof->Load("AliAnaVZEROQA.cxx++g");
-
-  gProof->UploadPackage("PWG4PartCorrBase.par");
-  gProof->EnablePackage("PWG4PartCorrBase");
-
-  gProof->UploadPackage("PWG4PartCorrDep.par");
-  gProof->EnablePackage("PWG4PartCorrDep");
-
-  //____________________________________________//
-  // Make the analysis manager
-  AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisQAManager");
-  AliESDInputHandler* esdH = new AliESDInputHandler();
-  mgr->SetInputEventHandler(esdH);  
-  mgr->SetDebugLevel(10);
-
-  // Wagons
-  gROOT->LoadMacro("AddTaskQAsym.C");
-  AddTaskQAsym(runNumber);
-
-  gROOT->LoadMacro("AddTaskVZEROQA.C");
-  AddTaskVZEROQA(runNumber);
-
-  gROOT->LoadMacro("AddTaskCalorimeterQA.C");
-  AliAnalysisTaskParticleCorrelation *taskQAcalo = AddTaskCalorimeterQA("ESD", kFALSE, kTRUE);
-
-  
-  if (!mgr->InitAnalysis()) return;
-  mgr->PrintStatus();
-  mgr->StartAnalysis("proof",
-                    Form("/ALIREC/aliprod/run%d",runNumber));
-
-  timer.Stop();
-  timer.Print();
-}