]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - TRD/AliTRD.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (non-incremental) | history | HEAD
Config.C modified, protection in AliMuonForwardTrackFinder.cxx in case of zero MUON...
[u/mrichter/AliRoot.git] / TRD / AliTRD.cxx
... / ...
CommitLineData
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/* $Id$ */
17
18///////////////////////////////////////////////////////////////////////////////
19// //
20// Transition Radiation Detector //
21// This class contains the basic functions for the Transition Radiation //
22// Detector. //
23// //
24///////////////////////////////////////////////////////////////////////////////
25
26#include <TGeoGlobalMagField.h>
27#include <TVirtualMC.h>
28#include <TGeoManager.h>
29
30#include "AliMC.h"
31#include "AliMagF.h"
32#include "AliRun.h"
33#include "AliLoader.h"
34
35#include "AliTRD.h"
36#include "AliTRDdigitizer.h"
37#include "AliTRDdigitsManager.h"
38#include "AliTRDgeometry.h"
39#include "AliTRDhit.h"
40#include "AliTRDrawData.h"
41#include "AliTRDCommonParam.h"
42
43ClassImp(AliTRD)
44
45//_____________________________________________________________________________
46AliTRD::AliTRD()
47 :AliDetector()
48 ,fGeometry(0)
49 ,fGasDensity(0)
50 ,fFoilDensity(0)
51 ,fGasNobleFraction(0)
52 ,fPrimaryIonisation(0)
53{
54 //
55 // Default constructor
56 //
57
58}
59
60//_____________________________________________________________________________
61AliTRD::AliTRD(const char *name, const char *title)
62 :AliDetector(name,title)
63 ,fGeometry(0)
64 ,fGasDensity(0)
65 ,fFoilDensity(0)
66 ,fGasNobleFraction(0)
67 ,fPrimaryIonisation(0)
68{
69 //
70 // Standard constructor for the TRD
71 //
72
73 // Check that FRAME is there otherwise we have no place where to put TRD
74 AliModule *frame = gAlice->GetModule("FRAME");
75 if (!frame) {
76 AliError("TRD needs FRAME to be present\n");
77 exit(1);
78 }
79
80 // Define the TRD geometry
81 if ((frame->IsVersion() == 0) ||
82 (frame->IsVersion() == 1)) {
83 fGeometry = new AliTRDgeometry();
84 }
85 else {
86 AliError("Could not find valid FRAME version\n");
87 exit(1);
88 }
89
90 // Allocate the hit array
91 fHits = new TClonesArray("AliTRDhit",405);
92 gAlice->GetMCApp()->AddHitList(fHits);
93
94}
95
96//_____________________________________________________________________________
97AliTRD::~AliTRD()
98{
99 //
100 // TRD destructor
101 //
102
103 if (fGeometry) {
104 delete fGeometry;
105 fGeometry = 0;
106 }
107
108 if (fHits) {
109 delete fHits;
110 fHits = 0;
111 }
112
113}
114
115//_____________________________________________________________________________
116void AliTRD::Hits2Digits()
117{
118 //
119 // Create digits
120 //
121
122 AliTRDdigitizer digitizer("TRDdigitizer","TRD digitizer class");
123 AliLog::SetClassDebugLevel("TRDdigitizer",AliDebugLevel());
124
125 // Initialization
126 digitizer.InitDetector();
127
128 if (!fLoader->TreeH()) {
129 fLoader->LoadHits("read");
130 }
131 fLoader->LoadDigits("recreate");
132
133 AliRunLoader *runLoader = fLoader->GetRunLoader();
134
135 for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
136 runLoader->GetEvent(iEvent);
137 digitizer.Open(runLoader,iEvent);
138 digitizer.MakeDigits();
139 digitizer.WriteDigits();
140 }
141
142 fLoader->UnloadHits();
143 fLoader->UnloadDigits();
144
145}
146
147//_____________________________________________________________________________
148void AliTRD::Hits2SDigits()
149{
150 //
151 // Create summable digits
152 //
153
154 AliTRDdigitizer digitizer("TRDdigitizer","TRD digitizer class");
155 // For the summable digits
156 digitizer.SetSDigits(kTRUE);
157 AliLog::SetClassDebugLevel("TRDdigitizer",AliDebugLevel());
158
159 // Initialization
160 digitizer.InitDetector();
161
162 if (!fLoader->TreeH()) {
163 fLoader->LoadHits("read");
164 }
165 fLoader->LoadSDigits("recreate");
166
167 AliRunLoader *runLoader = fLoader->GetRunLoader();
168
169 for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
170 runLoader->GetEvent(iEvent);
171 digitizer.Open(runLoader,iEvent);
172 digitizer.MakeDigits();
173 digitizer.WriteDigits();
174 }
175
176 fLoader->UnloadHits();
177 fLoader->UnloadSDigits();
178
179}
180
181//_____________________________________________________________________________
182AliDigitizer *AliTRD::CreateDigitizer(AliDigitizationInput* digInput) const
183{
184 //
185 // Creates a new digitizer object
186 //
187
188 return new AliTRDdigitizer(digInput);
189
190}
191
192//_____________________________________________________________________________
193void AliTRD::SDigits2Digits()
194{
195 //
196 // Create final digits from summable digits
197 //
198
199 // Create the TRD digitizer
200 AliTRDdigitizer digitizer("TRDdigitizer","TRD digitizer class");
201 AliLog::SetClassDebugLevel("TRDdigitizer",AliDebugLevel());
202
203 // Set the parameter
204 digitizer.SetEvent(gAlice->GetEvNumber());
205
206 // Initialization
207 digitizer.InitDetector();
208
209 // Read the s-digits via digits manager
210 AliTRDdigitsManager sdigitsManager;
211
212 AliLog::SetClassDebugLevel("TRDdigitisManager",AliDebugLevel());
213 sdigitsManager.SetSDigits(kTRUE);
214 sdigitsManager.CreateArrays();
215
216 if (!fLoader->TreeS()) {
217 if (fLoader->LoadSDigits("read")) {
218 return;
219 }
220 }
221 if (!fLoader->TreeS()) {
222 AliError(Form("Error while reading SDigits for event %d",gAlice->GetEvNumber()));
223 return;
224 }
225
226 sdigitsManager.ReadDigits(fLoader->TreeS());
227
228 // Add the s-digits to the input list
229 digitizer.AddSDigitsManager(&sdigitsManager);
230
231 // Convert the s-digits to normal digits
232 digitizer.SDigits2Digits();
233
234 // Store the digits
235 if (!fLoader->TreeD()) {
236 fLoader->MakeTree("D");
237 }
238 if (digitizer.MakeBranch(fLoader->TreeD())){
239 digitizer.WriteDigits();
240 }
241
242}
243
244//_____________________________________________________________________________
245void AliTRD::Digits2Raw()
246{
247 //
248 // Convert digits of the current event to raw data
249 //
250
251 fLoader->LoadDigits();
252 TTree *digits = fLoader->TreeD();
253 if (!digits) {
254 AliError("No digits tree");
255 return;
256 }
257
258 AliTRDrawData rawWriter;
259 if (!rawWriter.Digits2Raw(digits)) {
260 AliError("The raw writer could not load the digits tree");
261 }
262
263 fLoader->UnloadDigits();
264
265}
266
267//_____________________________________________________________________________
268void AliTRD::AddHit(Int_t track, Int_t det, Float_t *hits, Int_t q
269 , Float_t time, Bool_t inDrift)
270{
271 //
272 // Add a hit for the TRD
273 //
274
275 TClonesArray &lhits = *fHits;
276 AliTRDhit *hit = new(lhits[fNhits++]) AliTRDhit(fIshunt
277 ,track
278 ,det
279 ,hits
280 ,q
281 ,time);
282
283 if (inDrift) {
284 hit->SetDrift();
285 }
286 else {
287 hit->SetAmplification();
288 }
289
290 if (q < 0) {
291 hit->SetTRphoton();
292 }
293
294}
295
296//_____________________________________________________________________________
297void AliTRD::CreateGeometry()
298{
299 //
300 // Creates the volumes for the TRD chambers
301 //
302
303 // Check that FRAME is there otherwise we have no place where to put the TRD
304 AliModule *frame = gAlice->GetModule("FRAME");
305 if (!frame) {
306 AliFatal("The TRD needs the FRAME to be defined first");
307 }
308
309 fGeometry->CreateGeometry(fIdtmed->GetArray() - 1299);
310
311}
312
313//_____________________________________________________________________________
314void AliTRD::CreateMaterials()
315{
316 //
317 // Create the materials for the TRD
318 //
319
320 Int_t isxfld = ((AliMagF *) TGeoGlobalMagField::Instance()->GetField())->Integ();
321 Float_t sxmgmx = ((AliMagF *) TGeoGlobalMagField::Instance()->GetField())->Max();
322
323 //////////////////////////////////////////////////////////////////////////
324 // Define Materials
325 //////////////////////////////////////////////////////////////////////////
326
327 // Aluminum
328 AliMaterial( 1,"Al", 26.98, 13.0, 2.7, 8.9, 37.2);
329 // Copper
330 AliMaterial( 2,"Cu", 63.54, 29.0, 8.96, 1.43, 14.8);
331 // Carbon
332 AliMaterial( 3,"C" , 12.01, 6.0, 2.265, 18.8, 74.4);
333 // Carbon for fiber mats
334 AliMaterial( 4,"C2", 12.01, 6.0, 1.75, 18.8, 74.4);
335 // Zinc
336 AliMaterial( 5,"Sn", 118.71, 50.0, 7.31, 1.21, 14.8);
337 // Silicon
338 AliMaterial( 6,"Si", 28.09, 14.0, 2.33, 9.36, 37.2);
339 // Iron
340 AliMaterial( 7,"Fe", 55.85, 26.0, 7.87, 1.76, 14.8);
341
342 // Air
343 Float_t aAir[4] = { 12.011 , 14.0 , 15.9994 , 36.0 };
344 Float_t zAir[4] = { 6.0 , 7.0 , 8.0 , 18.0 };
345 Float_t wAir[4] = { 0.000124, 0.755267, 0.231781, 0.012827 };
346 Float_t dAir = 1.20479e-03;
347 AliMixture(51,"Air", aAir, zAir, dAir, 4, wAir );
348 // Polyethilene (CH2)
349 Float_t ape[2] = { 12.011 , 1.0079 };
350 Float_t zpe[2] = { 6.0 , 1.0 };
351 Float_t wpe[2] = { 1.0 , 2.0 };
352 Float_t dpe = 0.95;
353 AliMixture(52,"Polyethilene", ape, zpe, dpe, -2, wpe );
354 // Gas mixtures
355 // Xe/CO2-gas-mixture (85% / 15%)
356 Float_t aXeCO2[3] = { 131.29 , 12.0107 , 15.9994 };
357 Float_t zXeCO2[3] = { 54.0 , 6.0 , 8.0 };
358 Float_t wXeCO2[3] = { 8.5 , 1.5 , 3.0 };
359 Float_t fxc = 0.85;
360 Float_t dxe = 0.00549; // at 20C
361 Float_t dco = 0.00186; // at 20C
362 Float_t dgmXe = fxc * dxe + (1.0 - fxc) * dco;
363 // Ar/CO2-gas-mixture
364 Float_t aArCO2[3] = { 39.948 , 12.0107 , 15.9994 };
365 Float_t zArCO2[3] = { 18.0 , 6.0 , 8.0 };
366 Float_t wArCO2[3] = { 8.2 , 1.8 , 3.6 };
367 Float_t fac = 0.82;
368 Float_t dar = 0.00166; // at 20C
369 Float_t dgmAr = fac * dar + (1.0 - fac) * dco;
370 if (AliTRDCommonParam::Instance()->IsXenon()) {
371 AliMixture(53,"XeCO2", aXeCO2, zXeCO2, dgmXe, -3, wXeCO2);
372 }
373 else if (AliTRDCommonParam::Instance()->IsArgon()) {
374 AliInfo("Gas mixture: Ar C02 (80/20)");
375 AliMixture(53,"ArCO2", aArCO2, zArCO2, dgmAr, -3, wArCO2);
376 }
377 else {
378 AliFatal("Wrong gas mixture");
379 exit(1);
380 }
381 // G10
382 Float_t aG10[4] = { 1.0079, 12.011 , 15.9994, 28.086 };
383 Float_t zG10[4] = { 1.0 , 6.0 , 8.0 , 14.0 };
384 Float_t wG10[4] = { 0.023 , 0.194 , 0.443 , 0.340 };
385 Float_t dG10 = 2.0;
386 AliMixture(54,"G10", aG10, zG10, dG10, 4,wG10 );
387 // Water
388 Float_t awa[2] = { 1.0079, 15.9994 };
389 Float_t zwa[2] = { 1.0 , 8.0 };
390 Float_t wwa[2] = { 2.0 , 1.0 };
391 Float_t dwa = 1.0;
392 AliMixture(55,"Water", awa, zwa, dwa, -2,wwa );
393 // Rohacell (C5H8O2), X0 = 535.005cm
394 Float_t arh[3] = { 12.011 , 1.0079, 15.9994 };
395 Float_t zrh[3] = { 6.0 , 1.0 , 8.0 };
396 Float_t wrh[3] = { 5.0 , 8.0 , 2.0 };
397 Float_t drh = 0.075;
398 AliMixture(56,"Rohacell", arh, zrh, drh, -3,wrh );
399 // Epoxy (C18H19O3)
400 Float_t aEpoxy[3] = { 15.9994, 1.0079, 12.011 };
401 Float_t zEpoxy[3] = { 8.0 , 1.0 , 6.0 };
402 Float_t wEpoxy[3] = { 3.0 , 19.0 , 18.0 };
403 Float_t dEpoxy = 1.8 ;
404 AliMixture(57,"Epoxy", aEpoxy,zEpoxy,dEpoxy,-3,wEpoxy);
405 // Araldite, low density epoxy (C18H19O3)
406 Float_t aAral[3] = { 15.9994, 1.0079, 12.011 };
407 Float_t zAral[3] = { 8.0 , 1.0 , 6.0 };
408 Float_t wAral[3] = { 3.0 , 19.0 , 18.0 };
409 Float_t dAral = 1.12; // Hardener: 1.15, epoxy: 1.1, mixture: 1/2
410 AliMixture(58,"Araldite", aAral, zAral, dAral, -3,wAral );
411 // Mylar
412 Float_t aMy[3] = { 12.011 , 1.0 , 15.9994 };
413 Float_t zMy[3] = { 6.0 , 1.0 , 8.0 };
414 Float_t wMy[3] = { 5.0 , 4.0 , 2.0 };
415 Float_t dMy = 1.39;
416 AliMixture(59,"Mylar", aMy, zMy, dMy, -3,wMy );
417 // Polypropylene (C3H6) for radiator fibers
418 Float_t app[2] = { 12.011 , 1.0079 };
419 Float_t zpp[2] = { 6.0 , 1.0 };
420 Float_t wpp[2] = { 3.0 , 6.0 };
421 Float_t dpp = 0.068;
422 AliMixture(60,"Polypropylene",app, zpp, dpp, -2,wpp );
423 // Aramide for honeycomb
424 Float_t aAra[4] = { 1.0079, 12.011 , 15.9994, 14.0067 };
425 Float_t zAra[4] = { 1.0 , 6.0 , 8.0 , 7.0 };
426 Float_t wAra[4] = { 3.0 , 1.0 , 1.0 , 1.0 };
427 Float_t dAra = 0.032;
428 AliMixture(61,"Aramide", aAra, zAra, dAra, -4,wAra );
429 // GFK for Wacosit (Epoxy + Si)
430 Float_t aGFK[4] = { 1.0079, 12.011 , 15.9994, 28.086 };
431 Float_t zGFK[4] = { 1.0 , 6.0 , 8.0 , 14.0 };
432 Float_t wGFK[4] = { 0.0445, 0.5031, 0.1118, 0.340 };
433 Float_t dGFK = 2.0;
434 AliMixture(62,"GFK", aGFK, zGFK, dGFK, 4,wGFK );
435
436 //////////////////////////////////////////////////////////////////////////
437 // Tracking Media Parameters
438 //////////////////////////////////////////////////////////////////////////
439
440 // General tracking parameter
441 Float_t tmaxfd = -10.0;
442 Float_t stemax = -1.0e10;
443 Float_t deemax = -0.1;
444 Float_t epsil = 1.0e-4;
445 Float_t stmin = -0.001;
446
447 // Al Frame
448 AliMedium( 1,"Al Frame" , 1,0,isxfld,sxmgmx
449 ,tmaxfd,stemax,deemax,epsil,stmin);
450 // Air
451 AliMedium( 2,"Air" ,51,0,isxfld,sxmgmx
452 ,tmaxfd,stemax,deemax,epsil,stmin);
453 // Wires
454 AliMedium( 3,"Wires" , 2,0,isxfld,sxmgmx
455 ,tmaxfd,stemax,deemax,epsil,stmin);
456 // All other ROB materials (caps, etc.)
457 AliMedium( 4,"ROB Other" , 2,0,isxfld,sxmgmx
458 ,tmaxfd,stemax,deemax,epsil,stmin);
459 // Cu pads
460 AliMedium( 5,"Padplane" , 2,1,isxfld,sxmgmx
461 ,tmaxfd,stemax,deemax,epsil,stmin);
462 // Fee + cables
463 AliMedium( 6,"Readout" , 2,0,isxfld,sxmgmx
464 ,tmaxfd,stemax,deemax,epsil,stmin);
465 // C frame (Wacosit)
466 AliMedium( 7,"Wacosit" ,62,0,isxfld,sxmgmx
467 ,tmaxfd,stemax,deemax,epsil,stmin);
468 // INOX of cooling bus bars
469 AliMedium( 8,"Cooling bus", 7,0,isxfld,sxmgmx
470 ,tmaxfd,stemax,deemax,epsil,stmin);
471 // Gas-mixture (Xe/CO2)
472 AliMedium( 9,"Gas-mix" ,53,1,isxfld,sxmgmx
473 ,tmaxfd,stemax,deemax,epsil,stmin);
474 // Honeycomb
475 AliMedium(10,"Honeycomb" ,61,0,isxfld,sxmgmx
476 ,tmaxfd,stemax,deemax,epsil,stmin);
477 // Araldite glue
478 AliMedium(11,"Glue" ,58,0,isxfld,sxmgmx
479 ,tmaxfd,stemax,deemax,epsil,stmin);
480 // G10-plates
481 AliMedium(13,"G10-plates" ,54,0,isxfld,sxmgmx
482 ,tmaxfd,stemax,deemax,epsil,stmin);
483 // Cooling water
484 AliMedium(14,"Water" ,55,0,isxfld,sxmgmx
485 ,tmaxfd,stemax,deemax,epsil,stmin);
486 // Rohacell for the radiator
487 AliMedium(15,"Rohacell" ,56,0,isxfld,sxmgmx
488 ,tmaxfd,stemax,deemax,epsil,stmin);
489 // Al layer in MCMs
490 AliMedium(16,"MCM-Al" , 1,0,isxfld,sxmgmx
491 ,tmaxfd,stemax,deemax,epsil,stmin);
492 // Sn layer in MCMs
493 AliMedium(17,"MCM-Sn" , 5,0,isxfld,sxmgmx
494 ,tmaxfd,stemax,deemax,epsil,stmin);
495 // Cu layer in MCMs
496 AliMedium(18,"MCM-Cu" , 2,0,isxfld,sxmgmx
497 ,tmaxfd,stemax,deemax,epsil,stmin);
498 // G10 layer in MCMs
499 AliMedium(19,"MCM-G10" ,54,0,isxfld,sxmgmx
500 ,tmaxfd,stemax,deemax,epsil,stmin);
501 // Si in readout chips
502 AliMedium(20,"Chip-Si" , 6,0,isxfld,sxmgmx
503 ,tmaxfd,stemax,deemax,epsil,stmin);
504 // Epoxy in readout chips
505 AliMedium(21,"Chip-Ep" ,57,0,isxfld,sxmgmx
506 ,tmaxfd,stemax,deemax,epsil,stmin);
507 // PE in connectors
508 AliMedium(22,"Conn-PE" ,52,0,isxfld,sxmgmx
509 ,tmaxfd,stemax,deemax,epsil,stmin);
510 // Cu in connectors
511 AliMedium(23,"Chip-Cu" , 2,0,isxfld,sxmgmx
512 ,tmaxfd,stemax,deemax,epsil,stmin);
513 // Al of cooling pipes
514 AliMedium(24,"Cooling" , 1,0,isxfld,sxmgmx
515 ,tmaxfd,stemax,deemax,epsil,stmin);
516 // Cu in services
517 AliMedium(25,"Serv-Cu" , 2,0,isxfld,sxmgmx
518 ,tmaxfd,stemax,deemax,epsil,stmin);
519 // Carbon fiber mat
520 AliMedium(26,"Carbon" , 4,0,isxfld,sxmgmx
521 ,tmaxfd,stemax,deemax,epsil,stmin);
522 // Mylar foil
523 AliMedium(27,"Mylar" ,59,0,isxfld,sxmgmx
524 ,tmaxfd,stemax,deemax,epsil,stmin);
525 // Polypropylene fibers
526 AliMedium(28,"Fiber" ,60,0,isxfld,sxmgmx
527 ,tmaxfd,stemax,deemax,epsil,stmin);
528
529 // Save the density values for the TRD absorbtion
530 Float_t dmy = 1.39;
531 fFoilDensity = dmy;
532 if (AliTRDCommonParam::Instance()->IsXenon()) {
533 fGasDensity = dgmXe;
534 fGasNobleFraction = fxc;
535 }
536 else if (AliTRDCommonParam::Instance()->IsArgon()) {
537 fGasDensity = dgmAr;
538 fGasNobleFraction = fac;
539 }
540
541}
542
543//_____________________________________________________________________________
544void AliTRD::Init()
545{
546 //
547 // Initialize the TRD detector after the geometry has been created
548 //
549
550 AliDebug(1,"++++++++++++++++++++++++++++++++++++++++++++++");
551
552 if (fGeometry->IsVersion() != 1) {
553 AliError("Not a valid geometry");
554 }
555
556}
557
558//_____________________________________________________________________________
559void AliTRD::ResetDigits()
560{
561 //
562 // Reset number of digits and the digits array for this detector
563 //
564
565 fNdigits = 0;
566
567 if (fDigits) {
568 fDigits->Clear();
569 }
570
571}
572
573//_____________________________________________________________________________
574void AliTRD::SetTreeAddress()
575{
576 //
577 // Set the branch addresses for the trees.
578 //
579
580 if (fLoader->TreeH() &&
581 (fHits == 0x0)) {
582 fHits = new TClonesArray("AliTRDhit",405);
583 }
584 AliDetector::SetTreeAddress();
585
586}
587
588//_____________________________________________________________________________
589Bool_t AliTRD::Raw2SDigits(AliRawReader *rawReader)
590{
591 //
592 // Converts RAW data to SDigits
593 //
594
595 AliLoader *loader = fRunLoader->GetLoader("TRDLoader");
596 if (!loader) {
597 AliError("Can not get TRD loader from Run Loader");
598 return kFALSE;
599 }
600
601 TTree *tree = 0;
602 tree = loader->TreeS();
603 if (!tree) {
604 loader->MakeTree("S");
605 tree = loader->TreeS();
606 }
607
608 AliTRDdigitizer digitizer("TRDdigitizer","TRD digitizer class");
609
610 AliTRDrawData *rawdata = new AliTRDrawData();
611 AliTRDdigitsManager *digitsManager = rawdata->Raw2Digits(rawReader);
612
613 // Create the s-digits manager
614 AliTRDdigitsManager *sdigitsManager = new AliTRDdigitsManager();
615
616 if (sdigitsManager) {
617
618 sdigitsManager->SetSDigits(kTRUE);
619 sdigitsManager->CreateArrays();
620
621 // Convert the digits into s-digits
622 digitizer.Digits2SDigits(digitsManager,sdigitsManager);
623
624 sdigitsManager->MakeBranch(tree);
625 sdigitsManager->WriteDigits();
626
627 delete digitsManager;
628
629 return kTRUE;
630
631 }
632 else {
633
634 return kFALSE;
635
636 }
637
638}
639
640//_____________________________________________________________________________
641AliLoader *AliTRD::MakeLoader(const Char_t *topfoldername)
642{
643 //
644 // Create a loader for the TRD tracklets
645 //
646
647 fLoader = new AliLoader(GetName(),topfoldername);
648
649 AliInfo("Adding Tracklets-loader");
650
651 AliDataLoader *dl = new AliDataLoader("TRD.Tracklets.root"
652 ,"tracklets"
653 ,"tracklets");
654 fLoader->AddDataLoader(dl);
655
656 dl = new AliDataLoader("TRD.GtuTracks.root"
657 ,"gtutracks"
658 ,"gtutracks");
659 fLoader->AddDataLoader(dl);
660
661 return fLoader;
662
663}
664
665//_____________________________________________________________________________
666AliTRD &AliTRD::operator=(const AliTRD &trd)
667{
668 //
669 // Assignment operator
670 //
671
672 if (this != &trd) {
673 ((AliTRD &) trd).Copy(*this);
674 }
675
676 return *this;
677
678}