Stupid bug fix in new superlight mode (from Zurich airport)
[u/mrichter/AliRoot.git] / ITS / AliITSOnlineSDD.cxx
CommitLineData
348f80b7 1/**************************************************************************
2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16///////////////////////////////////////////////////////////////////
17// //
18// Implementation of the base class for SDD detector algorithms //
19// Origin: F.Prino, Torino, prino@to.infn.it //
20// //
21///////////////////////////////////////////////////////////////////
22
23#include "AliITSOnlineSDD.h"
9aa65e3e 24#include "TH2F.h"
348f80b7 25
26ClassImp(AliITSOnlineSDD)
27//______________________________________________________________________
979b5a5f 28 AliITSOnlineSDD::AliITSOnlineSDD():TObject(),fDDL(0),fCarlos(0),fSide(0),fFirstGoodTB(0),fLastGoodTB(0)
348f80b7 29{
30 // default constructor
750296dd 31 SetFirstGoodTB();
32 SetLastGoodTB();
348f80b7 33}
34//______________________________________________________________________
979b5a5f 35 AliITSOnlineSDD::AliITSOnlineSDD(Int_t nddl, Int_t ncarlos, Int_t sid):TObject(),fDDL(0),fCarlos(0),fSide(0),fFirstGoodTB(0),fLastGoodTB(0)
348f80b7 36{
37 // standard constructor
979b5a5f 38 SetDDL(nddl);
39 SetCarlos(ncarlos);
348f80b7 40 SetDetectorSide(sid);
750296dd 41 SetFirstGoodTB();
42 SetLastGoodTB();
348f80b7 43}
9aa65e3e 44//______________________________________________________________________
45TH2F* AliITSOnlineSDD::ApplyZeroSuppression(TH2F* hRaw, Float_t basl, Int_t tL, Int_t tH){
46 // apply zero suppression
47
48 Int_t nx=hRaw->GetNbinsX();
49 Float_t xmin=hRaw->GetXaxis()->GetXmin();
50 Float_t xmax=hRaw->GetXaxis()->GetXmax();
51 Int_t ny=hRaw->GetNbinsY();
52 Float_t ymin=hRaw->GetYaxis()->GetXmin();
53 Float_t ymax=hRaw->GetYaxis()->GetXmax();
54 TH2F* hZeroSupp=new TH2F(Form("%s_ZeroSupp",hRaw->GetName()),"",nx,xmin,xmax,ny,ymin,ymax);
55 for(Int_t ix=1; ix<=nx;ix++){
56 for(Int_t iy=1; iy<=ny;iy++){
57 // N
58 // Get "quintuple": WCE
59 // S
60 Float_t cC=hRaw->GetBinContent(ix,iy)-basl;
61 Int_t nLow=0, nHigh=0;
62 if(cC<tL){
63 hZeroSupp->SetBinContent(ix,iy,0.);
64 continue;
65 }
66 nLow++; // cC is greater than tL
67 if(cC>tH) nHigh++;
68 Float_t wW=0.;
69 if(iy>0) wW=hRaw->GetBinContent(ix,iy-1)-basl;
70 if(wW>tL) nLow++;
71 if(wW>tH) nHigh++;
72 Float_t eE=0.;
73 if(iy<ny) eE=hRaw->GetBinContent(ix,iy+1)-basl;
74 if(eE>tL) nLow++;
75 if(eE>tH) nHigh++;
76 Float_t nN=0.;
77 if(ix<nx) nN=hRaw->GetBinContent(ix+1,iy)-basl;
78 if(nN>tL) nLow++;
79 if(nN>tH) nHigh++;
80 Float_t sS=0.;
81 if(ix>0) sS=hRaw->GetBinContent(ix-1,iy)-basl;
82 if(sS>tL) nLow++;
83 if(sS>tH) nHigh++;
84 if(nLow>=2 && nHigh>=1){
85 hZeroSupp->SetBinContent(ix,iy,cC);
86 }else{
87 hZeroSupp->SetBinContent(ix,iy,0.);
88 }
89 }
90 }
91 return hZeroSupp;
92}