+#include <TClonesArray.h>
#include "AliITSClusterFinderV2SDD.h"
#include "AliITSRecPoint.h"
#include "AliITSDetTypeRec.h"
#include "AliITSCalibrationSDD.h"
#include "AliITSDetTypeRec.h"
#include "AliITSsegmentationSDD.h"
-#include <TClonesArray.h>
#include "AliITSdigitSDD.h"
#include "AliITSgeomTGeo.h"
ClassImp(AliITSClusterFinderV2SDD)
-AliITSClusterFinderV2SDD::AliITSClusterFinderV2SDD(AliITSDetTypeRec* dettyp):AliITSClusterFinderV2(dettyp),
-fNySDD(256),
-fNzSDD(256),
-fZpitchSDD(0.02940),
-fHwSDD(3.5085),
-fHlSDD(3.7632),
-fTimeOffsetSDD(0.){
+AliITSClusterFinderV2SDD::AliITSClusterFinderV2SDD(AliITSDetTypeRec* dettyp):AliITSClusterFinderV2(dettyp)
+{
//Default constructor
- SetTimeOffset();
}
//------------------------------------------------------------
// Actual SDD cluster finder
//------------------------------------------------------------
- Int_t kNzBins = fNzSDD + 2;
- const Int_t kMAXBIN=kNzBins*(fNySDD+2);
+ Int_t nAnodes = GetSeg()->NpzHalf();
+ Int_t nzBins = nAnodes+2;
+ Int_t nTimeBins = GetSeg()->Npx();
+ Int_t nxBins = nTimeBins+2;
+ const Int_t kMaxBin=nzBins*nxBins;
AliBin *bins[2];
- bins[0]=new AliBin[kMAXBIN];
- bins[1]=new AliBin[kMAXBIN];
+ bins[0]=new AliBin[kMaxBin];
+ bins[1]=new AliBin[kMaxBin];
AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(fModule);
+ if(cal==0){
+ AliError(Form("Calibration object not present for SDD module %d\n",fModule));
+ return;
+ }
AliITSresponseSDD* res = (AliITSresponseSDD*)cal->GetResponse();
const char *option=res->ZeroSuppOption();
AliITSdigitSDD *d=0;
Int_t i, ndigits=digits->GetEntriesFast();
for (i=0; i<ndigits; i++) {
d=(AliITSdigitSDD*)digits->UncheckedAt(i);
- Int_t y=d->GetCoord2()+1; //y
- Int_t z=d->GetCoord1()+1; //z
- Int_t q=d->GetSignal();
- if(!((strstr(option,"1D")) || (strstr(option,"2D")))){
- Float_t baseline = cal->GetBaseline(d->GetCoord1());
- if(q>baseline) q-=(Int_t)baseline;
- else q=0;
- }
- if(q<cal->GetThresholdAnode(d->GetCoord1())) continue;
- //if (q<3) continue;
-
- if (z <= fNzSDD) {
- bins[0][y*kNzBins+z].SetQ(q);
- bins[0][y*kNzBins+z].SetMask(1);
- bins[0][y*kNzBins+z].SetIndex(i);
- } else {
- z-=fNzSDD;
- bins[1][y*kNzBins+z].SetQ(q);
- bins[1][y*kNzBins+z].SetMask(1);
- bins[1][y*kNzBins+z].SetIndex(i);
+ Int_t ian=d->GetCoord1();
+ Int_t itb=d->GetCoord2();
+ Int_t iSide=0;
+ if (ian >= nAnodes) iSide=1;
+ Float_t gain=cal->GetChannelGain(ian);
+ Float_t charge=d->GetSignal();
+ if(strstr(option,"ZS")) charge+=(Float_t)cal->GetZSLowThreshold(iSide);
+ Float_t baseline = cal->GetBaseline(ian);
+ if(charge>baseline) charge-=baseline;
+ else charge=0;
+
+ if(gain>0.){ // Bad channels have gain=0.
+ charge/=gain;
+ if(charge<cal->GetThresholdAnode(ian)) continue;
+ Int_t q=(Int_t)(charge+0.5);
+ Int_t y=itb+1;
+ Int_t z=ian+1;
+ if (z <= nAnodes){
+ bins[0][y*nzBins+z].SetQ(q);
+ bins[0][y*nzBins+z].SetMask(1);
+ bins[0][y*nzBins+z].SetIndex(i);
+ } else {
+ z-=nAnodes;
+ bins[1][y*nzBins+z].SetQ(q);
+ bins[1][y*nzBins+z].SetMask(1);
+ bins[1][y*nzBins+z].SetIndex(i);
+ }
}
}
- FindClustersSDD(bins, kMAXBIN, kNzBins, digits);
+ FindClustersSDD(bins, kMaxBin, nzBins, digits);
delete[] bins[0];
delete[] bins[1];
const TGeoHMatrix *mT2L=AliITSgeomTGeo::GetTracking2LocalMatrix(fModule);
AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(fModule);
+ if(cal==0){
+ AliError(Form("Calibration object not present for SDD module %d\n",fModule));
+ return;
+ }
Int_t ncl=0;
TClonesArray &cl=*clusters;
for (Int_t s=0; s<2; s++)
for (k=0; k<npeaks-1; k++){//mark adjacent peaks
if (idx[k] < 0) continue; //this peak is already removed
for (l=k+1; l<npeaks; l++) {
- if (idx[l] < 0) continue; //this peak is already removed
- Int_t ki=idx[k]/nzBins, kj=idx[k] - ki*nzBins;
- Int_t li=idx[l]/nzBins, lj=idx[l] - li*nzBins;
- Int_t di=TMath::Abs(ki - li);
- Int_t dj=TMath::Abs(kj - lj);
- if (di>1 || dj>1) continue;
- if (bins[s][idx[k]].GetQ() > bins[s][idx[l]].GetQ()) {
- msk[l]=msk[k];
- idx[l]*=-1;
- } else {
- msk[k]=msk[l];
- idx[k]*=-1;
- break;
- }
+ if (idx[l] < 0) continue; //this peak is already removed
+ Int_t ki=idx[k]/nzBins, kj=idx[k] - ki*nzBins;
+ Int_t li=idx[l]/nzBins, lj=idx[l] - li*nzBins;
+ Int_t di=TMath::Abs(ki - li);
+ Int_t dj=TMath::Abs(kj - lj);
+ if (di>1 || dj>1) continue;
+ if (bins[s][idx[k]].GetQ() > bins[s][idx[l]].GetQ()) {
+ msk[l]=msk[k];
+ idx[l]*=-1;
+ } else {
+ msk[k]=msk[l];
+ idx[k]*=-1;
+ break;
+ }
}
}
}
for (k=0; k<npeaks; k++) {
- if (idx[k] < 0) continue; //removed peak
- AliITSRecPoint c;
- MakeCluster(idx[k], nzBins, bins[s], msk[k], c);
- //mi change
- Int_t milab[10];
- for (Int_t ilab=0;ilab<10;ilab++){
- milab[ilab]=-2;
- }
- Int_t maxi=0,mini=0,maxj=0,minj=0;
- //AliBin *bmax=&bins[s][idx[k]];
- //Float_t max = TMath::Max(TMath::Abs(bmax->GetQ())/5.,3.);
-
- for (Int_t di=-2; di<=2;di++){
- for (Int_t dj=-3;dj<=3;dj++){
- Int_t index = idx[k]+di+dj*nzBins;
- if (index<0) continue;
- if (index>=nMaxBin) continue;
- AliBin *b=&bins[s][index];
- Int_t nAnode=index%nzBins-1;
- Int_t adcSignal=b->GetQ();
- if(adcSignal>cal->GetThresholdAnode(nAnode)){
- if (di>maxi) maxi=di;
- if (di<mini) mini=di;
- if (dj>maxj) maxj=dj;
- if (dj<minj) minj=dj;
- }
- //
- if(digits) {
- if (TMath::Abs(di)<2&&TMath::Abs(dj)<2){
- AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex());
- for (Int_t itrack=0;itrack<10;itrack++){
- Int_t track = (d->GetTracks())[itrack];
- if (track>=0) {
- AddLabel(milab, track);
- }
- }
- }
- }
- }
- }
-
-
- Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ();
- y/=q; z/=q;
-
-
- Float_t yyyy = y;
-
- z=(z-0.5)*fZpitchSDD;
- z-=fHlSDD;
- Float_t zdet = z;
- Float_t timebin = GetSeg()->Dpx(0);
- Float_t driftTime = (yyyy-0.5)*timebin - fTimeOffsetSDD;
- Float_t xdet = cal->GetDriftPath(driftTime,0);
- xdet=xdet/10000.-fHwSDD;
- if (s) xdet=-xdet;
-
- CorrectPosition(zdet,xdet);
-
- {
- Double_t loc[3]={xdet,0.,zdet},trk[3]={0.,0.,0.};
- mT2L->MasterToLocal(loc,trk);
- y=trk[1];
- z=trk[2];
- }
- q/=5.243; //to have MPV 1 MIP = 86.4 KeV --> this must go to calibr.
- Float_t hit[5] = {y, z, 0.0030*0.0030, 0.0020*0.0020, q};
- Int_t info[3] = {maxj-minj+1, maxi-mini+1, fNlayer[fModule]};
- if (digits) {
- // AliBin *b=&bins[s][idx[k]];
- // AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex());
- {
- //Int_t lab[3];
- //lab[0]=(d->GetTracks())[0];
- //lab[1]=(d->GetTracks())[1];
- //lab[2]=(d->GetTracks())[2];
- //CheckLabels(lab);
- CheckLabels2(milab);
- }
- }
- milab[3]=fNdet[fModule];
-
- AliITSRecPoint cc(milab,hit,info);
- cc.SetType(npeaks);
-
- if(clusters) new (cl[ncl]) AliITSRecPoint(cc);
- else {
- fDetTypeRec->AddRecPoint(cc);
- }
- ncl++;
+ if (idx[k] < 0) continue; //removed peak
+ AliITSRecPoint c;
+ MakeCluster(idx[k], nzBins, bins[s], msk[k], c);
+ //mi change
+ Int_t milab[10];
+ for (Int_t ilab=0;ilab<10;ilab++){
+ milab[ilab]=-2;
+ }
+ Int_t maxi=0,mini=0,maxj=0,minj=0;
+ //AliBin *bmax=&bins[s][idx[k]];
+ //Float_t max = TMath::Max(TMath::Abs(bmax->GetQ())/5.,3.);
+
+ for (Int_t di=-2; di<=2;di++){
+ for (Int_t dj=-3;dj<=3;dj++){
+ Int_t index = idx[k]+di+dj*nzBins;
+ if (index<0) continue;
+ if (index>=nMaxBin) continue;
+ AliBin *b=&bins[s][index];
+ Int_t nAnode=index%nzBins-1;
+ Int_t adcSignal=b->GetQ();
+ if(adcSignal>cal->GetThresholdAnode(nAnode)){
+ if (di>maxi) maxi=di;
+ if (di<mini) mini=di;
+ if (dj>maxj) maxj=dj;
+ if (dj<minj) minj=dj;
+ }
+ //
+ if(digits) {
+ if (TMath::Abs(di)<2&&TMath::Abs(dj)<2){
+ AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex());
+ for (Int_t itrack=0;itrack<10;itrack++){
+ Int_t track = (d->GetTracks())[itrack];
+ if (track>=0) {
+ AddLabel(milab, track);
+ }
+ }
+ }
+ }
+ }
+ }
+
+
+ Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ();
+ y/=q; z/=q;
+ Float_t zAnode=z-0.5; // to have anode in range 0.-255. and centered on the mid of the pitch
+ Float_t timebin=y-0.5; // to have time bin in range 0.-255. amd centered on the mid of the bin
+ if(s==1) zAnode += GetSeg()->NpzHalf(); // right side has anodes from 256. to 511.
+ Float_t zdet = GetSeg()->GetLocalZFromAnode(zAnode);
+ Float_t driftTime = GetSeg()->GetDriftTimeFromTb(timebin) - cal->GetTimeOffset();
+ Float_t driftPathMicron = cal->GetDriftPath(driftTime,zAnode);
+ const Double_t kMicronTocm = 1.0e-4;
+ Float_t xdet=(driftPathMicron-GetSeg()->Dx())*kMicronTocm; // xdet is negative
+ if (s==0) xdet=-xdet; // left side has positive local x
+
+ Float_t corrx=0, corrz=0;
+ cal->GetCorrections(zdet,xdet,corrz,corrx,GetSeg());
+ zdet+=corrz;
+ xdet+=corrx;
+
+ Double_t loc[3]={xdet,0.,zdet},trk[3]={0.,0.,0.};
+ mT2L->MasterToLocal(loc,trk);
+ y=trk[1];
+ z=trk[2];
+
+ q/=cal->GetADC2keV(); //to have MPV 1 MIP = 86.4 KeV
+ Float_t hit[5] = {y, z, 0.0030*0.0030, 0.0020*0.0020, q};
+ Int_t info[3] = {maxj-minj+1, maxi-mini+1, fNlayer[fModule]};
+ if (digits) {
+ // AliBin *b=&bins[s][idx[k]];
+ // AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex());
+ {
+ //Int_t lab[3];
+ //lab[0]=(d->GetTracks())[0];
+ //lab[1]=(d->GetTracks())[1];
+ //lab[2]=(d->GetTracks())[2];
+ //CheckLabels(lab);
+ CheckLabels2(milab);
+ }
+ }
+ milab[3]=fNdet[fModule];
+
+ AliITSRecPoint cc(milab,hit,info);
+ cc.SetType(npeaks);
+ if(clusters) new (cl[ncl]) AliITSRecPoint(cc);
+ else {
+ fDetTypeRec->AddRecPoint(cc);
+ }
+ ncl++;
}
}
}
-
-
-
+//______________________________________________________________________
void AliITSClusterFinderV2SDD::RawdataToClusters(AliRawReader* rawReader,TClonesArray** clusters){
//------------------------------------------------------------
// This function creates ITS clusters from raw data
//------------------------------------------------------------
rawReader->Reset();
AliITSRawStreamSDD inputSDD(rawReader);
+ AliITSDDLModuleMapSDD *ddlmap=(AliITSDDLModuleMapSDD*)fDetTypeRec->GetDDLModuleMapSDD();
+ inputSDD.SetDDLModuleMap(ddlmap);
FindClustersSDD(&inputSDD,clusters);
}
// Actual SDD cluster finder for raw data
//------------------------------------------------------------
Int_t nClustersSDD = 0;
- Int_t kNzBins = fNzSDD + 2;
- Int_t kMaxBin = kNzBins * (fNySDD+2);
- AliBin *binsSDDInit = new AliBin[kMaxBin];
- AliBin *binsSDD1 = new AliBin[kMaxBin];
- AliBin *binsSDD2 = new AliBin[kMaxBin];
- AliBin* bins[2] = {NULL, NULL};
-
+ Int_t nAnodes = GetSeg()->NpzHalf();
+ Int_t nzBins = nAnodes+2;
+ Int_t nTimeBins = GetSeg()->Npx();
+ Int_t nxBins = nTimeBins+2;
+ const Int_t kMaxBin=nzBins*nxBins;
+ AliBin *bins[2];
+ AliBin *ddlbins[kHybridsPerDDL]; // 12 modules (=24 hybrids) of 1 DDL read "in parallel"
+ for(Int_t iHyb=0;iHyb<kHybridsPerDDL;iHyb++) ddlbins[iHyb]=new AliBin[kMaxBin];
// read raw data input stream
- while (kTRUE) {
- Bool_t next = input->Next();
- if (!next || input->IsCompletedModule()) {
+ while (input->Next()) {
+ Int_t iModule = input->GetModuleID();
+ if(iModule<0){
+ AliWarning(Form("Invalid SDD module number %d\n", iModule));
+ continue;
+ }
+
+ Int_t iCarlos =((AliITSRawStreamSDD*)input)->GetCarlosId();
+ Int_t iSide = ((AliITSRawStreamSDD*)input)->GetChannel();
+ Int_t iHybrid=iCarlos*2+iSide;
+ if (input->IsCompletedModule()) {
// when all data from a module was read, search for clusters
- Int_t iModule = input->GetModuleID();
- if (bins[0]) {
- clusters[iModule] = new TClonesArray("AliITSRecPoint");
- fModule = iModule;
- FindClustersSDD(bins, kMaxBin, kNzBins, NULL, clusters[iModule]);
- Int_t nClusters = clusters[iModule]->GetEntriesFast();
- nClustersSDD += nClusters;
- bins[0] = bins[1] = NULL;
+ if(iCarlos<0){
+ AliWarning(Form("Invalid SDD carlos number %d on module %d\n", iCarlos,iModule));
+ continue;
+ }
+ clusters[iModule] = new TClonesArray("AliITSRecPoint");
+ fModule = iModule;
+ bins[0]=ddlbins[iCarlos*2]; // first hybrid of the completed module
+ bins[1]=ddlbins[iCarlos*2+1]; // second hybrid of the completed module
+ FindClustersSDD(bins, kMaxBin, nzBins, NULL, clusters[iModule]);
+ Int_t nClusters = clusters[iModule]->GetEntriesFast();
+ nClustersSDD += nClusters;
+ for(Int_t iBin=0;iBin<kMaxBin; iBin++){
+ ddlbins[iCarlos*2][iBin].Reset();
+ ddlbins[iCarlos*2+1][iBin].Reset();
}
-
- if (!next) break;
- bins[0]=binsSDD1;
- bins[1]=binsSDD2;
- memcpy(binsSDD1,binsSDDInit,sizeof(AliBin)*kMaxBin);
- memcpy(binsSDD2,binsSDDInit,sizeof(AliBin)*kMaxBin);
}else{
// fill the current digit into the bins array
- AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(input->GetModuleID());
+ if(iHybrid<0 || iHybrid>=kHybridsPerDDL){
+ AliWarning(Form("Invalid SDD hybrid number %d on module %d\n", iHybrid,iModule));
+ continue;
+ }
+ AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(iModule);
+ if(cal==0){
+ AliError(Form("Calibration object not present for SDD module %d\n",iModule));
+ continue;
+ }
AliITSresponseSDD* res = (AliITSresponseSDD*)cal->GetResponse();
const char *option=res->ZeroSuppOption();
- Int_t q=input->GetSignal();
- if(!((strstr(option,"1D")) || (strstr(option,"2D")))){
- Float_t baseline = cal->GetBaseline(input->GetCoord1());
- if(q>baseline) q-=(Int_t)baseline;
- else q=0;
- }
- if(q>=cal->GetThresholdAnode(input->GetCoord1())) {
- Int_t iz = input->GetCoord1()+1;
- //Int_t side = ((iz <= fNzSDD) ? 0 : 1);
- Int_t side = ((AliITSRawStreamSDD*)input)->GetChannel();
- // iz -= side*fNzSDD;
- Int_t index = (input->GetCoord2()+1) * kNzBins + iz;
- bins[side][index].SetQ(q);
- bins[side][index].SetMask(1);
- bins[side][index].SetIndex(index);
+ Float_t charge=input->GetSignal();
+ Int_t chan=input->GetCoord1()+nAnodes*iSide;
+ Float_t gain=cal->GetChannelGain(chan);
+ if(strstr(option,"ZS")) charge+=(Float_t)cal->GetZSLowThreshold(iSide);
+ Float_t baseline = cal->GetBaseline(chan);
+ if(charge>baseline) charge-=baseline;
+ else charge=0;
+ if(gain>0.){ // Bad channels have gain=0
+ charge/=gain;
+ if(charge>=cal->GetThresholdAnode(chan)) {
+ Int_t q=(Int_t)(charge+0.5);
+ Int_t iz = input->GetCoord1();
+ Int_t itb = input->GetCoord2();
+ Int_t index = (itb+1) * nzBins + (iz+1);
+ if(index<kMaxBin){
+ ddlbins[iHybrid][index].SetQ(q);
+ ddlbins[iHybrid][index].SetMask(1);
+ ddlbins[iHybrid][index].SetIndex(index);
+ }else{
+ AliWarning(Form("Invalid SDD cell: Anode=%d TimeBin=%d",iz,itb));
+ }
+ }
}
}
}
- delete[] binsSDD1;
- delete[] binsSDD2;
- delete[] binsSDDInit;
-
+ for(Int_t iHyb=0;iHyb<kHybridsPerDDL;iHyb++) delete [] ddlbins[iHyb];
Info("FindClustersSDD", "found clusters in ITS SDD: %d", nClustersSDD);
}
-//_________________________________________________________________________
-void AliITSClusterFinderV2SDD::CorrectPosition(Float_t &z, Float_t&y){
-
- //correction of coordinates using the maps stored in the DB
- AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(fModule);
- static const Int_t knbint = cal->GetMapTimeNBin();
- static const Int_t knbina = cal->Chips()*cal->Channels();
- Float_t stepa = (GetSeg()->Dpz(0))/10000.; //anode pitch in cm
- Float_t stept = (GetSeg()->Dx()/cal->GetMapTimeNBin()/2.)/10.;
-
- Int_t bint = TMath::Abs((Int_t)(y/stept));
- if(y>=0) bint+=(Int_t)(knbint/2.);
- if(bint>knbint) AliError("Wrong bin number!");
-
- Int_t bina = TMath::Abs((Int_t)(z/stepa));
- if(z>=0) bina+=(Int_t)(knbina/2.);
- if(bina>knbina) AliError("Wrong bin number!");
-
- Float_t devz = cal->GetMapACell(bina,bint)/10000.;
- Float_t devx = cal->GetMapTCell(bina,bint)/10000.;
- z+=devz;
- y+=devx;
-
-
-}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff