From 0beb4580e6b23467f1d621bbae452113ebd85889 Mon Sep 17 00:00:00 2001 From: dainese Date: Mon, 7 Feb 2011 11:51:07 +0000 Subject: [PATCH] Setter for centrality in PbPb (Chiara) --- .../macros/makeTFile4CutsD0toKpi.C | 297 +++++++----------- 1 file changed, 122 insertions(+), 175 deletions(-) diff --git a/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C b/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C index 2c5a7e306a4..31f3ca37e8c 100644 --- a/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C +++ b/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C @@ -35,30 +35,31 @@ void makeInputAliAnalysisTaskSED0Mass(){ AliESDtrackCuts::kAny); // default is kBoth, otherwise kAny esdTrackCuts->SetMinDCAToVertexXY(0.); - esdTrackCuts->SetPtRange(0.3,1.e10); + esdTrackCuts->SetPtRange(0.8,1.e10); RDHFD0toKpi->AddTrackCuts(esdTrackCuts); const Int_t nvars=9; - const Int_t nptbins=11; + const Int_t nptbins=13; Float_t* ptbins; ptbins=new Float_t[nptbins+1]; ptbins[0]=0.; - ptbins[1]=1.; - ptbins[2]=2.; - ptbins[3]=3.; - ptbins[4]=4.; - ptbins[5]=5.; - ptbins[6]=6.; - ptbins[7]=8.; - ptbins[8]=12.; - ptbins[9]=16.; - ptbins[10]=20.; - ptbins[11]=24.; - //ptbins[12]=99999.; - + ptbins[1]=0.5; + ptbins[2]=1.; + ptbins[3]=2.; + ptbins[4]=3.; + ptbins[5]=4.; + ptbins[6]=5.; + ptbins[7]=6.; + ptbins[8]=8.; + ptbins[9]=12.; + ptbins[10]=16.; + ptbins[11]=20.; + ptbins[12]=24.; + ptbins[13]=9999.; + RDHFD0toKpi->SetPtBins(nptbins+1,ptbins); @@ -80,125 +81,56 @@ void makeInputAliAnalysisTaskSED0Mass(){ // printf(" d0d0 [cm^2] < %f\n",fD0toKpiCuts[7]); // printf(" cosThetaPoint > %f\n",fD0toKpiCuts[8]); - /* - //setting PPR cut values - rdcutsvalPPR[0][0]=0.7; - rdcutsvalPPR[1][0]=0.04; - rdcutsvalPPR[2][0]=0.8; - rdcutsvalPPR[3][0]=0.5; - rdcutsvalPPR[4][0]=0.5; - rdcutsvalPPR[5][0]=0.05; - rdcutsvalPPR[6][0]=0.05; - rdcutsvalPPR[7][0]=-0.0002; - rdcutsvalPPR[8][0]=0.5; - - rdcutsvalPPR[0][1]=rdcutsvalPPR[0][2]=0.7; - rdcutsvalPPR[1][1]=rdcutsvalPPR[1][2]=0.02; - rdcutsvalPPR[2][1]=rdcutsvalPPR[2][2]=0.8; - rdcutsvalPPR[3][1]=rdcutsvalPPR[3][2]=0.7; - rdcutsvalPPR[4][1]=rdcutsvalPPR[4][2]=0.7; - rdcutsvalPPR[5][1]=rdcutsvalPPR[5][2]=0.05; - rdcutsvalPPR[6][1]=rdcutsvalPPR[6][2]=0.05; - rdcutsvalPPR[7][1]=rdcutsvalPPR[7][2]=-0.0002; - rdcutsvalPPR[8][1]=rdcutsvalPPR[8][2]=0.6; - - rdcutsvalPPR[0][3]=0.7; - rdcutsvalPPR[1][3]=0.02; - rdcutsvalPPR[2][3]=0.8; - rdcutsvalPPR[3][3]=0.7; - rdcutsvalPPR[4][3]=0.7; - rdcutsvalPPR[5][3]=0.05; - rdcutsvalPPR[6][3]=0.05; - rdcutsvalPPR[7][3]=-0.0001; - rdcutsvalPPR[8][3]=0.8; - - rdcutsvalPPR[0][4]=0.7; - rdcutsvalPPR[1][4]=0.02; - rdcutsvalPPR[2][4]=0.8; - rdcutsvalPPR[3][4]=0.7; - rdcutsvalPPR[4][4]=0.7; - rdcutsvalPPR[5][4]=0.05; - rdcutsvalPPR[6][4]=0.05; - rdcutsvalPPR[7][4]=-0.00005; - rdcutsvalPPR[8][4]=0.8; - */ - Double_t arrcuts[9]={0.2,0.03,0.8,0.3,0.3,0.1,0.1,-0.0004,0.7}; - //cout<<"here"<SetCuts(nvars,nptbins,rdcutsvalmine); @@ -227,12 +159,27 @@ void makeInputAliAnalysisTaskSED0Mass(){ RDHFD0toKpi->SetUseDefaultPID(kFALSE); //to use the AliAODPidHF - //activate pileup rejection - RDHFD0toKpi->SetOptPileup(AliRDHFCuts::kRejectPileupEvent); + //activate pileup rejection (for pp) + //RDHFD0toKpi->SetOptPileup(AliRDHFCuts::kRejectPileupEvent); + + //Do not recalculate the vertex + RDHFD0toKpi->SetRemoveDaughtersFromPrim(kFALSE); //activate for pp + + TString cent=""; + //centrality selection (Pb-Pb) + Float_t minc=20,maxc=80; + RDHFD0toKpi->SetMinCentrality(minc); + RDHFD0toKpi->SetMaxCentrality(maxc); + cent=Form("%.0f%.0f",minc,maxc); + RDHFD0toKpi->SetUseCentrality(AliRDHFCuts::kCentV0M); //kCentOff,kCentV0M,kCentTRK,kCentTKL,kCentCL1,kCentInvalid + + //temporary + RDHFD0toKpi->SetFixRefs(); cout<<"This is the odject I'm going to save:"<PrintAll(); - TFile* fout=new TFile("D0toKpiCutsStdpileup.root","recreate"); //set this!! + TFile* fout=new TFile(Form("D0toKpiCuts%s%s%sRecVtx%sPileupRej.root", RDHFD0toKpi->GetUseCentrality()==0 ? "pp" : "PbPb",cent.Data(),RDHFD0toKpi->GetIsPrimaryWithoutDaughters() ? "" : "No",RDHFD0toKpi->GetOptPileUp() ? "" : "No"),"recreate"); //set this!! + fout->cd(); RDHFD0toKpi->Write(); fout->Close(); @@ -258,13 +205,13 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){ esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); esdTrackCuts->SetMinDCAToVertexXY(0.); esdTrackCuts->SetEtaRange(-0.8,0.8); - esdTrackCuts->SetPtRange(0.1,1.e10); + esdTrackCuts->SetPtRange(0.8,1.e10); RDHFD0toKpi->AddTrackCuts(esdTrackCuts); const Int_t nvars=9; - const Int_t nptbins=14; //change this when adding pt bins! + const Int_t nptbins=13; //change this when adding pt bins! Float_t ptbins[nptbins+1]; ptbins[0]=0.; ptbins[1]=0.5; @@ -275,15 +222,13 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){ ptbins[6]=5.; ptbins[7]=6.; ptbins[8]=8.; - ptbins[9]=10.; - ptbins[10]=12.; - ptbins[11]=14.; - ptbins[12]=16.; - ptbins[13]=24.; - ptbins[14]=9999.; - + ptbins[9]=12.; + ptbins[10]=16.; + ptbins[11]=20.; + ptbins[12]=24.; + ptbins[13]=9999.; + RDHFD0toKpi->SetPtBins(nptbins+1,ptbins); - Float_t** rdcutsvalmine; rdcutsvalmine=new Float_t*[nvars]; @@ -304,29 +249,20 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){ // printf(" cosThetaPoint > %f\n",fD0toKpiCuts[8]); - - Float_t cutsMatrixD0toKpiStand[nptbins][nvars]={{0.2,400.*1E-4,0.7,0.3,0.3,1000.*1E-4,1000.*1E-4,-0.00005,0.8},/* pt<0.5*/ - {0.2,400.*1E-4,0.7,0.3,0.3,1000.*1E-4,1000.*1E-4,-0.00005,0.8},/* 0.524*/ - - // Float_t cutsMatrixD0toKpiStand[nptbins][nvars]={{0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.6},/* 624 */ + Float_t cutsMatrixD0toKpiStand[nptbins][nvars]={{0.3,400.*1E-4,0.8,0.8,0.8,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 024*/ - //CREATE TRANSPOSE MATRIX...REVERSE INDICES as required by AliRDHFCuts Float_t **cutsMatrixTransposeStand=new Float_t*[nvars]; for(Int_t iv=0;iv"<(name.Data(),tighterval[ival][jpt]); } } + Bool_t flagPID=kTRUE; RDHFD0toKpi->SetUsePID(flagPID); + RDHFD0toKpi->PrintAll(); printf("Use PID? %s\n",flagPID ? "yes" : "no"); - //pid settings AliAODPidHF* pidObj=new AliAODPidHF(); //pidObj->SetName("pid4D0"); @@ -490,11 +424,24 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){ RDHFD0toKpi->SetUseDefaultPID(kFALSE); //to use the AliAODPidHF - //activate pileup rejection - RDHFD0toKpi->SetOptPileup(AliRDHFCuts::kRejectPileupEvent); + //activate pileup rejection (for pp) + //RDHFD0toKpi->SetOptPileup(AliRDHFCuts::kRejectPileupEvent); + + //Do not recalculate the vertex + RDHFD0toKpi->SetRemoveDaughtersFromPrim(kFALSE); //activate for pp + + TString cent=""; + //centrality selection (Pb-Pb) + Float_t minc=20,maxc=80; + RDHFD0toKpi->SetMinCentrality(minc); + RDHFD0toKpi->SetMaxCentrality(maxc); + cent=Form("%.0f%.0f",minc,maxc); + RDHFD0toKpi->SetUseCentrality(AliRDHFCuts::kCentV0M); //kCentOff,kCentV0M,kCentTRK,kCentTKL,kCentCL1,kCentInvalid + //temporary + RDHFD0toKpi->SetFixRefs(); - TFile* fout=new TFile("cuts4SignifMaxim.root","recreate"); //set this!! + TFile* fout=new TFile(Form("cuts4SignifMaxim%s%s%sRecVtx%sPileupRej.root", RDHFD0toKpi->GetUseCentrality()==0 ? "pp" : "PbPb",cent.Data(),RDHFD0toKpi->GetIsPrimaryWithoutDaughters() ? "" : "No",RDHFD0toKpi->GetOptPileUp() ? "" : "No"),"recreate"); //set this!! fout->cd(); RDHFD0toKpi->Write(); -- 2.43.0