From c216e461c0dfad88fe24d4ee71f547b3de85e8aa Mon Sep 17 00:00:00 2001 From: rbailhac Date: Tue, 8 May 2012 16:50:20 +0000 Subject: [PATCH] Update config --- .../configs/PbPb/ConfigHFE_FLOW_TOFTPC.C | 46 ++++++++++--------- 1 file changed, 25 insertions(+), 21 deletions(-) diff --git a/PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TOFTPC.C b/PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TOFTPC.C index 6bde87d89b6..8f2597baa4f 100644 --- a/PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TOFTPC.C +++ b/PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TOFTPC.C @@ -1,32 +1,41 @@ -AliAnalysisTaskHFEFlow* ConfigHFE_FLOW_TOFTPC(Bool_t useMC, Int_t tpcCls, Double_t tpcClsr,Int_t tpcClspid,Int_t itsCls, Double_t dcaxy, - Double_t dcaz, Double_t tofsig, Double_t tpcdedx0, Double_t tpcdedx1, Double_t tpcdedx2, Double_t tpcdedx3, Double_t tpcdedx4, Int_t vzero) +AliAnalysisTaskHFEFlow* ConfigHFE_FLOW_TOFTPC(Bool_t useMC, Int_t tpcCls, Double_t tpcClsr,Int_t tpcClspid, Double_t tpcsharedfraction, Int_t itsCls, Double_t chi2peritscl, Int_t pixellayer, Double_t dcaxy, Double_t dcaz, Double_t tofsig, Double_t tpcdedx0, Double_t tpcdedx1, Double_t tpcdedx2, Double_t tpcdedx3, Double_t tpcdedx4, Int_t vzero, Int_t debuglevel) { // // HFE flow task // printf("Summary settings flow task\n"); printf("TPC number of tracking clusters %d\n",tpcCls); - printf("TPC ratio clusters %f\n",tpcClsr*0.1); + printf("TPC ratio clusters %f\n",tpcClsr*0.01); printf("TPC number of pid clusters %d\n",tpcClspid); + printf("Maximal fraction of TPC shared cluster %f\n",tpcsharedfraction*0.01); printf("ITS number of clusters %d\n",itsCls); - printf("dcaxy %f\n",dcaxy); - printf("dcaz %f\n",dcaz); + printf("Maximal chi2 per ITS cluster %f\n",chi2peritscl); + printf("Requirement on the pixel layer %d\n",pixellayer); + printf("dcaxy %f\n",dcaxy*0.01); + printf("dcaz %f\n",dcaz*0.01); printf("TOF sigma %f\n",tofsig*0.1); + printf("TPC min sigma cut 0: %f\n",tpcdedx0*0.01); + printf("TPC min sigma cut 1: %f\n",tpcdedx1*0.01); + printf("TPC min sigma cut 2: %f\n",tpcdedx2*0.01); + printf("TPC min sigma cut 3: %f\n",tpcdedx3*0.01); + printf("TPC min sigma cut 4: %f\n",tpcdedx4*0.01); printf("VZERO event plane %d\n",vzero); - + printf("Debug level %d\n",debuglevel); + // Cut HFE AliHFEcuts *hfecuts = new AliHFEcuts("hfeCuts","HFE Standard Cuts"); hfecuts->CreateStandardCuts(); hfecuts->SetMinNClustersTPC(tpcCls); hfecuts->SetMinNClustersTPCPID(tpcClspid); - hfecuts->SetMinRatioTPCclusters(tpcClsr*0.1); + hfecuts->SetMinRatioTPCclusters(tpcClsr*0.01); hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable); + hfecuts->SetFractionOfSharedTPCClusters(tpcsharedfraction*0.01); hfecuts->SetMinNClustersITS(itsCls); - hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst); + hfecuts->SetCutITSpixel(pixellayer); hfecuts->SetCheckITSLayerStatus(kFALSE); - hfecuts->SetMaxImpactParam(dcaxy,dcaz); + hfecuts->SetMaxImpactParam(dcaxy*0.01,dcaz*0.01); //hfecuts->UnsetVertexRequirement(); hfecuts->SetVertexRange(10.); @@ -34,7 +43,7 @@ AliAnalysisTaskHFEFlow* ConfigHFE_FLOW_TOFTPC(Bool_t useMC, Int_t tpcCls, Double hfecuts->SetTOFPIDStep(kTRUE); // Name - TString appendix(TString::Format("TPC%dTPCr%dITS%dDCAr%dz%dVZERO%d",tpcCls,(Int_t)tpcClsr, itsCls,(Int_t)dcaxy,(Int_t)dcaz,vzero)); + TString appendix(TString::Format("TPC%dTPCr%dTPCpid%dTPCShared%dITScl%dChi2perITS%dPixelLayer%dDCAr%dz%dTOFsig%dTPCmindedx0%dTPCmindedx1%dTPCmindedx2%dTPCmindedx3%dTPCmindedx4%dVZERO%dDebugLevel%decorr%d",tpcCls,(Int_t)tpcClsr,tpcClspid,(Int_t) tpcsharedfraction,itsCls,(Int_t) chi2peritscl,(Int_t) pixellayer,(Int_t)dcaxy,(Int_t)dcaz,(Int_t)tofsig,(Int_t)tpcdedx0,(Int_t)tpcdedx1,(Int_t)tpcdedx2,(Int_t)tpcdedx3,(Int_t)tpcdedx4,vzero,debuglevel,(Int_t)withetacorrection)); printf("appendix %s\n", appendix.Data()); // The task @@ -52,19 +61,14 @@ AliAnalysisTaskHFEFlow* ConfigHFE_FLOW_TOFTPC(Bool_t useMC, Int_t tpcCls, Double if(vzero>=1) task->SetVZEROEventPlane(kTRUE); if(vzero==2) task->SetVZEROEventPlaneA(kTRUE); if(vzero==3) task->SetVZEROEventPlaneC(kTRUE); - + task->SetDebugLevel(debuglevel); // Define PID AliHFEpid *pid = task->GetPID(); if(useMC) pid->SetHasMCData(kTRUE); pid->AddDetector("TOF", 0); pid->AddDetector("TPC", 1); - - //pid->AddDetector("BAYES", 0); - //pid->ConfigureBayesDetectorMask(AliPIDResponse::kDetTPC+AliPIDResponse::kDetTOF+AliPIDResponse::kDetTRD); - //pid->ConfigureBayesPIDThreshold(0.9); - TString datatype=gSystem->Getenv("CONFIG_FILE"); if(!useMC) { @@ -74,11 +78,11 @@ AliAnalysisTaskHFEFlow* ConfigHFE_FLOW_TOFTPC(Bool_t useMC, Int_t tpcCls, Double Double_t params_centr_10_20[1]; Double_t params_centr_20_30[1]; Double_t params_centr_per[1]; - params_centr_0_5[0]=tpcdedx0; // cut tuned for 0-10% - params_centr_5_10[0]=tpcdedx1; // cut tuned for 0-10% - params_centr_10_20[0]=tpcdedx2; - params_centr_20_30[0]=tpcdedx3; - params_centr_per[0]=tpcdedx4; + params_centr_0_5[0]=tpcdedx0*0.01; // cut tuned for 0-10% + params_centr_5_10[0]=tpcdedx1*0.01; // cut tuned for 0-10% + params_centr_10_20[0]=tpcdedx2*0.01; + params_centr_20_30[0]=tpcdedx3*0.01; + params_centr_per[0]=tpcdedx4*0.01; char *cutmodel; cutmodel="pol0"; -- 2.39.3