AliAnalysisVertexingHF* ConfigVertexingHF() {
printf("Call to AliAnalysisVertexingHF parameters setting :\n");
vHF = new AliAnalysisVertexingHF();
//--- switch-off candidates finding (default: all on)
//vHF->SetD0toKpiOff();
vHF->SetJPSItoEleOff();
//vHF->Set3ProngOff();
vHF->SetLikeSignOn(); // like-sign pairs and triplets
// vHF->SetLikeSign3prongOff();
vHF->Set4ProngOff();
//vHF->SetDstarOff();
vHF->SetFindVertexForDstar(kFALSE);
//--- secondary vertex with KF?
//vHF->SetSecVtxWithKF();
vHF->SetCascadesOff();
vHF->SetFindVertexForCascades(kFALSE);
vHF->SetMassCutBeforeVertexing(kTRUE); // PbPb
vHF->SetV0TypeForCascadeVertex(AliRDHFCuts::kAllV0s);
//--- set cuts for single-track selection
// displaced tracks
AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
esdTrackCuts->SetRequireTPCRefit(kTRUE);
esdTrackCuts->SetMinNClustersTPC(70);
esdTrackCuts->SetRequireITSRefit(kTRUE);
//esdTrackCuts->SetMinNClustersITS(4);
esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
AliESDtrackCuts::kAny);
// |d0|>75 micron for pt<2GeV, no cut above 2
esdTrackCuts->SetMinDCAToVertexXYPtDep("0.0075*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))");
esdTrackCuts->SetMaxDCAToVertexXY(1.);
esdTrackCuts->SetMaxDCAToVertexZ(1.);
esdTrackCuts->SetPtRange(0.6,1.e10);
esdTrackCuts->SetEtaRange(-0.8,+0.8);
AliAnalysisFilter *trkFilter = new AliAnalysisFilter("trackFilter");
trkFilter->AddCuts(esdTrackCuts);
vHF->SetTrackFilter(trkFilter);
// D* soft pion tracks
AliESDtrackCuts *esdTrackCutsSoftPi = new AliESDtrackCuts("AliESDtrackCuts","default");
esdTrackCutsSoftPi->SetMinNClustersITS(4);
esdTrackCutsSoftPi->SetMaxDCAToVertexXY(1.);
esdTrackCutsSoftPi->SetMaxDCAToVertexZ(1.);
esdTrackCutsSoftPi->SetPtRange(0.2,1.e10);
esdTrackCutsSoftPi->SetEtaRange(-0.8,+0.8);
AliAnalysisFilter *trkFilterSoftPi = new AliAnalysisFilter("trackFilterSoftPi");
trkFilterSoftPi->AddCuts(esdTrackCutsSoftPi);
vHF->SetTrackFilterSoftPi(trkFilterSoftPi);
//--- set cuts for candidates selection
Int_t nptbins=2; Float_t ptlimits[2]={0.,1000000.};
AliRDHFCutsD0toKpi *cutsD0toKpi = new AliRDHFCutsD0toKpi("CutsD0toKpi");
cutsD0toKpi->SetStandardCutsPbPb2010();
cutsD0toKpi->SetMinCentrality(-10);
cutsD0toKpi->SetMaxCentrality(110);
cutsD0toKpi->SetUseSpecialCuts(kFALSE);
cutsD0toKpi->SetMinPtCandidate(0.);
cutsD0toKpi->SetUsePID(kFALSE);
cutsD0toKpi->SetUsePhysicsSelection(kFALSE);
cutsD0toKpi->SetMaxVtxZ(1.e6);
cutsD0toKpi->SetTriggerClass("");
Float_t cutsArrayD0toKpi[11]={0.4,999999.,1.1,0.,0.,999999.,999999.,0.,0.5,-1,0.};
cutsD0toKpi->SetPtBins(nptbins,ptlimits);
cutsD0toKpi->SetCuts(11,cutsArrayD0toKpi);
cutsD0toKpi->AddTrackCuts(esdTrackCuts);
vHF->SetCutsD0toKpi(cutsD0toKpi);
AliRDHFCutsJpsitoee *cutsJpsitoee = new AliRDHFCutsJpsitoee("CutsJpsitoee");
Float_t cutsArrayJpsitoee[9]={0.350,100000.,1.1,0.,0.,100000.,100000.,100000000.,-1.1};
cutsJpsitoee->SetCuts(9,cutsArrayJpsitoee);
cutsJpsitoee->AddTrackCuts(esdTrackCuts);
vHF->SetCutsJpsitoee(cutsJpsitoee);
AliRDHFCutsDplustoKpipi *cutsDplustoKpipi = new AliRDHFCutsDplustoKpipi("CutsDplustoKpipi");
cutsDplustoKpipi->SetStandardCutsPbPb2010();
cutsDplustoKpipi->SetUsePID(kFALSE);
Float_t cutsArrayDplustoKpipi[14]={0.25,0.3,0.3,0.,0.,0.01,0.05,0.05,0.,0.88,0.,10000000000.,0.,-1.};
cutsDplustoKpipi->SetPtBins(nptbins,ptlimits);
cutsDplustoKpipi->SetCuts(14,cutsArrayDplustoKpipi);
cutsDplustoKpipi->AddTrackCuts(esdTrackCuts);
cutsDplustoKpipi->SetMinPtCandidate(2.);
vHF->SetCutsDplustoKpipi(cutsDplustoKpipi);
AliRDHFCutsDstoKKpi *cutsDstoKKpi = new AliRDHFCutsDstoKKpi("CutsDstoKKpi");
cutsDstoKKpi->SetStandardCutsPbPb2010();
cutsDstoKKpi->SetUsePID(kFALSE);
Float_t cutsArrayDstoKKpi[20]={0.35,0.3,0.3,0.,0.,0.005,0.06,0.,0.,0.9,0.,100000.,0.035,0.0001,-1.,1.,0.,0.,0.,-1.};
cutsDstoKKpi->SetPtBins(nptbins,ptlimits);
cutsDstoKKpi->SetCuts(20,cutsArrayDstoKKpi);
cutsDstoKKpi->AddTrackCuts(esdTrackCuts);
cutsDstoKKpi->SetMinPtCandidate(2.);
vHF->SetCutsDstoKKpi(cutsDstoKKpi);
AliRDHFCutsLctopKpi *cutsLctopKpi = new AliRDHFCutsLctopKpi("CutsLctopKpi");
cutsLctopKpi->SetStandardCutsPbPb2010();
cutsLctopKpi->SetUsePID(kFALSE);
Float_t cutsArrayLctopKpi[13]={0.13,0.9,1.,0.,0.,0.01,0.04,0.006,1.,0.5,0.,0.05,0.4};
cutsLctopKpi->SetPtBins(nptbins,ptlimits);
cutsLctopKpi->SetCuts(13,cutsArrayLctopKpi);
cutsLctopKpi->AddTrackCuts(esdTrackCuts);
cutsLctopKpi->SetMinPtCandidate(2.);
vHF->SetCutsLctopKpi(cutsLctopKpi);
AliRDHFCutsD0toKpipipi *cutsD0toKpipipi = new AliRDHFCutsD0toKpipipi("CutsD0toKpipipi");
Float_t cutsArrayD0toKpipipi[9]={0.2,0.04,0.00,0.01,0.02,0.8,0.,0.1,0.};
cutsD0toKpipipi->SetCuts(9,cutsArrayD0toKpipipi);
cutsD0toKpipipi->AddTrackCuts(esdTrackCuts);
vHF->SetCutsD0toKpipipi(cutsD0toKpipipi);
// D* pt dependent cuts ------------------------------------------
AliRDHFCutsDStartoKpipi *cutsDStartoKpipi = new AliRDHFCutsDStartoKpipi("CutsDStartoKpipi");
cutsDStartoKpipi->SetUsePID(kFALSE);
const Int_t nvars=16;
const Int_t nptbins=2;
Float_t* ptbins;
ptbins=new Float_t[nptbins+1];
ptbins[0]=0.;
ptbins[1]=5.;
ptbins[2]=999.;
cutsDStartoKpipi->SetPtBins(nptbins+1,ptbins);
Float_t** rdcutsvalmine;
rdcutsvalmine=new Float_t*[nvars];
for(Int_t iv=0;iv<nvars;iv++){
rdcutsvalmine[iv]=new Float_t[nptbins];
}
//0-5
rdcutsvalmine[0][0]=0.10; //D0 inv mass window
rdcutsvalmine[1][0]=0.06; // dca
rdcutsvalmine[2][0]=0.9; // thetastar
rdcutsvalmine[3][0]=0.5; // pt Pion
rdcutsvalmine[4][0]=0.5; // Pt Kaon
rdcutsvalmine[5][0]=0.1; // d0K
rdcutsvalmine[6][0]=0.1; // d0Pi
rdcutsvalmine[7][0]=0.0001; // d0xd0
rdcutsvalmine[8][0]=0.8; // costhetapoint
rdcutsvalmine[9][0]=0.15; // Dstar inv mass window
rdcutsvalmine[10][0]=0.03; // half width of (M_Kpipi-M_D0)
rdcutsvalmine[11][0]=0.1; // Pt min of Pi soft
rdcutsvalmine[12][0]=100.; // Pt max of pi soft
rdcutsvalmine[13][0]=9999.; // theta
rdcutsvalmine[14][0]=0.9; // |cosThetaPointXY|
rdcutsvalmine[15][0]=1.; // NormDecayLenghtXY
//5-999
rdcutsvalmine[0][1]=0.10; //D0 inv mass window
rdcutsvalmine[1][1]=0.06; // dca
rdcutsvalmine[2][1]=0.9; // thetastar
rdcutsvalmine[3][1]=0.5; // pt Pion
rdcutsvalmine[4][1]=0.5; // Pt Kaon
rdcutsvalmine[5][1]=0.1; // d0K
rdcutsvalmine[6][1]=0.1; // d0Pi
rdcutsvalmine[7][1]=0.0001; // d0xd0
rdcutsvalmine[8][1]=0.7; // costhetapoint
rdcutsvalmine[9][1]=0.15; // Dstar inv mass window
rdcutsvalmine[10][1]=0.03; // half width of (M_Kpipi-M_D0)
rdcutsvalmine[11][1]=0.1; // Pt min of Pi soft
rdcutsvalmine[12][1]=100.; // Pt max of pi soft
rdcutsvalmine[13][1]=9999.; // theta
rdcutsvalmine[14][1]=0.8; // |cosThetaPointXY|
rdcutsvalmine[15][1]=0.; // NormDecayLenghtXY
cutsDStartoKpipi->SetCuts(nvars,nptbins,rdcutsvalmine);
cutsDStartoKpipi->AddTrackCuts(esdTrackCuts);
cutsDStartoKpipi->AddTrackCutsSoftPi(esdTrackCutsSoftPi);
cutsDStartoKpipi->SetMinPtCandidate(2.);
vHF->SetCutsDStartoKpipi(cutsDStartoKpipi);
//--------------------------------------------------------
AliRDHFCutsLctoV0 *cutsLctoV0 = new AliRDHFCutsLctoV0("CutsLctoV0");
Float_t cutsArrayLctoV0[17]={1.0,1.0,0.05,0.05,0.0,0.0,0.0,1000.,1000.,0.99,3.,1000.,0.,0.,0.,0.,0.0};
cutsLctoV0->SetCuts(17,cutsArrayLctoV0);
cutsLctoV0->AddTrackCuts(esdTrackCuts);
vHF->SetCutsLctoV0(cutsLctoV0);
//
//--- set this if you want to reconstruct primary vertex candidate by
// candidate using other tracks in the event (for pp, broad
// interaction region)
//vHF->SetRecoPrimVtxSkippingTrks();
//--- OR set this if you want to remove the candidate daughters from
// the primary vertex, without recostructing it from scratch
//vHF->SetRmTrksFromPrimVtx();
//--- check the settings
vHF->PrintStatus();
//--- verbose
// AliLog::SetClassDebugLevel("AliAnalysisVertexingHF",2);
return vHF;
}
AliCFSingleTrackEfficiencyTask *AddSingleTrackEfficiencyTaskAnalysisQA(TString collSyst="pp",
const Bool_t readAOD = 0, // Flag to read AOD:1 or ESD:0
TString suffix="default", // suffix for the output directory
AliPID::EParticleType specie=AliPID::kPion,
Int_t pdgcode=0, //particle specie
Bool_t useMCtruthForKine=kFALSE,
Double_t maxRadiusForPrimaries=999.,
ULong64_t triggerMask=AliVEvent::kAnyINT,
TString centralityEstimator = "V0M",
Int_t fBit=0,
Bool_t TPCRefit = kTRUE,
Int_t minclustersTPC = 0,
Bool_t ITSRefit = kTRUE,
Int_t spdHits=AliESDtrackCuts::kAny,
Int_t minclustersITS = 0,
Bool_t useRecoEvSelCutsForKine=kFALSE,
Int_t configuration=AliCFSingleTrackEfficiencyTask::kFast,
Int_t usageOfBayesianPID=AliSingleTrackEffCuts::kNoBayesianPID)
{
Info("AliCFSingleTrackEfficiencyTask","SETUP CONTAINER");
//
// Setting up the container
//
// Variables
const Int_t nvar = 7; // number of variables on the grid: pt, y, phi, theta, zvtx, multiplicity, centrality
UInt_t nstep = 8; // number of container steps
const UInt_t ipt = 0;
const UInt_t ieta = 1;
const UInt_t iphi = 2;
const UInt_t itheta = 3;
const UInt_t izvtx = 4;
const UInt_t imult = 5;
const UInt_t icent = 6;
//
// Containter bining
const Int_t nbinpt=30;
// A2. Bins variation by hand for other variables
const Int_t nbin2 = configuration==AliCFSingleTrackEfficiencyTask::kSlow ? 8 : 8; //bins in eta
const Int_t nbin3 = configuration==AliCFSingleTrackEfficiencyTask::kSlow ? 9 : 18; //bins in phi
const Int_t nbin4 = configuration==AliCFSingleTrackEfficiencyTask::kSlow ? 9 : 1; //bins in theta
const Int_t nbin5 = configuration==AliCFSingleTrackEfficiencyTask::kSlow ? 10 : 4; //bins in zvtx
const Int_t nbinmult=8;
const Int_t nbincent=12;
//arrays for the number of bins in each dimension
Int_t iBin[nvar];
iBin[0]=nbinpt;
iBin[1]=nbin2;
iBin[2]=nbin3;
iBin[3]=nbin4;
iBin[4]=nbin5;
iBin[5]=nbinmult;
iBin[6]=nbincent;
//arrays for lower bounds :
Double_t binLimpT[nbinpt+1] = {0.00,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.50,0.60,
0.70,0.80,0.90,1.00,1.25,1.50,1.75,2.00,2.50,3.00,
3.50,4.00,5.00,6.00,7.00,8.00,10.0,12.0,16.0,20.0,30.0};
Double_t *binLim2 = new Double_t[iBin[1]+1];
Double_t *binLim3 = new Double_t[iBin[2]+1];
Double_t *binLim4 = new Double_t[iBin[3]+1];
Double_t *binLim5 = new Double_t[iBin[4]+1];
Double_t binLimmult[nbinmult+1] = {0.,5.,10.,20.,30.,40.,50.,80.,200.};
Double_t binLimcent[nbincent+1] = {0.,2.5,5.0,7.5,10.,15.,20.,30.,40.,50.,60.,80.,100.};
if(collSyst=="Pb-Pb" || collSyst=="PbPb" || collSyst=="Xe-Xe" || collSyst=="XeXe"){
binLimmult[1]=100.;
binLimmult[2]=500.;
binLimmult[3]=1000.;
binLimmult[4]=2000.;
binLimmult[5]=3000.;
binLimmult[6]=4000.;
binLimmult[7]=5000.;
binLimmult[8]=10000.;
}
// Other Variables
for(Int_t i=0; i<=nbin2; i++) binLim2[i]=(Double_t)etamin + (etamax-etamin)/nbin2*(Double_t)i ;
for(Int_t i=0; i<=nbin3; i++) binLim3[i]=(Double_t)phimin + (phimax-phimin)/nbin3*(Double_t)i ;
for(Int_t i=0; i<=nbin4; i++) binLim4[i]=(Double_t)thetamin + (thetamax-thetamin)/nbin4*(Double_t)i ;
for(Int_t i=0; i<=nbin5; i++) binLim5[i]=(Double_t)zvtxmin + (zvtxmax-zvtxmin)/nbin5*(Double_t)i ;
// Container
AliCFContainer* container = new AliCFContainer(Form("container%s",suffix.Data()),"container for tracks",nstep,nvar,iBin);
container -> SetBinLimits(ipt,binLimpT); // pt
container -> SetBinLimits(ieta,binLim2); // eta
container -> SetBinLimits(iphi,binLim3); // phi
container -> SetBinLimits(itheta,binLim4); // theta
container -> SetBinLimits(izvtx,binLim5); // Zvtx
container -> SetBinLimits(imult,binLimmult);// multiplicity
container -> SetBinLimits(icent,binLimcent);// centrality
// Variable Titles
container -> SetVarTitle(ipt,"pt");
container -> SetVarTitle(ieta, "eta");
container -> SetVarTitle(iphi,"phi");
container -> SetVarTitle(itheta, "theta");
container -> SetVarTitle(izvtx, "Zvtx");
container -> SetVarTitle(imult, "Multiplicity");
container -> SetVarTitle(icent, "Centrality");
// Step Titles
container -> SetStepTitle(0, " MC Particle with Generated Cuts");
container -> SetStepTitle(1, " MC Particle with Kine Acceptance Cuts");
container -> SetStepTitle(2, " MC Particle with Track Ref Acceptance Cuts");
container -> SetStepTitle(3, " Total Reconstructed Particle ");
container -> SetStepTitle(4, " Reco Particle With Kine Acceptance Cuts");
container -> SetStepTitle(5, " Reco Particle to MC True pt particles ");
container -> SetStepTitle(6, " Reco Particle With Quality Cuts");
container -> SetStepTitle(7, " Reco PID With Quality Cuts");
// SET TLIST FOR QA HISTOS
TList* qaList = new TList();
TObjArray* emptyList = new TObjArray(0);
//CREATE THE INTERFACE TO CORRECTION FRAMEWORK USED IN THE TASK
printf("CREATE INTERFACE AND CUTS\n");
AliCFManager* man = new AliCFManager();
man->SetNStepEvent(2);
man->SetEventContainer(container);
man->SetEventCutsList(0,emptyList);//evtmcList);
man->SetEventCutsList(1,emptyList);//evtrecoList);
man->SetParticleContainer(container);
man->SetParticleCutsList(0,emptyList);//mcGenList);
man->SetParticleCutsList(1,emptyList);//mcKineList);
man->SetParticleCutsList(2,emptyList);//mcaccList);
man->SetParticleCutsList(3,emptyList);//evtrecoPureList);
man->SetParticleCutsList(4,emptyList);//recKineList);
man->SetParticleCutsList(5,emptyList);
man->SetParticleCutsList(6,emptyList);
man->SetParticleCutsList(7,emptyList);
// Simulated particle & event cuts
AliSingleTrackEffCuts* cuts = new AliSingleTrackEffCuts();
cuts->SetPtRange(ptmin,ptmax);
cuts->SetEtaRange(etamin,etamax);
cuts->SetIsCharged(ischarged);
cuts->SetMinVtxContr(1);
cuts->SetMaxVtxZ(zvtxmax);
cuts->SetNumberOfClusters(mintrackrefsITS,mintrackrefsTPC,mintrackrefsTOF,mintrackrefsMUON);
cuts->SetTriggerMask(triggerMask);
cuts->SetIsAOD(readAOD);
cuts->SetMaxRadiusOfParticleOrigin(maxRadiusForPrimaries);
//
// Pid selection here
//
if(pdgcode>0){
cuts->SetUsePid(true);
cuts->SetParticleSpecie(specie);
cuts->SetPdgCode(pdgcode);
//
const Int_t nlims=1;
Float_t plims[nlims+1]={0.,999.}; //TPC limits in momentum [GeV/c]
Float_t sigmas[nlims]={3.};
cuts->SetUseTPCPid();
cuts->SetTPCSigmaPtBins(nlims,plims,sigmas);
cuts->SetMaximumPTPC(4.);
//
const Int_t nlims2=1;
Float_t plims2[nlims2+1]={0.,999.}; //TPC limits in momentum [GeV/c]
Float_t sigmas2[nlims2]={3.};
cuts->SetUseTOFPid();
cuts->SetTOFSigmaPtBins(nlims2,plims2,sigmas2);
cuts->SetMaximumPTOF(4.);
if(usageOfBayesianPID>0) {
cuts->SetUseCombinPID(usageOfBayesianPID);
if(usageOfBayesianPID==AliSingleTrackEffCuts::kThresholdBayesianProb)
cuts->SetPIDThreshold(thresholdPID);
}
}
//
// Track Quality cuts via ESD track cuts
//
AliESDtrackCuts* QualityCuts = new AliESDtrackCuts();
QualityCuts->SetRequireSigmaToVertex(kFALSE);
QualityCuts->SetMinNClustersTPC(minclustersTPC);
QualityCuts->SetMinNClustersITS(minclustersITS);
QualityCuts->SetRequireTPCRefit(TPCRefit);
QualityCuts->SetRequireITSRefit(ITSRefit);
QualityCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,spdHits);
QualityCuts->SetMinDCAToVertexXY(0.);
QualityCuts->SetEtaRange(etamin,etamax);
QualityCuts->SetPtRange(ptmin,ptmax);
//CREATE THE TASK
printf("CREATE CF Single track task\n");
AliCFSingleTrackEfficiencyTask *task = new AliCFSingleTrackEfficiencyTask("AliCFSingleTrackEfficiencyTask",QualityCuts,cuts);
if(readAOD && fBit>=0){
task->SetFilterBit(kTRUE);
task->SetFilterType(fBit);
}else{
task->SetFilterBit(kFALSE);
}
// task->SelectCollisionCandidates(triggerMask);//AliVEvent::kMB);
if(centralityEstimator != "") task->SetUseCentrality(kTRUE,centralityEstimator);
task->SetConfiguration(configuration);
task->SetUseGeneratedKine(useMCtruthForKine);
if(useRecoEvSelCutsForKine) task->ApplyRecoEventSelectionsToFillMCKine();
task->SetCFManager(man); //here is set the CF manager
//
// Get the pointer to the existing analysis manager via the static access method.
//==============================================================================
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
if (!mgr) {
::Error("AddTask", "No analysis manager to connect to.");
return NULL;
}
// This task requires an ESD or AOD input handler and an AOD output handler.
// Check this using the analysis manager.
//===============================================================================
TString type = mgr->GetInputEventHandler()->GetDataType();
if (!type.Contains("ESD") && !type.Contains("AOD")) {
::Error("AddSingleTrackEfficiencyTaskPbPb", "AliCFSingleTrackEfficiency task needs the manager to have an ESD or AOD input handler.");
return NULL;
}
mgr->AddTask(task);
printf(" Create the output container\n");
//
// Create and connect containers for input/output
//
// ----- output data -----
TString outputfile = AliAnalysisManager::GetCommonFileName();
TString input1name="cchain0";
TString output2name="HistEventsProcessed", output3name="container",output4name="list",output5name="ESDtrackCuts",output6name="MCtrackCuts";
outputfile += ":PWGPP_CFSingleTrack";
// outputfile += suffix;
output2name += suffix;
output3name += suffix;
output4name += suffix;
output5name += suffix;
output6name += suffix;
// ------ input data ------
AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
// ----- output data -----
// output TH1I for event counting
AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(output2name, TH1I::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
// output Correction Framework Container (for acceptance & efficiency calculations)
AliAnalysisDataContainer *coutput2 = mgr->CreateContainer(output3name, AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
// output QA histograms
AliAnalysisDataContainer *coutput3 = mgr->CreateContainer(output4name, TList::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
// output ESD track cuts for book keeping
AliAnalysisDataContainer *coutput4 = mgr->CreateContainer(output5name, AliESDtrackCuts::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
// output event and particle selection cuts for book keeping
AliAnalysisDataContainer *coutput5 = mgr->CreateContainer(output6name, AliSingleTrackEffCuts::Class(),AliAnalysisManager::kOutputContainer,outputfile.Data());
mgr->ConnectInput(task,0,mgr->GetCommonInputContainer());
mgr->ConnectOutput(task,1,coutput1);
mgr->ConnectOutput(task,2,coutput2);
mgr->ConnectOutput(task,3,coutput3);
mgr->ConnectOutput(task,4,coutput4);
mgr->ConnectOutput(task,5,coutput5);
return task;
}
AliAnalysisTaskEmcalHFCJQA* AddTaskEmcalHFCJQA(
const char *ntracks = "Tracks",
const char *nclusters = "CaloClusters",
const char *njets = "Jets",
const char *nrho = "Rho",
Int_t nCentBins = 1,
Double_t jetradius = 0.4,
Double_t jetptcut = 0,
Double_t jetareacut = 0.6,
const char *type = "TPC",
Int_t leadhadtype = 0,
const char *taskname = "AliAnalysisTaskEmcalHFCJQA",
TString cutfile ="HFCJCuts.root",
UInt_t triggerMask =-1,//AliVEvent::kEMC1 | AliVEvent::kEMC7 | AliVEvent::kEMC8,/*kMB kEMC7 (kEMC8) kEMCEJE kEMCEGA*/
Bool_t isMC = kFALSE,
Float_t minC = 0.,
Float_t maxC = 7.5
)
{
// Get the pointer to the existing analysis manager via the static access method.
//==============================================================================
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
if (!mgr)
{
::Error("AddTaskEmcalHFCJQA", "No analysis manager to connect to.");
return NULL;
}
// Check the analysis type using the event handlers connected to the analysis manager.
//==============================================================================
if (!mgr->GetInputEventHandler())
{
::Error("AddTaskEmcalHFCJQA", "This task requires an input event handler");
return NULL;
}
//-------------------------------------------------------
// Init the task and do settings
//-------------------------------------------------------
TString name(taskname);
if (strcmp(njets,"")) {
name += "_";
name += njets;
}
if (strcmp(nrho,"")) {
name += "_";
name += nrho;
}
if (strcmp(type,"")) {
name += "_";
name += type;
}
Printf("name: %s",name.Data());
AliAnalysisTaskEmcalHFCJQA* jetTask = new AliAnalysisTaskEmcalHFCJQA(name);
jetTask->SetReadMC(isMC);
jetTask->SetDebug(-1);
jetTask->SetFilterBit(AliAODTrack::kTrkGlobalNoDCA);
//Defaut parameters
AliParticleContainer *trackCont = jetTask->AddParticleContainer(ntracks);
trackCont->SetClassName("AliVTrack");
AliClusterContainer *clusterCont = jetTask->AddClusterContainer(nclusters);
TString strType(type);
AliJetContainer *jetCont = jetTask->AddJetContainer(njets,strType,jetradius);
if(jetCont) {
jetCont->SetRhoName(nrho);
jetCont->ConnectParticleContainer(trackCont);
jetCont->ConnectClusterContainer(clusterCont);
jetCont->SetZLeadingCut(0.98,0.98);
jetCont->SetPercAreaCut(0.6);
jetCont->SetJetPtCut(jetptcut);
jetCont->SetLeadingHadronType(leadhadtype);
}
//=========================CUTS=========================
AliRDHFJetsCuts *cuts;
bool kFileExists=kFALSE;
//if(!gSystem->AccessPathName(cutfile.Data(),kFileExists)){
if(gSystem->AccessPathName(cutfile.Data(),kFileExists))
{
Printf("\n==CutObject not Defined in .root File. Using Standard Cuts==\n");
// possible (!not standard!) selection for pp2012 data triggered with EMCAL
cuts=new AliRDHFJetsCuts();
hfTask = new AliAnalysisTaskEmcalHFCJQA("AliAnalysisTaskEmcalHFCJQA");
// AliRDHFJetsCutsVertex *cuts2=new AliRDHFJetsCutsVertex("jetCuts");
//cuts for jets
// cuts->SetJetRadius(0.4);
// cuts->SetMaxEtaJet(0.5);//0.9-R
// cuts->SetMinPtJet(0);
// cuts->SetMaxPtJet(200);
// cuts->ResetMaskAndEnableMBTrigger();
//cuts->SetUseAnyTrigger();
//cuts->SetTriggerMask(0);
cuts->SetTriggerMask(AliVEvent::kEMC1 | AliVEvent::kEMC7 | AliVEvent::kEMC8);
cuts->SetTriggerClass("");
AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
esdTrackCuts->SetMinNClustersTPC(70);
esdTrackCuts->SetMaxChi2PerClusterTPC(4);
esdTrackCuts->SetRequireTPCRefit(kTRUE);
esdTrackCuts->SetRequireITSRefit(kTRUE);
esdTrackCuts->SetMinNClustersITS(2);
esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
esdTrackCuts->SetMinDCAToVertexXY(0.);
esdTrackCuts->SetEtaRange(-0.8,0.8);
esdTrackCuts->SetPtRange(1,1.e10);
if(minC>0&&minC<maxC)
{
// Pb-Pb
cuts->SetTriggerClass("");
cuts->ResetMaskAndEnableMBTrigger();
cuts->EnableCentralTrigger();
cuts->EnableSemiCentralTrigger();
cuts->SetUseCentrality(AliRDHFCuts::kCentV0M);
cuts->SetMinCentrality(minC);
cuts->SetMaxCentrality(maxC);
}
cuts->AddTrackCuts(esdTrackCuts);
//cuts for vertexing
// ::Error("AddTaskSEHFJets","No Cut Object");
}
else
{
TFile *f=TFile::Open(cutfile.Data());
//cuts= (AliRDHFCutsD0toKpi*)f->Get("EventTrackCuts");
cuts= (AliRDHFJetsCuts*)f->Get("HFCJCuts");
cout<<"\n==========================================\n Cutfile used:\n"<<cutfile.Data()<<endl;
//cuts->PrintAll();
//jetTask->SetJetCuts(cuts);
}
if(triggerMask>0) cuts->SetTriggerMask(triggerMask);
jetTask->SetJetCuts(cuts);
delete cuts;
//==========================================================
//-------------------------------------------------------
// Final settings, pass to manager and set the containers
//-------------------------------------------------------
mgr->AddTask(jetTask);
// Create containers for input/output
AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer() ;
TString contname(name);
contname += "_histos";
AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(contname.Data(),
TList::Class(),AliAnalysisManager::kOutputContainer,
Form("%s", AliAnalysisManager::GetCommonFileName()));
mgr->ConnectInput (jetTask, 0, cinput1 );
mgr->ConnectOutput (jetTask, 1, coutput1 );
return jetTask;
}
//________________________________________________________________
AliAnalysisCuts* SetupTrackCuts(Int_t cutInstance)
{
std::cout << "SetupTrackCuts()" <<std::endl;
//AliAnalysisCuts* trackCuts=0x0;
AliESDtrackCuts *fesdTrackCuts = new AliESDtrackCuts();
//global
fesdTrackCuts->SetPtRange( 0.2 , 100. );
fesdTrackCuts->SetEtaRange( -0.8 , 0.8 );
fesdTrackCuts->SetAcceptKinkDaughters(kFALSE);
fesdTrackCuts->SetRequireSigmaToVertex(kFALSE);
fesdTrackCuts->SetDCAToVertex2D(kFALSE);
fesdTrackCuts->SetMaxDCAToVertexZ(3.);
fesdTrackCuts->SetMaxDCAToVertexXY(1.);
fesdTrackCuts->SetRequireTPCRefit(kTRUE);
fesdTrackCuts->SetRequireITSRefit(kTRUE);
if(cutInstance == 0){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 1){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(130);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 2){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 3){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(130);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 4){
fesdTrackCuts->SetMinNClustersITS(6);
fesdTrackCuts->SetMaxChi2PerClusterITS(2.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 5){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 6){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 7){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 8){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(80);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 9){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 10){
fesdTrackCuts->SetMinNClustersITS(5);
// fesdTrackCuts->SetMaxChi2PerClusterITS(6);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(120);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 11){
fesdTrackCuts->SetMinNClustersITS(5);
// fesdTrackCuts->SetMaxChi2PerClusterITS(100);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(120);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 12){
fesdTrackCuts->SetMinNClustersITS(6);
// fesdTrackCuts->SetMaxChi2PerClusterITS(100);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 13){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 14){
fesdTrackCuts->SetMinNClustersITS(6);
fesdTrackCuts->SetMaxChi2PerClusterITS(2.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.7);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 15){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(4.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 16){
fesdTrackCuts->SetMinNClustersITS(4);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
// fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 17){
fesdTrackCuts->SetMinNClustersITS(3);
// fesdTrackCuts->SetMaxChi2PerClusterITS(4);
// fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kBoth);
fesdTrackCuts->SetMinNClustersTPC(100);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(4);
}
if(cutInstance == 18){
fesdTrackCuts->SetMinNClustersITS(4);
// fesdTrackCuts->SetMaxChi2PerClusterITS(6);
// fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(80);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.5);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
if(cutInstance == 19){
fesdTrackCuts->SetMinNClustersITS(5);
fesdTrackCuts->SetMaxChi2PerClusterITS(3.5);
fesdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kFirst);
fesdTrackCuts->SetMinNClustersTPC(120);
fesdTrackCuts->SetMinNCrossedRowsTPC(100);
fesdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.9);
fesdTrackCuts->SetMaxChi2PerClusterTPC(3);
}
return fesdTrackCuts;
}