hfecuts->SetPtRange(2, 50);\r
hfecuts->SetMaxImpactParam(1,2);\r
\r
- AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFE v2");\r
+ AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFEanalysisEMCal");\r
printf("task ------------------------ %p\n ", task);\r
task->SetHFECuts(hfecuts);\r
task->SetInvariantMassCut(0.05);\r
pid->AddDetector("EMCAL", 1);\r
// change E/p cuts\r
AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);\r
- emcpid->SetEoPMax(1.2);\r
+ emcpid->SetEoPMax(1.3);\r
emcpid->SetEoPMim(0.9);\r
\r
Double_t params[4];\r
if(useMC){\r
// Monte-Carlo needs modelling of the falling mean with momentum at low momentum\r
// for high momentum it is consistent with a flat -0.94\r
- cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";\r
- Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};\r
- for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];\r
+ //cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";\r
+ //Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};\r
+ //for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];\r
+ cutmodel = "pol0";\r
+ params[0] = -1.0; //sigma min\r
} else {\r
// Data is consistent with a flat 0.12\r
cutmodel = "pol0";\r