1 // Author: Uli Frankenfeld <mailto:franken@fi.uib.no>
2 //*-- Copyright © Uli
6 #include "AliL3Logging.h"
7 #include "AliL3InterMerger.h"
8 #include "AliL3Track.h"
9 #include "AliL3TrackSegmentData.h"
10 #include "AliL3Transform.h"
11 #include "AliL3TrackArray.h"
13 //_____________________________________________________________
16 // The L3 track segment merger
19 ClassImp(AliL3InterMerger)
21 AliL3InterMerger::AliL3InterMerger():AliL3Merger(1){
24 // SetParameter(2,2,.3,.3,.3);
25 SetParameter(1,0.5,0.0005,0.05,0.1);
26 fRowMax = fRowMin = 0;
30 AliL3InterMerger::~AliL3InterMerger(){
35 void AliL3InterMerger::SlowMerge(){
36 Int_t nrow= fRowMax-fRowMin+1;
37 void *ntuple=GetNtuple();
38 AliL3TrackArray * tracks = GetInTracks(0);
39 const Int_t kNIn =tracks->GetNTracks();
42 for(Int_t in=0;in<kNIn;in++)
43 tracks->GetCheckedTrack(in)->CalculateHelix();
45 Int_t inmin=-1,outmin=-1;
47 for(Int_t out=0;out<kNIn;out++){
48 AliL3Track *outertrack=tracks->GetCheckedTrack(out);
49 if(!outertrack) continue;
50 for(Int_t in=0;in<kNIn;in++){
52 AliL3Track *innertrack=tracks->GetCheckedTrack(in);
53 if(!innertrack) continue;
54 if(outertrack->GetNHits()+innertrack->GetNHits()>nrow) continue;
56 Double_t diff = TrackDiff(innertrack,outertrack);
58 if(diff>=0&&diff<min){
65 if(inmin>=0&&outmin>=0){
66 AliL3Track *outertrack=tracks->GetTrack(outmin);
67 AliL3Track *innertrack=tracks->GetTrack(inmin);
71 MultiMerge(tracks,tr,2);
72 outertrack->CalculatePoint(tr[0]->GetLastPointX());
73 innertrack->CalculatePoint(tr[0]->GetLastPointX());
74 PrintDiff(innertrack,outertrack);
75 FillNtuple(ntuple,innertrack,outertrack);
76 tracks->Remove(outmin);
77 tracks->Remove(inmin);
81 LOG(AliL3Log::kInformational,"AliL3InterMerger::SlowMerge","Result")
82 <<AliL3Log::kDec<<"Merged Tracks: "<<tracks->GetNTracks()-kNIn<<ENDLOG;
85 sprintf(name,"ntuple_i_%d.root",fPatch);
86 WriteNtuple(name,ntuple);
89 void AliL3InterMerger::MMerge(){
91 GetOutTracks()->AddTracks(GetInTracks(0));
94 Int_t AliL3InterMerger::Merge(){
95 Int_t nrow= fRowMax-fRowMin+1;
96 Double_t xval =AliL3Transform::Row2X((fRowMax+fRowMin)/2);
97 AliL3TrackArray * tracks = GetInTracks(0);
98 const Int_t kNIn =tracks->GetNTracks();
100 for(Int_t in=0;in<kNIn;in++){
101 AliL3Track *t = tracks->GetCheckedTrack(in);
104 t->CalculatePoint(xval);
107 for(Int_t out=0;out<kNIn;out++){
108 AliL3Track *outertrack=tracks->GetCheckedTrack(out);
109 if(!outertrack) continue;
110 for(Int_t in=0;in<kNIn;in++){
111 if(in==out) continue;
112 AliL3Track *innertrack=tracks->GetCheckedTrack(in);
113 if(!innertrack) continue;
114 if(outertrack->GetNHits()+innertrack->GetNHits()>nrow) continue;
116 if(IsTrack(innertrack,outertrack)){
119 SortGlobalTracks(tr,2);
121 Double_t r0 = pow(tr[0]->GetLastPointX(),2)+
122 pow(tr[0]->GetLastPointY(),2);
123 Double_t r1 = pow(tr[1]->GetFirstPointX(),2)+
124 pow(tr[1]->GetFirstPointY(),2);
126 MultiMerge(tracks,tr,2);
134 Int_t nmerged = tracks->GetNTracks()-kNIn;
135 LOG(AliL3Log::kInformational,"AliL3InterMerger::Merge","Result")
136 <<AliL3Log::kDec<<"Merged Tracks: "<<nmerged<<ENDLOG;
138 // GetOutTracks()->AddTracks(GetInTracks(0));