3 // Author: Uli Frankenfeld <mailto:franken@fi.uib.no>
4 //*-- Copyright © ALICE HLT Group
6 #include "AliL3StandardIncludes.h"
8 #include "AliL3Logging.h"
9 #include "AliL3InterMerger.h"
10 #include "AliL3Track.h"
11 #include "AliL3TrackSegmentData.h"
12 #include "AliL3Transform.h"
13 #include "AliL3TrackArray.h"
15 /** \class AliL3InterMerger
17 //_____________________________________________________________
20 // The L3 track segment merger
29 ClassImp(AliL3InterMerger)
31 AliL3InterMerger::AliL3InterMerger()
36 // SetParameter(2,2,.3,.3,.3);
37 SetParameter(1,0.5,0.0005,0.05,0.1);
38 fRowMax = fRowMin = 0;
42 AliL3InterMerger::~AliL3InterMerger(){
47 void AliL3InterMerger::SlowMerge(){
48 Int_t nrow= fRowMax-fRowMin+1;
49 void *ntuple=GetNtuple();
50 AliL3TrackArray * tracks = GetInTracks(0);
51 const Int_t kNIn =tracks->GetNTracks();
54 for(Int_t in=0;in<kNIn;in++)
55 tracks->GetCheckedTrack(in)->CalculateHelix();
57 Int_t inmin=-1,outmin=-1;
59 for(Int_t out=0;out<kNIn;out++){
60 AliL3Track *outertrack=tracks->GetCheckedTrack(out);
61 if(!outertrack) continue;
62 for(Int_t in=0;in<kNIn;in++){
64 AliL3Track *innertrack=tracks->GetCheckedTrack(in);
65 if(!innertrack) continue;
66 if(outertrack->GetNHits()+innertrack->GetNHits()>nrow) continue;
68 Double_t diff = TrackDiff(innertrack,outertrack);
70 if(diff>=0&&diff<min){
77 if(inmin>=0&&outmin>=0){
78 AliL3Track *outertrack=tracks->GetTrack(outmin);
79 AliL3Track *innertrack=tracks->GetTrack(inmin);
83 MultiMerge(tracks,tr,2);
84 outertrack->CalculatePoint(tr[0]->GetLastPointX());
85 innertrack->CalculatePoint(tr[0]->GetLastPointX());
86 PrintDiff(innertrack,outertrack);
87 FillNtuple(ntuple,innertrack,outertrack);
88 tracks->Remove(outmin);
89 tracks->Remove(inmin);
93 LOG(AliL3Log::kInformational,"AliL3InterMerger::SlowMerge","Result")
94 <<AliL3Log::kDec<<"Merged Tracks: "<<tracks->GetNTracks()-kNIn<<ENDLOG;
97 sprintf(name,"ntuple_i_%d.root",fPatch);
98 WriteNtuple(name,ntuple);
101 void AliL3InterMerger::MMerge(){
103 GetOutTracks()->AddTracks(GetInTracks(0));
106 Int_t AliL3InterMerger::Merge(){
107 Int_t nrow= fRowMax-fRowMin+1;
108 Double_t xval =AliL3Transform::Row2X((fRowMax+fRowMin)/2);
109 AliL3TrackArray * tracks = GetInTracks(0);
110 const Int_t kNIn =tracks->GetNTracks();
112 for(Int_t in=0;in<kNIn;in++){
113 AliL3Track *t = tracks->GetCheckedTrack(in);
116 t->CalculatePoint(xval);
119 for(Int_t out=0;out<kNIn;out++){
120 AliL3Track *outertrack=tracks->GetCheckedTrack(out);
121 if(!outertrack) continue;
122 for(Int_t in=0;in<kNIn;in++){
123 if(in==out) continue;
124 AliL3Track *innertrack=tracks->GetCheckedTrack(in);
125 if(!innertrack) continue;
126 if(outertrack->GetNHits()+innertrack->GetNHits()>nrow) continue;
128 if(IsTrack(innertrack,outertrack)){
131 SortGlobalTracks(tr,2);
133 Double_t r0 = pow(tr[0]->GetLastPointX(),2)+
134 pow(tr[0]->GetLastPointY(),2);
135 Double_t r1 = pow(tr[1]->GetFirstPointX(),2)+
136 pow(tr[1]->GetFirstPointY(),2);
138 MultiMerge(tracks,tr,2);
146 Int_t nmerged = tracks->GetNTracks()-kNIn;
147 LOG(AliL3Log::kInformational,"AliL3InterMerger::Merge","Result")
148 <<AliL3Log::kDec<<"Merged Tracks: "<<nmerged<<ENDLOG;
150 // GetOutTracks()->AddTracks(GetInTracks(0));