Fixes for the TPC PID electron selection and the contamination
authorsma <sma@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 27 Sep 2010 10:32:06 +0000 (10:32 +0000)
committersma <sma@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 27 Sep 2010 10:32:06 +0000 (10:32 +0000)
correction.

PWG3/hfe/AddTaskHFE.C
PWG3/hfe/AliHFEpid.cxx

index 6ac0a8b383f7cf1adbbd6c00080aa5fe80854190..d99b932508bb97a76ad6564fb6e322044506f70c 100644 (file)
@@ -9,16 +9,16 @@ AliAnalysisTask *AddTaskHFE(){
     ::Error("AddTaskHFE", "This task requires an input event handler");
     return NULL;
   }  
-  Bool_t MCthere=kTRUE;
-  AliMCEventHandler *mcH = dynamic_cast<AliMCEventHandler*>(mgr->GetMCtruthEventHandler());
-  if(!mcH){
-    MCthere=kFALSE;
-  }
   TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
   if (type=="AOD"){
     ::Error("AddTaskHFE", "The tasks exits because AODs are in input");
     return NULL;
   }
+  Bool_t MCthere=kTRUE;
+  AliMCEventHandler *mcH = dynamic_cast<AliMCEventHandler*>(mgr->GetMCtruthE    ventHandler());
+  if(!mcH){
+    MCthere=kFALSE;
+  }
 
   //============= Set Task Name ===================
   //TString taskName=("AliAnalysisTaskHFE.cxx+");
@@ -46,10 +46,10 @@ AliAnalysisTask *AddTaskHFE(){
   task->SetHFECuts(hfecuts);
   if(!MCthere){
     TF1 *hBackground = new TF1("hadronicBackgroundFunction", "[0]+[1]*TMath::Erf([2]*x+[3])", 0, 20);
-    hBackground->SetParameter(0, 0.3336);
-    hBackground->SetParameter(1, 0.3306);
-    hBackground->SetParameter(2, 1.1255);
-    hBackground->SetParameter(3, -3.3747);
+    hBackground->SetParameter(0, 0.2194);
+    hBackground->SetParameter(1, 0.217);
+    hBackground->SetParameter(2, 0.6829);
+    hBackground->SetParameter(3, -2.697);
     task->SetBackGroundFactorsFunction(hBackground);
   }
 
index 2a918a0af33adaef16cc4462d5d1d283c52e061c..ecdcc7be7ca23b5fdb2f6039e5c586f2ce10500c 100644 (file)
@@ -532,8 +532,9 @@ void AliHFEpid::InitStrategy6(){
   //upperCut->SetParameter(0, 2.7);
   //upperCut->SetParameter(1, -0.4357);
   lowerCut->SetParameter(0, -2.65);
-  lowerCut->SetParameter(1, -0.6757);
-  if(HasMCData()) lowerCut->SetParameter(2, -1);
+  lowerCut->SetParameter(1, -0.8757);
+//  lowerCut->SetParameter(2, -1);
+  if(HasMCData()) lowerCut->SetParameter(2, -0.9);
   else lowerCut->SetParameter(2, 0);
   tpcpid->SetUpperSigmaCut(upperCut);
   tpcpid->SetLowerSigmaCut(lowerCut);