+++ /dev/null
-AliAnalysisTask *AddTask_GammaConvdPhi_PbPb(TString v0Cut = "1090002002092970023220000000",
- TString pionCut = "01522045009000",
- TString photoncut = "",
- Bool_t pbpb = kTRUE) {
-
- ////////////////////CURRENTLY NOT WORKING ///////////////////////////////
-
- // standard with task
- printf("========================================================================================\n");
- printf("dPhiAnalysis: Initialising AliAnalysisTaskdPhi\n");
- printf("========================================================================================\n");
-
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- if (!mgr) {
- Error("AddTask_slindal_dPhi", "No analysis manager found.");
- return 0;
- }
-
- AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
- Bool_t isAOD=man->GetInputEventHandler()->IsA()==AliAODInputHandler::Class();
-
- AliAnalysisTaskdPhi *task = new AliAnalysisTaskdPhi((TString("slindalTask_dPhi")+"_" + v0Cut));
-
- ///Axes for histrograms
- Double_t cptbins[6] = {0.5, 1.0, 1.5, 3.0, 5.0, 10};
- task->GetAxiscPt().Set(5, cptbins);
-
- Double_t tptbins[3] = {3.0, 5.0, 10.0};
- task->GetAxistPt().Set(2, tptbins);
-
- if(pbpb) {
- Double_t centBins[7] = {0, 5, 10, 30, 50, 60, 90};
- task->GetAxisCent().Set(6, centBins);
- } else {
- Double_t centBins[2] = {-9999, 9999};
- task->GetAxisCent().Set(1, centBins);
- }
-
- Double_t zbins[6] = { -10, -5, -1.5, 1.5, 5, 10};
- task->GetAxisZ().Set(5, zbins);
-
- Double_t mbins[17] = {0.07, 0.09, 0.1, 0.11, 0.12, 0.125, 0.1275, 0.13, 0.14, 0.1425, 0.145, 0.15, 0.16, 0.18, 0.2, 0.24, 0.26};
- task->GetAxisPiMass().Set(16, mbins);
-
- AliESDtrackCuts * trackCuts = new AliESDtrackCuts();
- TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
- trackCuts->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep, 100);
- trackCuts->SetMaxChi2PerClusterTPC(4);
- trackCuts->SetRequireTPCStandAlone(kTRUE);
- trackCuts->SetAcceptKinkDaughters(kFALSE);
- trackCuts->SetRequireTPCRefit(kTRUE);
- trackCuts->SetMaxFractionSharedTPCClusters(0.4);
-
- trackCuts->SetMaxDCAToVertexXY(2.4);
- trackCuts->SetMaxDCAToVertexZ(3.2);
- trackCuts->SetDCAToVertex2D(kTRUE);
-
- trackCuts->SetMaxChi2PerClusterITS(36);
- trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
-
- trackCuts->SetRequireSigmaToVertex(kFALSE);
-
- trackCuts->SetEtaRange(-0.9, 0.9);
- trackCuts->SetPtRange(cptbins[0], 1000000.0);
-
- trackCuts->SetRequireITSRefit(kFALSE);
- task->SetTrackFilter(trackCuts);
-
- ///open side nclusters graph
- AliESDtrackCuts * trackCuts2 = new AliESDtrackCuts();
- TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","65.+30./20.*x");
- trackCuts2->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep, 100);
- trackCuts2->SetMaxChi2PerClusterTPC(4);
- trackCuts2->SetRequireTPCStandAlone(kTRUE);
- trackCuts2->SetAcceptKinkDaughters(kFALSE);
- trackCuts2->SetRequireTPCRefit(kTRUE);
- trackCuts2->SetMaxFractionSharedTPCClusters(0.4);
-
- trackCuts2->SetMaxDCAToVertexXY(2.4);
- trackCuts2->SetMaxDCAToVertexZ(3.2);
- trackCuts2->SetDCAToVertex2D(kTRUE);
-
- trackCuts2->SetMaxChi2PerClusterITS(36);
- trackCuts2->SetMaxChi2TPCConstrainedGlobal(36);
-
- trackCuts2->SetRequireSigmaToVertex(kFALSE);
-
- trackCuts2->SetEtaRange(-0.9, 0.9);
- trackCuts2->SetPtRange(cptbins[0], 1000000.0);
-
- trackCuts2->SetRequireITSRefit(kFALSE);
- task->AddTrackFilter(trackCuts2, kTRUE);
-
- ///Tight side nclusters graphs
- AliESDtrackCuts * trackCuts3 = new AliESDtrackCuts();
- TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","75.+30./20.*x");
- trackCuts3->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep, 100);
- trackCuts3->SetMaxChi2PerClusterTPC(4);
- trackCuts3->SetRequireTPCStandAlone(kTRUE);
- trackCuts3->SetAcceptKinkDaughters(kFALSE);
- trackCuts3->SetRequireTPCRefit(kTRUE);
- trackCuts3->SetMaxFractionSharedTPCClusters(0.4);
-
- trackCuts3->SetMaxDCAToVertexXY(2.4);
- trackCuts3->SetMaxDCAToVertexZ(3.2);
- trackCuts3->SetDCAToVertex2D(kTRUE);
-
- trackCuts3->SetMaxChi2PerClusterITS(36);
- trackCuts3->SetMaxChi2TPCConstrainedGlobal(36);
-
- trackCuts3->SetRequireSigmaToVertex(kFALSE);
-
- trackCuts3->SetEtaRange(-0.9, 0.9);
- trackCuts3->SetPtRange(cptbins[0], 1000000.0);
-
- trackCuts3->SetRequireITSRefit(kFALSE);
- task->AddTrackFilter(trackCuts3, kFALSE);
-
-
- ////Tight xyz track cuts
- AliESDtrackCuts4 * trackCuts4 = new AliESDtrackCuts();
- TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
- trackCuts4->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep, 100);
- trackCuts4->SetMaxChi2PerClusterTPC(4);
- trackCuts4->SetRequireTPCStandAlone(kTRUE);
- trackCuts4->SetAcceptKinkDaughters(kFALSE);
- trackCuts4->SetRequireTPCRefit(kTRUE);
- trackCuts4->SetMaxFractionSharedTPCClusters(0.4);
-
- trackCuts4->SetMaxDCAToVertexXY(2.2);
- trackCuts4->SetMaxDCAToVertexZ(3.0);
- trackCuts4->SetDCAToVertex2D(kTRUE);
-
- trackCuts4->SetMaxChi2PerClusterITS(36);
- trackCuts4->SetMaxChi2TPCConstrainedGlobal(36);
-
- trackCuts4->SetRequireSigmaToVertex(kFALSE);
-
- trackCuts4->SetEtaRange(-0.9, 0.9);
- trackCuts4->SetPtRange(cptbins[0], 1000000.0);
-
- trackCuts4->SetRequireITSRefit(kFALSE);
- task->AddTrackFilter(trackCuts4, kFALSE);
-
-
- ////lose xyz track cuts
- AliESDtrackCuts5 * trackCuts5 = new AliESDtrackCuts();
- TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
- trackCuts5->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep, 100);
- trackCuts5->SetMaxChi2PerClusterTPC(4);
- trackCuts5->SetRequireTPCStandAlone(kTRUE);
- trackCuts5->SetAcceptKinkDaughters(kFALSE);
- trackCuts5->SetRequireTPCRefit(kTRUE);
- trackCuts5->SetMaxFractionSharedTPCClusters(0.4);
-
- trackCuts5->SetMaxDCAToVertexXY(2.6);
- trackCuts5->SetMaxDCAToVertexZ(3.6);
- trackCuts5->SetDCAToVertex2D(kTRUE);
-
- trackCuts5->SetMaxChi2PerClusterITS(36);
- trackCuts5->SetMaxChi2TPCConstrainedGlobal(36);
-
- trackCuts5->SetRequireSigmaToVertex(kFALSE);
-
- trackCuts5->SetEtaRange(-0.9, 0.9);
- trackCuts5->SetPtRange(cptbins[0], 1000000.0);
-
- trackCuts5->SetRequireITSRefit(kFALSE);
- task->AddTrackFilter(trackCuts5);
-
-
-
- ///open side shared tpt clusters
- AliESDtrackCuts * trackCuts2 = new AliESDtrackCuts();
- TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
- trackCuts2->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep, 100);
- trackCuts2->SetMaxChi2PerClusterTPC(4);
- trackCuts2->SetRequireTPCStandAlone(kTRUE);
- trackCuts2->SetAcceptKinkDaughters(kFALSE);
- trackCuts2->SetRequireTPCRefit(kTRUE);
- trackCuts2->SetMaxFractionSharedTPCClusters(0.44);
-
- trackCuts2->SetMaxDCAToVertexXY(2.4);
- trackCuts2->SetMaxDCAToVertexZ(3.2);
- trackCuts2->SetDCAToVertex2D(kTRUE);
-
- trackCuts2->SetMaxChi2PerClusterITS(36);
- trackCuts2->SetMaxChi2TPCConstrainedGlobal(36);
-
- trackCuts2->SetRequireSigmaToVertex(kFALSE);
-
- trackCuts2->SetEtaRange(-0.9, 0.9);
- trackCuts2->SetPtRange(cptbins[0], 1000000.0);
-
- trackCuts2->SetRequireITSRefit(kFALSE);
- task->AddTrackFilter(trackCuts2, kTRUE);
-
-
- ///tight side shared tpt clusters
- AliESDtrackCuts * trackCuts2 = new AliESDtrackCuts();
- TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
- trackCuts2->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep, 100);
- trackCuts2->SetMaxChi2PerClusterTPC(4);
- trackCuts2->SetRequireTPCStandAlone(kTRUE);
- trackCuts2->SetAcceptKinkDaughters(kFALSE);
- trackCuts2->SetRequireTPCRefit(kTRUE);
- trackCuts2->SetMaxFractionSharedTPCClusters(0.36);
-
- trackCuts2->SetMaxDCAToVertexXY(2.4);
- trackCuts2->SetMaxDCAToVertexZ(3.2);
- trackCuts2->SetDCAToVertex2D(kTRUE);
-
- trackCuts2->SetMaxChi2PerClusterITS(36);
- trackCuts2->SetMaxChi2TPCConstrainedGlobal(36);
-
- trackCuts2->SetRequireSigmaToVertex(kFALSE);
-
- trackCuts2->SetEtaRange(-0.9, 0.9);
- trackCuts2->SetPtRange(cptbins[0], 1000000.0);
-
- trackCuts2->SetRequireITSRefit(kFALSE);
- task->AddTrackFilter(trackCuts2, kTRUE);
-
-
-
- ///Pizero cuts
- pionString = "01630031009000";
- TString * picutarr[2][4];
-
- picutarr[0][
-
-
-
-
-
- AliConversionMesonCuts * picuts = new AliConversionMesonCuts();
- picuts->InitializeCutsFromCutString(pionString);
- task->SetMesonFilter(picuts);
-
-
- AliConversionMesonCuts * picuts2 = new AliConversionMesonCuts();
- picuts2->InitializeCutsFromCutString(pionString2);
- task->AddMesonFilter(picuts2, kFALSE);
-
-
- AliConversionMesonCuts * picuts3 = new AliConversionMesonCuts();
- picuts3->InitializeCutsFromCutString(pionString3);
- task->AddMesonFilter(picuts3);
-
-
-
-
- ///Pion cuts
- AliConversionMesonCuts * picuts = new AliConversionMesonCuts("dphi_pioncuts");
- picuts->InitializeCutsFromCutString(pionCut);
- task->SetMesonFilter(picuts);
-
- if(photoncut.Length() > 0) {
- ///V0 analysis cuts (applied before pion analysis)
- AliConversionCuts * gcuts = new AliConversionCuts();
- gcuts->InitializeCutsFromCutString(photoncut);
- task->SetPhotonFilter(gcuts);
- }
-
- //================================================
- // data containers
- //================================================
- // find input container
- //below the trunk version
- AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
- TString baseString("slindal_dPhi_");
- baseString += v0Cut + "_" + pionCut;
-
- //dumm output container
- AliAnalysisDataContainer *coutput0 =
- mgr->CreateContainer(baseString+"tree",
- TTree::Class(),
- AliAnalysisManager::kExchangeContainer,
- "slindal_default");
-
- //define output containers, please use 'username'_'somename'
- AliAnalysisDataContainer *coutput1 =
- mgr->CreateContainer(baseString+"me", TList::Class(),
- AliAnalysisManager::kOutputContainer,baseString+".root");
-
- //define output containers, please use 'username'_'somename'
-
-
- //========= Add PID Reponse to ANALYSIS manager ====
- Bool_t isMC = kFALSE;
- if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
- gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
- AddTaskPIDResponse(isMC);
- }
-
- //========= Set Cutnumber for V0Reader ================================
- // TString cutnumber = "100000000008400100150000000";
- //========= Add V0 Reader to ANALYSIS manager if not yet existent =====
- AliV0ReaderV1 * fV0ReaderV1 = (AliV0ReaderV1*)mgr->GetTask("V0ReaderV1");
- if( !fV0reader ){
- AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
-
- fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
- fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
- fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
-
- if (!mgr) {
- Error("AddTask_V0ReaderV1", "No analysis manager found.");
- return;
- }
-
- // Set AnalysisCut Number
- AliConversionCuts *fCuts=NULL;
- if(v0Cut!=""){
- fCuts= new AliConversionCuts(v0Cut.Data(),v0Cut.Data());
- if(pbpb) {
- fCuts->SetPreSelectionCutFlag(kFALSE);
- } else {
- fCuts->SetPreSelectionCutFlag(kTRUE);
- }
- if(fCuts->InitializeCutsFromCutString(v0Cut)){
- fV0ReaderV1->SetConversionCuts(fCuts);
- //fCuts->SetFillCutHistograms("",kTRUE);
- }
-
-
- }
-
- fV0ReaderV1->Init();
- mgr->AddTask(fV0ReaderV1);
- mgr->ConnectInput(fV0ReaderV1,0,cinput);
- task->SetV0Reader(fV0ReaderV1);
- task->SaveReaderHists();
- } else {
- ///V0 analysis cuts (applied before pion analysis)
- AliConversionCuts * v0cuts = new AliConversionCuts();
- v0cuts->InitializeCutsFromCutString(v0Cut);
- task->SetV0Reader(fV0ReaderV1);
- task->SetV0Filter(v0cuts);
- }
-
- task->SelectCollisionCandidates(AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral);
-
- mgr->AddTask(task);
- mgr->ConnectInput (task, 0, cinput );
- mgr->ConnectOutput (task, 0, coutput0);
- mgr->ConnectOutput (task, 1, coutput1);
-
- return task;
-}
-
-
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff