Update config
authorrbailhac <rbailhac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 8 May 2012 16:50:20 +0000 (16:50 +0000)
committerrbailhac <rbailhac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 8 May 2012 16:50:20 +0000 (16:50 +0000)
PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TOFTPC.C

index 6bde87d..8f2597b 100644 (file)
@@ -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";