//DEFINITION OF A FEW CONSTANTS
-const Float_t ptmin = 2.0 ;
-const Float_t ptmax = 50.0 ;
const Float_t phimin = 0.;
const Float_t phimax = 2.*TMath::Pi();
const Float_t etamin = -0.9;
const Float_t etamax = 0.9;
-const Float_t dcarmin = -0.2;
-const Float_t dcarmax = 0.2;
-const Float_t chi2TPCmin = 0.0;
-const Float_t chi2TPCmax = 4.;
-
-const Float_t ptmin1 = ptmin ;
-const Float_t ptmax1 = 10.0 ;
-const Float_t ptmin2 = ptmax1 ;
-const Float_t ptmax2 = 20.0 ;
-const Float_t ptmin3 = ptmax2 ;
-const Float_t ptmax3 = ptmax ;
const Int_t mintrackrefsTPC = 1;
const Int_t mintrackrefsITS = 1;
-const Int_t charge = 1;
-AliPWG4HighPtSpectra* AddTaskPWG4HighPtSpectra()
+AliPWG4HighPtSpectra* AddTaskPWG4HighPtSpectra(char *prodType = "LHC10e14",Int_t trackType = 0)
{
// Creates HighPtSpectra analysis task and adds it to the analysis manager.
UInt_t ipt = 0;
UInt_t iphi = 1;
UInt_t ieta = 2;
- UInt_t idcar = 3;
- UInt_t ichi2TPC = 4;
//Setting up the container grid...
- UInt_t nstep = 5; //Steps/Modes for containers
+ UInt_t nstep = 6; //Steps/Modes for containers
Int_t kStepReconstructed = 0;
Int_t kStepReconstructedTPCOnly = 1;
Int_t kStepSecondaries = 2;
Int_t kStepReconstructedMC = 3;
Int_t kStepMCAcceptance = 4;
-
- const Int_t nvar = 5; //number of variables on the grid: pt:phi:eta:DCAR:Chi2/NClusTPC
- const Int_t nbin11 = (int)(ptmax1-ptmin1);
- const Int_t nbin12 = (int)((ptmax2-ptmin2)/2.)+nbin11;
- const Int_t nbin13 = (int)((ptmax3-ptmin3)/5.)+nbin12;
+ Int_t kStepReconstructedTPCOnlyMC = 5;
+
+ //redefine pt ranges in case of Jet-Jet production
+ Float_t ptBinEdges[2][2];
+ Float_t ptmin = 2.0 ;
+ Float_t ptmax = 50.0 ;
+ Float_t binWidth3 = 5.;
+ if(!strcmp(prodType, "LHC10e14")) {
+ ptmin = 0.0 ;
+ ptmax = 500.0 ;
+
+ ptBinEdges[0][0] = 100.;
+ ptBinEdges[0][1] = 5.;
+ ptBinEdges[1][0] = 300.;
+ ptBinEdges[1][1] = 10.;
+ binWidth3 = 20.;
+ } else {
+ ptmin = 2.0 ;
+ ptmax = 50.0 ;
+
+ ptBinEdges[0][0] = 10.;
+ ptBinEdges[0][1] = 0.5;
+ ptBinEdges[1][0] = 60.;
+ ptBinEdges[1][1] = 2.;
+ binWidth3 = 5.;
+ }
+
+ const Int_t nvar = 3; //number of variables on the grid: pt:phi:eta
+ const Int_t nbin11 = (int)((ptBinEdges[0][0]-ptmin)/ptBinEdges[0][1]);
+ const Int_t nbin12 = (int)((ptBinEdges[1][0]-ptBinEdges[0][0])/ptBinEdges[1][1])+nbin11;
+ const Int_t nbin13 = (int)((ptmax-ptBinEdges[1][0])/binWidth3)+nbin12;
const Int_t nbin1 = nbin13; //bins in pt
- const Int_t nbin2 = 18;//36; //bins in phi
- const Int_t nbin3 = 8; //bins in eta
- const Int_t nbin4 = 40; //bins in DCAR
- const Int_t nbin5 = 20; //bins in Chi2/#NclusTPC
+ const Int_t nbin2 = 18; //bins in phi
+ const Int_t nbin3 = 2; //bins in eta
//arrays for the number of bins in each dimension
Int_t iBin[nvar];
iBin[0]=nbin1;
iBin[1]=nbin2;
iBin[2]=nbin3;
- iBin[3]=nbin4;
- iBin[4]=nbin5;
-
+
//arrays for lower bounds :
Double_t *binLim1=new Double_t[nbin1+1];
Double_t *binLim2=new Double_t[nbin2+1];
Double_t *binLim3=new Double_t[nbin3+1];
- Double_t *binLim4=new Double_t[nbin4+1];
- Double_t *binLim5=new Double_t[nbin5+1];
-
+
//values for bin lower bounds
for(Int_t i=0; i<=nbin1; i++) {
- if(i<=nbin11) binLim1[i]=(Double_t)ptmin1 + (ptmax1-ptmin1)/nbin11*(Double_t)i ;
- if(i<=nbin12 && i>nbin11) binLim1[i]=(Double_t)ptmin2 + (ptmax2-ptmin2)/(nbin12-nbin11)*((Double_t)i-(Double_t)nbin11) ;
- if(i<=nbin13 && i>nbin12) binLim1[i]=(Double_t)ptmin3 + (ptmax3-ptmin3)/(nbin13-nbin12)*((Double_t)i-(Double_t)nbin12) ;
+ if(i<=nbin11) binLim1[i]=(Double_t)ptmin + (ptBinEdges[0][0]-ptmin)/nbin11*(Double_t)i ;
+ if(i<=nbin12 && i>nbin11) binLim1[i]=(Double_t)ptBinEdges[0][0] + (ptBinEdges[1][0]-ptBinEdges[0][0])/(nbin12-nbin11)*((Double_t)i-(Double_t)nbin11) ;
+ if(i<=nbin13 && i>nbin12) binLim1[i]=(Double_t)ptBinEdges[1][0] + (ptmax-ptBinEdges[1][0])/(nbin13-nbin12)*((Double_t)i-(Double_t)nbin12) ;
}
- // for(Int_t i=0; i<=nbin1; i++) binLim1[i]=(Double_t)ptmin + (ptmax-ptmin)/nbin1*(Double_t)i ;
for(Int_t i=0; i<=nbin2; i++) binLim2[i]=(Double_t)phimin + (phimax-phimin)/nbin2*(Double_t)i ;
for(Int_t i=0; i<=nbin3; i++) binLim3[i]=(Double_t)etamin + (etamax-etamin)/nbin3*(Double_t)i ;
- for(Int_t i=0; i<=nbin4; i++) binLim4[i]=(Double_t)dcarmin + (dcarmax-dcarmin)/nbin4*(Double_t)i ;
- for(Int_t i=0; i<=nbin5; i++) binLim5[i]=(Double_t)chi2TPCmin + (chi2TPCmax-chi2TPCmin)/nbin5*(Double_t)i ;
AliCFContainer* containerPos = new AliCFContainer("containerPos","container for positive tracks",nstep,nvar,iBin);
containerPos -> SetBinLimits(ipt,binLim1);
containerPos -> SetBinLimits(iphi,binLim2);
containerPos -> SetBinLimits(ieta,binLim3);
- containerPos -> SetBinLimits(idcar,binLim4);
- containerPos -> SetBinLimits(ichi2TPC,binLim5);
AliCFContainer* containerNeg = new AliCFContainer("containerNeg","container for negative tracks",nstep,nvar,iBin);
//setting the bin limits
containerNeg -> SetBinLimits(ipt,binLim1);
containerNeg -> SetBinLimits(iphi,binLim2);
containerNeg -> SetBinLimits(ieta,binLim3);
- containerNeg -> SetBinLimits(idcar,binLim4);
- containerNeg -> SetBinLimits(ichi2TPC,binLim5);
//CREATE THE CUTS -----------------------------------------------
//Use AliESDtrackCuts
AliESDtrackCuts *trackCuts = new AliESDtrackCuts("AliESDtrackCuts","Standard Cuts");
+ if(trackType==0) {
+ trackCuts = trackCuts->GetStandardITSTPCTrackCuts2010(kTRUE);//Primary Track Quality Selection for Global tracks
+ }
+ else if(trackType==1) {
+ //Cuts global tracks with ITSrefit requirement
+ // TPC
+ trackCuts->SetMinNClustersTPC(70);
+ trackCuts->SetMaxChi2PerClusterTPC(4);
+ trackCuts->SetAcceptKinkDaughters(kFALSE);
+ trackCuts->SetRequireTPCRefit(kTRUE);
+ // ITS
+ trackCuts->SetRequireITSRefit(kTRUE);
+
+ trackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01");
+ trackCuts->SetMaxDCAToVertexZ(2);
+ trackCuts->SetDCAToVertex2D(kFALSE);
+ trackCuts->SetRequireSigmaToVertex(kFALSE);
+ }
trackCuts->SetEtaRange(-0.9,0.9);
trackCuts->SetPtRange(0.15, 1e10);
- trackCuts = trackCuts->GetStandardITSTPCTrackCuts2009(kTRUE);//Primary Track Quality Selection for Global tracks
AliESDtrackCuts *trackCutsTPConly = new AliESDtrackCuts("AliESDtrackCuts","Standard Cuts TPC only tracks");
+ trackCutsTPConly = trackCutsTPConly->GetStandardTPCOnlyTrackCuts();//TPC only Track Quality Cuts
+ trackCutsTPConly->SetMinNClustersTPC(70);
trackCutsTPConly->SetEtaRange(-0.9,0.9);
trackCutsTPConly->SetPtRange(0.15, 1e10);
- trackCutsTPConly = trackCutsTPConly->GetStandardTPCOnlyTrackCuts();//TPC only Track Quality Cuts
-
- // Gen-Level kinematic cuts
+ // Gen-Level kinematic cuts
AliCFTrackKineCuts *mcKineCuts = new AliCFTrackKineCuts("mcKineCuts","MC-level kinematic cuts");
mcKineCuts->SetPtRange(0.15,1e10);
mcKineCuts->SetRapidityRange(-0.9,0.9);//-0.5,0.5);
TObjArray* recTPConlyList = new TObjArray(0);
TObjArray* secList = new TObjArray(0) ;
TObjArray* recMCList = new TObjArray(0);
+ TObjArray* recTPConlyMCList = new TObjArray(0);
printf("CREATE MC KINE CUTS\n");
TObjArray* mcList = new TObjArray(0) ;
manPos->SetParticleCutsList(kStepSecondaries,secList);
manPos->SetParticleCutsList(kStepReconstructedMC,recMCList);
manPos->SetParticleCutsList(kStepMCAcceptance,mcList);
-
+ manPos->SetParticleCutsList(kStepReconstructedTPCOnlyMC,recTPConlyMCList);
AliCFManager* manNeg = new AliCFManager("manNeg","Manager for Negative tracks") ;
manNeg->SetParticleContainer(containerNeg);
manNeg->SetParticleCutsList(kStepSecondaries,secList);
manNeg->SetParticleCutsList(kStepReconstructedMC,recMCList);
manNeg->SetParticleCutsList(kStepMCAcceptance,mcList);
-
+ manNeg->SetParticleCutsList(kStepReconstructedTPCOnlyMC,recTPConlyMCList);
printf("Create task AliPWG4HighPtSpectra\n");
- AliPWG4HighPtSpectra *taskPWG4HighPtSpectra = new AliPWG4HighPtSpectra("taskPWG4HighPtSpectra");
+ AliPWG4HighPtSpectra *taskPWG4HighPtSpectra = new AliPWG4HighPtSpectra(Form("AliPWG4HighPtSpectra%d",trackType));
+ taskPWG4HighPtSpectra->SetTrackType(trackType);
taskPWG4HighPtSpectra->SetCuts(trackCuts);
taskPWG4HighPtSpectra->SetCutsTPConly(trackCutsTPConly);
taskPWG4HighPtSpectra->SetCFManagerPos(manPos); //here is set the CF manager +
//------ output containers ------
TString outputfile = AliAnalysisManager::GetCommonFileName();
- outputfile += ":PWG4_HighPtSpectra";
+ outputfile += Form(":PWG4_HighPtSpectra%d",trackType);
- AliAnalysisDataContainer *coutput0 = mgr->CreateContainer("chist0HighPtSpectra", TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
- AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("ccontainer0HighPtSpectra", AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
- AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("ccontainer1HighPtSpectra", AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
- AliAnalysisDataContainer *cout_cuts0 = mgr->CreateContainer("qa_trackCuts", AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
- AliAnalysisDataContainer *cout_cuts1 = mgr->CreateContainer("qa_trackCutsTPConly", AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
+ AliAnalysisDataContainer *coutput0 = mgr->CreateContainer(Form("chist0HighPtSpectra%d",trackType), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
+ AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(Form("ccontainer0HighPtSpectra%d",trackType), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
+ AliAnalysisDataContainer *coutput2 = mgr->CreateContainer(Form("ccontainer1HighPtSpectra%d",trackType), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
+ AliAnalysisDataContainer *cout_cuts0 = mgr->CreateContainer(Form("qa_trackCuts%d",trackType), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
+ AliAnalysisDataContainer *cout_cuts1 = mgr->CreateContainer(Form("qa_trackCutsTPConly%d",trackType), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
mgr->AddTask(taskPWG4HighPtSpectra);