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 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////////////////
20 // Class for TRD reconstruction //
22 ///////////////////////////////////////////////////////////////////////////////
26 #include "AliRunLoader.h"
27 #include "AliRawReader.h"
29 #include "AliESDTrdTrack.h"
32 #include "AliTRDReconstructor.h"
33 #include "AliTRDclusterizerV1.h"
34 #include "AliTRDtracker.h"
35 #include "AliTRDpidESD.h"
36 #include "AliTRDtrigger.h"
37 #include "AliTRDtrigParam.h"
38 #include "AliTRDgtuTrack.h"
39 #include "AliTRDrawData.h"
40 #include "AliTRDdigitsManager.h"
42 ClassImp(AliTRDReconstructor)
44 Bool_t AliTRDReconstructor::fgkSeedingOn = kFALSE;
45 Int_t AliTRDReconstructor::fgStreamLevel = 0; // Stream (debug) level
47 //_____________________________________________________________________________
48 void AliTRDReconstructor::ConvertDigits(AliRawReader *rawReader
49 , TTree *digitsTree) const
52 // Convert raw data digits into digit objects in a root tree
55 AliTRDrawData rawData;
56 AliTRDdigitsManager *manager = rawData.Raw2Digits(rawReader);
57 manager->MakeBranch(digitsTree);
58 manager->WriteDigits();
62 //_____________________________________________________________________________
63 void AliTRDReconstructor::Reconstruct(AliRunLoader *runLoader
64 , AliRawReader *rawReader) const
67 // Reconstruct clusters
70 AliInfo("Reconstruct TRD clusters from RAW data");
72 AliLoader *loader = runLoader->GetLoader("TRDLoader");
73 loader->LoadRecPoints("recreate");
75 runLoader->CdGAFile();
76 Int_t nEvents = runLoader->GetNumberOfEvents();
78 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
79 if (!rawReader->NextEvent()) break;
80 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
81 clusterer.Open(runLoader->GetFileName(),iEvent);
82 clusterer.ReadDigits(rawReader);
83 clusterer.MakeClusters();
84 clusterer.WriteClusters(-1);
87 loader->UnloadRecPoints();
90 // Trigger (tracklets, LTU)
92 // loader->LoadTracks();
93 // if (loader->TreeT()) {
94 // AliError("Tracklets already exist");
97 // AliInfo("Trigger tracklets will be produced");
99 // AliTRDtrigger trdTrigger("Trigger","Trigger class");
101 // AliTRDtrigParam *trigp = new AliTRDtrigParam("TRDtrigParam"
102 // ,"TRD Trigger parameters");
104 // Float_t field = AliTracker::GetBz() * 0.1; // Tesla
105 // AliInfo(Form("Trigger set for magnetic field = %f Tesla \n",field));
106 // trigp->SetField(field);
108 // trdTrigger.SetParameter(trigp);
110 // rawReader->RewindEvents();
112 // for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
113 // if (!rawReader->NextEvent()) break;
114 // trdTrigger.Open(runLoader->GetFileName(),iEvent);
115 // trdTrigger.ReadDigits(rawReader);
116 // trdTrigger.MakeTracklets();
117 // trdTrigger.WriteTracklets(-1);
120 // loader->UnloadTracks();
124 //_____________________________________________________________________________
125 void AliTRDReconstructor::Reconstruct(AliRawReader *rawReader
126 , TTree *clusterTree) const
129 // Reconstruct clusters
132 AliInfo("Reconstruct TRD clusters from RAW data");
134 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
135 clusterer.OpenOutput(clusterTree);
136 clusterer.ReadDigits(rawReader);
137 clusterer.MakeClusters();
140 // No trigger, since we don't have the output tree for tracklets
145 //_____________________________________________________________________________
146 void AliTRDReconstructor::Reconstruct(TTree *digitsTree
147 , TTree *clusterTree) const
150 // Reconstruct clusters
153 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
154 clusterer.OpenOutput(clusterTree);
155 clusterer.ReadDigits(digitsTree);
156 clusterer.MakeClusters();
159 // No trigger, since we don't have the output tree for tracklets
164 //_____________________________________________________________________________
165 void AliTRDReconstructor::Reconstruct(AliRunLoader *runLoader) const
168 // Reconstruct clusters
171 AliLoader *loader = runLoader->GetLoader("TRDLoader");
172 loader->LoadRecPoints("recreate");
174 runLoader->CdGAFile();
175 Int_t nEvents = runLoader->GetNumberOfEvents();
177 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
178 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
179 clusterer.Open(runLoader->GetFileName(),iEvent);
180 clusterer.ReadDigits();
181 clusterer.MakeClusters();
182 clusterer.WriteClusters(-1);
185 loader->UnloadRecPoints();
188 // Trigger (tracklets, LTU)
190 loader->LoadTracks("RECREATE");
191 AliInfo("Trigger tracklets will be produced");
193 AliTRDtrigger trdTrigger("Trigger","Trigger class");
195 AliTRDtrigParam *trigp = new AliTRDtrigParam("TRDtrigParam"
196 ,"TRD Trigger parameters");
198 Float_t field = AliTracker::GetBz() * 0.1; // Tesla
199 AliInfo(Form("Trigger set for magnetic field = %f Tesla \n",field));
200 trigp->SetField(field);
202 trdTrigger.SetParameter(trigp);
204 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
205 trdTrigger.Open(runLoader->GetFileName(),iEvent);
206 trdTrigger.ReadDigits();
207 trdTrigger.MakeTracklets();
208 trdTrigger.WriteTracklets(-1);
211 loader->UnloadTracks();
215 //_____________________________________________________________________________
216 AliTracker *AliTRDReconstructor::CreateTracker(AliRunLoader *runLoader) const
219 // Create a TRD tracker
222 runLoader->CdGAFile();
224 return new AliTRDtracker(gFile);
228 //_____________________________________________________________________________
229 void AliTRDReconstructor::FillESD(AliRunLoader* /*runLoader*/
230 , AliRawReader* /*rawReader*/
241 // No trigger, since we don't have the output tree for tracklets
246 //_____________________________________________________________________________
247 void AliTRDReconstructor::FillESD(AliRawReader* /*rawReader*/
248 , TTree* /*clusterTree*/
259 // No trigger, since we don't have the output tree for tracklets
264 //_____________________________________________________________________________
265 void AliTRDReconstructor::FillESD(TTree* /*digitsTree*/
266 , TTree* /*clusterTree*/
277 // No trigger, since we don't have the output tree for tracklets
282 //_____________________________________________________________________________
283 void AliTRDReconstructor::FillESD(AliRunLoader *runLoader
294 // Trigger (tracks, GTU)
296 // AliTRDtrigger trdTrigger("Trigger","Trigger class");
298 // AliTRDtrigParam *trigp = new AliTRDtrigParam("TRDtrigParam"
299 // ,"TRD Trigger parameters");
301 // Float_t field = AliTracker::GetBz() * 0.1; // Tesla
302 // AliInfo(Form("Trigger set for magnetic field = %f Tesla \n",field));
303 // trigp->SetField(field);
306 // trdTrigger.SetParameter(trigp);
307 // trdTrigger.SetRunLoader(runLoader);
308 // trdTrigger.Init();
310 // Int_t iEvent = runLoader->GetEventNumber();
311 // runLoader->GetEvent(iEvent);
312 // trdTrigger.ReadTracklets(runLoader);
314 // AliESDTrdTrack *TrdTrack = new AliESDTrdTrack();
315 // AliTRDgtuTrack *GtuTrack;
317 // Int_t nTracks = trdTrigger.GetNumberOfTracks();
318 // for (Int_t iTrack = 0; iTrack < nTracks; iTrack++) {
320 // GtuTrack = trdTrigger.GetTrack(iTrack);
322 // TrdTrack->SetYproj(GtuTrack->GetYproj());
323 // TrdTrack->SetZproj(GtuTrack->GetZproj());
324 // TrdTrack->SetSlope(GtuTrack->GetSlope());
325 // TrdTrack->SetDetector(GtuTrack->GetDetector());
326 // TrdTrack->SetTracklets(GtuTrack->GetTracklets());
327 // TrdTrack->SetPlanes(GtuTrack->GetPlanes());
328 // TrdTrack->SetClusters(GtuTrack->GetClusters());
329 // TrdTrack->SetPt(GtuTrack->GetPt());
330 // TrdTrack->SetPhi(GtuTrack->GetPhi());
331 // TrdTrack->SetEta(GtuTrack->GetEta());
332 // TrdTrack->SetLabel(GtuTrack->GetLabel());
333 // TrdTrack->SetPID(GtuTrack->GetPID());
334 // TrdTrack->SetIsElectron(GtuTrack->IsElectron());
336 // esd->AddTrdTrack(TrdTrack);