1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 ////////////////////////////////////////////////////////////////////////////
20 // Author: J. Klein (Jochen.Klein@cern.ch) //
22 ////////////////////////////////////////////////////////////////////////////
24 /* $Id: AliTRDtrackGTU.cxx 27566 2008-07-24 15:31:08Z cblume $ */
27 #include "TObjArray.h"
31 #include "AliESDTrdTrack.h"
33 #include "AliTRDgtuParam.h"
34 #include "AliTRDtrackGTU.h"
35 #include "AliTRDtrackletGTU.h"
36 #include "AliTRDtrackletMCM.h"
37 #include "AliESDTrdTrack.h"
39 ClassImp(AliTRDtrackGTU)
41 AliTRDtrackGTU::AliTRDtrackGTU() :
60 fTracklets = new TClonesArray("AliTRDtrackletGTU", 6);
61 for (Int_t iTracklet = 0; iTracklet < 6; iTracklet++)
62 new ((*fTracklets)[iTracklet]) AliTRDtrackletGTU();
63 // fTracklets->BypassStreamer(kFALSE);
66 AliTRDtrackGTU::~AliTRDtrackGTU()
74 void AliTRDtrackGTU::AddTracklet(const AliTRDtrackletGTU * const tracklet, Int_t layer)
76 // add a tracklet to this track
78 if ( (fTrackletMask & (1 << layer)) != 0 ) {
79 AliError(Form("Only one tracklet per layer (%i) possible! Mask: 0x%02x", layer, fTrackletMask));
83 new ((*fTracklets)[layer]) AliTRDtrackletGTU(*tracklet);
85 fTrackletMask |= (1 << layer);
88 AliTRDtrackletGTU* AliTRDtrackGTU::GetTracklet(Int_t layer)
90 // get a pointer to the tracklet in the layer specified
92 if (IsTrackletInLayer(layer))
93 return ((AliTRDtrackletGTU*) (*fTracklets)[layer]);
98 Int_t AliTRDtrackGTU::GetNTracklets() const
100 // returns the number of tracklets in this track
105 Bool_t AliTRDtrackGTU::IsTrackletInLayer(Int_t layer) const
107 // checks for a tracklet in the given layer
109 if ( (GetTrackletMask() & (1 << layer)) != 0)
115 void AliTRDtrackGTU::SetFitParams(Float_t a, Float_t b, Float_t c)
117 // set the fit parameters
124 Int_t AliTRDtrackGTU::GetZSubChannel()
126 // returns the z-subchannel
128 if (fZSubChannel < 0) {
129 for (Int_t layer = 0; layer < AliTRDgtuParam::GetNLayers(); layer++)
131 if (IsTrackletInLayer(layer))
132 fZSubChannel = ((AliTRDtrackletGTU*) (*fTracklets)[layer])->GetSubChannel(GetZChannel());
138 Int_t AliTRDtrackGTU::GetYapprox()
140 // returns an approximated y-position for the track
142 for (Int_t layer = 0; layer < AliTRDgtuParam::GetNLayers(); layer++)
144 if (IsTrackletInLayer(layer))
145 return ((AliTRDtrackletGTU*) (*fTracklets)[layer])->GetYProj();
150 AliESDTrdTrack* AliTRDtrackGTU::CreateTrdTrack() const
152 // creates an AliESDTrdTrack to be added to the ESD
154 AliESDTrdTrack *trk = new AliESDTrdTrack();
155 trk->SetPt(1./128. * fPt);
157 trk->SetDetector((Char_t) (fSector * 5 + fStack));
159 trk->SetLabel(fLabel);
160 AliInfo(Form("setting detector to: %i (sector: %i, stack: %i), readback: %i", fSector * 30 + fStack * 6, fSector, fStack, trk->GetDetector()));
164 Bool_t AliTRDtrackGTU::CookLabel()
166 TH1F *h = new TH1F("trkref", "trkref", 100000, 0, 100000);
167 for (Int_t iTracklet = 0; iTracklet < 6; iTracklet++) {
168 if (!IsTrackletInLayer(iTracklet))
170 h->Fill( ((AliTRDtrackletGTU*) (*fTracklets)[iTracklet])->GetLabel());
172 if (h->GetEntries() > 0)
173 fLabel = h->GetMaximumBin() - 1;
177 return (fLabel >= 0);