// Implementation of the ITS clusterer V2 class //
// //
// Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch //
+// Unfolding switch from AliITSRecoParam: D. Elia, INFN Bari //
// //
-///////////////////////////////////////////////////////////////////////////
+////////////////////////////////////////////////////////////////////////////
-#include "AliRun.h"
+#include "AliITSCalibrationSPD.h"
#include "AliITSClusterFinderV2SPD.h"
-#include "AliITSclusterV2.h"
+#include "AliITSRecPoint.h"
+#include "AliITSgeomTGeo.h"
+#include "AliITSDetTypeRec.h"
+#include "AliITSReconstructor.h"
#include "AliRawReader.h"
#include "AliITSRawStreamSPD.h"
-
#include <TClonesArray.h>
-#include "AliITS.h"
-#include "AliITSgeom.h"
#include "AliITSdigitSPD.h"
+#include "AliITSFOSignalsSPD.h"
ClassImp(AliITSClusterFinderV2SPD)
-extern AliRun *gAlice;
-
-AliITSClusterFinderV2SPD::AliITSClusterFinderV2SPD():AliITSClusterFinderV2(){
+//__________________________________________________________________________
+AliITSClusterFinderV2SPD::AliITSClusterFinderV2SPD(AliITSDetTypeRec* dettyp):AliITSClusterFinder(dettyp),
+fLastSPD1(AliITSgeomTGeo::GetModuleIndex(2,1,1)-1),
+fNySPD(256),
+fNzSPD(160),
+fYpitchSPD(0.0050),
+fZ1pitchSPD(0.0425),
+fZ2pitchSPD(0.0625),
+fHwSPD(0.64),
+fHlSPD(3.48){
//Default constructor
- AliITSgeom *geom=(AliITSgeom*)fITS->GetITSgeom();
-
- fLastSPD1=geom->GetModuleIndex(2,1,1)-1;
- fNySPD=256; fNzSPD=160;
- fYpitchSPD=0.0050;
- fZ1pitchSPD=0.0425; fZ2pitchSPD=0.0625;
- fHwSPD=0.64; fHlSPD=3.48;
fYSPD[0]=0.5*fYpitchSPD;
for (Int_t m=1; m<fNySPD; m++) fYSPD[m]=fYSPD[m-1]+fYpitchSPD;
fZSPD[0]=fZ1pitchSPD;
}
}
-
-
+//__________________________________________________________________________
void AliITSClusterFinderV2SPD::FindRawClusters(Int_t mod){
-
//Find clusters V2
SetModule(mod);
FindClustersSPD(fDigits);
}
-
+//__________________________________________________________________________
void AliITSClusterFinderV2SPD::RawdataToClusters(AliRawReader* rawReader, TClonesArray** clusters){
- //------------------------------------------------------------
+ //------------------------------------------------------------
// This function creates ITS clusters from raw data
//------------------------------------------------------------
rawReader->Reset();
FindClustersSPD(&inputSPD, clusters);
}
-
+//__________________________________________________________________________
Int_t AliITSClusterFinderV2SPD::ClustersSPD(AliBin* bins, TClonesArray* digits,TClonesArray* clusters,Int_t maxBins,Int_t nzbins,Int_t iModule,Bool_t rawdata){
//Cluster finder for SPD (from digits and from rawdata)
+ 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(iModule);
+
+ if (repa->GetSPDRemoveNoisyFlag()) {
+ // Loop on noisy pixels and reset them
+ AliITSCalibrationSPD *cal =
+ (AliITSCalibrationSPD*) fDetTypeRec->GetCalibrationModel(iModule);
+ for(Int_t ipix = 0; ipix<cal->GetNrBad(); ipix++){
+ Int_t row, col;
+ cal->GetBadPixel(ipix,row,col);
+ printf(" module %d row %d col %d \n",iModule,row,col);
+ Int_t index = (row+1) * nzbins + (col+1);
+
+ bins[index].SetQ(0);
+ bins[index].SetMask(0xFFFFFFFE);
+ }
+ }
+
+ if (repa->GetSPDRemoveDeadFlag()) {
+ // Loop on dead pixels and reset them
+ AliITSCalibrationSPD *cal =
+ (AliITSCalibrationSPD*) fDetTypeRec->GetSPDDeadModel(iModule);
+ if (cal->IsBad()) return 0; // if all ladder is dead, return to save time
+ for(Int_t ipix = 0; ipix<cal->GetNrBad(); ipix++){
+ Int_t row, col;
+ cal->GetBadPixel(ipix,row,col);
+ Int_t index = (row+1) * nzbins + (col+1);
+ bins[index].SetQ(0);
+ bins[index].SetMask(0xFFFFFFFE);
+ }
+ }
+
Int_t nclu=0;
for(Int_t iBin =0; iBin < maxBins;iBin++){
if(bins[iBin].IsUsed()) continue;
zmin=dig->GetCoord1();
zmax=zmin;
}
+ if(iModule == 24 || iModule == 25) printf("\n");
for (Int_t idx = 0; idx < nBins; idx++) {
Int_t iy;
Int_t iz;
AliITSdigitSPD* dig = (AliITSdigitSPD*)digits->UncheckedAt(idxBins[idx]);
iy = dig->GetCoord2();
iz = dig->GetCoord1();
+ //if(iModule == 24 || iModule == 25) printf(" || iy %d iz %d in Module %d \n",iy,iz,iModule);
}
if (ymin > iy) ymin = iy;
if (ymax < iy) ymax = iy;
if((iModule <= fLastSPD1) &&idy<3) idy=3;
if((iModule > fLastSPD1) &&idy<4) idy=4;
Int_t idz=3;
+
+ // Switch the unfolding OFF/ON
+ if(!repa->GetUseUnfoldingInClusterFinderSPD()) {
+ idy=ymax-ymin+1;
+ idz=zmax-zmin+1;
+ }
+
for(Int_t iiz=zmin; iiz<=zmax;iiz+=idz){
for(Int_t iiy=ymin;iiy<=ymax;iiy+=idy){
if(TMath::Abs(iz-iiz)>0.75*idz) continue;
}
ndigits++;
- Float_t qBin;
+ Float_t qBin=0.;
if(rawdata) qBin = bins[idxBins[idx]].GetQ();
if(!rawdata){
AliITSdigitSPD* dig = (AliITSdigitSPD*)digits->UncheckedAt(idxBins[idx]);
z /= q;
y -= fHwSPD;
z -= fHlSPD;
+
Float_t hit[5]; //y,z,sigma(y)^2, sigma(z)^2, charge
- hit[0] = -(-y+fYshift[iModule]);
- if(iModule <= fLastSPD1) hit[0] = -hit[0];
- hit[1] = -z+fZshift[iModule];
+ {
+ Double_t loc[3]={y,0.,z},trk[3]={0.,0.,0.};
+ mT2L->MasterToLocal(loc,trk);
+ hit[0]=trk[1];
+ hit[1]=trk[2];
+ }
hit[2] = fYpitchSPD*fYpitchSPD/12.;
hit[3] = fZ1pitchSPD*fZ1pitchSPD/12.;
- hit[4] = (zmax-zmin+1)*100 + (ymax-ymin+1);
+ hit[4] = 1.;
+
if(!rawdata) milab[3]=fNdet[iModule];
Int_t info[3] = {ymax-ymin+1,zmax-zmin+1,fNlayer[iModule]};
if(!rawdata){
- AliITSclusterV2 cl(milab,hit,info);
- fITS->AddClusterV2(cl);
+ AliITSRecPoint cl(milab,hit,info);
+ fDetTypeRec->AddRecPoint(cl);
}
else{
Int_t label[4]={milab[0],milab[1],milab[2],milab[3]};
new (clusters->AddrAt(nclu))
- AliITSclusterV2(label, hit,info);
+ AliITSRecPoint(label, hit,info);
}
nclu++;
}// for iiy
return nclu;
}
-
-
-
-
-void AliITSClusterFinderV2SPD::FindClustersSPD(AliITSRawStream* input,
+//__________________________________________________________________________
+void AliITSClusterFinderV2SPD::FindClustersSPD(AliITSRawStreamSPD* input,
TClonesArray** clusters)
{
//------------------------------------------------------------
- // Actual SPD cluster finder for raw data
+ // SPD cluster finder for raw data (this method is called once per event)
+ // Now also fills fast-or fired map
//------------------------------------------------------------
+
Int_t nClustersSPD = 0;
Int_t kNzBins = fNzSPD + 2;
Int_t kNyBins = fNySPD + 2;
Int_t kMaxBin = kNzBins * kNyBins;
+ AliBin *binsSPD = new AliBin[kMaxBin];
+ AliBin *binsSPDInit = new AliBin[kMaxBin];
AliBin* bins = NULL;
// read raw data input stream
// when all data from a module was read, search for clusters
if (bins) {
- clusters[iModule] = new TClonesArray("AliITSclusterV2");
+ clusters[iModule] = new TClonesArray("AliITSRecPoint");
Int_t nClusters = ClustersSPD(bins,0,clusters[iModule],kMaxBin,kNzBins,iModule,kTRUE);
nClustersSPD += nClusters;
- delete bins;
+ bins = NULL;
}
if (!next) break;
- bins = new AliBin[kMaxBin];
+ bins = binsSPD;
+ memcpy(binsSPD,binsSPDInit,sizeof(AliBin)*kMaxBin);
}
- // fill the current digit into the bins array
- Int_t index = (input->GetCoord2()+1) * kNzBins + (input->GetCoord1()+1);
- bins[index].SetIndex(index);
- bins[index].SetMask(1);
- bins[index].SetQ(1);
+ if (next) {
+ // fill the current digit into the bins array
+ Int_t index = (input->GetCoord2()+1) * kNzBins + (input->GetCoord1()+1);
+ bins[index].SetIndex(index);
+ bins[index].SetQ(1);
+ bins[index].SetMask(1);
+ }
}
+ delete [] binsSPDInit;
+ delete [] binsSPD;
+
+ // AliDebug(1,Form("found clusters in ITS SPD: %d", nClustersSPD));
Info("FindClustersSPD", "found clusters in ITS SPD: %d", nClustersSPD);
+
+ // Fill the FastOr fired map
+ for (UInt_t eq=0; eq<20; eq++) {
+ for (UInt_t hs=0; hs<6; hs++) {
+ for (UInt_t chip=0; chip<10; chip++) {
+ if (input->GetFastOrSignal(eq,hs,chip)) {
+ fDetTypeRec->SetFastOrFiredMapOnline(eq,hs,chip);
+ }
+ }
+ }
+ }
+
}
-
-
-
+//__________________________________________________________________________
void AliITSClusterFinderV2SPD::FindClustersSPD(TClonesArray *digits) {
//------------------------------------------------------------
- // Actual SPD cluster finder
+ // SPD cluster finder for digits (this method is called for each module)
+ // Now also fills the fast-or fired map
//------------------------------------------------------------
Int_t ndigits=digits->GetEntriesFast();
AliBin *bins=new AliBin[kMAXBIN];
- Int_t k;
- AliITSdigitSPD *d=0;
- for (k=0; k<ndigits; k++) {
- d=(AliITSdigitSPD*)digits->UncheckedAt(k);
- Int_t i=d->GetCoord2()+1; //y
- Int_t j=d->GetCoord1()+1;
- Int_t index=i*kNzBins+j;
- bins[index].SetIndex(k);
- bins[index].SetMask(1);
+ Int_t idig;
+ AliITSdigitSPD *digit=0;
+ for (idig=0; idig<ndigits; idig++) {
+ digit=(AliITSdigitSPD*)digits->UncheckedAt(idig);
+ Int_t i=digit->GetCoord2()+1; //y
+ Int_t j=digit->GetCoord1()+1;
+ Int_t index=i*kNzBins+j;
+
+ bins[index].SetIndex(idig);
+ bins[index].SetQ(1);
+ bins[index].SetMask(1);
}
+
- ClustersSPD(bins,digits,0,kMAXBIN,kNzBins,fModule,kFALSE);
+ Int_t nClustersSPD = ClustersSPD(bins,digits,0,kMAXBIN,kNzBins,fModule,kFALSE);
delete [] bins;
-}
-
+ AliDebug(1,Form("found clusters in ITS SPD: %d", nClustersSPD));
+
+ // Fill the FastOr fired map
+ AliITSFOSignalsSPD* foSignals = fDetTypeRec->GetFOSignals();
+ if (foSignals) {
+ UInt_t eq = AliITSRawStreamSPD::GetOnlineEqIdFromOffline(fModule);
+ UInt_t hs = AliITSRawStreamSPD::GetOnlineHSFromOffline(fModule);
+ for(UInt_t ch=0; ch<5; ch++) {
+ UInt_t chip = AliITSRawStreamSPD::GetOnlineChipFromOffline(fModule,ch*32);
+ if (foSignals->GetSignal(eq,hs,chip)) {
+ fDetTypeRec->SetFastOrFiredMapOnline(eq,hs,chip);
+ }
+ }
+ }
+
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff