+++ /dev/null
-#!/bin/bash
-
-source /afs/cern.ch/alice/caf/caf-lxplus.sh -alien v4-18-12-AN
-
- [ -d qa_pp ] || mkdir qa_pp
-
- [ -z $1 ] && { echo "Usage: qa_pp.sh <run_number>"; exit 1; }
-
- cd qa_pp
- root.exe qa_pp.C\($1\)
-
- rm plot_macros/QAsym.proof.root
- ln -s /home/alishift/acrcaf/qa_pp/run$1.root plot_macros/QAsym.proof.root
-
+++ /dev/null
-AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kUseKinematics = kFALSE, Bool_t kPrintSettings = kFALSE)\r
-{\r
- // Creates a PartCorr task for calorimeters performance studies, configures it and adds it to the analysis manager.\r
- \r
- // Get the pointer to the existing analysis manager via the static access method.\r
- //==============================================================================\r
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
- if (!mgr) {\r
- ::Error("AddTaskPartCorr", "No analysis manager to connect to.");\r
- return NULL;\r
- } \r
- \r
- // Check the analysis type using the event handlers connected to the analysis manager.\r
- //==============================================================================\r
- if (!mgr->GetInputEventHandler()) {\r
- ::Error("AddTaskPartCorr", "This task requires an input event handler");\r
- return NULL;\r
- }\r
- TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
- \r
- // Configure analysis\r
- //===========================================================================\r
- \r
- //Reader\r
- //For this particular analysis few things done by the reader.\r
- //Nothing else needs to be set.\r
- AliCaloTrackReader * reader = 0x0;\r
- if(data=="AOD") reader = new AliCaloTrackAODReader();\r
- else if(data=="ESD") reader = new AliCaloTrackESDReader();\r
- //reader->SetDebug(10);//10 for lots of messages\r
- if(kPrintSettings) reader->Print("");\r
- \r
- if(kUseKinematics){\r
- if(inputDataType == "ESD"){\r
- reader->SwitchOnStack(); \r
- reader->SwitchOffAODMCParticles(); \r
- }\r
- else if(inputDataType == "AOD"){\r
- reader->SwitchOffStack(); \r
- reader->SwitchOnAODMCParticles(); \r
- }\r
- }\r
- \r
- reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches.\r
-\r
- \r
- // ##### Analysis algorithm settings ####\r
- //Only needed now for MC data\r
- //AliFiducialCut * fidCut = new AliFiducialCut();\r
- //fidCut->DoCTSFiducialCut(kFALSE) ;\r
- //fidCut->DoEMCALFiducialCut(kTRUE) ;\r
- //fidCut->DoPHOSFiducialCut(kTRUE) ;\r
- \r
- AliAnaCalorimeterQA *emcalQA = new AliAnaCalorimeterQA();\r
- //emcalQA->SetDebug(2); //10 for lots of messages\r
- emcalQA->SetCalorimeter("EMCAL");\r
- if(kUseKinematics) emcalQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.\r
- else emcalQA->SwitchOffDataMC() ;\r
- emcalQA->AddToHistogramsName("EMCAL_"); //Begining of histograms name\r
- //emcalQA->SetFiducialCut(fidCut);\r
- emcalQA->SwitchOffFidutialCut();\r
- if(kPrintSettings) emcalQA->Print(""); \r
- //emcalQA->SwithOnPlotsMaking();\r
- //emcalQA->GetMCAnalysisUtils()->SetDebug(10);\r
- \r
- AliAnaCalorimeterQA *phosQA = new AliAnaCalorimeterQA();\r
- //phosQA->SetDebug(2); //10 for lots of messages\r
- phosQA->SetCalorimeter("PHOS");\r
- if(kUseKinematics) phosQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.\r
- else phosQA->SwitchOffDataMC() ; \r
- phosQA->AddToHistogramsName("PHOS_");//Begining of histograms name\r
- //phosQA->SetFiducialCut(fidCut);\r
- phosQA->SwitchOffFidutialCut();\r
- if(kPrintSettings)phosQA->Print(""); \r
- //phosQA->GetMCAnalysisUtils()->SetDebug(10);\r
- //phosQA->SwithOnPlotsMaking();\r
- \r
- // #### Configure Maker ####\r
- AliAnaPartCorrMaker * maker = new AliAnaPartCorrMaker();\r
- maker->SetReader(reader);//pointer to reader\r
- maker->AddAnalysis(emcalQA,0);\r
- maker->AddAnalysis(phosQA,1);\r
- maker->SetAnaDebug(0) ;\r
- maker->SwitchOnHistogramsMaker() ;\r
- maker->SwitchOffAODsMaker() ;\r
- if(kPrintSettings) maker->Print("");\r
- \r
-\r
- printf("======================== \n");\r
- printf(" End Configuration of Calorimeter QA \n");\r
- printf("======================== \n");\r
- \r
- // Create task\r
- //===========================================================================\r
- AliAnalysisTaskParticleCorrelation * task = new AliAnalysisTaskParticleCorrelation ("CalorimeterPerformance");\r
- task->SetConfigFileName(""); //Don't configure the analysis via configuration file.\r
- task->SetDebugLevel(0);\r
- task->SetAnalysisMaker(maker); \r
- mgr->AddTask(task);\r
- \r
- AliAnalysisDataContainer *cout_pc = mgr->CreateContainer("Calo.Performance",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, "Calo.Performance.root");\r
- \r
- // Create ONLY the output containers for the data produced by the task.\r
- // Get and connect other common input/output containers via the manager as below\r
- //==============================================================================\r
- mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer());\r
- mgr->ConnectOutput (task, 1, cout_pc);\r
- \r
- return task;\r
-}\r
-\r
-\r
+++ /dev/null
-AliAnalysisTaskQASym * AddTaskQAsym(Int_t runNumber)\r
-\r
-{\r
- // Creates a QA task exploiting simple symmetries phi, eta +/-, charge ...\r
- \r
- // Get the pointer to the existing analysis manager via the static access method.\r
- //==============================================================================\r
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
- if (!mgr) {\r
- ::Error("AddTaskQAsym", "No analysis manager to connect to.");\r
- return NULL;\r
- } \r
- \r
- // Check the analysis type using the event handlers connected to the analysis manager.\r
- //==============================================================================\r
- if (!mgr->GetInputEventHandler()) {\r
- ::Error("AddTasQAsym", "This task requires an input event handler");\r
- return NULL;\r
- }\r
- TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
- \r
- // Configure analysis\r
- //===========================================================================\r
- \r
- \r
- //Task for global tracks\r
- AliAnalysisTaskQASym *task0 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_Global");\r
- task0->SetTrackType(0);\r
- task0->SelectCollisionCandidates();\r
- //Task for ITS tracks \r
- AliAnalysisTaskQASym *task1 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_ITS");\r
- task1->SetTrackType(1);\r
- task1->SelectCollisionCandidates();\r
- //Task for TPC tracks \r
- AliAnalysisTaskQASym *task2 = new AliAnalysisTaskQASym("AliAnalysisTaskQASym_TPC");\r
- task2->SetTrackType(2);\r
- task2->SelectCollisionCandidates();\r
-\r
- //cuts for global tracks\r
- AliESDtrackCuts* esdTrackCutsL0 = new AliESDtrackCuts("AliESDtrackCuts0","Global");\r
- esdTrackCutsL0->SetMinNClustersTPC(70);\r
- esdTrackCutsL0->SetRequireTPCRefit(kTRUE);\r
- esdTrackCutsL0->SetMaxDCAToVertexXY(3.);\r
- esdTrackCutsL0->SetMaxDCAToVertexZ(3.);\r
- esdTrackCutsL0->SetAcceptKinkDaughters(kFALSE);\r
-\r
- AliESDtrackCuts* esdTrackCutsL1 = new AliESDtrackCuts("AliESDtrackCuts1","ITS");\r
- esdTrackCutsL1->SetMaxDCAToVertexXY(3.);\r
- esdTrackCutsL1->SetMaxDCAToVertexZ(3.);\r
- esdTrackCutsL1->SetAcceptKinkDaughters(kFALSE);\r
- esdTrackCutsL1->SetRequireITSRefit(kTRUE);\r
- esdTrackCutsL1->SetRequireITSStandAlone(kTRUE);\r
- \r
- //cuts for TPC tracks\r
- AliESDtrackCuts* esdTrackCutsL2 = new AliESDtrackCuts("AliESDtrackCuts2","TPC");\r
- esdTrackCutsL2->SetRequireTPCRefit(kFALSE);\r
- esdTrackCutsL2->SetAcceptKinkDaughters(kFALSE);\r
- //jacek's cuts:\r
- esdTrackCutsL2->SetMinNClustersTPC(70);\r
- // cut on max ncl=160 in Task\r
- esdTrackCutsL2->SetMaxDCAToVertexXY(3.);\r
- esdTrackCutsL2->SetMaxDCAToVertexZ(3.);\r
- esdTrackCutsL2->SetMaxChi2PerClusterTPC(3.999);\r
- //cut minChi=0 in task\r
- //esdTrackCutsL2->SetPRange(0.15,16); // not needed for QA\r
- //esdTrackCutsL2->SetEtaRange(-0.8, 0.7999); // not needed for QA\r
- \r
-\r
- task0->SetCuts(esdTrackCutsL0);\r
- task1->SetCuts(esdTrackCutsL1);\r
- task2->SetCuts(esdTrackCutsL2);\r
-\r
- mgr->AddTask(task0);\r
- mgr->AddTask(task1);\r
- mgr->AddTask(task2);\r
- \r
- AliAnalysisDataContainer *cout0 = 0;\r
- AliAnalysisDataContainer *cout1 = 0;\r
- AliAnalysisDataContainer *cout2 = 0;\r
- \r
- if(runNumber>0){ \r
- cout0 = mgr->CreateContainer("QAsymHists_Global",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber));\r
- cout1 = mgr->CreateContainer("QAsymHists_ITS",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber));\r
- cout2 = mgr->CreateContainer("QAsymHists_TPC",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, Form("run%d.root",runNumber));\r
- }\r
- \r
- else{\r
- cout0 = mgr->CreateContainer("QAsymHists_Global",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, \r
- Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName()));\r
- cout1 = mgr->CreateContainer("QAsymHists_ITS",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, \r
- Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName()));\r
- cout2 = mgr->CreateContainer("QAsymHists_TPC",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, \r
- Form("%s:PWG1_QAsymHists",AliAnalysisManager::GetCommonFileName()));\r
- }\r
-\r
-\r
- mgr->ConnectInput (task0, 0, mgr->GetCommonInputContainer());\r
- mgr->ConnectInput (task1, 0, mgr->GetCommonInputContainer());\r
- mgr->ConnectInput (task2, 0, mgr->GetCommonInputContainer());\r
- \r
- mgr->ConnectOutput (task0, 0, mgr->GetCommonOutputContainer());\r
- mgr->ConnectOutput (task1, 0, mgr->GetCommonOutputContainer());\r
- mgr->ConnectOutput (task2, 0, mgr->GetCommonOutputContainer());\r
-\r
- mgr->ConnectOutput (task0, 1, cout0);\r
- mgr->ConnectOutput (task1, 1, cout1);\r
- mgr->ConnectOutput (task2, 1, cout2);\r
-\r
- \r
- return task0;\r
-\r
-}\r
-\r
-\r
+++ /dev/null
-AliAnalysisTaskSE* AddTaskVZEROQA(Int_t runNumber)\r
-{\r
- // Creates a QA task exploiting simple symmetries phi, eta +/-, charge ...\r
- \r
- // Get the pointer to the existing analysis manager via the static access method.\r
- //==============================================================================\r
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
- if (!mgr) {\r
- ::Error("AddTaskQAsym", "No analysis manager to connect to.");\r
- return NULL;\r
- } \r
- \r
- // Check the analysis type using the event handlers connected to the analysis manager.\r
- //==============================================================================\r
- if (!mgr->GetInputEventHandler()) {\r
- ::Error("AddTasQAsym", "This task requires an input event handler");\r
- return NULL;\r
- }\r
- TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
- \r
- // Configure analysis\r
- //===========================================================================\r
- \r
- \r
- \r
- AliAnaVZEROQA* task = new AliAnaVZEROQA("AliAnaVZEROQA");\r
- mgr->AddTask(task);\r
- \r
- AliAnalysisDataContainer *cout = mgr->CreateContainer("QAVZEROHists",TList::Class(),\r
- AliAnalysisManager::kOutputContainer, "VZERO.Performance.root");\r
-\r
- mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer());\r
- mgr->ConnectOutput (task, 1, cout);\r
-\r
- return task;\r
- \r
- \r
-}\r
-\r
-\r
+++ /dev/null
-#include "TChain.h"
-#include "TTree.h"
-#include "TH1F.h"
-#include "TH2F.h"
-#include "TH3F.h"
-#include "TCanvas.h"
-#include "TList.h"
-#include "TParticle.h"
-#include "TParticlePDG.h"
-#include "TProfile.h"
-#include "TNtuple.h"
-#include "TFile.h"
-
-#include "AliAnalysisTask.h"
-#include "AliAnalysisManager.h"
-
-#include "AliESDEvent.h"
-#include "AliLog.h"
-#include "AliESDVertex.h"
-#include "AliESDInputHandler.h"
-#include "AliESDtrackCuts.h"
-#include "AliMultiplicity.h"
-
-#include "AliAnalysisTaskQASym.h"
-#include "AliExternalTrackParam.h"
-#include "AliTrackReference.h"
-#include "AliHeader.h"
-#include "AliGenEventHeader.h"
-#include "AliGenDPMjetEventHeader.h"
-
-// Analysis Task for basic QA on the ESD
-
-// Authors: Jan Fiete Grosse-Oetringhaus, Christian Klein-Boesing,
-// Andreas Morsch, Eva Sicking
-
-ClassImp(AliAnalysisTaskQASym)
-
- //________________________________________________________________________
- AliAnalysisTaskQASym::AliAnalysisTaskQASym(const char *name)
- : AliAnalysisTaskSE(name)
- ,fTrackType(0)
- ,fFieldOn(kTRUE)
-
- ,fHists(0)
-
- ,fHistRECpt(0)
- ,fEta(0)
- ,fEtaPhi(0)
- ,fEtaPt(0)
- ,fQPt(0)
- ,fDca(0)
- ,fqRec(0)
- ,fsigmaPt(0)
-
- ,fRecPtPos(0)
- ,fRecPtNeg(0)
- ,fRecPhiPos(0)
- ,fRecPhiNeg(0)
- ,fRecEtaPos(0)
- ,fRecEtaNeg(0)
- ,fRecEtaPtPos(0)
- ,fRecEtaPtNeg(0)
- ,fRecDcaPos(0)
- ,fRecDcaNeg(0)
- ,fRecDcaNegInv(0)
- ,fRecDPos(0)
- ,fRecDNeg(0)
-
-
- ,fRecQPtPosEta(0)
- ,fRecQPtNegEta(0)
- ,fRecPtPosEta(0)
- ,fRecPtNegEta(0)
- ,fRecPhiPosEta(0)
- ,fRecPhiNegEta(0)
- ,fRecDcaPosEta(0)
- ,fRecDcaNegEta(0)
- ,fRecDPosEta(0)
- ,fRecDNegEta(0)
-
- ,fRecPtPosVz(0)
- ,fRecPtNegVz(0)
- ,fRecEtaPosVz(0)
- ,fRecEtaNegVz(0)
- ,fRecPhiPosVz(0)
- ,fRecPhiNegVz(0)
- ,fSignedDcaPosVz(0)
- ,fSignedDcaNegVz(0)
- ,fRecQPtPosEtaVz(0)
- ,fRecQPtNegEtaVz(0)
- ,fRecEtaPtPosVz(0)
- ,fRecEtaPtNegVz(0)
-
-
- ,fDeltaPhiAll(0)
- ,fDeltaPhiLeading(0)
- ,fDiffDcaD(0)
-
- ,fPhiRec(0)
- ,fThetaRec(0)
- ,fNumber(0)
- ,fVx(0)
- ,fVy(0)
- ,fVz(0)
- ,test(0)
-
- ,fRecDcaPosPhi(0)
- ,fRecDcaNegPhi(0)
- ,fRecPtPosPhi(0)
- ,fRecPtNegPhi(0)
- ,fRecEtaPosPhi(0)
- ,fRecEtaNegPhi(0)
- ,fRecQPtPhi(0)
- ,fRecEtaPtPosPhi(0)
- ,fRecEtaPtNegPhi(0)
-
- ,fRecPtPosEtaPos(0)
- ,fRecPtNegEtaPos(0)
- ,fRecPtPosEtaNeg(0)
- ,fRecPtNegEtaNeg(0)
-
- ,fRec1PtPosEtaPos(0)
- ,fRec1PtNegEtaPos(0)
- ,fRec1PtPosEtaNeg(0)
- ,fRec1PtNegEtaNeg(0)
-
- ,fRecPhiPosEtaPos(0)
- ,fRecPhiNegEtaPos(0)
- ,fRecPhiPosEtaNeg(0)
- ,fRecPhiNegEtaNeg(0)
-
- ,fRecDcaPosPhiEtaPos(0)
- ,fRecDcaNegPhiEtaPos(0)
- ,fRecDcaPosPhiEtaNeg(0)
- ,fRecDcaNegPhiEtaNeg(0)
- ,fRecPtPosPhiEtaPos(0)
- ,fRecPtNegPhiEtaPos(0)
- ,fRecPtPosPhiEtaNeg(0)
- ,fRecPtNegPhiEtaNeg(0)
-
-
- ,fRecDcaPhiPtPosEtaPos(0)
- ,fRecDcaPhiPtNegEtaPos(0)
- ,fRecDcaPhiPtPosEtaNeg(0)
- ,fRecDcaPhiPtNegEtaNeg(0)
- ,fEtavPt(0)
-
-
-
- ,sdca(0)
- ,xy(0)
- ,z(0)
- ,xvertexcor(0)
- ,yvertexcor(0)
-
- ,fCuts(0)
-
-{
- // Constructor
- for(Int_t i = 0;i<18;++i){
- fRecPtTpcSector[i] = 0;
- fRecEtaTpcSector[i] = 0;
- fSignedDcaTpcSector[i] = 0;
- fRecQPtTpcSector[i] = 0;
- fRecEtaPtTpcSector[i] = 0;
- }
-
- for(Int_t i = 0;i< 7;++i){
- fRecPtPosLadder[i] = 0;
- fRecPtNegLadder[i] = 0;
- fRecPhiPosLadder[i] = 0;
- fRecPhiNegLadder[i] = 0;
- fRecEtaPosLadder[i] = 0;
- fRecEtaNegLadder[i] = 0;
- fSignDcaPos[i] = 0;
- fSignDcaNeg[i] = 0;
- fSignDcaNegInv[i] = 0;
- fPtSigmaPos[i] =0;
- fPtSigmaNeg[i] =0;
- fqPtRec[i] =0;
- fDcaSigmaPos[i] =0;
- fDcaSigmaNeg[i] =0;
- }
-
- DefineOutput(1, TList::Class());
-
-
-
-}
-
-
-//________________________________________________________________________
-void AliAnalysisTaskQASym::UserCreateOutputObjects()
-{
- // Create histograms
- // Called once
-
- Bool_t oldStatus = TH1::AddDirectoryStatus();
- TH1::AddDirectory(kFALSE);
-
- Double_t range = 0.2;
- Double_t pt = 20.;
-
- fHists = new TList();
- // test = new TNtuple("test","test",
- // "pt:phi:theta:x:y:z:charge");
- fHistRECpt = new TH1F("fHistRECpt",
- " p_{T}",
- 100, 0., 0.6);
- fEta = new TH1F("fEta",
- " #eta",
- 200, -2., 2.);
- fEtavPt = new TH2F("fEtavPt",
- " #eta -p_{T}",
- 200, -2., 2.,
- 100, 0, 1.5);
- fEtaPhi = new TH2F("fEtaPhi",
- " #eta - #phi",
- 200, -2., 2., 128, 0., 2. * TMath::Pi());
-
- fThetaRec = new TH1F("fThetaRec",
- " #theta",
- 180, 0., TMath::Pi());
- fPhiRec = new TH1F("fPhiRec",
- " #phi",
- 180, 0., 2*TMath::Pi());
- fNumber = new TH1F("fNumber",
- "number of tracks per event",
- 50, 0.5, 49.5);
- fVx = new TH1F("fVx",
- "X of vertex",
- 100, -1., 1.);
- fVy = new TH1F("fVy",
- "Y of vertex",
- 100, -1., 1.);
- fVz = new TH1F("fVz",
- "Z of vertex",
- 200, -50., 50.);
-
- fEtaPt = new TH1F("fEtaPt",
- " #eta/p_{T} ",
- 100, -1., 1.);
-
- fQPt = new TH1F("fQPt",
- " charge/p_{T} ",
- 100, -1., 1.);
-
- fDca = new TH1F("fDca",
- " dca ",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
-
-
- fqRec = new TH1F("fqRec",
- " charge all reconstructed particle",
- 21, -9.5, 10.5);
-
- fsigmaPt = new TH1F("fsigmaPt",
- "Log_{10}(#sigma_{p_{T}})",
- 200, -4., 8.);
-
-
-
-
- //------------
- for(Int_t ITSlayer_case=0;ITSlayer_case<7;ITSlayer_case++){
-
- fSignDcaPos[ITSlayer_case] = new TH1F(Form("fSignDcaPos%d", ITSlayer_case),
- " Signed dca",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fSignDcaPos[ITSlayer_case]->GetXaxis()->SetTitle("dca");
- fSignDcaPos[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-
- fSignDcaNeg[ITSlayer_case] = new TH1F(Form("fSignDcaNeg%d", ITSlayer_case),
- " Signed dcas",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fSignDcaNeg[ITSlayer_case]->GetXaxis()->SetTitle("dca");
- fSignDcaNeg[ITSlayer_case]->GetYaxis()->SetTitle("");
-
- fSignDcaNegInv[ITSlayer_case] = new TH1F(Form("fSignDcaNegInv%d", ITSlayer_case),
- " inverse Signed dca ",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fSignDcaNegInv[ITSlayer_case]->GetXaxis()->SetTitle("-dca");
- fSignDcaNegInv[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-
-
-
- fPtSigmaPos[ITSlayer_case] = new TH1F(Form("fPtSigmaPos%d", ITSlayer_case),
- " #sigma_{pT} ",
- 208, -4., 8.);
- fPtSigmaPos[ITSlayer_case]->GetXaxis()->SetTitle("Log_{10}(#sigma_{pT})");
- fPtSigmaPos[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-
- fPtSigmaNeg[ITSlayer_case] = new TH1F(Form("fPtSigmaNeg%d",ITSlayer_case),
- " #sigma_{pT}",
- 208, -4., 8.);
- fPtSigmaNeg[ITSlayer_case]->GetXaxis()->SetTitle("Log_{10}(#sigma_{pT})");
- fPtSigmaNeg[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-
-
-
-
- fqPtRec[ITSlayer_case] = new TH1F(Form("fqPtRec%d",ITSlayer_case),
- "q/ p_{T}",
- 200, -100., 100.);
- fqPtRec[ITSlayer_case]->GetXaxis()->SetTitle("q_{tr}/p_{T, tr} (GeV/c)");
- fqPtRec[ITSlayer_case]->GetYaxis()->SetTitle("");
-
-
-
-
-
- fDcaSigmaPos[ITSlayer_case] = new TH2F(Form("fDcaSigmaPos%d", ITSlayer_case),
- " p_{T} shift vs #sigma_{pT} ",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200, -4., 4. );
- fDcaSigmaPos[ITSlayer_case]->GetXaxis()->SetTitle("signed DCA)");
- fDcaSigmaPos[ITSlayer_case]->GetYaxis()->SetTitle("log_{10}(#sigma_{pT})");
-
-
- fDcaSigmaNeg[ITSlayer_case] = new TH2F(Form("fDcaSigmaNeg%d", ITSlayer_case),
- " p_{T} shift vs #sigma_{pT} ",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200, -4., 4. );
- fDcaSigmaNeg[ITSlayer_case]->GetXaxis()->SetTitle("signed DCA");
- fDcaSigmaNeg[ITSlayer_case]->GetYaxis()->SetTitle("log_{10}(#sigma_{pT})");
-
-
- }
-
-
-
- // YIELDs---------- positive and negative particles
-
- fRecPtPos = new TH1F("fRecPtPos",
- " p_{T}",
- 100, 0.,pt);
- fRecPtPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
- fRecPtNeg = new TH1F("fRecPtNeg",
- " p_{T} ",
- 100, 0., pt);
- fRecPtNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
- fRecPhiPos = new TH1F("fRecPhiPos",
- "#phi",
- 361, 0., 360.);
- fRecPhiPos->GetXaxis()->SetTitle("#phi (deg)");
-
- fRecPhiNeg = new TH1F("fRecPhiNeg",
- "#phi ",
- 361, 0., 360.);
- fRecPhiNeg->GetXaxis()->SetTitle("#phi (deg)");
-
- fRecEtaPos = new TH1F("fRecEtaPos",
- "#eta",
- 200, -2., 2.);
- fRecEtaPos->GetXaxis()->SetTitle("#eta");
-
- fRecEtaNeg = new TH1F("fRecEtaNeg",
- "#eta",
- 200, -2., 2.);
- fRecEtaNeg->GetXaxis()->SetTitle("#eta");
-
- fRecEtaPtPos = new TH1F("fRecEtaPtPos",
- "#eta/p_{T}",
- 200, -0.1, .1);
- fRecEtaPtPos->GetXaxis()->SetTitle("#eta/p_{T}");
-
- fRecEtaPtNeg = new TH1F("fRecEtaPtNeg",
- "#eta/p_{T}",
- 200, -.1, .1);
- fRecEtaPtNeg->GetXaxis()->SetTitle("#eta/p_{T}");
-
- fRecDcaPos = new TH1F("fRecDcaPos",
- " dca",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaPos->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaNeg = new TH1F("fRecDcaNeg",
- " dca",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaNeg->GetXaxis()->SetTitle("dca (cm)");
-
- fRecDcaNegInv = new TH1F("fRecDcaNegInv",
- " dca",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaNegInv->GetXaxis()->SetTitle("dca (cm)");
-
-
- fRecDPos = new TH1F("fRecDPos",
- " d",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDPos->GetXaxis()->SetTitle("d (cm)");
- fRecDNeg = new TH1F("fRecDNeg",
- "d",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDNeg->GetXaxis()->SetTitle("d (cm)");
-
-
- // YIELDs ---------------- positive and negative eta
-
-
- fRecQPtPosEta = new TH1F("fRecQPtPosEta",
- "q/p_{T}",
- 200, -0.5, 0.5);
- fRecQPtPosEta->GetXaxis()->SetTitle("q/p_{T} ");
-
- fRecQPtNegEta = new TH1F("fRecQPtNegEta",
- "q/p_{T}",
- 200, -0.5, 0.5);
- fRecQPtNegEta->GetXaxis()->SetTitle("q/p_{T}");
-
- fRecPtPosEta = new TH1F("fRecPtPosEta",
- " p_{T} ",
- 100, 0., pt);
- fRecPtPosEta->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecPtNegEta = new TH1F("fRecPtNegEta",
- " p_{T}",
- 100, 0., pt);
- fRecPtNegEta->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecPhiPosEta = new TH1F("fRecPhiPosEta",
- "#phi",
- 361, 0., 360);
- fRecPhiPosEta->GetXaxis()->SetTitle("#phi (deg)");
-
- fRecPhiNegEta = new TH1F("fRecPhiNegEta",
- "#phi ",
- 361, 0, 360);
- fRecPhiNegEta->GetXaxis()->SetTitle("#phi (deg)");
-
- fRecDcaPosEta = new TH1F("fRecDcaPosEta",
- " dca ",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaPosEta->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaNegEta = new TH1F("fRecDcaNegEta",
- " dca",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaNegEta->GetXaxis()->SetTitle("dca (cm)");
-
- fRecDPosEta = new TH1F("fRecDPosEta",
- " d",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDPosEta->GetXaxis()->SetTitle("d (cm)");
- fRecDNegEta = new TH1F("fRecDNegEta",
- "d",
- 100, -5., 5.);
- fRecDNegEta->GetXaxis()->SetTitle("d (cm)");
-
- fRecDcaPosPhi = new TH2F("fRecDcaPosPhi",
- " dca vs. phi",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
- fRecDcaPosPhi->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
- fRecDcaNegPhi = new TH2F("fRecDcaNegPhi",
- " dca vs. phi",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
- fRecDcaNegPhi->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
- fRecPtPosPhi = new TH2F("fRecPtPosPhi",
- " log(p_T) vs. phi",
- 100, -2.5, 2., 180, 0, TMath::Pi()*2);
- fRecPtPosPhi->GetXaxis()->SetTitle("log_{10}(p_{T})");
- fRecPtPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
- fRecPtNegPhi = new TH2F("fRecPtNegPhi",
- " log(p_T) vs. phi",
- 100,-2.5 , 2., 180, 0, TMath::Pi()*2);
- fRecPtNegPhi->GetXaxis()->SetTitle("log_{10}(p_{T})");
- fRecPtNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
- fRecEtaPosPhi = new TH2F("fRecEtaPosPhi",
- " eta vs. phi",
- 100, -1.5, 1.5, 180, 0, TMath::Pi()*2);
- fRecEtaPosPhi->GetXaxis()->SetTitle("#eta");
- fRecEtaPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
- fRecEtaNegPhi = new TH2F("fRecEtaNegPhi",
- " eta vs. phi",
- 100, -1.5, 1.5, 180, 0, TMath::Pi()*2);
- fRecEtaNegPhi->GetXaxis()->SetTitle("#eta");
- fRecEtaNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
- fRecQPtPhi = new TH2F("fRecQPtPhi",
- " charge/p_T vs. phi",
- 100,-1. , 1., 180, 0, TMath::Pi()*2);
- fRecQPtPhi->GetXaxis()->SetTitle("charge/p_{T}");
- fRecQPtPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
- fRecEtaPtPosPhi = new TH2F("fRecEtaPtPosPhi",
- " eta/p_T vs. phi",
- 100, -5, 5., 180, 0, TMath::Pi()*2);
- fRecEtaPtPosPhi->GetXaxis()->SetTitle("#eta/p_{T}");
- fRecEtaPtPosPhi->GetYaxis()->SetTitle("#phi (rad.)");
- fRecEtaPtNegPhi = new TH2F("fRecEtaPtNegPhi",
- " eta/p_T vs. phi",
- 100,-5 , 5., 180, 0, TMath::Pi()*2);
- fRecEtaPtNegPhi->GetXaxis()->SetTitle("#eta/p_{T}");
- fRecEtaPtNegPhi->GetYaxis()->SetTitle("#phi (rad.)");
-
-
-
-
-
- fRecDcaPosPhiEtaPos = new TH2F("fRecDcaPosPhiEtaPos",
- " dca vs. phi",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
- fRecDcaPosPhiEtaPos->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaPosPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
- fRecDcaNegPhiEtaPos = new TH2F("fRecDcaNegPhiEtaPos",
- " dca vs. phi",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
- fRecDcaNegPhiEtaPos->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaNegPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
-
- fRecPtPosPhiEtaPos = new TH2F("fRecPtPosPhiEtaPos",
- " log(p_T) vs. phi",
- 100, -2.5, 2., 180, 0, TMath::Pi()*2);
- fRecPtPosPhiEtaPos->GetXaxis()->SetTitle("log_{10}(p_{T})");
- fRecPtPosPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
- fRecPtNegPhiEtaPos = new TH2F("fRecPtNegPhiEtaPos",
- " log(p_T) vs. phi",
- 100,-2.5 , 2., 180, 0, TMath::Pi()*2);
- fRecPtNegPhiEtaPos->GetXaxis()->SetTitle("log_{10}(p_{T})");
- fRecPtNegPhiEtaPos->GetYaxis()->SetTitle("#phi (rad.)");
-
-
- fRecDcaPosPhiEtaNeg = new TH2F("fRecDcaPosPhiEtaNeg",
- " dca vs. phi",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
- fRecDcaPosPhiEtaNeg->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaPosPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
- fRecDcaNegPhiEtaNeg = new TH2F("fRecDcaNegPhiEtaNeg",
- " dca vs. phi",
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 180, 0, TMath::Pi()*2);
- fRecDcaNegPhiEtaNeg->GetXaxis()->SetTitle("dca (cm)");
- fRecDcaNegPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
-
- fRecPtPosPhiEtaNeg = new TH2F("fRecPtPosPhiEtaNeg",
- " log(p_T) vs. phi",
- 100, -2.5, 2., 180, 0, TMath::Pi()*2);
- fRecPtPosPhiEtaNeg->GetXaxis()->SetTitle("log_{10}(p_{T})");
- fRecPtPosPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
- fRecPtNegPhiEtaNeg = new TH2F("fRecPtNegPhiEtaNeg",
- " log(p_T) vs. phi",
- 100,-2.5 , 2., 180, 0, TMath::Pi()*2);
- fRecPtNegPhiEtaNeg->GetXaxis()->SetTitle("log_{10}(p_{T})");
- fRecPtNegPhiEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
-
-
-
-
- // YIELDs ---------------- for TPC sectors
- for(Int_t sector=0; sector<18;sector++){
-
-
- fRecPtTpcSector[sector] = new TH1F(Form("fRecPtTpcSector%02d",sector),
- Form("p_{T} distribution: TPC sector %d",
- sector),100, 0., pt);
- fRecPtTpcSector[sector]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecEtaTpcSector[sector] = new TH1F(Form("fRecEtaTpcSector%02d",sector),
- Form("#eta distribution: TPC sector %d",
- sector),200, -2., 2.);
- fRecEtaTpcSector[sector]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
- fSignedDcaTpcSector[sector] = new TH1F(Form("fSignedDcaTpcSector%02d",sector),
- Form("dca distribution: TPC sector %d",
- sector),200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9) );
- fSignedDcaTpcSector[sector]->GetXaxis()->SetTitle("dca");
-
- fRecQPtTpcSector[sector] = new TH1F(Form("fRecQPtTpcSector%02d",sector),
- Form("Q/ p_{T} distribution: TPC sector %d",
- sector),100, -1., 1.);
- fRecQPtTpcSector[sector]->GetXaxis()->SetTitle("Q/p_{T} (GeV/c)");
-
- fRecEtaPtTpcSector[sector] = new TH1F(Form("fRecEtaPtTpcSector%02d",sector),
- Form("#eta/ p_{T} distribution: TPC sector %d",
- sector),100, -1., 1.);
- fRecEtaPtTpcSector[sector]->GetXaxis()->SetTitle("#eta/p_{T} (GeV/c)");
-
- }
- // YIELDS ITS ladder
- for(Int_t i=0;i<7;i++){
- fRecPtPosLadder[i] = new TH1F(Form("fRecPtPosLadder%d", i),
- " p_{T} distribution",
- 100, 0., pt);
- fRecPtPosLadder[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
- fRecPtNegLadder[i] = new TH1F(Form("fRecPtNegLadder%d",i),
- " p_{T} distribution ",
- 100, 0., pt);
- fRecPtNegLadder[i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
- fRecPhiPosLadder[i] = new TH1F(Form("fRecPhiPosLadder%d",i),
- "#phi distribution: all pos eta",
- 361, 0., 360);
- fRecPhiPosLadder[i]->GetXaxis()->SetTitle("#phi (deg)");
-
- fRecPhiNegLadder[i] = new TH1F(Form("fRecPhiNegLadder%d", i),
- "#phi distribution: all neg eta",
- 361, 0, 360);
- fRecPhiNegLadder[i]->GetXaxis()->SetTitle("#phi (deg)");
-
-
-
- fRecEtaPosLadder[i] = new TH1F(Form("fRecEtaPosLadder%d",i),
- "#eta distribution",
- 200, -2., 2.);
- fRecEtaPosLadder[i]->GetXaxis()->SetTitle("#eta)");
-
- fRecEtaNegLadder[i] = new TH1F(Form("fRecEtaNegLadder%d", i),
- "#eta distribution",
- 200, -2., 2.);
- fRecEtaNegLadder[i]->GetXaxis()->SetTitle("#eta");
- }
-
- Double_t vzmax = 15;
-
- fRecPtPosVz = new TH2F("fRecPtPosVz",
- "p_{T} distribution vs Vz()",
- 100, -1., 2., 200,-vzmax,vzmax);
- fRecPtPosVz->GetXaxis()->SetTitle("log_{10}(p_{T})");
-
- fRecPtNegVz = new TH2F("fRecPtNegVz",
- "p_{T} distribution vs Vz()",
- 100, -1., 2.,200,-vzmax,vzmax);
- fRecPtNegVz->GetXaxis()->SetTitle("Log_{10}(p_{T})");
-
-
- fRecEtaPosVz= new TH2F("fRecEtaPosVz",
- "#eta distribution vs Vz()",
- 100, -2., 2., 200,-vzmax,vzmax);
- fRecEtaPosVz->GetXaxis()->SetTitle("#eta");
- fRecEtaNegVz = new TH2F("fRecEtaNegVz",
- "#eta distribution vs Vz()",
- 100, -2., 2.,200,-vzmax,vzmax);
- fRecEtaNegVz->GetXaxis()->SetTitle("#eta");
-
- fRecPhiPosVz= new TH2F("fRecPhiPosVz",
- "#eta distribution vs Vz()",
- 361, 0., 360., 200,-vzmax,vzmax);
- fRecPhiPosVz->GetXaxis()->SetTitle("#phi (deg)");
- fRecPhiNegVz = new TH2F("fRecPhiNegVz",
- "dca vs Vz()",
- 361, 0., 360.,200,-vzmax,vzmax);
- fRecPhiNegVz->GetXaxis()->SetTitle("#phi (deg)");
-
- fSignedDcaPosVz= new TH2F("fSignedDcaPosVz",
- "#eta distribution vs Vz()",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9), 200,-vzmax,vzmax);
- fSignedDcaPosVz->GetXaxis()->SetTitle("dca (cm)");
- fSignedDcaNegVz = new TH2F("fSignedDcaNegVz",
- "dca vs Vz()",
- 200, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9),200,-vzmax,vzmax);
- fSignedDcaNegVz->GetXaxis()->SetTitle("dca (cm)");
-
- fRecQPtPosEtaVz= new TH2F("fRecQPtPosEtaVz",
- " Q/p_{T} distribution vs Vz()",
- 100, -1., 1., 200,-vzmax,vzmax);
- fRecQPtPosEtaVz->GetXaxis()->SetTitle("Q/p_{T}");
- fRecQPtNegEtaVz = new TH2F("fRecQPtNegEtaVz",
- " Q/p_{T} distribution vs Vz()",
- 100, -1., 1.,200,-vzmax,vzmax);
- fRecQPtNegEtaVz->GetXaxis()->SetTitle("Q/p_{T}");
-
-
- fRecEtaPtPosVz= new TH2F("fRecEtaPtPosVz",
- " #eta/p_{T} distribution vs Vz()",
- 100, -1., 1., 200,-vzmax,vzmax);
- fRecEtaPtPosVz->GetXaxis()->SetTitle("#eta/p_{T");
- fRecEtaPtNegVz = new TH2F("fRecEtaPtNegVz",
- " #eta/p_{T} distribution vs Vz()",
- 100, -1., 1.,200,-vzmax,vzmax);
- fRecEtaPtNegVz->GetXaxis()->SetTitle("#eta/p_{T}");
-
- //new
- fDeltaPhiAll = new TH1F("fDeltaPhiAll",
- " #Delta #phi",200,-360,360);
- fDeltaPhiAll->GetXaxis()->SetTitle("#Delta #phi");
-
-
- fDeltaPhiLeading = new TH2F("fDeltaPhiLeading",
- " #Delta #phi",361,-360,360, 361,0, 360);
- fDeltaPhiLeading->GetXaxis()->SetTitle("#Delta #phi (deg.)");
- fDeltaPhiLeading->GetYaxis()->SetTitle("#phi_{leading particle} (deg.)");
-
- fDiffDcaD = new TH1F("fDiffDcaD",
- "dca-d",
- 200, -1., 1.);
-
-
- fRecPtPosEtaPos = new TH1F("fRecPtPosEtaPos",
- "p_{T} distribution",100,0,pt);
- fRecPtPosEtaPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecPtNegEtaPos = new TH1F("fRecPtNegEtaPos",
- "p_{T} distribution",100,0,pt);
- fRecPtNegEtaPos->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecPtPosEtaNeg = new TH1F("fRecPtPosEtaNeg",
- "p_{T} distribution",100,0,pt);
- fRecPtPosEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecPtNegEtaNeg = new TH1F("fRecPtNegEtaNeg",
- "p_{T} distribution",100,0,pt);
- fRecPtNegEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
-
- fRec1PtPosEtaPos = new TH1F("fRec1PtPosEtaPos",
- "1/p_{T} distribution",100,0,0.5);
- fRec1PtPosEtaPos->GetXaxis()->SetTitle("p_{T} (c/GeV)");
-
- fRec1PtNegEtaPos = new TH1F("fRec1PtNegEtaPos",
- "1/p_{T} distribution",100,0,0.5);
- fRec1PtNegEtaPos->GetXaxis()->SetTitle("p_{T} (c/GeV)");
-
- fRec1PtPosEtaNeg = new TH1F("fRec1PtPosEtaNeg",
- "1/p_{T} distribution",100,0,0.5);
- fRec1PtPosEtaNeg->GetXaxis()->SetTitle("p_{T} (c/GeV)");
-
- fRec1PtNegEtaNeg = new TH1F("fRec1PtNegEtaNeg",
- "1/p_{T} distribution",100,0,0.5);
- fRec1PtNegEtaNeg->GetXaxis()->SetTitle("1/p_{T} (c/GeV)");
-
-
-
- fRecPhiPosEtaPos = new TH1F("fRecPhiPosEtaPos",
- "#phi",180,0,2*TMath::Pi());
- fRecPhiPosEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
-
- fRecPhiNegEtaPos = new TH1F("fRecPhiNegEtaPos",
- "#phi",180,0,2*TMath::Pi());
- fRecPhiNegEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
-
- fRecPhiPosEtaNeg = new TH1F("fRecPhiPosEtaNeg",
- "#phi",180,0,2*TMath::Pi());
- fRecPhiPosEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
-
- fRecPhiNegEtaNeg = new TH1F("fRecPhiNegEtaNeg",
- "#phi",180,0,2*TMath::Pi());
- fRecPhiNegEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
-
-
-
- fRecDcaPhiPtPosEtaPos = new TH3F("fRecDcaPhiPtPosEtaPos",
- "#phi- p_{T} - DCA",
- 180,0,2*TMath::Pi(),
- 100,0,pt,
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaPhiPtPosEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
- fRecDcaPhiPtPosEtaPos->GetYaxis()->SetTitle("p_{T} (GeV/c)");
- fRecDcaPhiPtPosEtaPos->GetZaxis()->SetTitle("dca (cm)");
-
- fRecDcaPhiPtPosEtaNeg = new TH3F("fRecDcaPhiPtPosEtaNeg",
- "#phi- p_{T} - DCA",
- 180,0,2*TMath::Pi(),
- 100,0,pt,
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaPhiPtPosEtaNeg->GetZaxis()->SetTitle("dca (cm)");
- fRecDcaPhiPtPosEtaNeg->GetXaxis()->SetTitle("#phi (rad.)");
- fRecDcaPhiPtPosEtaNeg->GetYaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecDcaPhiPtNegEtaPos = new TH3F("fRecDcaPhiPtNegEtaPos",
- "#phi- p_{T} - DCA",
- 180,0,2*TMath::Pi(),
- 100,0,pt,
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaPhiPtNegEtaPos->GetZaxis()->SetTitle("dca (cm)");
- fRecDcaPhiPtNegEtaPos->GetXaxis()->SetTitle("#phi (rad.)");
- fRecDcaPhiPtNegEtaPos->GetYaxis()->SetTitle("p_{T} (GeV/c)");
-
- fRecDcaPhiPtNegEtaNeg = new TH3F("fRecDcaPhiPtNegEtaNeg",
- "#phi- p_{T} - DCA",
- 180,0,2*TMath::Pi(),
- 100,0,pt,
- 100, -range*(1+Int_t(fTrackType/2)*9), range*(1+Int_t(fTrackType/2)*9));
- fRecDcaPhiPtNegEtaNeg->GetZaxis()->SetTitle("dca (cm)");
- fRecDcaPhiPtNegEtaNeg->GetYaxis()->SetTitle("#phi (rad.)");
- fRecDcaPhiPtNegEtaNeg->GetXaxis()->SetTitle("p_{T} (GeV/c)");
-
-
-
-
- fHists->SetOwner();
-
- fHists->Add(fHistRECpt);
- fHists->Add(fEta);
- fHists->Add(fEtavPt);
- fHists->Add(fEtaPhi);
- fHists->Add(fThetaRec);
- fHists->Add(fPhiRec);
- fHists->Add(fNumber);
- fHists->Add(fVx);
- fHists->Add(fVy);
- fHists->Add(fVz);
-
- fHists->Add(fEtaPt);
- fHists->Add(fQPt);
- fHists->Add(fDca);
-
- fHists->Add(fDeltaPhiAll);
- fHists->Add(fDeltaPhiLeading);
- fHists->Add(fDiffDcaD);
-
- fHists->Add(fqRec);
- fHists->Add(fsigmaPt);
-
- fHists->Add(fRecPtPos);
- fHists->Add(fRecPtNeg);
- fHists->Add(fRecPhiPos);
- fHists->Add(fRecPhiNeg);
- fHists->Add(fRecEtaPos);
- fHists->Add(fRecEtaNeg);
- fHists->Add(fRecEtaPtPos);
- fHists->Add(fRecEtaPtNeg);
- fHists->Add(fRecDcaPos);
- fHists->Add(fRecDcaNeg);
- fHists->Add(fRecDcaNegInv);
- fHists->Add(fRecDPos);
- fHists->Add(fRecDNeg);
-
-
- fHists->Add(fRecQPtPosEta);
- fHists->Add(fRecQPtNegEta);
- fHists->Add(fRecPtPosEta);
- fHists->Add(fRecPtNegEta);
- fHists->Add(fRecPhiPosEta);
- fHists->Add(fRecPhiNegEta);
- fHists->Add(fRecDcaPosEta);
- fHists->Add(fRecDcaNegEta);
- fHists->Add(fRecDPosEta);
- fHists->Add(fRecDNegEta);
-
-
- for(Int_t i=0;i<18;i++){
- fHists->Add(fRecPtTpcSector[i]);
- fHists->Add(fRecEtaTpcSector[i]);
- fHists->Add(fSignedDcaTpcSector[i]);
- fHists->Add(fRecQPtTpcSector[i]);
- fHists->Add(fRecEtaPtTpcSector[i]);
- }
-
- for(Int_t i=0;i<7;i++){
- fHists->Add(fRecPtPosLadder[i]);
- fHists->Add(fRecPtNegLadder[i]);
- fHists->Add(fRecPhiPosLadder[i]);
- fHists->Add(fRecPhiNegLadder[i]);
- fHists->Add(fRecEtaPosLadder[i]);
- fHists->Add(fRecEtaNegLadder[i]);
- }
-
- fHists->Add(fRecPtPosVz);
- fHists->Add(fRecPtNegVz);
- fHists->Add(fRecEtaPosVz);
- fHists->Add(fRecEtaNegVz);
- fHists->Add(fRecPhiPosVz);
- fHists->Add(fRecPhiNegVz);
- fHists->Add(fSignedDcaPosVz);
- fHists->Add(fSignedDcaNegVz);
- fHists->Add(fRecQPtPosEtaVz);
- fHists->Add(fRecQPtNegEtaVz);
- fHists->Add(fRecEtaPtPosVz);
- fHists->Add(fRecEtaPtNegVz);
-
-
- for(Int_t i=0;i<7;i++){
- fHists->Add(fSignDcaPos[i]);
- fHists->Add(fSignDcaNeg[i]);
- fHists->Add(fSignDcaNegInv[i]);
-
- fHists->Add(fPtSigmaPos[i]);
- fHists->Add(fPtSigmaNeg[i]);
- fHists->Add(fqPtRec[i]);
-
- fHists->Add(fDcaSigmaPos[i]);
- fHists->Add(fDcaSigmaNeg[i]);
-
-
- }
-
- fHists->Add(fRecDcaPosPhi);
- fHists->Add(fRecDcaNegPhi);
- fHists->Add(fRecPtPosPhi);
- fHists->Add(fRecPtNegPhi);
- fHists->Add(fRecEtaPosPhi);
- fHists->Add(fRecEtaNegPhi);
- fHists->Add(fRecQPtPhi);
- fHists->Add(fRecEtaPtPosPhi);
- fHists->Add(fRecEtaPtNegPhi);
-
- fHists->Add(fRecPtPosEtaPos);
- fHists->Add(fRecPtNegEtaPos);
- fHists->Add(fRecPtPosEtaNeg);
- fHists->Add(fRecPtNegEtaNeg);
-
- fHists->Add(fRec1PtPosEtaPos);
- fHists->Add(fRec1PtNegEtaPos);
- fHists->Add(fRec1PtPosEtaNeg);
- fHists->Add(fRec1PtNegEtaNeg);
-
-
- fHists->Add(fRecPhiPosEtaPos);
- fHists->Add(fRecPhiNegEtaPos);
- fHists->Add(fRecPhiPosEtaNeg);
- fHists->Add(fRecPhiNegEtaNeg);
-
- fHists->Add(fRecDcaPosPhiEtaPos);
- fHists->Add(fRecDcaNegPhiEtaPos);
- fHists->Add(fRecPtPosPhiEtaPos);
- fHists->Add(fRecPtNegPhiEtaPos);
- fHists->Add(fRecDcaPosPhiEtaNeg);
- fHists->Add(fRecDcaNegPhiEtaNeg);
- fHists->Add(fRecPtPosPhiEtaNeg);
- fHists->Add(fRecPtNegPhiEtaNeg);
-
- fHists->Add(fRecDcaPhiPtPosEtaPos);
- fHists->Add(fRecDcaPhiPtPosEtaNeg);
- fHists->Add(fRecDcaPhiPtNegEtaPos);
- fHists->Add(fRecDcaPhiPtNegEtaNeg);
-
-
-
-
-
-// for (Int_t i=0; i<fHists->GetEntries(); ++i) {
-// TH1 *h1 = dynamic_cast<TH1*>(fHists->At(i));
-// if (h1){
-// // Printf("%s ",h1->GetName());
-// h1->Sumw2();
-// }
-// }
- // BKC
-
- TH1::AddDirectory(oldStatus);
-}
-
-//__________________________________________________________
-
-void AliAnalysisTaskQASym::UserExec(Option_t *)
-{
- AliVEvent *event = InputEvent();
- if (!event) {
- Printf("ERROR: Could not retrieve event");
- return;
- }
-
-
- if(Entry()==0){
- AliESDEvent* esd = dynamic_cast<AliESDEvent*>(event);
- if(esd){
- Printf("We are reading from ESD");
- }
-
- }
-
-
-
- if(fDebug>1)Printf("There are %d tracks in this event", event->GetNumberOfTracks());
-
-
- Int_t leadingTrack = 0;
- Float_t leadingEnergy = -20.;
- Float_t leadingPhi = 0;//TMath::Pi();
-
-
- if(event->GetNumberOfTracks()!=0) fNumber->Fill(event->GetNumberOfTracks());
-
- const AliVVertex* vertex = event->GetPrimaryVertex();
- Float_t vz = vertex->GetZ();
- if (TMath::Abs(vz) > 10.) return;
-
- for (Int_t iTrack = 0; iTrack < event->GetNumberOfTracks(); iTrack++) {
-
-
- AliVParticle *track = event->GetTrack(iTrack);
- AliESDtrack *esdtrack = dynamic_cast<AliESDtrack*>(track);
- esdtrack->PropagateToDCA(event->GetPrimaryVertex(),
- event->GetMagneticField(), 10000.);
-
- if (!track) {
- Printf("ERROR: Could not receive track %d", iTrack);
- continue;
- }
- //__________
- // run Task for global tracks or ITS tracks or TPC tracks
- AliESDtrack *tpcP = 0x0;
- const AliExternalTrackParam *tpcPin = 0x0;
- Double_t phiIn=0;
- if(fTrackType==0){
- //Fill all histograms with global tracks
- tpcP = esdtrack;
- phiIn = tpcP->Phi();
- if (!tpcP) continue;
- if (!fCuts->AcceptTrack(tpcP)) continue;
- }
- else if(fTrackType==1){
- //Fill all histograms with ITS tracks
- tpcP = esdtrack;
- phiIn = tpcP->Phi();
- if (!tpcP) continue;
- if (!fCuts->AcceptTrack(tpcP)) continue;
- }
- else if(fTrackType==2){
- //Fill all histograms with TPC track information
- tpcPin = esdtrack->GetInnerParam();
- if (!tpcPin) continue;
- phiIn=tpcPin->Phi();
-
- tpcP = AliESDtrackCuts::GetTPCOnlyTrack(dynamic_cast<AliESDEvent*>(event),esdtrack->GetID());
- if (!tpcP) continue;
- if (!fCuts->AcceptTrack(tpcP)) continue;
- if(tpcP->GetNcls(1)>160)continue;//jacek's track cut
- if(tpcP->GetConstrainedChi2TPC()<0)continue; // jacek's track cut
- }
- else{
- Printf("ERROR: wrong track type \n");
- continue;
- }
- //___________
- //
-
-
- if(tpcP->E()>leadingEnergy){
- leadingTrack=iTrack;
- leadingEnergy=tpcP->E();
- leadingPhi=phiIn;
- }
-
-
- fqRec->Fill(tpcP->Charge());
-
-
- Double_t sigmapt = tpcP->GetSigma1Pt2();
- sigmapt= sqrt(sigmapt);
- sigmapt= sigmapt *(tpcP->Pt()*tpcP->Pt());
-
- if(sigmapt == 0.)continue;
- fsigmaPt->Fill(TMath::Log10(sigmapt));
-
-
- // hits in ITS layer
- Int_t cas=-1;
- if(tpcP->HasPointOnITSLayer(0))
- cas=0;
- else if(!tpcP->HasPointOnITSLayer(0)
- && tpcP->HasPointOnITSLayer(1))
- cas=1;
- else if(!tpcP->HasPointOnITSLayer(0)
- && !tpcP->HasPointOnITSLayer(1)
- && tpcP->HasPointOnITSLayer(2))
- cas=2;
- else if(!tpcP->HasPointOnITSLayer(0)
- && !tpcP->HasPointOnITSLayer(1)
- && !tpcP->HasPointOnITSLayer(2)
- && tpcP->HasPointOnITSLayer(3))
- cas=3;
- else if(!tpcP->HasPointOnITSLayer(0)
- && !tpcP->HasPointOnITSLayer(1)
- && !tpcP->HasPointOnITSLayer(2)
- && !tpcP->HasPointOnITSLayer(3)
- && tpcP->HasPointOnITSLayer(4))
- cas=4;
- else if( !tpcP->HasPointOnITSLayer(0)
- && !tpcP->HasPointOnITSLayer(1)
- && !tpcP->HasPointOnITSLayer(2)
- && !tpcP->HasPointOnITSLayer(3)
- && !tpcP->HasPointOnITSLayer(4)
- && tpcP->HasPointOnITSLayer(5))
- cas=5;
- else
- cas=6;
-
-
-
- //-------------------
-
- xvertexcor = tpcP->Xv() - vertex->GetX(); // coordinate corrected for vertex position
- yvertexcor = tpcP->Yv() - vertex->GetY(); // "
- sdca = (tpcP->Py()*xvertexcor - tpcP->Px()*yvertexcor)/tpcP->Pt();
-
-
- fqPtRec[cas]->Fill(tpcP->Charge()/tpcP->Pt());
-
-
-
- fHistRECpt->Fill(tpcP->Pt());
- fEta->Fill(tpcP->Eta());
- fEtavPt->Fill(tpcP->Eta(), tpcP->Pt());
- fEtaPhi->Fill(tpcP->Eta(), phiIn);
- fThetaRec->Fill(tpcP->Theta());
- fPhiRec->Fill(phiIn);
- fVx->Fill(tpcP->Xv());
- fVy->Fill(tpcP->Yv());
- fVz->Fill(tpcP->Zv());
-
-
- fEtaPt->Fill(tpcP->Eta()/tpcP->Pt());
- fQPt->Fill(tpcP->Charge()/tpcP->Pt());
- fDca->Fill(sdca);
- fRecQPtPhi->Fill(tpcP->Charge()/tpcP->Pt(), phiIn);
-
-
- tpcP->GetImpactParameters(xy,z);
- fDiffDcaD->Fill(sdca+xy);
-
-
-
- //for positive particles
-
- if(tpcP->Charge()>0){
- fRecPtPos->Fill(tpcP->Pt());
- fRecPtPosLadder[cas]->Fill(tpcP->Pt());
- fRecPtPosVz->Fill(TMath::Log10(tpcP->Pt()),tpcP->Zv());
- fRecPhiPos->Fill(TMath::RadToDeg()*phiIn);
-
-
- fRecPhiPosLadder[cas]->Fill(TMath::RadToDeg()*phiIn);
- fRecPhiPosVz->Fill(TMath::RadToDeg()*phiIn,tpcP->Zv());
- fSignedDcaPosVz->Fill(sdca,tpcP->Zv());
-
- fRecEtaPos->Fill(tpcP->Eta());
- fRecEtaPosLadder[cas]->Fill(tpcP->Eta());
- fRecEtaPtPos->Fill(tpcP->Eta()/tpcP->Pt());
- fRecEtaPosVz->Fill(tpcP->Eta(),tpcP->Zv());
- fRecEtaPtPosVz->Fill(tpcP->Eta()/tpcP->Pt(),tpcP->Zv());
-
- fRecDcaPos->Fill(sdca);
- fRecDcaPosPhi->Fill(sdca, phiIn);
- fRecPtPosPhi->Fill(TMath::Log10(tpcP->Pt()), phiIn);
- fRecEtaPtPosPhi->Fill(tpcP->Eta()/tpcP->Pt(), phiIn);
- fRecEtaPosPhi->Fill(tpcP->Eta(), phiIn);
- fRecDPos->Fill(xy);
- fSignDcaPos[cas]->Fill(sdca);
-
-
- fDcaSigmaPos[cas]->Fill(sdca, TMath::Log10(sigmapt));
-
- fPtSigmaPos[cas]->Fill(TMath::Log10(sigmapt));
- //pos eta
- if(tpcP->Eta()>0){
- fRecPtPosEtaPos->Fill(tpcP->Pt());
- fRec1PtPosEtaPos->Fill(1/tpcP->Pt());
- fRecPhiPosEtaPos->Fill(phiIn);
- fRecDcaPosPhiEtaPos->Fill(sdca, phiIn);
- fRecPtPosPhiEtaPos->Fill(TMath::Log10(tpcP->Pt()), phiIn);
- fRecDcaPhiPtPosEtaPos->Fill(phiIn, tpcP->Pt(), sdca);
- }
- //neg eta
- else{
- fRecPtPosEtaNeg->Fill(tpcP->Pt());
- fRec1PtPosEtaNeg->Fill(1/tpcP->Pt());
- fRecPhiPosEtaNeg->Fill(phiIn);
- fRecDcaPosPhiEtaNeg->Fill(sdca, phiIn);
- fRecPtPosPhiEtaNeg->Fill(TMath::Log10(tpcP->Pt()), phiIn);
- fRecDcaPhiPtPosEtaNeg->Fill(phiIn, tpcP->Pt(), sdca);
- }
-
- }
- //and negative particles
- else {
- fRecPtNeg->Fill(tpcP->Pt());
- fRecPtNegLadder[cas]->Fill(tpcP->Pt());
- fRecPtNegVz->Fill(TMath::Log10(tpcP->Pt()),tpcP->Zv());
-
- fRecPhiNeg->Fill(TMath::RadToDeg()*phiIn);
- fRecPhiNegLadder[cas]->Fill(TMath::RadToDeg()*phiIn);
- fRecPhiNegVz->Fill(TMath::RadToDeg()*phiIn,tpcP->Zv());
- fSignedDcaNegVz->Fill(sdca,tpcP->Zv());
- fRecEtaPtNegVz->Fill(tpcP->Eta()/tpcP->Pt(),tpcP->Zv());
-
- fRecEtaNeg->Fill(tpcP->Eta());
- fRecEtaNegLadder[cas]->Fill(tpcP->Eta());
- fRecEtaPtNeg->Fill(tpcP->Eta()/tpcP->Pt());
- fRecEtaNegVz->Fill(tpcP->Eta(),tpcP->Zv());
-
- fRecDcaNeg->Fill(sdca);
- fRecDcaNegInv->Fill(-sdca);
- fRecDcaNegPhi->Fill(sdca, phiIn);
- fRecPtNegPhi->Fill(TMath::Log10(tpcP->Pt()), phiIn);
- fRecEtaNegPhi->Fill(tpcP->Eta(), phiIn);
- fRecEtaPtNegPhi->Fill(tpcP->Eta()/tpcP->Pt(), phiIn);
- fRecDNeg->Fill(xy);
- fSignDcaNeg[cas]->Fill(sdca);
- fSignDcaNegInv[cas]->Fill(-sdca);
-
-
- fDcaSigmaNeg[cas]->Fill(sdca,TMath::Log10(sigmapt));
-
- fPtSigmaNeg[cas]->Fill(TMath::Log10(sigmapt));
-
- //pos eta
- if(tpcP->Eta()>0){
- fRecPtNegEtaPos->Fill(tpcP->Pt());
- fRec1PtNegEtaPos->Fill(1/tpcP->Pt());
- fRecPhiNegEtaPos->Fill(phiIn);
- fRecDcaNegPhiEtaPos->Fill(sdca, phiIn);
- fRecPtNegPhiEtaPos->Fill(TMath::Log10(tpcP->Pt()), phiIn);
- fRecDcaPhiPtNegEtaPos->Fill(phiIn, tpcP->Pt(), sdca);
- }
- //neg eta
- else{
- fRecPtNegEtaNeg->Fill(tpcP->Pt());
- fRec1PtNegEtaNeg->Fill(1/tpcP->Pt());
- fRecPhiNegEtaNeg->Fill(phiIn);
- fRecDcaNegPhiEtaNeg->Fill(sdca, phiIn);
- fRecPtNegPhiEtaNeg->Fill(TMath::Log10(tpcP->Pt()), phiIn);
- fRecDcaPhiPtNegEtaNeg->Fill(phiIn, tpcP->Pt(), sdca);
- }
-
- }
-
-
-
- //all particles with positive eta
- if(tpcP->Eta()>0){
- fRecQPtPosEta->Fill(tpcP->Charge()/tpcP->Pt());
- fRecPtPosEta->Fill(tpcP->Pt());
- fRecPhiPosEta->Fill(TMath::RadToDeg()*phiIn);
- fRecQPtPosEtaVz->Fill(tpcP->Charge()/tpcP->Pt(),tpcP->Zv());
- fRecDcaPosEta->Fill(sdca);
- fRecDPosEta->Fill(xy);
- }
- //all particles with negative eta (and eta==0)
- else{
- fRecQPtNegEta->Fill(tpcP->Charge()/tpcP->Pt());
- fRecPtNegEta->Fill(tpcP->Pt());
- fRecPhiNegEta->Fill(TMath::RadToDeg()*phiIn);
- fRecQPtNegEtaVz->Fill(tpcP->Charge()/tpcP->Pt(),tpcP->Zv());
- fRecDcaNegEta->Fill(sdca);
- fRecDNegEta->Fill(xy);
-
- }
-
-
- fRecPtTpcSector[Int_t(phiIn*
- TMath::RadToDeg()/20)]->Fill(tpcP->Pt());
- fRecEtaTpcSector[Int_t(phiIn*
- TMath::RadToDeg()/20)]->Fill(tpcP->Eta());
- fSignedDcaTpcSector[Int_t(phiIn*
- TMath::RadToDeg()/20)]->Fill(sdca);
- fRecQPtTpcSector[Int_t(phiIn*
- TMath::RadToDeg()/20)]->Fill(tpcP->Charge()/tpcP->Pt());
- fRecEtaPtTpcSector[Int_t(phiIn*
- TMath::RadToDeg()/20)]->Fill(tpcP->Eta()/tpcP->Pt());
-
-
-
-// // another track loop
-// for (Int_t iTrack2 = 0; iTrack2 < event->GetNumberOfTracks(); iTrack2++) {
-
-// if(LeadingTrack==iTrack2) continue;
-
-// AliVParticle *track2 = event->GetTrack(iTrack2);
-// AliESDtrack* esdtrack2 = dynamic_cast<AliESDtrack*>(track2);
-// if (!track2) {
-// Printf("ERROR: Could not receive track %d", iTrack);
-// continue;
-// }
-// if (!fCuts->AcceptTrack(esdtrack2)) continue;
-// //propagate to dca
-// esdtrack2->PropagateToDCA(event->GetPrimaryVertex(),
-// event->GetMagneticField(), 10000.);
-
-// fDeltaPhiLeading->Fill((LeadingPhi-esdtrack2->Phi())*TMath::RadToDeg(),
-// LeadingPhi*TMath::RadToDeg() );
-
-
-
-// }//second track loop
- }//first track loop
-
-
-
-
- // Post output data.
- // PostData(1, fHistPt);
- PostData(1, fHists);
-}
-
-
-
-
-
-//________________________________________________________________________
-void AliAnalysisTaskQASym::Terminate(Option_t *)
-{
-
-
-}
-
-
-
-
-
+++ /dev/null
-#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
+++ /dev/null
-{
-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();
-}
-
-
-
-}
+++ /dev/null
-{
-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");
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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()");
-
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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}");
-
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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()");
-
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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()");
-
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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()");
-
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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()");
-
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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("");
-
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-{
-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();
-
-}
-
-
-
-}
+++ /dev/null
-//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);
-
- }
-
-
-}
-
-
-
+++ /dev/null
-//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);
-
-
-
- }
-
-
-
-}
+++ /dev/null
-
-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();
-}