- // Trigger (tracklets, LTU)
-
- loader->LoadTracks();
- if (loader->TreeT()) {
- Info("Reconstruct","Tracklets already exist");
- return;
- }
- Info("Reconstruct","Trigger tracklets will be produced");
-
- AliTRDtrigger trdTrigger("Trigger","Trigger class");
-
- AliTRDtrigParam *trigp = new AliTRDtrigParam("TRDtrigParam","TRD Trigger parameters");
-
- if (runLoader->GetAliRun() == 0x0) runLoader->LoadgAlice();
- gAlice = runLoader->GetAliRun();
- Double_t x[3] = { 0.0, 0.0, 0.0 };
- Double_t b[3];
- gAlice->Field(x,b); // b[] is in kilo Gauss
- Float_t field = b[2] * 0.1; // Tesla
- Info("Reconstruct","Trigger set for magnetic field = %f Tesla \n",field);
-
- trigp->SetField(field);
- trigp->Init();
- trdTrigger.SetParameter(trigp);
-
- rawReader->RewindEvents();
-
- for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
- if (!rawReader->NextEvent()) break;
- trdTrigger.Open(runLoader->GetFileName(), iEvent);
- trdTrigger.ReadDigits(rawReader);
- trdTrigger.MakeTracklets();
- trdTrigger.WriteTracklets(-1);
- }
-
- loader->UnloadTracks();
-