Adding definition of Dmax in constants and semantic modifications for the generation...
[u/mrichter/AliRoot.git] / MUON / AliMUONConstants.cxx
CommitLineData
f665c1ea 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
88cb7938 16/* $Id$ */
f665c1ea 17
18#include "AliMUONConstants.h"
19
fe126bce 20#include "AliLog.h"
21
9d5f6a64 22#include "TMath.h"
fe126bce 23#include "TClass.h"
9d5f6a64 24
25///
26/// This class holds various constants to be used in many places,
27/// such as the number of tracking and trigger chambers,
28/// some geometrical constants (to build the initial geometry for simulation)
29/// and mathieson distribution default values.
30/// Those constants should as much as possible replace hard-coded values
31/// which are to be considered strictly illegal in the MUON code (or any code,
32/// by the way).
33///
f665c1ea 34
37c0cd40 35Int_t AliMUONConstants::fgNCh = 14;
f665c1ea 36Int_t AliMUONConstants::fgNTrackingCh = 10;
37Int_t AliMUONConstants::fgNTriggerCh = 4;
38Int_t AliMUONConstants::fgNTriggerCircuit = 234;
edb96c6a 39Int_t AliMUONConstants::fgNDetElem = 228;
24c9eaa0 40
89b0b88d 41
42
e516b01d 43Float_t AliMUONConstants::fgDefaultChamberZ[14] =
89b0b88d 44 {-526.16, -545.24, -676.4, -695.4, // St12
45 -967.5, -998.5, -1276.5, -1307.5, -1406.6, -1437.6,// updated 08/05, EDMS id 335328 (A. Tournaire)
e516b01d 46 -1603.5, -1620.5, -1703.5, -1720.5}; // M1 & M2
47
43adb63d 48Float_t AliMUONConstants::fgDefaultRatioTriggerChamber[4] =
49{1., 1.01060, 1.06236, 1.07296};
50
51
2c9844e7 52Float_t AliMUONConstants::fgSt345inclination = 0.794; // in degrees, rotation axis is X axis
53
e516b01d 54Float_t AliMUONConstants::fgDzCh = 15.5/2.;
55Float_t AliMUONConstants::fgDzSlat = 8.5/2.;
56
1c334adf 57Float_t AliMUONConstants::fgSqrtKx3 = 0.7131;
58Float_t AliMUONConstants::fgSqrtKy3 = 0.7642;
7e4a628d 59
1c334adf 60Float_t AliMUONConstants::fgSqrtKx3St1 = 0.7000;
61Float_t AliMUONConstants::fgSqrtKy3St1 = 0.7550;
7e4a628d 62
1c334adf 63Float_t AliMUONConstants::fgChargeCorrel = 0.11;
f4a7360f 64Float_t AliMUONConstants::fgChargeCorrelSt1 = 1.0; //???
1c334adf 65Float_t AliMUONConstants::fgPitch = 0.25;
66Float_t AliMUONConstants::fgPitchSt1 = 0.21;
c47b11f0 67
f4a7360f 68// From Alain TOURNAIRE
69// ALICE / ALICE Engineering baseline / Dimuonspectrometer (DIS) v7-1
70// EDMS Id 335328 for "search in EDMS
124dd59f 71// These are the diameter (Dmin == innner and Dmax - outner) values of the active surface
72// In the case of Dmax, the value corresponds to the maximum diameter of the active surface with 2pi coverture in phi
73Float_t AliMUONConstants::fgDmin[7] = { 36.4, 46.2, 63.0, 79.0, 79.0, 99.0, 100.0}; // cm
74Float_t AliMUONConstants::fgDmax[7] = { 176.6, 229.0, 308.84, 418.2, 522.0, 850.0, 900.0}; // cm
f4a7360f 75
e118b27e 76Int_t AliMUONConstants::fgMaxZoom = 20;
ecfa008b 77
63ed9c6b 78ClassImp(AliMUONConstants)
79
24c9eaa0 80//______________________________________________________________________________
29fc2c86 81Int_t AliMUONConstants::ChamberNumber(Float_t z)
82{
83 // return chamber number according z position of hit. Should be taken from geometry ?
84
85 Float_t dMaxChamber = DzSlat() + DzCh() + 0.25; // cm st 3 &4 & 5
fe126bce 86 dMaxChamber += 3.00; // factor for inclination of chamber
87 // dMaxChamber += Rmax(4) * TMath::Sin(fgSt345inclination*TMath::Pi()/360.);
88 // factor for inclination of chamber
29fc2c86 89 if ( z > (DefaultChamberZ(4)+50.)) dMaxChamber = 7.; // cm stations 1 & 2
90 Int_t iChamber;
91
92 for (iChamber = 0; iChamber < 10; iChamber++) {
29fc2c86 93 if (TMath::Abs(z-DefaultChamberZ(iChamber)) < dMaxChamber) {
94 return iChamber;
95 }
96 }
fe126bce 97
98 if ( z > DefaultChamberZ(NTrackingCh()-1) ) {
99 AliWarningClass(Form("No chamber number found for z = %f",z));
100 // for (iChamber = 0; iChamber < 10; iChamber++) {
101 // cout << iChamber << " zpos: " << DefaultChamberZ(iChamber)
102 // << " from " << DefaultChamberZ(iChamber) + dMaxChamber
103 // << " to " << DefaultChamberZ(iChamber) - dMaxChamber
104 // << " delta " << dMaxChamber << endl;
105 //}
106 }
29fc2c86 107 return -1;
108}