New plots for trending injector efficiencies (Melinda)
[u/mrichter/AliRoot.git] / ITS / AliITSMapSDD.cxx
CommitLineData
efb451bf 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/* $Id$ */
17
18///////////////////////////////////////////////////////////////////
19// //
20// Implementation of the base class for SDD map corrections //
21// Origin: F.Prino, Torino, prino@to.infn.it //
22// //
23///////////////////////////////////////////////////////////////////
24
25#include "TH1F.h"
26#include "TH2F.h"
27#include "AliITSMapSDD.h"
54af1add 28#include "AliITSCorrMap1DSDD.h"
ba0a07bf 29
efb451bf 30ClassImp(AliITSMapSDD)
31//______________________________________________________________________
54af1add 32AliITSMapSDD::AliITSMapSDD():TNamed("defaultmap","")
efb451bf 33{
34 // default constructor
54af1add 35 for(Int_t iAn=0;iAn<fgkNAnodPts; iAn++){
36 for(Int_t iDr=0;iDr<fgkNDrifPts; iDr++){
37 fMap[iAn][iDr]=0;
38 }
39 }
efb451bf 40}
41//______________________________________________________________________
54af1add 42AliITSMapSDD::AliITSMapSDD(Char_t *mapname):TNamed(mapname,"")
efb451bf 43{
44 // standard constructor
54af1add 45 for(Int_t iAn=0;iAn<fgkNAnodPts; iAn++){
46 for(Int_t iDr=0;iDr<fgkNDrifPts; iDr++){
47 fMap[iAn][iDr]=0;
48 }
49 }
50}
51
52//______________________________________________________________________
53void AliITSMapSDD::SetMap(TH2F* hmap){
54 // Fill map staring from 2D histo
55 // with anodes on x axis and drift dist. on y axis
56 for(Int_t iAn=0;iAn<fgkNAnodPts; iAn++){
57 for(Int_t iDr=0;iDr<fgkNDrifPts; iDr++){
58 SetCellContent(iAn,iDr,hmap->GetBinContent(iAn+1,iDr+1));
59 }
60 }
efb451bf 61}
62//______________________________________________________________________
c9a38d3d 63Float_t AliITSMapSDD::GetCorrection(Float_t z, Float_t x, AliITSsegmentationSDD *seg){
64 // returns correction in cm starting from local coordinates on the module
65 const Double_t kMicronTocm = 1.0e-4;
66 Int_t nAnodes=seg->Npz();
67 Int_t nAnodesHybrid=seg->NpzHalf();
68 Int_t bina =(Int_t) seg->GetAnodeFromLocal(x,z);
69 if(bina>nAnodes) AliError("Wrong anode anumber!");
70 if(bina>=nAnodesHybrid) bina-=nAnodesHybrid;
54af1add 71 Float_t stept = seg->Dx()*kMicronTocm/(Float_t)fgkNDrifPts;
c9a38d3d 72 Int_t bint = TMath::Abs((Int_t)(x/stept));
54af1add 73 if(bint==fgkNDrifPts) bint-=1;
74 if(bint>=fgkNDrifPts) AliError("Wrong bin number along drift direction!");
c9a38d3d 75 return kMicronTocm*GetCellContent(bina,bint);
76}
77//______________________________________________________________________
efb451bf 78TH2F* AliITSMapSDD::GetMapHisto() const{
79 // Returns a TH2F histogram with map of residuals
2c4e6a6a 80 TString hname;
81 hname.Form("h%s",GetName());
82 TH2F* hmap=new TH2F(hname.Data(),"",fgkNAnodPts,-0.5,255.5,fgkNDrifPts,0.,35.);
54af1add 83 for(Int_t iAn=0;iAn<fgkNAnodPts; iAn++){
84 for(Int_t iDr=0;iDr<fgkNDrifPts; iDr++){
ab516493 85 hmap->SetBinContent(iAn+1,iDr+1,GetCellContent(iAn,iDr));
efb451bf 86 }
87 }
88 return hmap;
89}
90//______________________________________________________________________
91TH1F* AliITSMapSDD::GetResidualDistr(Float_t dmin, Float_t dmax) const{
92 // Returns a TH1F histogram with distribution of residual
2c4e6a6a 93 TString hname;
94 hname.Form("hd%s",GetName());
95 TH1F* hd=new TH1F(hname.Data(),"",100,dmin,dmax);
54af1add 96 for(Int_t iAn=0;iAn<fgkNAnodPts; iAn++){
97 for(Int_t iDr=0;iDr<fgkNDrifPts; iDr++){
ab516493 98 hd->Fill(GetCellContent(iAn,iDr));
efb451bf 99 }
100 }
101 return hd;
102}
54af1add 103//______________________________________________________________________
104AliITSCorrMapSDD* AliITSMapSDD::ConvertToNewFormat() const{
105 // convert correction map to new format
106 Char_t* name=(Char_t*)GetName();
107 AliITSCorrMapSDD* newmap=new AliITSCorrMap1DSDD(name,fgkNDrifPts);
108 for(Int_t i=0; i<fgkNDrifPts; i++){
109 newmap->SetCellContent(0,i,GetCellContent(0,i));
110 }
111 return newmap;
112}
113