#include "AliTRDhit.h"
#include "AliTRDcluster.h"
#include "AliTRDseedV1.h"
+#include "AliTRDtrackletMCM.h"
+#include "AliTRDmcmSim.h"
#include "AliTRDtrackV1.h"
#include "AliTRDtrackerV1.h"
#include "AliTRDpadPlane.h"
#include "AliTRDdigitsManager.h"
-#include "AliTRDdataArrayDigits.h"
#include "AliTRDarrayADC.h"
#include "AliTRDSignalIndex.h"
#include "AliTRDgeometry.h"
for (Int_t ic = 0; ic < ncols; ic++) {
dy = pp->GetColSize(ic);
for (Int_t it = 0; it < ntbs; it++) {
- q = fData.GetDataUnchecked(ir, ic, it);
+ q = fData.GetData(ir, ic, it);
if (q < threshold) continue;
Double_t x[6] = {0., 0., Double_t(q), 0., 0., 0.};
- Int_t roc[3] = {ir, ic, 0};
- Bool_t out = kTRUE;
- transform.Transform(&x[0], &roc[0], UInt_t(it), out, 0);
+ //Int_t roc[3] = {ir, ic, 0};
+ //Bool_t out = kTRUE;
+ //transform.Transform(&x[0], &roc[0], UInt_t(it), out, 0);
scale = q < 512 ? q/512. : 1.;
color = 50+int(scale*50.);
} // end time loop
} // end col loop
} // end row loop
- fData.Compress(1);
+ fData.Compress();
// rotate to global coordinates
//RefitPlex();
}
adc = data->GetData(row, col, time);
if(adc <= 1) continue;
- fData.SetDataUnchecked(row, col, time, adc);
+ fData.SetData(row, col, time, adc);
//fIndex->AddIndexTBin(row,col,time);
//printf("\tr[%d] c[%d] t[%d] ADC[%d]\n", row, col, time, adc);
}
}
- fData.Compress(1);
+ fData.Compress();
}
Float_t tilt = trklt->GetTilt();
Float_t g[3];
AliTRDcluster *c = 0x0;
- for(Int_t ic=0; ic<AliTRDseed::knTimebins; ic++){
+ for(Int_t ic=0; ic<AliTRDseedV1::kNclusters; ic++){
if(!(c = trklt->GetClusters(ic))) continue;
if(!fClusters) AddElement(fClusters = new AliEveTRDClusters());
dx = x0 - c->GetX();
Int_t sec = AliTRDgeometry::GetSector(trklt->GetDetector());
Double_t alpha = AliTRDgeometry::GetAlpha() * (sec<9 ? sec + .5 : sec - 17.5);
- trklt->Fit(kTRUE);
+ //trklt->Fit(kTRUE);
y0 = trklt->GetYfit(0);
dydx = trklt->GetYfit(1);
Double_t xg = x0 * TMath::Cos(alpha) - y0 * TMath::Sin(alpha);
//SetPoint(1, x0, y0, z0);
}
+AliEveTRDTracklet::AliEveTRDTracklet(AliTRDtrackletMCM *tracklet) :
+ TEveLine(),
+ fClusters(0x0)
+{
+ SetName("TRD tracklet");
+ AliTRDtrackletMCM *trkl = new AliTRDtrackletMCM(*tracklet);
+ SetUserData(trkl);
+
+ SetTitle(Form("Det: %i, ROB: %i, MCM: %i, Label: %i\n0x%08x",
+ trkl->GetDetector(), trkl->GetROB(), trkl->GetMCM(), trkl->GetLabel(),
+ trkl->GetTrackletWord()));
+ SetLineColor(kRed);
+
+ AliTRDgeometry *geo = new AliTRDgeometry();
+
+ Float_t length = 10;
+ Double_t x[3];
+ Double_t p[3];
+ x[0] = trkl->GetX();
+ x[1] = trkl->GetY();
+ x[2] = trkl->GetZ();
+
+ geo->RotateBack(trkl->GetDetector(), x, p);
+ SetPoint(0, p[0], p[1], p[2]);
+
+ x[0] -= length;
+ x[1] += length * trkl->GetdYdX();
+ x[2] *= x[0] / (x[0] + length);
+ geo->RotateBack(trkl->GetDetector(), x, p);
+ SetPoint(1, p[0], p[1], p[2]);
+}
+
+AliEveTRDTracklet::~AliEveTRDTracklet()
+{
+ AliTRDtrackletMCM *trkl = dynamic_cast<AliTRDtrackletMCM*> ((AliTRDtrackletBase*) GetUserData());
+ delete trkl;
+}
+
+void AliEveTRDTracklet::ShowMCM(Option_t *opt) const
+{
+ AliTRDtrackletMCM *trkl = dynamic_cast<AliTRDtrackletMCM*> ((AliTRDtrackletBase*) GetUserData());
+ if (!trkl)
+ return;
+ printf("Det: %3i, ROB: %i, MCM: %2i\n", trkl->GetDetector(), trkl->GetROB(), trkl->GetMCM());
+ AliTRDmcmSim *mcm = new AliTRDmcmSim();
+ AliRunLoader *rl = AliEveEventManager::AssertRunLoader();
+ mcm->LoadMCM(rl, trkl->GetDetector(), trkl->GetROB(), trkl->GetMCM());
+ mcm->Tracklet();
+ mcm->Draw(opt);
+}
+
//______________________________________________________________________________
void AliEveTRDTracklet::Print(Option_t *o) const
{
AliTRDtrackerV1::SetNTimeBins(24);
-
fRim = new AliRieman(trk->GetNumberOfClusters());
AliTRDseedV1 *tracklet = 0x0;
for(Int_t il=0; il<AliTRDgeometry::kNlayer; il++){
AddElement(new AliEveTRDTracklet(tracklet));
AliTRDcluster *c = 0x0;
- tracklet->ResetClusterIter(kFALSE);
- while((c = tracklet->PrevCluster())){
+// tracklet->ResetClusterIter(kFALSE);
+// while((c = tracklet->PrevCluster())){
+ for(Int_t ic=AliTRDseedV1::kNtb; ic--;){
+ if(!(c=tracklet->GetClusters(ic))) continue;
Float_t xc = c->GetX();
Float_t yc = c->GetY();
Float_t zc = c->GetZ();
Float_t zt = tracklet->GetZref(0) - (tracklet->GetX0()-xc)*tracklet->GetZref(1);
yc -= tracklet->GetTilt()*(zc-zt);
- fRim->AddPoint(xc, yc, zc, 1, 10);
+ fRim->AddPoint(xc, yc, zc, .05, 2.3);
}
}
- fRim->Update();
+ if(trk->GetNumberOfTracklets()>1) fRim->Update();
SetStatus(fTrackState);
}
if(!fPoints){
fPoints = new AliTrackPoint[nc];
- AliTRDcluster *c = trk->GetCluster(0);
- Double_t x = c->GetX();
- Int_t sec = c->GetDetector()/30;
- fAlpha = AliTRDgeometry::GetAlpha() * (sec<9 ? sec + .5 : sec - 17.5);
+ // define the radial span of the track in the TRD
+ Double_t xmin = -1., xmax = -1.;
+ Int_t det = 0;
+ AliTRDseedV1 *trklt = 0x0;
+ for(Int_t ily=AliTRDgeometry::kNlayer; ily--;){
+ if(!(trklt = trk->GetTracklet(ily))) continue;
+ if(xmin<0.) xmin = trklt->GetX0() - AliTRDgeometry::CamHght() - AliTRDgeometry::CdrHght();
+ if(xmax<0.) xmax = trklt->GetX0();
+ det = trklt->GetDetector();
+ }
+ Int_t sec = det/AliTRDgeometry::kNdets;
+ fAlpha = AliTRDgeometry::GetAlpha() * (sec<9 ? sec + .5 : sec - 17.5); //trk->GetAlpha()
- Double_t dx = (trk->GetCluster(trk->GetNumberOfClusters()-1)->GetX()-x)/nc;
+ Double_t dx =(xmax - xmin)/nc;
for(Int_t ip=0; ip<nc; ip++){
- fPoints[ip].SetXYZ(x, 0., 0.);
- x+=dx;
+ fPoints[ip].SetXYZ(xmin, 0., 0.);
+ xmin+=dx;
}
BUILD = kTRUE;
}
} else {
//printf("Rieman track\n");
// if(trk->GetNumberOfTracklets() >=4) AliTRDtrackerV1::FitRiemanTilt(trk, 0x0, kTRUE, nc, fPoints);
- Float_t x = 0.;
+
+/* Float_t x = 0.;
for(Int_t ip = nc; ip--;){
x = fPoints[ip].GetX();
fPoints[ip].SetXYZ(x, fRim->GetYat(x), fRim->GetZat(x));
- }
+ }*/
}
}
Float_t global[3];
- for(Int_t ip=0; ip<nc; ip++){
+ for(Int_t ip=0; ip<0/*nc*/; ip++){
fPoints[ip].Rotate(-fAlpha).GetXYZ(global);
SetPoint(ip, global[0], global[1], global[2]);
}
SetLineColor(kMagenta);
}
} else {
- if(TESTBIT(s, kPID) == AliTRDReconstructor::kLQPID){
+ if(TESTBIT(s, kPID) == AliTRDpidUtil::kLQ){
//printf("PID color kLQPID\n");
//trk->GetReconstructor()->SetOption("!nn");
} else {