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 ///////////////////////////////////////////////////////
19 // TRD main trigger class for L1 //
22 ///////////////////////////////////////////////////////
27 #include "AliTriggerInput.h"
29 #include "AliTRDTriggerL1.h"
30 #include "AliTRDtrigParam.h"
31 #include "AliTRDtrigger.h"
32 #include "AliTRDgtuTrack.h"
33 #include "AliTRDgeometry.h"
35 ClassImp(AliTRDTriggerL1)
37 //_____________________________________________________________________________
38 AliTRDTriggerL1::AliTRDTriggerL1():AliTriggerDetector()
45 //_____________________________________________________________________________
46 void AliTRDTriggerL1::CreateInputs()
49 fInputs.AddLast(new AliTriggerInput( "TRD_HadrLPt_L1", "Single hadron low pt ", 0x01 ));
50 fInputs.AddLast(new AliTriggerInput( "TRD_HadrHPt_L1", "Single hadron high pt", 0x02 ));
51 fInputs.AddLast(new AliTriggerInput( "TRD_Unlike_EPair_L1", "Unlike electron pair", 0x04 ));
52 fInputs.AddLast(new AliTriggerInput( "TRD_Unlike_EPair_HPt_L1", "Unlike electron pair high pt", 0x08 ));
53 fInputs.AddLast(new AliTriggerInput( "TRD_Like_EPair_L1", "Like electron pair", 0x10 ));
54 fInputs.AddLast(new AliTriggerInput( "TRD_Like_EPair_HPt_L1", "Like electron pair high pt", 0x20 ));
55 fInputs.AddLast(new AliTriggerInput( "TRD_Electron_L1", "Single electron", 0x40 ));
56 fInputs.AddLast(new AliTriggerInput( "TRD_Electron_HPt_L1", "Single electron high pt", 0x80 ));
60 //_____________________________________________________________________________
61 void AliTRDTriggerL1::Trigger()
64 AliRunLoader* runLoader = gAlice->GetRunLoader();
66 AliLoader *loader=runLoader->GetLoader("TRDLoader");
68 Int_t nEvents = runLoader->GetNumberOfEvents();
70 // Trigger (tracklets, LTU)
72 AliTRDgeometry *geo = AliTRDgeometry::GetGeometry(runLoader);
74 AliTRDtrigger trdTrigger("Trigger","Trigger class");
76 AliTRDtrigParam *trigp = new AliTRDtrigParam("TRDtrigParam","TRD Trigger parameters");
78 gAlice = runLoader->GetAliRun();
79 Double_t x[3] = { 0.0, 0.0, 0.0 };
81 gAlice->Field(x,b); // b[] is in kilo Gauss
82 Float_t field = b[2] * 0.1; // Tesla
83 Info("Trigger","Trigger set for magnetic field = %f Tesla \n",field);
85 trigp->SetField(field);
88 trdTrigger.SetParameter(trigp);
89 trdTrigger.SetRunLoader(runLoader);
92 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
94 trdTrigger.Open(runLoader->GetFileName(), iEvent);
95 trdTrigger.ReadDigits();
96 trdTrigger.MakeTracklets(kTRUE);
97 trdTrigger.WriteTracklets(-1);
101 // Trigger (tracks, GTU)
103 Float_t highPt = trigp->GetHighPt();
109 const Int_t maxEle = 1000;
111 Int_t electronPlus, electronMinus;
112 Int_t sectorElePlus[maxEle], sectorEleMinus[maxEle];
113 Float_t ptElePlus[maxEle], ptEleMinus[maxEle];
114 Int_t hadronLowPt, hadronHighPt;
119 AliTRDgtuTrack *gtuTrack;
120 Int_t nTracks = trdTrigger.GetNumberOfTracks();
121 for (Int_t iTrack = 0; iTrack < nTracks; iTrack++) {
123 gtuTrack = trdTrigger.GetTrack(iTrack);
125 pid = gtuTrack->GetPID();
126 isElectron = gtuTrack->IsElectron();
127 pt = gtuTrack->GetPt();
128 det = gtuTrack->GetDetector();
130 sec = geo->GetSector(det);
135 sectorEleMinus[electronMinus] = sec;
136 ptEleMinus[electronMinus] = pt;
139 sectorElePlus[electronPlus] = sec;
140 ptElePlus[electronPlus] = pt;
146 if (TMath::Abs(pt) < highPt) {
156 loader->UnloadTracks();
160 if (hadronLowPt) SetInput("TRD_Hadr_LPt_L1");
161 if (hadronHighPt) SetInput("TRD_Hadr_HPt_L1");
163 // electron-positron pairs (open angle > 80 deg)
165 Int_t secPlus, secMinus, secDiff;
166 Bool_t electronUnlikePair = kFALSE;
167 Bool_t electronUnlikePairHPt = kFALSE;
169 if (electronMinus > 0 && electronPlus > 0) {
170 for (Int_t iPlus = 0; iPlus < electronPlus; iPlus++) {
171 secPlus = sectorElePlus[iPlus];
172 for (Int_t iMinus = 0; iMinus < electronMinus; iMinus++) {
173 secMinus = sectorEleMinus[iMinus];
174 secDiff = TMath::Abs(secPlus-secMinus);
175 if (secDiff > 9) secDiff = 18 - secDiff;
177 electronUnlikePair = kTRUE;
178 if (TMath::Abs(ptElePlus[iPlus]) > highPt && TMath::Abs(ptEleMinus[iMinus]) > highPt) {
179 electronUnlikePairHPt = kTRUE;
186 if (electronUnlikePair) SetInput("TRD_Unlike_EPair_L1");
187 if (electronUnlikePairHPt) SetInput("TRD_Unlike_EPair_HPt_L1");
189 // like electron/positron pairs
191 Bool_t ele1, ele1HPt;
192 Bool_t ele2, ele2HPt;
200 if (electronPlus > 1) {
201 for (Int_t iPlus = 0; iPlus < electronPlus; iPlus++) {
207 if (TMath::Abs(ptElePlus[iPlus]) > highPt) {
210 } else if (!ele2HPt) {
217 if (ele1 && ele2 ) SetInput("TRD_Like_EPair_L1");
218 if (ele1HPt && ele2HPt) SetInput("TRD_Like_EPair_HPt_L1");
226 if (electronMinus > 1) {
227 for (Int_t iMinus = 0; iMinus < electronMinus; iMinus++) {
233 if (TMath::Abs(ptEleMinus[iMinus]) > highPt) {
236 } else if (!ele2HPt) {
243 if (ele1 && ele2 ) SetInput("TRD_Like_EPair_L1");
244 if (ele1HPt && ele2HPt) SetInput("TRD_Like_EPair_HPt_L1");
246 // single electron/positron
248 if (electronPlus > 0 || electronMinus > 0) {
249 SetInput("TRD_Electron_L1");
250 for (Int_t iPlus = 0; iPlus < electronPlus; iPlus++) {
251 if (TMath::Abs(ptElePlus[iPlus]) > highPt) SetInput("TRD_Electron_HPt_L1");
254 for (Int_t iMinus = 0; iMinus < electronMinus; iMinus++) {
255 if (TMath::Abs(ptEleMinus[iMinus]) > highPt) SetInput("TRD_Electron_HPt_L1");