]>
Commit | Line | Data |
---|---|---|
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 | ||
16 | ||
17 | /* $Id$ */ | |
18 | ||
19 | //_________________________________________________________________________ | |
20 | // Base Class for JetFinder Algorithms | |
21 | // -- | |
22 | //*-- Author: Mark Horner (LBL/UCT) | |
23 | // -- | |
24 | ||
25 | ||
26 | ||
27 | #include "AliEMCALJetFinderAlgo.h" | |
28 | #include "AliRun.h" | |
29 | #include "AliEMCALGeometry.h" | |
30 | #include "AliMagF.h" | |
31 | ClassImp(AliEMCALJetFinderAlgo) | |
32 | ||
33 | AliEMCALJetFinderAlgo::AliEMCALJetFinderAlgo() | |
34 | { | |
35 | fDebug =0; | |
36 | fOutputPointer=0; | |
37 | fInputPointer=0; | |
38 | fOutputAllocated=kFALSE; | |
39 | } | |
40 | ||
41 | AliEMCALJetFinderAlgo::~AliEMCALJetFinderAlgo() | |
42 | { | |
43 | if (fOutputAllocated) | |
44 | delete fOutputPointer; | |
45 | } | |
46 | ||
47 | void AliEMCALJetFinderAlgo::SetOutput(AliEMCALJetFinderOutput* output) { | |
48 | if (fOutputAllocated) | |
49 | delete fOutputPointer; | |
50 | fOutputPointer=output; | |
51 | fOutputAllocated=kFALSE; | |
52 | } | |
53 | ||
54 | void AliEMCALJetFinderAlgo::InitInput(AliEMCALJetFinderInput* input) | |
55 | { | |
56 | // Take input data | |
57 | if (fDebug>1) Info("InitInput","Beginning InitInput"); | |
58 | fInputPointer = input; | |
59 | if (fOutputPointer==0) { | |
60 | ||
61 | if (fDebug>1) Info("InitInput","Allocating output object"); | |
62 | fOutputPointer=new AliEMCALJetFinderOutput(); | |
63 | fOutputAllocated=kTRUE; | |
64 | } | |
65 | fOutputPointer->Reset(kResetAll); | |
66 | // automatically copy parton and particle info to output object | |
67 | ||
68 | for (Int_t counter = 0 ; counter < fInputPointer->GetNPartons();counter++) | |
69 | { | |
70 | fOutputPointer->AddParton(fInputPointer->GetParton(counter)); | |
71 | } | |
72 | for (Int_t counter = 0 ; counter < fInputPointer->GetNParticles();counter++) | |
73 | { | |
74 | fOutputPointer->AddParticle(fInputPointer->GetParticle(counter)); | |
75 | } | |
76 | } | |
77 | ||
78 | Float_t AliEMCALJetFinderAlgo::PropagatePhi(Float_t pt, Float_t charge, Bool_t& curls) | |
79 | { | |
80 | // Propagates phi angle to EMCAL radius | |
81 | // // | |
82 | Float_t b = 0.0, rEMCAL = -1.0; | |
83 | if(rEMCAL<0) | |
84 | { | |
85 | b = gAlice->Field()->SolenoidField(); | |
86 | rEMCAL = AliEMCALGeometry::GetInstance()->GetIPDistance(); | |
87 | } | |
88 | Float_t dPhi = 0.; | |
89 | Float_t rB = 3335.6 * pt / b; // [cm] (case of |charge|=1) | |
90 | if (2.*rB < rEMCAL) | |
91 | { | |
92 | curls = kTRUE; | |
93 | return dPhi; | |
94 | } | |
95 | Float_t phi = TMath::ACos(1.-rEMCAL*rEMCAL/(2.*rB*rB)); | |
96 | dPhi = TMath::ATan2(1.-TMath::Cos(phi), TMath::Sin(phi)); | |
97 | dPhi = -TMath::Sign(dPhi, charge); | |
98 | return dPhi; | |
99 | } | |
100 |