#include "AliITSDetTypeRec.h"
#include "AliRawReader.h"
#include "AliITSRawStreamSDD.h"
+#include "AliITSRawStreamSDDCompressed.h"
#include "AliITSCalibrationSDD.h"
+#include "AliITSresponseSDD.h"
#include "AliITSDetTypeRec.h"
+#include "AliITSReconstructor.h"
#include "AliITSsegmentationSDD.h"
#include "AliITSdigitSDD.h"
#include "AliITSgeomTGeo.h"
+
ClassImp(AliITSClusterFinderV2SDD)
AliITSClusterFinderV2SDD::AliITSClusterFinderV2SDD(AliITSDetTypeRec* dettyp):AliITSClusterFinderV2(dettyp)
bins[0]=new AliBin[kMaxBin];
bins[1]=new AliBin[kMaxBin];
AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(fModule);
- AliITSresponseSDD* res = (AliITSresponseSDD*)cal->GetResponse();
- const char *option=res->ZeroSuppOption();
+ if(cal==0){
+ AliError(Form("Calibration object not present for SDD module %d\n",fModule));
+ return;
+ }
+
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
- Float_t gain=cal->GetChannelGain(d->GetCoord1());
- Float_t charge=d->GetSignal();
-
- if(!((strstr(option,"1D")) || (strstr(option,"2D")))){
- Float_t baseline = cal->GetBaseline(d->GetCoord1());
- if(charge>baseline) charge-=baseline;
- else charge=0;
- }
-
- if(gain>0) charge/=gain;
- if(charge<cal->GetThresholdAnode(d->GetCoord1())) continue;
- Int_t q=(Int_t)(charge+0.5);
- 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);
+ 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(); // returns expanded signal
+ // (10 bit, low threshold already added)
+ 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);
+ }
}
}
// Actual SDD cluster finder
//------------------------------------------------------------
+ static AliITSRecoParam *repa = NULL;
+ if(!repa){
+ repa = (AliITSRecoParam*) AliITSReconstructor::GetRecoParam();
+ if(!repa){
+ repa = AliITSRecoParam::GetHighFluxParam();
+ AliWarning("Using default AliITSRecoParam class");
+ }
+ }
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 (Int_t i=0; i<nMaxBin; i++) {
+ NoiseSuppress(i,s,nzBins,bins[s],cal);
if (bins[s][i].IsUsed()) continue;
Int_t idx[32]; UInt_t msk[32]; Int_t npeaks=0;
FindPeaks(i, nzBins, bins[s], idx, msk, npeaks);
}
for (k=0; k<npeaks; k++) {
+ if(repa->GetUseUnfoldingInClusterFinderSDD()==kFALSE) msk[k]=msk[0];
MarkPeak(TMath::Abs(idx[k]), nzBins, bins[s], msk[k]);
}
}
- Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ();
+ AliITSresponseSDD* rsdd = fDetTypeRec->GetResponseSDD();
+ 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 driftTime = GetSeg()->GetDriftTimeFromTb(timebin) - rsdd->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
-
- CorrectPosition(zdet,xdet);
+
+ 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
+ q/=rsdd->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) {
// This function creates ITS clusters from raw data
//------------------------------------------------------------
rawReader->Reset();
- AliITSRawStreamSDD inputSDD(rawReader);
- AliITSDDLModuleMapSDD *ddlmap=(AliITSDDLModuleMapSDD*)fDetTypeRec->GetDDLModuleMapSDD();
- inputSDD.SetDDLModuleMap(ddlmap);
- FindClustersSDD(&inputSDD,clusters);
+ AliITSRawStream* inputSDD;
+ if(fDetTypeRec->IsHLTmodeC()==kTRUE){
+ inputSDD=new AliITSRawStreamSDDCompressed(rawReader);
+ }else{
+ inputSDD=new AliITSRawStreamSDD(rawReader);
+ }
+ AliITSDDLModuleMapSDD *ddlmap=(AliITSDDLModuleMapSDD*)fDetTypeRec->GetDDLModuleMapSDD();
+ inputSDD->SetDDLModuleMap(ddlmap);
+ for(Int_t iddl=0; iddl<AliITSDDLModuleMapSDD::GetNDDLs(); iddl++){
+ for(Int_t icar=0; icar<AliITSDDLModuleMapSDD::GetNModPerDDL();icar++){
+ Int_t iMod=ddlmap->GetModuleNumber(iddl,icar);
+ if(iMod==-1) continue;
+ AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(iMod);
+ if(cal==0){
+ AliError(Form("Calibration object not present for SDD module %d\n",iMod));
+ continue;
+ }
+ Bool_t isZeroSupp=cal->GetZeroSupp();
+ if(isZeroSupp){
+ for(Int_t iSid=0; iSid<2; iSid++) inputSDD->SetZeroSuppLowThreshold(iMod-240,iSid,cal->GetZSLowThreshold(iSid));
+ }else{
+ for(Int_t iSid=0; iSid<2; iSid++) inputSDD->SetZeroSuppLowThreshold(iMod-240,iSid,0);
+ }
+ }
+ }
+ FindClustersSDD(inputSDD,clusters);
+ delete inputSDD;
}
void AliITSClusterFinderV2SDD::FindClustersSDD(AliITSRawStream* input,
continue;
}
- Int_t iCarlos =((AliITSRawStreamSDD*)input)->GetCarlosId();
- Int_t iSide = ((AliITSRawStreamSDD*)input)->GetChannel();
+ Int_t iCarlos =input->GetCarlosId();
+ Int_t iSide = input->GetChannel();
Int_t iHybrid=iCarlos*2+iSide;
if (input->IsCompletedModule()) {
// when all data from a module was read, search for clusters
continue;
}
AliITSCalibrationSDD* cal = (AliITSCalibrationSDD*)GetResp(iModule);
- AliITSresponseSDD* res = (AliITSresponseSDD*)cal->GetResponse();
- const char *option=res->ZeroSuppOption();
+ if(cal==0){
+ AliError(Form("Calibration object not present for SDD module %d\n",iModule));
+ continue;
+ }
Float_t charge=input->GetSignal();
Int_t chan=input->GetCoord1()+nAnodes*iSide;
Float_t gain=cal->GetChannelGain(chan);
- if(!((strstr(option,"1D")) || (strstr(option,"2D")))){
- Float_t baseline = cal->GetBaseline(chan);
- if(charge>baseline) charge-=baseline;
- else charge=0;
- }
- if(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));
+ 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));
+ }
}
}
}
Info("FindClustersSDD", "found clusters in ITS SDD: %d", nClustersSDD);
}
+//______________________________________________________________________
+void AliITSClusterFinderV2SDD::NoiseSuppress(Int_t k, Int_t sid,Int_t nzBins, AliBin* bins, AliITSCalibrationSDD* cal) const {
+ // applies zero suppression using the measured noise of each anode
+ // threshold values from ALICE-INT-1999-28 V10
+ Float_t xfactL=2.2;
+ Float_t xfactH=4.0;
+ //
+
+ Int_t iAn=(k%nzBins)-1;
+ if(iAn<0 || iAn>255) return;
+ if(sid==1) iAn+=256;
+ Int_t nLow=0, nHigh=0;
+ Float_t noise=cal->GetNoiseAfterElectronics(iAn);
+ Float_t noisem1=noise;
+ if(iAn>1) noisem1=cal->GetNoiseAfterElectronics(iAn-1);
+ Float_t noisep1=noise;
+ if(iAn<511) noisep1=cal->GetNoiseAfterElectronics(iAn+1);
+ Float_t tL=noise*xfactL;
+ Float_t tH=noise*xfactH;
+ Float_t tLp1=noisep1*xfactL;
+ Float_t tHp1=noisep1*xfactH;
+ Float_t tLm1=noisem1*xfactL;
+ Float_t tHm1=noisem1*xfactH;
+ Float_t cC=bins[k].GetQ();
+ if(cC<=tL){
+ bins[k].SetQ(0);
+ bins[k].SetMask(0xFFFFFFFE);
+ return;
+ }
+ nLow++; // cC is greater than tL
+ if(cC>tH) nHigh++;
+ Int_t sS=bins[k-1].GetQ();
+ if(sS>tLm1) nLow++;
+ if(sS>tHm1) nHigh++;
+ Int_t nN=bins[k+1].GetQ();
+ if(nN>tLp1) nLow++;
+ if(nN>tHp1) nHigh++;
+ Int_t eE=bins[k-nzBins].GetQ();
+ if(eE>tL) nLow++;
+ if(eE>tH) nHigh++;
+ Int_t wW=bins[k+nzBins].GetQ();
+ if(wW>tL) nLow++;
+ if(wW>tH) nHigh++;
+ if(nLow<3 || nHigh<1){
+ bins[k].SetQ(0);
+ bins[k].SetMask(0xFFFFFFFE);
+ }
+}
-//_________________________________________________________________________
-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 knbina = cal->Wings()*cal->Chips()*cal->Channels();
- Int_t bina =(Int_t) GetSeg()->GetAnodeFromLocal(y,z);
- if(bina>knbina) AliError("Wrong bin number along anodes!");
-
- static const Int_t knbint = cal->GetMapTimeNBin();
- const Double_t kMicronTocm = 1.0e-4;
- Float_t stept = GetSeg()->Dx()*kMicronTocm/cal->GetMapTimeNBin();
- Int_t bint = TMath::Abs((Int_t)(y/stept));
- if(bint==knbint) bint-=1;
- if(bint>=knbint) AliError("Wrong bin number along drift direction!");
-
-
- Float_t devz = cal->GetMapACell(bina,bint)*kMicronTocm;
- Float_t devx = cal->GetMapTCell(bina,bint)*kMicronTocm;
- z+=devz;
- y+=devx;
-}