From 6b53f1d89211c732e327a5d1065048ffd6066fe0 Mon Sep 17 00:00:00 2001 From: morsch Date: Mon, 29 Mar 2010 07:45:06 +0000 Subject: [PATCH] Obsolete --- prod/acrcaf/qa_pp.sh | 14 - prod/acrcaf/qa_pp/AddTaskCalorimeterQA.C | 113 -- prod/acrcaf/qa_pp/AddTaskQAsym.C | 120 -- prod/acrcaf/qa_pp/AddTaskVZEROQA.C | 40 - prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx | 1298 ----------------- prod/acrcaf/qa_pp/AliAnalysisTaskQASym.h | 192 --- prod/acrcaf/qa_pp/plot_macros/plot.C | 78 - prod/acrcaf/qa_pp/plot_macros/plotD.C | 73 - prod/acrcaf/qa_pp/plot_macros/plotDca.C | 61 - prod/acrcaf/qa_pp/plot_macros/plotDcaVz.C | 42 - prod/acrcaf/qa_pp/plot_macros/plotDeltaPhi.C | 47 - .../qa_pp/plot_macros/plotEta-over-pt.C | 56 - prod/acrcaf/qa_pp/plot_macros/plotEta.C | 57 - prod/acrcaf/qa_pp/plot_macros/plotEtaPtVz.C | 47 - prod/acrcaf/qa_pp/plot_macros/plotEtaVz.C | 42 - prod/acrcaf/qa_pp/plot_macros/plotLadderDca.C | 81 - prod/acrcaf/qa_pp/plot_macros/plotLadderEta.C | 76 - prod/acrcaf/qa_pp/plot_macros/plotLadderPhi.C | 77 - prod/acrcaf/qa_pp/plot_macros/plotLadderPt.C | 79 - .../qa_pp/plot_macros/plotLadderSigma.C | 80 - prod/acrcaf/qa_pp/plot_macros/plotPhi.C | 56 - prod/acrcaf/qa_pp/plot_macros/plotPhiVz.C | 42 - prod/acrcaf/qa_pp/plot_macros/plotPhi_Eta.C | 56 - prod/acrcaf/qa_pp/plot_macros/plotPt.C | 56 - prod/acrcaf/qa_pp/plot_macros/plotPtVz.C | 42 - prod/acrcaf/qa_pp/plot_macros/plotPt_Eta.C | 56 - prod/acrcaf/qa_pp/plot_macros/plotQPtVz.C | 52 - prod/acrcaf/qa_pp/plot_macros/plotQPt_Eta.C | 56 - prod/acrcaf/qa_pp/plot_macros/plotSectorDca.C | 76 - .../qa_pp/plot_macros/plotSectorDcaRatio.C | 77 - prod/acrcaf/qa_pp/plot_macros/plotSectorEta.C | 78 - .../qa_pp/plot_macros/plotSectorEtaPt.C | 74 - .../qa_pp/plot_macros/plotSectorEtaPtRatio.C | 77 - .../qa_pp/plot_macros/plotSectorEtaRatio.C | 78 - prod/acrcaf/qa_pp/plot_macros/plotSectorPt.C | 78 - .../qa_pp/plot_macros/plotSectorPtRatio.C | 78 - prod/acrcaf/qa_pp/plot_macros/plotSectorQPt.C | 75 - .../qa_pp/plot_macros/plotSectorQPtRatio.C | 78 - prod/acrcaf/qa_pp/plot_macros/pos-NegInv.C | 74 - prod/acrcaf/qa_pp/plot_macros/sigma.C | 95 -- prod/acrcaf/qa_pp/qa_pp.C | 69 - 41 files changed, 4096 deletions(-) delete mode 100755 prod/acrcaf/qa_pp.sh delete mode 100644 prod/acrcaf/qa_pp/AddTaskCalorimeterQA.C delete mode 100644 prod/acrcaf/qa_pp/AddTaskQAsym.C delete mode 100644 prod/acrcaf/qa_pp/AddTaskVZEROQA.C delete mode 100644 prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx delete mode 100644 prod/acrcaf/qa_pp/AliAnalysisTaskQASym.h delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plot.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotD.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotDca.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotDcaVz.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotDeltaPhi.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotEta-over-pt.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotEta.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotEtaPtVz.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotEtaVz.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotLadderDca.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotLadderEta.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotLadderPhi.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotLadderPt.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotLadderSigma.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotPhi.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotPhiVz.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotPhi_Eta.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotPt.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotPtVz.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotPt_Eta.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotQPtVz.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotQPt_Eta.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorDca.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorDcaRatio.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorEta.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPt.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPtRatio.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorEtaRatio.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorPt.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorPtRatio.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorQPt.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/plotSectorQPtRatio.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/pos-NegInv.C delete mode 100644 prod/acrcaf/qa_pp/plot_macros/sigma.C delete mode 100644 prod/acrcaf/qa_pp/qa_pp.C diff --git a/prod/acrcaf/qa_pp.sh b/prod/acrcaf/qa_pp.sh deleted file mode 100755 index c8bb119cd05..00000000000 --- a/prod/acrcaf/qa_pp.sh +++ /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 "; 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 index ac70c8cf37d..00000000000 --- a/prod/acrcaf/qa_pp/AddTaskCalorimeterQA.C +++ /dev/null @@ -1,113 +0,0 @@ -AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kUseKinematics = kFALSE, Bool_t kPrintSettings = kFALSE) -{ - // Creates a PartCorr task for calorimeters performance studies, configures it and adds it to the analysis manager. - - // Get the pointer to the existing analysis manager via the static access method. - //============================================================================== - AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); - if (!mgr) { - ::Error("AddTaskPartCorr", "No analysis manager to connect to."); - return NULL; - } - - // Check the analysis type using the event handlers connected to the analysis manager. - //============================================================================== - if (!mgr->GetInputEventHandler()) { - ::Error("AddTaskPartCorr", "This task requires an input event handler"); - return NULL; - } - TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD" - - // Configure analysis - //=========================================================================== - - //Reader - //For this particular analysis few things done by the reader. - //Nothing else needs to be set. - AliCaloTrackReader * reader = 0x0; - if(data=="AOD") reader = new AliCaloTrackAODReader(); - else if(data=="ESD") reader = new AliCaloTrackESDReader(); - //reader->SetDebug(10);//10 for lots of messages - if(kPrintSettings) reader->Print(""); - - if(kUseKinematics){ - if(inputDataType == "ESD"){ - reader->SwitchOnStack(); - reader->SwitchOffAODMCParticles(); - } - else if(inputDataType == "AOD"){ - reader->SwitchOffStack(); - reader->SwitchOnAODMCParticles(); - } - } - - reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches. - - - // ##### Analysis algorithm settings #### - //Only needed now for MC data - //AliFiducialCut * fidCut = new AliFiducialCut(); - //fidCut->DoCTSFiducialCut(kFALSE) ; - //fidCut->DoEMCALFiducialCut(kTRUE) ; - //fidCut->DoPHOSFiducialCut(kTRUE) ; - - AliAnaCalorimeterQA *emcalQA = new AliAnaCalorimeterQA(); - //emcalQA->SetDebug(2); //10 for lots of messages - emcalQA->SetCalorimeter("EMCAL"); - if(kUseKinematics) emcalQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet. - else emcalQA->SwitchOffDataMC() ; - emcalQA->AddToHistogramsName("EMCAL_"); //Begining of histograms name - //emcalQA->SetFiducialCut(fidCut); - emcalQA->SwitchOffFidutialCut(); - if(kPrintSettings) emcalQA->Print(""); - //emcalQA->SwithOnPlotsMaking(); - //emcalQA->GetMCAnalysisUtils()->SetDebug(10); - - AliAnaCalorimeterQA *phosQA = new AliAnaCalorimeterQA(); - //phosQA->SetDebug(2); //10 for lots of messages - phosQA->SetCalorimeter("PHOS"); - if(kUseKinematics) phosQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet. - else phosQA->SwitchOffDataMC() ; - phosQA->AddToHistogramsName("PHOS_");//Begining of histograms name - //phosQA->SetFiducialCut(fidCut); - phosQA->SwitchOffFidutialCut(); - if(kPrintSettings)phosQA->Print(""); - //phosQA->GetMCAnalysisUtils()->SetDebug(10); - //phosQA->SwithOnPlotsMaking(); - - // #### Configure Maker #### - AliAnaPartCorrMaker * maker = new AliAnaPartCorrMaker(); - maker->SetReader(reader);//pointer to reader - maker->AddAnalysis(emcalQA,0); - maker->AddAnalysis(phosQA,1); - maker->SetAnaDebug(0) ; - maker->SwitchOnHistogramsMaker() ; - maker->SwitchOffAODsMaker() ; - if(kPrintSettings) maker->Print(""); - - - printf("======================== \n"); - printf(" End Configuration of Calorimeter QA \n"); - printf("======================== \n"); - - // Create task - //=========================================================================== - AliAnalysisTaskParticleCorrelation * task = new AliAnalysisTaskParticleCorrelation ("CalorimeterPerformance"); - task->SetConfigFileName(""); //Don't configure the analysis via configuration file. - task->SetDebugLevel(0); - task->SetAnalysisMaker(maker); - mgr->AddTask(task); - - AliAnalysisDataContainer *cout_pc = mgr->CreateContainer("Calo.Performance",TList::Class(), - AliAnalysisManager::kOutputContainer, "Calo.Performance.root"); - - // Create ONLY the output containers for the data produced by the task. - // Get and connect other common input/output containers via the manager as below - //============================================================================== - mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer()); - mgr->ConnectOutput (task, 1, cout_pc); - - return task; -} - - diff --git a/prod/acrcaf/qa_pp/AddTaskQAsym.C b/prod/acrcaf/qa_pp/AddTaskQAsym.C deleted file mode 100644 index b16cf4a7898..00000000000 --- a/prod/acrcaf/qa_pp/AddTaskQAsym.C +++ /dev/null @@ -1,120 +0,0 @@ -AliAnalysisTaskQASym * AddTaskQAsym(Int_t runNumber) - -{ - // Creates a QA task exploiting simple symmetries phi, eta +/-, charge ... - - // Get the pointer to the existing analysis manager via the static access method. - //============================================================================== - AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); - if (!mgr) { - ::Error("AddTaskQAsym", "No analysis manager to connect to."); - return NULL; - } - - // Check the analysis type using the event handlers connected to the analysis manager. - //============================================================================== - if (!mgr->GetInputEventHandler()) { - ::Error("AddTasQAsym", "This task requires an input event handler"); - return NULL; - } - TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD" - - // Configure analysis - //=========================================================================== - - - //Task for global tracks - AliAnalysisTaskQASym *task0 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_Global"); - task0->SetTrackType(0); - task0->SelectCollisionCandidates(); - //Task for ITS tracks - AliAnalysisTaskQASym *task1 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_ITS"); - task1->SetTrackType(1); - task1->SelectCollisionCandidates(); - //Task for TPC tracks - AliAnalysisTaskQASym *task2 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_TPC"); - task2->SetTrackType(2); - task2->SelectCollisionCandidates(); - - //cuts for global tracks - AliESDtrackCuts* esdTrackCutsL0 = new AliESDtrackCuts("AliESDtrackCuts0","Global"); - esdTrackCutsL0->SetMinNClustersTPC(70); - esdTrackCutsL0->SetRequireTPCRefit(kTRUE); - esdTrackCutsL0->SetMaxDCAToVertexXY(3.); - esdTrackCutsL0->SetMaxDCAToVertexZ(3.); - esdTrackCutsL0->SetAcceptKinkDaughters(kFALSE); - - AliESDtrackCuts* esdTrackCutsL1 = new AliESDtrackCuts("AliESDtrackCuts1","ITS"); - esdTrackCutsL1->SetMaxDCAToVertexXY(3.); - esdTrackCutsL1->SetMaxDCAToVertexZ(3.); - esdTrackCutsL1->SetAcceptKinkDaughters(kFALSE); - esdTrackCutsL1->SetRequireITSRefit(kTRUE); - esdTrackCutsL1->SetRequireITSStandAlone(kTRUE); - - //cuts for TPC tracks - AliESDtrackCuts* esdTrackCutsL2 = new AliESDtrackCuts("AliESDtrackCuts2","TPC"); - esdTrackCutsL2->SetRequireTPCRefit(kFALSE); - esdTrackCutsL2->SetAcceptKinkDaughters(kFALSE); - //jacek's cuts: - esdTrackCutsL2->SetMinNClustersTPC(70); - // cut on max ncl=160 in Task - esdTrackCutsL2->SetMaxDCAToVertexXY(3.); - esdTrackCutsL2->SetMaxDCAToVertexZ(3.); - esdTrackCutsL2->SetMaxChi2PerClusterTPC(3.999); - //cut minChi=0 in task - //esdTrackCutsL2->SetPRange(0.15,16); // not needed for QA - //esdTrackCutsL2->SetEtaRange(-0.8, 0.7999); // not needed for QA - - - task0->SetCuts(esdTrackCutsL0); - task1->SetCuts(esdTrackCutsL1); - task2->SetCuts(esdTrackCutsL2); - - mgr->AddTask(task0); - mgr->AddTask(task1); - mgr->AddTask(task2); - - AliAnalysisDataContainer *cout0 = 0; - AliAnalysisDataContainer *cout1 = 0; - AliAnalysisDataContainer *cout2 = 0; - - if(runNumber>0){ - cout0 = mgr->CreateContainer("QAsymHists_Global",TList::Class(), - AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber)); - cout1 = mgr->CreateContainer("QAsymHists_ITS",TList::Class(), - AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber)); - cout2 = mgr->CreateContainer("QAsymHists_TPC",TList::Class(), - AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber)); - } - - else{ - cout0 = mgr->CreateContainer("QAsymHists_Global",TList::Class(), - AliAnalysisManager::kOutputContainer, - Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName())); - cout1 = mgr->CreateContainer("QAsymHists_ITS",TList::Class(), - AliAnalysisManager::kOutputContainer, - Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName())); - cout2 = mgr->CreateContainer("QAsymHists_TPC",TList::Class(), - AliAnalysisManager::kOutputContainer, - Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName())); - } - - - mgr->ConnectInput (task0, 0, mgr->GetCommonInputContainer()); - mgr->ConnectInput (task1, 0, mgr->GetCommonInputContainer()); - mgr->ConnectInput (task2, 0, mgr->GetCommonInputContainer()); - - mgr->ConnectOutput (task0, 0, mgr->GetCommonOutputContainer()); - mgr->ConnectOutput (task1, 0, mgr->GetCommonOutputContainer()); - mgr->ConnectOutput (task2, 0, mgr->GetCommonOutputContainer()); - - mgr->ConnectOutput (task0, 1, cout0); - mgr->ConnectOutput (task1, 1, cout1); - mgr->ConnectOutput (task2, 1, cout2); - - - return task0; - -} - - diff --git a/prod/acrcaf/qa_pp/AddTaskVZEROQA.C b/prod/acrcaf/qa_pp/AddTaskVZEROQA.C deleted file mode 100644 index a8ed634efe2..00000000000 --- a/prod/acrcaf/qa_pp/AddTaskVZEROQA.C +++ /dev/null @@ -1,40 +0,0 @@ -AliAnalysisTaskSE* AddTaskVZEROQA(Int_t runNumber) -{ - // Creates a QA task exploiting simple symmetries phi, eta +/-, charge ... - - // Get the pointer to the existing analysis manager via the static access method. - //============================================================================== - AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); - if (!mgr) { - ::Error("AddTaskQAsym", "No analysis manager to connect to."); - return NULL; - } - - // Check the analysis type using the event handlers connected to the analysis manager. - //============================================================================== - if (!mgr->GetInputEventHandler()) { - ::Error("AddTasQAsym", "This task requires an input event handler"); - return NULL; - } - TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD" - - // Configure analysis - //=========================================================================== - - - - AliAnaVZEROQA* task = new AliAnaVZEROQA("AliAnaVZEROQA"); - mgr->AddTask(task); - - AliAnalysisDataContainer *cout = mgr->CreateContainer("QAVZEROHists",TList::Class(), - AliAnalysisManager::kOutputContainer, "VZERO.Performance.root"); - - mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer()); - mgr->ConnectOutput (task, 1, cout); - - return task; - - -} - - diff --git a/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx b/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx deleted file mode 100644 index 7d642346dc8..00000000000 --- a/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.cxx +++ /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; iGetEntries(); ++i) { -// TH1 *h1 = dynamic_cast(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(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(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(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(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 index d20b5b0f14f..00000000000 --- a/prod/acrcaf/qa_pp/AliAnalysisTaskQASym.h +++ /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 index eaf3caadfdc..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plot.C +++ /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 index 4a9430252dc..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotD.C +++ /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 index e5fcba2d923..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotDca.C +++ /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 index e0bd0cf5e59..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotDcaVz.C +++ /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 index e3c9bed5921..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotDeltaPhi.C +++ /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 index 64eabc964c1..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotEta-over-pt.C +++ /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 index 8a67c8b5aa5..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotEta.C +++ /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 index 48706561307..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotEtaPtVz.C +++ /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 index 116fd8bbba2..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotEtaVz.C +++ /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 index efa68454099..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotLadderDca.C +++ /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 index 3340bfe40fc..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotLadderEta.C +++ /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 index a8eca3f7910..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotLadderPhi.C +++ /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 index 9a44c3ca49a..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotLadderPt.C +++ /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 index 38cded49c62..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotLadderSigma.C +++ /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 index bec80143a75..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotPhi.C +++ /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 index 6f894531497..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotPhiVz.C +++ /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 index d1bb7c5084e..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotPhi_Eta.C +++ /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 index 989091c1ebd..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotPt.C +++ /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 index 6a5d9b44f72..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotPtVz.C +++ /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 index a7f3b7f7066..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotPt_Eta.C +++ /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 index ca7dc63a115..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotQPtVz.C +++ /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 index cb9d8e8a86d..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotQPt_Eta.C +++ /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 index 37ca2de076e..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorDca.C +++ /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 index ff3b955d9a5..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorDcaRatio.C +++ /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 index d4ab68814c3..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorEta.C +++ /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 index b08fe24e949..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPt.C +++ /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 index 4fd0dd5e048..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaPtRatio.C +++ /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 index 604be1c00a8..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorEtaRatio.C +++ /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 index c0a8fadd15b..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorPt.C +++ /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 index 7c61ad986d9..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorPtRatio.C +++ /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 index e9b231ba335..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorQPt.C +++ /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 index 1c8b823e639..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/plotSectorQPtRatio.C +++ /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 index d92089bb788..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/pos-NegInv.C +++ /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 index 82b4ab8b2e0..00000000000 --- a/prod/acrcaf/qa_pp/plot_macros/sigma.C +++ /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 index 488fa64646d..00000000000 --- a/prod/acrcaf/qa_pp/qa_pp.C +++ /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(); -} -- 2.43.0