]>
Commit | Line | Data |
---|---|---|
50c7d9f7 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
139cbd96 | 16 | /* $Id$ */ |
17 | ||
18 | //--------------------------------------------------------------------- | |
19 | // Class for input particles | |
20 | // manages the search for jets | |
21 | // Authors: Elena Bruna elena.bruna@yale.edu | |
22 | // | |
23 | // ** 2011 magali.estienne@subatech.in2p3.fr & alexandre.shabetai@cern.ch | |
24 | // Modified accordingly to reader/finder splitting and new handling of neutral information | |
25 | //--------------------------------------------------------------------- | |
26 | ||
50c7d9f7 | 27 | #include <Riostream.h> |
50c7d9f7 | 28 | |
29 | #include "AliJetHeader.h" | |
50c7d9f7 | 30 | #include "AliFastJetHeaderV1.h" |
50c7d9f7 | 31 | #include "AliFastJetInput.h" |
139cbd96 | 32 | #include "AliJetCalTrk.h" |
50c7d9f7 | 33 | |
34 | #include "fastjet/PseudoJet.hh" | |
50c7d9f7 | 35 | |
50c7d9f7 | 36 | #include<vector> |
50c7d9f7 | 37 | |
38 | using namespace std; | |
39 | ||
50c7d9f7 | 40 | ClassImp(AliFastJetInput) |
41 | ||
42 | //////////////////////////////////////////////////////////////////////// | |
43 | ||
44 | AliFastJetInput::AliFastJetInput(): | |
139cbd96 | 45 | fHeader(0x0), |
46 | fCalTrkEvent(0x0), | |
47 | fInputParticles(0), | |
48 | fInputParticlesCh(0) | |
50c7d9f7 | 49 | { |
50 | // Default constructor | |
51 | } | |
139cbd96 | 52 | |
53 | //______________________________________________________________________ | |
50c7d9f7 | 54 | AliFastJetInput::AliFastJetInput(const AliFastJetInput &input): |
139cbd96 | 55 | TObject(input), |
56 | fHeader(input.fHeader), | |
57 | fCalTrkEvent(input.fCalTrkEvent), | |
58 | fInputParticles(input.fInputParticles), | |
59 | fInputParticlesCh(input.fInputParticlesCh) | |
50c7d9f7 | 60 | { |
61 | // copy constructor | |
62 | } | |
139cbd96 | 63 | |
50c7d9f7 | 64 | //______________________________________________________________________ |
139cbd96 | 65 | AliFastJetInput& AliFastJetInput::operator=(const AliFastJetInput& source) |
66 | { | |
67 | // Assignment operator. | |
68 | if(this!=&source){ | |
69 | TObject::operator=(source); | |
70 | fHeader = source.fHeader; | |
71 | fCalTrkEvent = source.fCalTrkEvent; | |
72 | fInputParticles = source.fInputParticles; | |
73 | fInputParticlesCh = source.fInputParticlesCh; | |
74 | } | |
75 | ||
76 | return *this; | |
50c7d9f7 | 77 | |
78 | } | |
50c7d9f7 | 79 | |
139cbd96 | 80 | //___________________________________________________________ |
81 | void AliFastJetInput::FillInput() | |
82 | { | |
83 | // fills input particles for FASTJET based analysis | |
84 | ||
50c7d9f7 | 85 | AliFastJetHeaderV1 *header = (AliFastJetHeaderV1*)fHeader; |
139cbd96 | 86 | Int_t debug = header->GetDebug(); // debug option |
87 | ||
88 | if(debug>0) cout<<"-------- AliFastJetInput::FillInput() ----------------"<<endl; | |
89 | ||
50c7d9f7 | 90 | fInputParticles.clear(); |
91 | fInputParticlesCh.clear(); | |
92 | ||
50c7d9f7 | 93 | // RUN ALGORITHM |
94 | // read input particles ----------------------------- | |
95 | vector<fastjet::PseudoJet> inputParticles; | |
139cbd96 | 96 | |
97 | if(fCalTrkEvent == 0) { cout << "Could not get the CalTrk Event" << endl; return; } | |
98 | Int_t nIn = fCalTrkEvent->GetNCalTrkTracks() ; | |
99 | if(nIn == 0) { if (debug>0) cout << "entries = 0 ; Event empty !!!" << endl ; return; } | |
100 | ||
101 | // Information extracted from fCalTrkEvent | |
102 | // load input vectors and calculate total energy in array | |
103 | Float_t px = -999., py = -999., pz = -999., en = -999.; | |
104 | ||
105 | // Fill charged tracks | |
106 | for(Int_t i = 0; i < fCalTrkEvent->GetNCalTrkTracks(); i++) | |
107 | { // loop for all input particles | |
108 | if (fCalTrkEvent->GetCalTrkTrack(i)->GetCutFlag() != 1) continue; | |
109 | px = fCalTrkEvent->GetCalTrkTrack(i)->GetPx(); | |
110 | py = fCalTrkEvent->GetCalTrkTrack(i)->GetPy(); | |
111 | pz = fCalTrkEvent->GetCalTrkTrack(i)->GetPz(); | |
112 | en = fCalTrkEvent->GetCalTrkTrack(i)->GetP(); | |
113 | ||
114 | fastjet::PseudoJet inputPart(px,py,pz,en); // create PseudoJet object | |
115 | inputPart.set_user_index(i); //label the particle into Fastjet algortihm | |
116 | fInputParticles.push_back(inputPart); // back of the inputParticles vector | |
117 | ||
118 | // only for charged particles (TPC+ITS) | |
119 | fastjet::PseudoJet inputPartCh(px,py,pz,en); // create PseudoJet object | |
120 | inputPartCh.set_user_index(i); //label the particle into Fastjet algortihm | |
121 | fInputParticlesCh.push_back(inputPartCh); // back of the inputParticles vector | |
122 | } // End loop on CalTrk | |
50c7d9f7 | 123 | |
124 | } | |
125 | ||
126 | //_____________________________________________________________________ | |
139cbd96 | 127 | Double_t AliFastJetInput::Thermalspectrum(const Double_t *x, const Double_t *par) |
50c7d9f7 | 128 | { |
139cbd96 | 129 | // compute an exponential function |
50c7d9f7 | 130 | return x[0]*TMath::Exp(-x[0]/par[0]); |
131 | ||
132 | } |