#include "AliITSDetTypeRec.h"
#include "AliRawReader.h"
#include "AliITSRawStreamSDD.h"
+#include "AliITSRawStreamSDDCompressed.h"
#include "AliITSCalibrationSDD.h"
#include "AliITSresponseSDD.h"
#include "AliITSDetTypeRec.h"
ClassImp(AliITSClusterFinderV2SDD)
-AliITSClusterFinderV2SDD::AliITSClusterFinderV2SDD(AliITSDetTypeRec* dettyp):AliITSClusterFinderV2(dettyp)
+AliITSClusterFinderV2SDD::AliITSClusterFinderV2SDD(AliITSDetTypeRec* dettyp):AliITSClusterFinder(dettyp)
{
//Default constructor
Int_t itb=d->GetCoord2();
Int_t iSide=0;
if (ian >= nAnodes) iSide=1;
- Float_t gain=cal->GetChannelGain(ian);
+ Float_t gain=cal->GetChannelGain(ian)/fDetTypeRec->GetAverageGainSDD();
Float_t charge=d->GetSignal(); // returns expanded signal
// (10 bit, low threshold already added)
Float_t baseline = cal->GetBaseline(ian);
void AliITSClusterFinderV2SDD::
FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
- TClonesArray *digits, TClonesArray *clusters) {
+ TClonesArray *digits, TClonesArray *clusters, Int_t jitter) {
//------------------------------------------------------------
// Actual SDD cluster finder
//------------------------------------------------------------
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(NoiseSuppress(i,s,nzBins,bins[s],cal)) continue;
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);
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++){
+
+ for (Int_t di=-5; di<=5;di++){
+ for (Int_t dj=-10;dj<=10;dj++){
Int_t index = idx[k]+di+dj*nzBins;
if (index<0) continue;
if (index>=nMaxBin) continue;
}
}
}
-
+ Int_t clSizAnode=maxi-mini+1;
+ Int_t clSizTb=maxj-minj+1;
+ if(repa->GetUseSDDClusterSizeSelection()){
+ if(clSizTb==1) continue; // cut common mode noise spikes
+ if(clSizAnode>3) continue; // cut common mode noise spikes
+ if(clSizTb>10) continue; // cut clusters on noisy anodes
+ }
AliITSresponseSDD* rsdd = fDetTypeRec->GetResponseSDD();
Float_t y=c.GetY(),z=c.GetZ(), q=c.GetQ();
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) - rsdd->GetTimeOffset();
- Float_t driftPathMicron = cal->GetDriftPath(driftTime,zAnode);
+ Float_t driftTimeUncorr = GetSeg()->GetDriftTimeFromTb(timebin)+jitter*rsdd->GetCarlosRXClockPeriod();
+ Float_t driftTime=driftTimeUncorr-rsdd->GetTimeZero(fModule);
+ Float_t driftSpeed = cal->GetDriftSpeedAtAnode(zAnode) + rsdd->GetDeltaVDrift(fModule);
+ Float_t driftPathMicron = driftTime*driftSpeed;
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;
+ if(repa->GetUseSDDCorrectionMaps()){
+ 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);
z=trk[2];
q/=rsdd->GetADC2keV(); //to have MPV 1 MIP = 86.4 KeV
+ if(cal-> IsAMAt20MHz()) q*=2.; // account for 1/2 sampling freq.
+ if(q<repa->GetMinClusterChargeSDD()) continue; // remove noise clusters
+
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]};
+ Int_t info[3] = {clSizTb, clSizAnode, fNlayer[fModule]};
if (digits) {
// AliBin *b=&bins[s][idx[k]];
// AliITSdigitSDD* d=(AliITSdigitSDD*)digits->UncheckedAt(b->GetIndex());
}
}
milab[3]=fNdet[fModule];
-
AliITSRecPoint cc(milab,hit,info);
cc.SetType(npeaks);
+ cc.SetDriftTime(driftTimeUncorr);
if(clusters) new (cl[ncl]) AliITSRecPoint(cc);
else {
fDetTypeRec->AddRecPoint(cc);
// This function creates ITS clusters from raw data
//------------------------------------------------------------
rawReader->Reset();
- AliITSRawStreamSDD inputSDD(rawReader);
+ AliITSRawStream* inputSDD;
+ if(fDetTypeRec->IsHLTmodeC()==kTRUE){
+ inputSDD=new AliITSRawStreamSDDCompressed(rawReader);
+ }else{
+ inputSDD=new AliITSRawStreamSDD(rawReader);
+ }
+
AliITSDDLModuleMapSDD *ddlmap=(AliITSDDLModuleMapSDD*)fDetTypeRec->GetDDLModuleMapSDD();
- inputSDD.SetDDLModuleMap(ddlmap);
+ 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);
}
Bool_t isZeroSupp=cal->GetZeroSupp();
if(isZeroSupp){
- for(Int_t iSid=0; iSid<2; iSid++) inputSDD.SetZeroSuppLowThreshold(iMod-240,iSid,cal->GetZSLowThreshold(iSid));
+ 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);
+ for(Int_t iSid=0; iSid<2; iSid++) inputSDD->SetZeroSuppLowThreshold(iMod-240,iSid,0);
}
}
}
- FindClustersSDD(&inputSDD,clusters);
-
+ FindClustersSDD(inputSDD,clusters);
+ delete inputSDD;
}
void AliITSClusterFinderV2SDD::FindClustersSDD(AliITSRawStream* input,
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];
+ Int_t vectModId[kModulesPerDDL];
+ for(Int_t iMod=0; iMod<kModulesPerDDL; iMod++) vectModId[iMod]=-1;
+
// read raw data input stream
while (input->Next()) {
Int_t iModule = input->GetModuleID();
AliWarning(Form("Invalid SDD module number %d\n", iModule));
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
- 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();
+ // store the module number
+ vectModId[iCarlos]=iModule;
+ }
+ else if (input->IsCompletedDDL()) {
+ // when all data from a DDL was read, search for clusters
+ Int_t jitter=input->GetJitter();
+ for(Int_t iMod=0; iMod<kModulesPerDDL; iMod++){
+ if(vectModId[iMod]>=0){
+ fModule = vectModId[iMod];
+ clusters[fModule] = new TClonesArray("AliITSRecPoint");
+ bins[0]=ddlbins[iMod*2]; // first hybrid of the module
+ bins[1]=ddlbins[iMod*2+1]; // second hybrid of the module
+ FindClustersSDD(bins, kMaxBin, nzBins, NULL, clusters[fModule],jitter);
+ Int_t nClusters = clusters[fModule]->GetEntriesFast();
+ nClustersSDD += nClusters;
+ vectModId[iMod]=-1;
+ }
+ for(Int_t iBin=0;iBin<kMaxBin; iBin++){
+ ddlbins[iMod*2][iBin].Reset();
+ ddlbins[iMod*2+1][iBin].Reset();
+ }
}
}else{
// fill the current digit into the bins array
}
Float_t charge=input->GetSignal();
Int_t chan=input->GetCoord1()+nAnodes*iSide;
- Float_t gain=cal->GetChannelGain(chan);
+ Float_t gain=cal->GetChannelGain(chan)/fDetTypeRec->GetAverageGainSDD();;
Float_t baseline = cal->GetBaseline(chan);
if(charge>baseline) charge-=baseline;
else charge=0;
}
//______________________________________________________________________
-void AliITSClusterFinderV2SDD::NoiseSuppress(Int_t k, Int_t sid,Int_t nzBins, AliBin* bins, AliITSCalibrationSDD* cal) const {
+Bool_t 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
+ // returns kTRUE if the digit should eb noise suppressed, kFALSE if it should be kept
Float_t xfactL=2.2;
Float_t xfactH=4.0;
//
Int_t iAn=(k%nzBins)-1;
- if(iAn<0 || iAn>255) return;
+ if(iAn<0 || iAn>255) return kTRUE;
if(sid==1) iAn+=256;
Int_t nLow=0, nHigh=0;
Float_t noise=cal->GetNoiseAfterElectronics(iAn);
Float_t tHp1=noisep1*xfactH;
Float_t tLm1=noisem1*xfactL;
Float_t tHm1=noisem1*xfactH;
- Float_t cC=bins[k].GetQ();
+ Int_t cC=bins[k].GetQ();
if(cC<=tL){
bins[k].SetQ(0);
bins[k].SetMask(0xFFFFFFFE);
- return;
+ return kTRUE;;
}
nLow++; // cC is greater than tL
if(cC>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);
- }
+ if(nLow<3 || nHigh<1) return kTRUE;
+ else return kFALSE;
}