]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HLT/JET/AliHLTJETReader.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
- temporary disabled the weighting for LHC12f1a,f1b,i3
[u/mrichter/AliRoot.git] / HLT / JET / AliHLTJETReader.cxx
CommitLineData
0734d112 1//-*- Mode: C++ -*-
2// $Id: AliHLTJETReaderHeader.cxx $
3/**************************************************************************
4 * This file is property of and copyright by the ALICE HLT Project *
5 * ALICE Experiment at CERN, All rights reserved. *
6 * *
7 * Primary Authors: Jochen Thaeder <thaeder@kip.uni-heidelberg.de> *
8 * for The ALICE HLT Project. *
9 * *
10 * Permission to use, copy, modify and distribute this software and its *
11 * documentation strictly for non-commercial purposes is hereby granted *
12 * without fee, provided that the above copyright notice appears in all *
13 * copies and that both the copyright notice and this permission notice *
14 * appear in the supporting documentation. The authors make no claims *
15 * about the suitability of this software for any purpose. It is *
16 * provided "as is" without express or implied warranty. *
17 **************************************************************************/
18
19/** @file AliHLTJETReader.cxx
20 @author Jochen Thaeder
21 @date
6ce099ba 22 @brief Reader for jet finder
0734d112 23*/
24
25// see header file for class documentation
26// or
27// refer to README to build package
28// or
29// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
30
0734d112 31#include "TLorentzVector.h"
32#include "TParticle.h"
33#include "TParticlePDG.h"
1f9fec4a 34#include "TDatabasePDG.h"
0734d112 35
6ce099ba 36#include "AliHLTJETReader.h"
37
38#include "AliHLTJETConeJetCandidate.h"
39
a7f38894 40using namespace std;
41
0734d112 42/** ROOT macro for the implementation of ROOT specific class methods */
43ClassImp(AliHLTJETReader)
44
45/*
46 * ---------------------------------------------------------------------------------
47 * Constructor / Destructor
48 * ---------------------------------------------------------------------------------
49 */
50
51// #################################################################################
52AliHLTJETReader::AliHLTJETReader()
53 :
54 AliJetReader(),
55 fESD(NULL),
56 fMC(NULL),
1f9fec4a 57 fHLTMC(NULL),
6ce099ba 58 fAOD(NULL),
59#ifdef HAVE_FASTJET
1f9fec4a 60 fMomentumVector(NULL),
6ce099ba 61#endif
62 fGrid(NULL),
63 fNJetCandidates(0),
64 fJetCandidates(NULL),
65 fSeedCuts(NULL),
66 fTrackCuts(NULL) {
0734d112 67 // see header file for class documentation
68 // or
69 // refer to README to build package
70 // or
71 // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
72
73}
74
75// #################################################################################
76AliHLTJETReader::~AliHLTJETReader() {
77 // see header file for class documentation
78
6ce099ba 79#ifdef HAVE_FASTJET
80 if ( fMomentumVector )
81 delete fMomentumVector;
82 fMomentumVector = NULL;
83#endif
84
85 if ( fGrid )
86 delete fGrid;
87 fGrid = NULL;
88
89 if ( fJetCandidates ) {
90 fJetCandidates->Clear();
91 delete fJetCandidates;
92 }
93 fJetCandidates = NULL;
94
95}
96
97// #################################################################################
98Int_t AliHLTJETReader::Initialize() {
99 // see header file for class documentation
100
101 Int_t iResult = 0;
102 AliHLTJETReaderHeader* readerHeader = NULL;
103
eaca09ae 104 HLTInfo(" -= AliHLTJETReader =- ");
6ce099ba 105
106 // -- Initialize reader header
107 // -----------------------------
108 if ( fReaderHeader ) {
109 readerHeader = GetReaderHeader();
110
111 iResult = readerHeader->Initialize();
112 if ( iResult )
113 HLTError("Error initializing HLT jet reader header");
114 }
115 else {
116 HLTError("Reader Header not present");
117 iResult = -EINPROGRESS;
118 }
119
1f9fec4a 120 // -- Initialize Algorithms
121 // --------------------------
eaca09ae 122 if ( !iResult ) {
123 if ( readerHeader->GetJetAlgorithm() >= AliHLTJETBase::kFFSCSquareCell )
124 iResult = InitializeFFSC();
125 else {
1f9fec4a 126#ifdef HAVE_FASTJET
eaca09ae 127 iResult = InitializeFastjet();
1f9fec4a 128#else
eaca09ae 129 HLTError("Error FastJet not present.");
130 iResult = -EINPROGRESS;
1f9fec4a 131#endif
eaca09ae 132 }
6ce099ba 133 }
eaca09ae 134
1f9fec4a 135 // -- Get ptr to cuts from reader
136 // --------------------------------
eaca09ae 137 if ( !iResult ) {
6ce099ba 138 fTrackCuts = readerHeader->GetTrackCuts();
139 if ( ! fTrackCuts ) {
140 HLTError("Error getting ptr to track cuts.");
141 iResult = -EINPROGRESS;
142 }
143 }
144
145 // -- Final check
146 // ----------------
147 if ( iResult )
148 HLTError("Error initializing HLT jet reader");
149
150 return iResult;
151}
152
153//#################################################################################
154void AliHLTJETReader::ResetEvent() {
155 // see header file for class documentation
156
1f9fec4a 157 // -- Reset FFSC algorithms
158 // --------------------------
159 if ( GetReaderHeader()->GetJetAlgorithm() >= AliHLTJETBase::kFFSCSquareCell ) {
160 // -- clear grid
161 fGrid->Reset();
162
163 // -- clear jet candidates
164 fJetCandidates->Clear();
165
166 fNJetCandidates = 0;
167 }
6ce099ba 168
1f9fec4a 169 // -- Reset for FastJet algorithms
170 // ---------------------------------
171 else {
3c57fb71 172#ifdef HAVE_FASTJET
1f9fec4a 173 // -- Clear input vector
174 if ( fMomentumVector )
175 fMomentumVector->clear();
3c57fb71 176#endif
1f9fec4a 177 }
6ce099ba 178
179 return;
0734d112 180}
181
1f9fec4a 182/*
183 * ---------------------------------------------------------------------------------
184 * Setter
185 * ---------------------------------------------------------------------------------
186 */
187
188// #################################################################################
189void AliHLTJETReader::SetInputEvent(const TObject* esd, const TObject* aod, const TObject* mc) {
190 // see header file for class documentation
191
192 // Needs "useMC" flag for running in analysis task only
193
194 AliHLTJETReaderHeader* readerHeader = GetReaderHeader();
195
196 // -- Fill ESD
197 if ( esd && !readerHeader->GetUseMC() )
198 fESD = dynamic_cast<AliESDEvent*> (const_cast<TObject*>(esd));
199
200 // -- Fill AOD
201 else if ( aod && !readerHeader->GetUseMC() )
202 fAOD = dynamic_cast<AliAODEvent*> (const_cast<TObject*>(aod));
203
204 // -- Fill MC
205 else if ( mc && readerHeader->GetUseMC() ) {
206
207 // -- if off-line MC event,
208 if ( !strcmp (mc->ClassName(),"AliMCEvent") )
209 fMC = dynamic_cast<AliMCEvent*> (const_cast<TObject*>(mc));
210 // -- if on-line MC event
211 else
212 fHLTMC = dynamic_cast<AliHLTMCEvent*> (const_cast<TObject*>(mc));
213 }
214
215 return;
216}
217
0734d112 218/*
219 * ---------------------------------------------------------------------------------
6ce099ba 220 * Fastjet Reader functionality
0734d112 221 * ---------------------------------------------------------------------------------
222 */
6ce099ba 223#ifdef HAVE_FASTJET
0734d112 224// #################################################################################
1f9fec4a 225Bool_t AliHLTJETReader::FillVector() {
0734d112 226 // see header file for class documentation
227
228 Bool_t bResult = kFALSE;
229
230 if ( fESD )
1f9fec4a 231 bResult = FillVectorESD();
0734d112 232 else if ( fMC )
1f9fec4a 233 bResult = FillVectorMC();
234 else if ( fHLTMC )
235 bResult = FillVectorHLTMC();
0734d112 236 else if ( fAOD )
1f9fec4a 237 bResult = FillVectorAOD();
0734d112 238
239 return bResult;
240}
241
242// #################################################################################
1f9fec4a 243Bool_t AliHLTJETReader::FillVectorMC() {
0734d112 244 // see header file for class documentation
245
1f9fec4a 246 Bool_t bResult = kTRUE;
247
0734d112 248 if ( ! fMC ) {
6ce099ba 249 HLTError( "No MC Event present!" );
0734d112 250 return kFALSE;
251 }
1f9fec4a 252
253 // -- Reset Event
254 ResetEvent();
255
256 Int_t nTracks = 0;
257
258 AliStack* stack = fMC->Stack();
259
260 for (Int_t iterStack = 0; iterStack < stack->GetNtrack() && bResult; iterStack++) {
261
262 TParticle *particle = stack->Particle(iterStack);
263 if ( !particle) {
264 HLTError( "Error reading particle %i out of %i", iterStack, stack->GetNtrack() );
265 bResult = kFALSE;
266 continue;
267 }
268
269 // ------------------------------
270 // -- Basic cuts on MC particle --> To be done better XXX
271 // ------------------------------
272
273 // -- primary
274 if ( !(stack->IsPhysicalPrimary(iterStack)) )
275 continue;
276
277 // -- final state
278 if ( particle->GetNDaughters() != 0 )
279 continue;
280
281 // -- particle in DB
282 TParticlePDG * particlePDG = particle->GetPDG();
283 if ( ! particlePDG ) {
284 particlePDG = TDatabasePDG::Instance()->GetParticle( particle->GetPdgCode() );
285
286 if ( ! particlePDG ) {
287 HLTError("Particle %i not in PDG database", particle->GetPdgCode() );
288 bResult = kFALSE;
289 continue;
290 }
291 }
292
293 // ------------------------
294 // -- Standard track cuts
295 // ------------------------
296
297 // -- Apply track cuts
298 if ( ! fTrackCuts->IsSelected(particle) )
299 continue;
300
301 // -- Create PseudoJet object
302 fastjet::PseudoJet part( particle->Px(), particle->Py(),
303 particle->Pz(), particle->Energy() );
304
305 // -- label the particle into Fastjet algortihm
306 part.set_user_index( iterStack );
307
3c57fb71 308#ifdef HAVE_FASTJET
1f9fec4a 309 // -- Add to input_particles vector
310 fMomentumVector->push_back(part);
3c57fb71 311#endif
1f9fec4a 312
313 nTracks++;
314 } // for (Int_t iterStack = 0; iterStack < stack->GetNtrack() && !bResult; iterStack++) {
315
316 HLTDebug(" Number of selected tracks %d", nTracks);
317
318 return kTRUE;
319}
320
321// #################################################################################
322Bool_t AliHLTJETReader::FillVectorHLTMC() {
323 // see header file for class documentation
324
325 if ( ! fHLTMC ) {
326 HLTError( "No HLT MC Event present!" );
327 return kFALSE;
328 }
329
330 // -- Reset Event
331 ResetEvent();
0734d112 332
333 Int_t nTracks = 0;
334
0734d112 335 TParticle* particle = NULL;
6ce099ba 336
0734d112 337 // -- Loop over particles
338 // ------------------------
1f9fec4a 339 while ( (particle = fHLTMC->NextParticle() ) ) {
6ce099ba 340
1f9fec4a 341 // -- Apply track cuts
6ce099ba 342 if ( ! fTrackCuts->IsSelected(particle) )
0734d112 343 continue;
6ce099ba 344
345 // -- Create PseudoJet object
346 fastjet::PseudoJet part( particle->Px(), particle->Py(),
347 particle->Pz(), particle->Energy() );
348
349 // -- label the particle into Fastjet algortihm
1f9fec4a 350 part.set_user_index( fHLTMC->GetIndex() );
6ce099ba 351
3c57fb71 352#ifdef HAVE_FASTJET
6ce099ba 353 // -- Add to input_particles vector
354 fMomentumVector->push_back(part);
3c57fb71 355#endif
0734d112 356
0734d112 357 nTracks++;
358
1f9fec4a 359 } // while ( (particle = fHLTMC->NextParticle() ) ) {
6ce099ba 360
1f9fec4a 361 HLTInfo(" Number of selected tracks %d", nTracks);
0734d112 362
363 return kTRUE;
364}
365
366// #################################################################################
1f9fec4a 367Bool_t AliHLTJETReader::FillVectorESD() {
0734d112 368 // see header file for class documentation
369
1f9fec4a 370 Bool_t bResult = kTRUE;
371
372 if ( ! fESD ) {
373 HLTError( "No ESD Event present!" );
374 return kFALSE;
375 }
376
377 // -- Reset Event
378 ResetEvent();
379
380 Int_t nTracks = 0;
381
382 // -- Loop over particles
383 // ------------------------
384 for ( Int_t iter = 0; iter < fESD->GetNumberOfTracks() && bResult; iter++ ) {
385
386 AliESDtrack* esdTrack = fESD->GetTrack(iter);
387 if ( ! esdTrack ) {
388 HLTError("Could not read ESD track %d from %d", iter, fESD->GetNumberOfTracks() );
389 bResult = kFALSE;
390 continue;
391 }
392
393 // -- Apply track cuts
394 if ( ! fTrackCuts->IsSelected(esdTrack) )
395 continue;
396
397 // -- Create PseudoJet object
398 fastjet::PseudoJet part( esdTrack->Px(), esdTrack->Py(),
399 esdTrack->Pz(), esdTrack->E() );
400
401 // -- label the particle into Fastjet algortihm
402 part.set_user_index( iter );
403
3c57fb71 404#ifdef HAVE_FASTJET
1f9fec4a 405 // -- Add to input_particles vector
406 fMomentumVector->push_back(part);
3c57fb71 407#endif
1f9fec4a 408
409 nTracks++;
410
411 } // for ( Int_t iter = 0; iter < fESD->GetNumberOfTracks() && !iResult; iter++ ) {
412
413 HLTInfo(" Number of selected tracks %d", nTracks);
414
1f9fec4a 415 return bResult;
0734d112 416}
417
418// #################################################################################
1f9fec4a 419Bool_t AliHLTJETReader::FillVectorAOD() {
0734d112 420 // see header file for class documentation
421
422 return kFALSE;
423}
6ce099ba 424#endif
0734d112 425
426/*
427 * ---------------------------------------------------------------------------------
6ce099ba 428 * Grid functionality
0734d112 429 * ---------------------------------------------------------------------------------
430 */
431
6ce099ba 432// #################################################################################
433Bool_t AliHLTJETReader::FillGrid() {
434 // see header file for class documentation
435
436 Bool_t bResult = kFALSE;
437
438 if ( fESD )
439 bResult = FillGridESD();
440 else if ( fMC )
441 bResult = FillGridMC();
1f9fec4a 442 else if ( fHLTMC )
443 bResult = FillGridHLTMC();
6ce099ba 444 else if ( fAOD )
445 bResult = FillGridAOD();
1f9fec4a 446
6ce099ba 447 return bResult;
448}
449
450// #################################################################################
451Bool_t AliHLTJETReader::FillGridMC() {
452 // see header file for class documentation
453
1f9fec4a 454 Bool_t bResult = kTRUE;
455
6ce099ba 456 if ( ! fMC ) {
457 HLTError( "No MC Event present!" );
458 return kFALSE;
459 }
460
1f9fec4a 461 // -- Reset Event
462 ResetEvent();
463
464 AliHLTJETReaderHeader* readerHeader = GetReaderHeader();
465
466 Int_t nTracks = 0;
467
468 AliStack* stack = fMC->Stack();
469
470 for (Int_t iterStack = 0; iterStack < stack->GetNtrack() && bResult; iterStack++) {
471
472 TParticle *particle = stack->Particle(iterStack);
473 if ( !particle) {
474 HLTError( "Error reading particle %i out of %i", iterStack, stack->GetNtrack() );
475 bResult = kFALSE;
476 continue;
477 }
478
479 // ------------------------------
480 // -- Basic cuts on MC particle --> To be done better XXX
481 // ------------------------------
482
483 // -- primary
484 if ( !(stack->IsPhysicalPrimary(iterStack)) )
485 continue;
486
487 // -- final state
488 if ( particle->GetNDaughters() != 0 )
489 continue;
490
491 // -- particle in DB
492 TParticlePDG * particlePDG = particle->GetPDG();
493 if ( ! particlePDG ) {
494 particlePDG = TDatabasePDG::Instance()->GetParticle( particle->GetPdgCode() );
495
496 if ( ! particlePDG ) {
497 HLTError("Particle %i not in PDG database", particle->GetPdgCode() );
498 bResult = kFALSE;
499 continue;
500 }
501 }
502
503 // ------------------------
504 // -- Standard track cuts
505 // ------------------------
506
507 // -- Apply track cuts
508 if ( ! fTrackCuts->IsSelected(particle) )
509 continue;
510
2cb809a5 511 const Float_t aEtaPhi[] = { static_cast<Float_t>(particle->Eta()), static_cast<Float_t>(particle->Phi()), static_cast<Float_t>(particle->Pt()) };
1f9fec4a 512 Int_t aGridIdx[] = { -1, -1, -1, -1, -1 };
513
514 fGrid->FillTrack(particle, aEtaPhi, aGridIdx);
515
516 nTracks++;
517
518 // -- Apply seed cuts
519 if ( ! fSeedCuts->IsSelected(particle) )
520 continue;
521
522 // -- Add Seed
523 AddSeed(aEtaPhi, const_cast<const Int_t*> (aGridIdx), readerHeader->GetConeRadius());
524
525 } // for (Int_t iterStack = 0; iterStack < stack->GetNtrack() && !bResult; iterStack++) {
526
527 HLTInfo(" Number of selected tracks %d", nTracks);
528 HLTInfo(" Number of seeds %d", fNJetCandidates);
529
530 return kTRUE;
531}
532
533// #################################################################################
534Bool_t AliHLTJETReader::FillGridHLTMC() {
535 // see header file for class documentation
536
537 if ( ! fHLTMC ) {
538 HLTError( "No HLT MC Event present!" );
539 return kFALSE;
540 }
541
542 // -- Reset Event
543 ResetEvent();
544
6ce099ba 545 AliHLTJETReaderHeader* readerHeader = GetReaderHeader();
546
547 Int_t nTracks = 0;
548 TParticle* particle = NULL;
549
550 // -- Loop over particles
551 // ------------------------
1f9fec4a 552 while ( ( particle = fHLTMC->NextParticle() ) ) {
6ce099ba 553
1f9fec4a 554 // HLTError("=== nTracks %d ===",nTracks);
555
6ce099ba 556 // -- Apply track cuts
557 if ( ! fTrackCuts->IsSelected(particle) )
558 continue;
559
2cb809a5 560 const Float_t aEtaPhi[] = { static_cast<Float_t>(particle->Eta()), static_cast<Float_t>(particle->Phi()), static_cast<Float_t>(particle->Pt()) };
6ce099ba 561 Int_t aGridIdx[] = { -1, -1, -1, -1, -1 };
562
1f9fec4a 563 // -- Fill grid
6ce099ba 564 fGrid->FillTrack(particle, aEtaPhi, aGridIdx);
565
566 nTracks++;
567
568 // -- Apply seed cuts
569 if ( ! fSeedCuts->IsSelected(particle) )
570 continue;
571
572 // -- Add Seed
573 AddSeed(aEtaPhi, const_cast<const Int_t*> (aGridIdx),
574 readerHeader->GetConeRadius());
575
1f9fec4a 576 } // while ( (particle = fHLTMC->NextParticle() ) ) {
6ce099ba 577
1f9fec4a 578 HLTInfo(" Number of selected tracks %d", nTracks);
579 HLTInfo(" Number of seeds %d", fNJetCandidates);
6ce099ba 580
581 return kTRUE;
582}
583
584// #################################################################################
585Bool_t AliHLTJETReader::FillGridESD() {
586 // see header file for class documentation
587
588 Bool_t bResult = kTRUE;
589
590 if ( ! fESD ) {
591 HLTError( "No ESD Event present!" );
592 return kFALSE;
593 }
594
1f9fec4a 595 // -- Reset Event
596 ResetEvent();
597
6ce099ba 598 AliHLTJETReaderHeader* readerHeader = GetReaderHeader();
599
600 Int_t nTracks = 0;
601
602 // -- Loop over particles
603 // ------------------------
f5561f3a 604 // for ( Int_t iter = 0; iter < fESD->GetNumberOfTracks() && !bResult; iter++ ) { JMT Coverity
605 for ( Int_t iter = 0; iter < fESD->GetNumberOfTracks() && bResult; iter++ ) {
6ce099ba 606
607 AliESDtrack* esdTrack = fESD->GetTrack(iter);
608 if ( ! esdTrack ) {
1f9fec4a 609 HLTError("Could not read ESD track %d from %d", iter, fESD->GetNumberOfTracks() );
6ce099ba 610 bResult = kFALSE;
611 continue;
612 }
613
614 // -- Apply track cuts
615 if ( ! fTrackCuts->IsSelected(esdTrack) )
616 continue;
617
2cb809a5 618 const Float_t aEtaPhi[] = { static_cast<Float_t>(esdTrack->Eta()), static_cast<Float_t>(esdTrack->Phi()), static_cast<Float_t>(esdTrack->Pt()) };
6ce099ba 619 Int_t aGridIdx[] = { -1, -1, -1, -1, -1 };
620
621 // -- Fill grid
1f9fec4a 622 fGrid->FillTrack(esdTrack, aEtaPhi, aGridIdx);
6ce099ba 623
624 nTracks++;
625
626 // -- Apply seed cuts
627 if ( ! fSeedCuts->IsSelected(esdTrack) )
628 continue;
629
630 // -- Add Seed
631 AddSeed(aEtaPhi, const_cast<const Int_t*> (aGridIdx),
632 readerHeader->GetConeRadius());
633
634 } // for ( Int_t iter = 0; iter < fESD->GetNumberOfTracks() && !iResult; iter++ ) {
635
1f9fec4a 636 HLTInfo(" Number of selected tracks %d", nTracks);
637 HLTInfo(" Number of seeds %d", fNJetCandidates);
6ce099ba 638
639 return bResult;
640}
641
642// #################################################################################
643Bool_t AliHLTJETReader::FillGridAOD() {
644 // see header file for class documentation
645
646 return kFALSE;
647}
648
6ce099ba 649/*
650 * ---------------------------------------------------------------------------------
651 * Seeds
652 * ---------------------------------------------------------------------------------
653 */
654
655//#################################################################################
656void AliHLTJETReader::AddSeed( const Float_t* aEtaPhi, const Int_t* aGridIdx,
657 const Float_t coneRadius ) {
658 // see header file for class documentation
659
1f9fec4a 660 Bool_t useWholeCell = kTRUE;
6ce099ba 661
1f9fec4a 662 if ( GetReaderHeader()->GetJetAlgorithm() == AliHLTJETBase::kFFSCRadiusCell )
663 useWholeCell = kFALSE ;
664
665 // -- Add track / particle
6ce099ba 666 new( (*fJetCandidates) [fNJetCandidates] ) AliHLTJETConeJetCandidate( aEtaPhi,
667 aGridIdx,
668 coneRadius,
669 useWholeCell );
670 fNJetCandidates++;
671
672 HLTDebug("Added Seed Pt=%f, Eta=%f, Phi=%f", aEtaPhi[kIdxPt],
1f9fec4a 673 aEtaPhi[kIdxEta], aEtaPhi[kIdxPhi] );
6ce099ba 674
675 return;
676}
1f9fec4a 677
678/*
679 * ---------------------------------------------------------------------------------
680 * Initialize - private
681 * ---------------------------------------------------------------------------------
682 */
683
684// #################################################################################
685Int_t AliHLTJETReader::InitializeFFSC() {
686 // see header file for class documentation
687
688 Int_t iResult = 0;
689 AliHLTJETReaderHeader* readerHeader = GetReaderHeader();
690
691 // -- Initialize grid
692 // --------------------
693 if ( fGrid )
694 delete fGrid;
695
696 if ( ! (fGrid = new AliHLTJETConeGrid()) ) {
697 HLTError("Error instanciating grid.");
698 iResult = -EINPROGRESS;
699 }
700
701 if ( ! iResult ) {
702 fGrid->SetEtaRange( readerHeader->GetFiducialEtaMin(),
703 readerHeader->GetFiducialEtaMax(),
704 readerHeader->GetGridEtaRange() );
705
706 fGrid->SetPhiRange( readerHeader->GetFiducialPhiMin(),
707 readerHeader->GetFiducialPhiMax(),
708 readerHeader->GetGridPhiRange() );
709
710 fGrid->SetBinning( readerHeader->GetGridEtaBinning(),
711 readerHeader->GetGridEtaBinning() );
712
713 fGrid->SetConeRadius( readerHeader->GetConeRadius() );
714
715 iResult = fGrid->Initialize();
716 }
717
718 // -- Initialize jet candidates
719 // ------------------------------
720 if ( ! iResult ) {
721 fJetCandidates = new TClonesArray("AliHLTJETConeJetCandidate", 30);
722 if ( ! fJetCandidates) {
723 HLTError("Error instanciating jet candidates.");
724 iResult = -EINPROGRESS;
725 }
726 }
438118d4 727
728 // -- Get ptr to cuts from reader
729 // --------------------------------
730
731 // -- Seed cuts
732 if ( ! iResult ) {
733 fSeedCuts = readerHeader->GetSeedCuts();
734 if ( ! fSeedCuts ) {
735 HLTError("Error getting ptr to seed cuts.");
736 iResult = -EINPROGRESS;
737 }
738 }
1f9fec4a 739
740 return iResult;
741}
742
743#ifdef HAVE_FASTJET
744// #################################################################################
745Int_t AliHLTJETReader::InitializeFastjet() {
746 // see header file for class documentation
747
748 Int_t iResult = 0;
749
750 // -- Initialize Vector
751 // ----------------------
752 if ( fMomentumVector )
753 delete fMomentumVector;
754
755 if ( ! (fMomentumVector = new vector<fastjet::PseudoJet>) ) {
756 HLTError("Error instanciating momentum vector.");
757 iResult = -EINPROGRESS;
758 }
759
760 return iResult;
761}
762#endif