StRoot  1
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
digtransport.cxx
1 // //
3 // DIGTransport //
4 // //
5 // Contains charge transport models parameters //
6 // //
7 // //
8 // //
9 // //
10 // //
11 // //
13 #include <digtransport.h>
14 
15 
16 #include <TROOT.h> // for gROOT object
17 #include <TMath.h>
18 #include <TMatrixD.h>
19 #include <TCanvas.h>
20 #include <TGraph.h>
21 #include <TAxis.h>
22 #include <TRandom3.h>
23 #include <TFile.h>
24 #include <TTree.h>
25 #include <TBranch.h>
26 #include <TClonesArray.h>
27 
28 //include other classes.h:
29 
30 
31 using namespace std;
32 
33 //==============================================================================
34 ClassImp(DIGTransport)
35 //______________________________________________________________________________
36 //
38 {
39  //
40  // default constructor
41  //
42 }
43 //______________________________________________________________________________
44 //
45 DIGTransport::DIGTransport(Float_t myvar)
46 {
47  fMyvar = myvar;
48 }
49 //______________________________________________________________________________
50 //
51 DIGTransport::~DIGTransport() {
52  //
53  // virtual destructor
54  //
55 }
56 //______________________________________________________________________________
57 //
58 DIGTransport::DIGTransport(DIGTransport & adigtransport) : TObject()
59 {
60  fMyvar = adigtransport.GetMyvar();
61 
62 }
63 //______________________________________________________________________________
64 //
65 void DIGTransport::Clear(const Option_t *)
66 {
67  // delete pointers. fDIGParticleArray->Clear("C");
68 
69 }
70 //______________________________________________________________________________
71 //
72 void DIGTransport::SetLorentz2DModel_Cp0(Float_t Lorentz2DModel_Cp0){
73  fLorentz2DModel_Cp0 = Lorentz2DModel_Cp0;
74 }
75 //______________________________________________________________________________
76 //
77 void DIGTransport::SetLorentz2DModel_Cp1(Float_t Lorentz2DModel_Cp1){
78  fLorentz2DModel_Cp1 = Lorentz2DModel_Cp1;
79 }
80 //______________________________________________________________________________
81 //
82 void DIGTransport::SetRangeLimit_InPitchUnit(Float_t RangeLimit_InPitchUnit){
83  fRangeLimit_InPitchUnit = RangeLimit_InPitchUnit;
84 }
85 //______________________________________________________________________________
86 //
87 void DIGTransport::SetChargeModel(Int_t ChargeModel){
88  fChargeModel = ChargeModel;
89 }
90 //______________________________________________________________________________
91 //
92 void DIGTransport::SetGauss2DModel_sigma1_Cp0(Float_t Gauss2DModel_sigma1_Cp0){
93  fGauss2DModel_sigma1_Cp0 = Gauss2DModel_sigma1_Cp0;
94 }
95 //______________________________________________________________________________
96 //
97 void DIGTransport::SetGauss2DModel_sigma1_Cp1(Float_t Gauss2DModel_sigma1_Cp1){
98  fGauss2DModel_sigma1_Cp1 = Gauss2DModel_sigma1_Cp1;
99 }
100 //______________________________________________________________________________
101 //
102 void DIGTransport::SetGauss2DModel_sigma2_Cp0(Float_t Gauss2DModel_sigma2_Cp0){
103  fGauss2DModel_sigma2_Cp0 = Gauss2DModel_sigma2_Cp0;
104 }
105 //______________________________________________________________________________
106 //
107 void DIGTransport::SetGauss2DModel_sigma2_Cp1(Float_t Gauss2DModel_sigma2_Cp1){
108  fGauss2DModel_sigma2_Cp1 = Gauss2DModel_sigma2_Cp1;
109 }
110 //______________________________________________________________________________
111 //
112 void DIGTransport::SetGauss2DModel_weight(Float_t Gauss2DModel_weight){
113  fGauss2DModel_weight = Gauss2DModel_weight;
114 }
115 //______________________________________________________________________________
116 //
117 void DIGTransport::SetLorGaussModel_Norm1_Cp0(Float_t LorGaussModel_Norm1_Cp0){
118  fLorGaussModel_Norm1_Cp0=LorGaussModel_Norm1_Cp0 ;
119 }
120 //______________________________________________________________________________
121 //
122 void DIGTransport::SetLorGaussModel_Norm1_Cp1(Float_t LorGaussModel_Norm1_Cp1){
123  fLorGaussModel_Norm1_Cp1=LorGaussModel_Norm1_Cp1 ;
124 }
125 
126 //______________________________________________________________________________
127 //
128 void DIGTransport::SetLorGaussModel_Norm1_Cp2(Float_t LorGaussModel_Norm1_Cp2){
129  fLorGaussModel_Norm1_Cp2=LorGaussModel_Norm1_Cp2 ;
130 }
131 
132 //______________________________________________________________________________
133 //
134 void DIGTransport::SetLorGaussModel_sigma_Cp0(Float_t LorGaussModel_sigma_Cp0){
135  fLorGaussModel_sigma_Cp0=LorGaussModel_sigma_Cp0 ;
136 }
137 
138 //______________________________________________________________________________
139 //
140 void DIGTransport::SetLorGaussModel_sigma_Cp1(Float_t LorGaussModel_sigma_Cp1){
141  fLorGaussModel_sigma_Cp1=LorGaussModel_sigma_Cp1 ;
142 }
143 
144 //______________________________________________________________________________
145 //
146 void DIGTransport::SetLorGaussModel_C_Cp0(Float_t LorGaussModel_C_Cp0){
147  fLorGaussModel_C_Cp0=LorGaussModel_C_Cp0 ;
148 }
149 
150 //______________________________________________________________________________
151 //
152 void DIGTransport::SetLorGaussModel_C_Cp1(Float_t LorGaussModel_C_Cp1){
153  fLorGaussModel_C_Cp1=LorGaussModel_C_Cp1 ;
154 }
155 
156 //______________________________________________________________________________
157 //
158 void DIGTransport::SetLorGaussModel_Norm_Cp0(Float_t LorGaussModel_Norm_Cp0){
159  fLorGaussModel_Norm_Cp0=LorGaussModel_Norm_Cp0 ;
160 }
161 //______________________________________________________________________________
162 //
163  void DIGTransport::SetLorGaussModel_Norm_Cp1(Float_t LorGaussModel_Norm_Cp1){
164  fLorGaussModel_Norm_Cp1 = LorGaussModel_Norm_Cp1;
165 }
166 //______________________________________________________________________________
167 //
168 void DIGTransport::PrintInfo() {
169  std::cout<<"---------DIGTransport properties------------- "<<endl;
170  std::cout<<"fChargeModel"<<endl;
171  std::cout<<fChargeModel<<endl;
172  if(fChargeModel==1){
173  std::cout<<"fLorentz2DModel_Cp0"<<endl;
174  std::cout<<fLorentz2DModel_Cp0<<endl;
175  std::cout<<"fLorentz2DModel_Cp1"<<endl;
176  std::cout<<fLorentz2DModel_Cp1<<endl;
177  std::cout<<"fRangeLimit_InPitchUnit"<<endl;
178  std::cout<<fRangeLimit_InPitchUnit<<endl;
179  }else if(fChargeModel==2){
180 
181  std::cout<<"fGauss2DModel_sigma1_Cp0"<<endl;
182  std::cout<<fGauss2DModel_sigma1_Cp0 <<endl;
183  std::cout<<"fGauss2DModel_sigma1_Cp1"<<endl;
184  std::cout<< fGauss2DModel_sigma1_Cp1 <<endl;
185  std::cout<<"fGauss2DModel_sigma2_Cp0"<<endl;
186  std::cout<<fGauss2DModel_sigma2_Cp0 <<endl;
187  std::cout<<"fGauss2DModel_sigma2_Cp1"<<endl;
188  std::cout<< fGauss2DModel_sigma2_Cp1 <<endl;
189  std::cout<<"fGauss2DModel_weight"<<endl;
190  std::cout<< fGauss2DModel_weight <<endl;
191  }
192 
193 
194 
195  //std::cout<<""<<endl;
196  //std::cout<<<<endl;
197 }
198 
199 //______________________________________________________________________________
200 //
201 void DIGTransport::SetMyvar(Float_t Myvar){
202  fMyvar=Myvar;
203 }
204 //______________________________________________________________________________
205 //