#include "AliTRDtrackletMCM.h"
#include "AliESDEvent.h"
+using std::ifstream;
ClassImp(AliTRDgtuSim)
AliTRDgtuSim::AliTRDgtuSim(AliRunLoader *rl)
return kTRUE;
}
-Bool_t AliTRDgtuSim::RunGTU(AliLoader *loader, AliESDEvent *esd)
+Bool_t AliTRDgtuSim::RunGTU(AliLoader *loader, AliESDEvent *esd, Int_t label)
{
// run the GTU on tracklets taken from the loader
// if specified the GTU tracks are written to the ESD event
}
}
else {
- LoadTracklets(esd);
+ LoadTracklets(esd, label);
}
AliDebug(1, Form("running on %i tracklets", fTrackletArray->GetEntriesFast()));
return kTRUE;
}
-Bool_t AliTRDgtuSim::LoadTracklets(const AliESDEvent *const esd)
+Bool_t AliTRDgtuSim::LoadTracklets(const AliESDEvent *const esd, Int_t label)
{
AliDebug(1,"Loading tracklets from ESD event ...");
for (Int_t iTracklet = 0; iTracklet < esd->GetNumberOfTrdTracklets(); iTracklet++) {
AliESDTrdTracklet *trkl = esd->GetTrdTracklet(iTracklet);
- new ((*fTrackletArray)[fTrackletArray->GetEntries()]) AliTRDtrackletGTU(trkl);
+ if (label < -1) {
+ if (trkl->GetLabel() == label)
+ new ((*fTrackletArray)[fTrackletArray->GetEntries()]) AliTRDtrackletGTU(trkl);
+ }
+ else
+ if (trkl->GetLabel() >= -1)
+ new ((*fTrackletArray)[fTrackletArray->GetEntries()]) AliTRDtrackletGTU(trkl);
}
return kTRUE;
TIter next(listOfTracks);
while (AliTRDtrackGTU *trk = (AliTRDtrackGTU*) next()) {
AliESDTrdTrack *trdtrack = trk->CreateTrdTrack();
- if (trdtrack->GetLabel() < 0)
- trdtrack->SetLabel(-2);
esd->AddTrdTrack(trdtrack);
delete trdtrack;
}