]>
Commit | Line | Data |
---|---|---|
c6bca3f2 | 1 | AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC,Bool_t MassConst,Bool_t MassWidthCut,Bool_t MassCal,Bool_t MassNonlinear,Double_t asspTCut,Double_t angleCut,Double_t MassCut, Double_t NsigCut, Int_t fqa){\r |
4a7f4892 | 2 | //\r |
3 | // HFE standard task configuration\r | |
4 | //\r | |
5 | \r | |
6b2cee73 | 6 | printf("Setting in this analysis\n");\r |
7 | printf("MassConst = %d/n",MassConst); \r | |
8 | printf("MassCal = %d/n", MassCal); \r | |
9 | printf("MassWidth = %d/n",MassWidthCut); \r | |
10 | printf("asspTCut = %d/n",asspTCut); \r | |
11 | printf("angleCut = %d/n",angleCut); \r | |
12 | printf("MassCut = %d/n",MassCut); \r | |
13 | printf("NsigCut = %d/n", NsigCut); \r | |
c6bca3f2 | 14 | printf("qa = %d/n", fqa); \r |
6b2cee73 | 15 | \r |
4a7f4892 | 16 | Bool_t kAnalyseTaggedTracks = kTRUE;\r |
17 | \r | |
18 | AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");\r | |
19 | hfecuts->CreateStandardCuts();\r | |
20 | hfecuts->SetMinNClustersTPC(100);\r | |
21 | hfecuts->SetMinRatioTPCclusters(0.6);\r | |
22 | hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);\r | |
23 | hfecuts->SetMinNClustersITS(3);\r | |
ca1f9ac2 | 24 | hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny);\r |
4a7f4892 | 25 | hfecuts->SetCheckITSLayerStatus(kFALSE);\r |
26 | hfecuts->SetVertexRange(10.);\r | |
27 | hfecuts->SetTOFPIDStep(kFALSE);\r | |
28 | hfecuts->SetPtRange(2, 50);\r | |
ca1f9ac2 | 29 | hfecuts->SetMaxImpactParam(3.,3.);\r |
4a7f4892 | 30 | \r |
3c6739a4 | 31 | AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFEanalysisEMCal");\r |
4a7f4892 | 32 | printf("task ------------------------ %p\n ", task);\r |
33 | task->SetHFECuts(hfecuts);\r | |
a009053a | 34 | task->SetMassConstraint(MassConst);\r |
3c475a9d | 35 | task->SetMassWidthCut(MassWidthCut);\r |
2f8ef315 | 36 | task->SetMassNonlinear(MassNonlinear);\r |
4642ba9f | 37 | //Double_t masscut = 0.05;\r |
38 | //if(!MassConst)masscut = 0.1;\r | |
39 | Double_t masscut = MassCut;\r | |
a009053a | 40 | task->SetInvariantMassCut(masscut);\r |
4b76080c | 41 | task->SetOpeningAngleCut(angleCut);\r |
6451f693 | 42 | task->SetMimpTassCut(asspTCut);\r |
6b2cee73 | 43 | task->SetMimNsigassCut(NsigCut); \r |
6451f693 | 44 | task->SetMassCalMethod(MassCal);\r |
c6bca3f2 | 45 | task->SetQAHist(fqa);\r |
4a7f4892 | 46 | \r |
47 | // Define PID\r | |
48 | AliHFEpid *pid = task->GetPID();\r | |
49 | if(useMC) pid->SetHasMCData(kTRUE);\r | |
50 | pid->AddDetector("TPC", 0);\r | |
89f41a30 | 51 | //pid->AddDetector("EMCAL", 1);\r |
4a7f4892 | 52 | \r |
53 | Double_t params[4];\r | |
54 | char *cutmodel;\r | |
55 | if(useMC){\r | |
56 | // Monte-Carlo needs modelling of the falling mean with momentum at low momentum\r | |
57 | // for high momentum it is consistent with a flat -0.94\r | |
3c6739a4 | 58 | //cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";\r |
59 | //Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};\r | |
60 | //for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];\r | |
61 | cutmodel = "pol0";\r | |
62 | params[0] = -1.0; //sigma min\r | |
4a7f4892 | 63 | } else {\r |
64 | // Data is consistent with a flat 0.12\r | |
65 | cutmodel = "pol0";\r | |
66 | //params[0] = -0.0015;\r | |
67 | //params[0] = -3.0;\r | |
68 | //params[0] = -0.05; //sigma min\r | |
69 | params[0] = -1.0; //sigma min\r | |
70 | }\r | |
2f04fc90 | 71 | //pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); \r |
72 | for(Int_t a=0;a<11;a++)pid->ConfigureTPCcentralityCut(a,cutmodel,params,3.0);\r | |
73 | \r | |
4a7f4892 | 74 | \r |
bc9ab47c | 75 | // change E/p cuts\r |
76 | AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);\r | |
77 | emcpid->SetEoPMax(1.3);\r | |
78 | emcpid->SetEoPMim(0.9);\r | |
79 | \r | |
4a7f4892 | 80 | printf("*************************************\n");\r |
81 | printf("Configuring standard Task:\n");\r | |
82 | // task->PrintStatus();\r | |
83 | pid->PrintStatus();\r | |
84 | printf("*************************************\n"); \r | |
85 | return task;\r | |
86 | }\r |