Code from MUON-dev joined
[u/mrichter/AliRoot.git] / MUON / AliMUONSegmentationV04.cxx
CommitLineData
4c039060 1/**************************************************************************
2 * Copyright(c) 1998-1999, 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$Log$
a9e2aefa 18Revision 1.1.2.1 2000/06/09 21:38:15 morsch
19AliMUONSegmentationV04 code from AliMUONSegResV04.cxx
20
4c039060 21*/
22
a897a37a 23/////////////////////////////////////////////////////
24// Segmentation and Response classes version 04 //
25/////////////////////////////////////////////////////
26
a897a37a 27
a9e2aefa 28#include "AliMUONSegmentationV04.h"
29#include <TMath.h>
a897a37a 30
31//___________________________________________
a9e2aefa 32ClassImp(AliMUONSegmentationV04)
a897a37a 33
34
a9e2aefa 35void AliMUONSegmentationV04::Init(AliMUONChamber* Chamber)
a897a37a 36{
a9e2aefa 37 printf("\n Initialise segmentation v04 \n");
a897a37a 38//
39// Fill the arrays fCx (x-contour) and fNpxS (ix-contour) for each sector
40// These arrays help in converting from real to pad co-ordinates and
41// vice versa
a9e2aefa 42//
43// Segmentation is defined by rectangular modules approximating
44// concentric circles as shown below
a897a37a 45//
46// PCB module size in cm
a9e2aefa 47 const Float_t kDxPCB=40, kDyPCB=40;
a897a37a 48// PCB distribution (7 rows with 1+3 segmentation regions)
49 const Int_t kpcb[7][4] = {{1, 2, 2, 2},
50 {0, 3, 2, 2},
51 {0, 2, 2, 2},
52 {0, 0, 3, 3},
53 {0, 0, 2, 3},
54 {0, 0, 0, 4},
55 {0, 0, 0, 3}};
56
57
58//
59// 3 3 3 | 3 3 3
60// 3 3 3 3 | 3 3 3 3
61// 3 3 3 2 2 | 2 2 3 3 3
62// 3 3 3 2 2 2 | 2 2 2 3 3 3
63// 3 3 2 2 1 1 | 1 1 2 2 3 3
64// 3 3 2 2 1 1 1 | 1 1 1 2 2 3 3
65// 3 3 2 2 1 1 0 | 0 1 1 2 2 3 3
66// ------------------------------
67// 3 3 2 2 1 1 0 | 0 1 1 2 2 3 3
68// 3 3 2 2 1 1 1 | 1 1 1 2 2 3 3
69// 3 3 2 2 1 1 | 1 1 2 2 3 3
70// 3 3 3 2 2 2 | 2 2 2 3 3 3
71// 3 3 3 2 2 | 2 2 3 3 3
72// 3 3 3 3 | 3 3 3 3
73// 3 3 3 | 3 3 3
74//
75// number of pad rows per PCB
76//
a9e2aefa 77 Int_t nPyPCB=Int_t(kDyPCB/fDpy);
a897a37a 78//
79// maximum number of pad rows
a9e2aefa 80 fNpy=7*nPyPCB;
a897a37a 81//
82// Calculate padsize along x
83 fDpxD[fNsec-1]=fDpx;
84 if (fNsec > 1) {
85 for (Int_t i=fNsec-2; i>=0; i--){
86 fDpxD[i]=fDpxD[fNsec-1]/fNDiv[i];
a9e2aefa 87 printf("\n test ---dx %d %f \n",i,fDpxD[i]);
a897a37a 88 }
89 }
90//
91// fill the arrays defining the pad segmentation boundaries
92//
93// loop over pcb module rows
94 Int_t iy=0;
95 for (Int_t irow=0; irow<7; irow++) {
96//
97// loop over pads along the anode wires
a9e2aefa 98 for (Int_t i=0; i<=nPyPCB; i++) {
a897a37a 99// iy counts the padrow
100 iy++;
101// Loop over sectors (isec=0 is the dead space surounding the beam pipe)
102 for (Int_t isec=0; isec<4; isec++) {
103 if (isec==0) {
a9e2aefa 104 fNpxS[0][iy]=kpcb[irow][0]*Int_t(kDxPCB/fDpxD[0]);
105 fCx[0][iy]=kpcb[irow][0]*kDxPCB;
a897a37a 106 } else {
107 fNpxS[isec][iy]=fNpxS[isec-1][iy]
a9e2aefa 108 +kpcb[irow][isec]*Int_t(kDxPCB/fDpxD[isec]);
a897a37a 109
110 fCx[isec][iy]=fCx[isec-1][iy]
a9e2aefa 111 +kpcb[irow][isec]*kDxPCB;
a897a37a 112 }
113 } // sectors
114 } // pad raws in module
115 } // PCB rows
a897a37a 116}
117
a9e2aefa 118void AliMUONSegmentationV04::GiveTestPoints(Int_t &n, Float_t *x, Float_t *y)
a897a37a 119{
a9e2aefa 120// Returns test point on the pad plane.
121// Used during determination of the segmoid correction of the COG-method
a897a37a 122 n=3;
123 x[0]=(fCx[1][1]+fCx[0][1])/2/TMath::Sqrt(2.);
124 y[0]=x[0];
125 x[1]=(fCx[2][1]+fCx[1][1])/2/TMath::Sqrt(2.);
126 y[1]=x[1];
127 x[2]=(fCx[3][1]+fCx[2][1])/2/TMath::Sqrt(2.);
128 y[2]=x[2];
a897a37a 129}
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145