]>
Commit | Line | Data |
---|---|---|
1 | /////////////////////////////////////////////////////////////////////////////// | |
2 | // AddTaskJetFlowToyMC // | |
3 | // Author: Redmer A. Bertens, Utrecht University, 2013, rbertens@cern.ch // | |
4 | /////////////////////////////////////////////////////////////////////////////// | |
5 | ||
6 | /* AddTask macro for jet flow toy mc task | |
7 | * task uses an afterburner to tune vn in the pico track | |
8 | * selection which can be used by a jet finder | |
9 | * note that this task does not generate MC particles, it changes | |
10 | * the azimuthal distribution of already available tracks | |
11 | */ | |
12 | ||
13 | class AliAnalysisDataContainer; | |
14 | class AliAnalysisTaskJetFlowMC; | |
15 | ||
16 | AliAnalysisTaskJetFlowMC* AddTaskJetFlowMC( | |
17 | const char *outputTracks = "JetFlowToyMC", | |
18 | const char *inputTracks = "PicoTracks", | |
19 | const char *name ="AliAnalysisTaskJetFlowMC" | |
20 | ) | |
21 | { | |
22 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); | |
23 | if (!mgr) return 0x0; | |
24 | if (!mgr->GetInputEventHandler()) return 0x0; | |
25 | TString fileName = AliAnalysisManager::GetCommonFileName(); | |
26 | fileName += ":"; | |
27 | fileName += name; | |
28 | // create the task | |
29 | AliAnalysisTaskJetFlowMC *eTask = new AliAnalysisTaskJetFlowMC(name); | |
30 | eTask->SetTracksOutName(outputTracks); | |
31 | eTask->SetTracksInName(inputTracks); | |
32 | // connect input and output | |
33 | mgr->AddTask(eTask); | |
34 | mgr->ConnectInput (eTask, 0, mgr->GetCommonInputContainer()); | |
35 | mgr->ConnectOutput (eTask, 1, mgr->CreateContainer(Form("%s_container", fileName.Data()), TList::Class(), AliAnalysisManager::kOutputContainer, fileName.Data())); | |
36 | return eTask; | |
37 | } | |
38 | ||
39 | namespace TaskJetFlowMC{ | |
40 | TF1* GetSpectrum() { | |
41 | TF1* spectrum = new TF1("fspectrum", "[0]*(TMath::Power([1], 2)*x*TMath::Exp(-[1]*x))+(x>1)*[2]*(1.17676e-01*TMath::Sqrt(0.1396*0.1396+x*x)*TMath::Power(1.+1./[3]/8.21795e-01*TMath::Sqrt(0.1396*0.1396+x*x),-1.*[3]))*(1/(1 + TMath::Exp(-(x - [4])/[5])))", .2, 200.); | |
42 | fspectrum->SetParameters(2434401791.20528 ,2.98507 ,10069622.25117 ,5.50000 ,2.80000 ,0.20000 ); | |
43 | return fspectrum; | |
44 | } | |
45 | } |