FairRoot/PandaRoot
Functions | Variables
poormantracks.C File Reference

Go to the source code of this file.

Functions

TVector3 fVertex (0., 0., 0.)
 
int poormantracks (int nevt=250, int laut=0, int seed=-1)
 
TVector3 RollVertexBox (double widx, double widy, double widz)
 
void SmearVertex (TVector3 &vertex)
 
TVector3 RollMomentumBox (const TVector3 &vtx, double dtheta, double ptmin, double ptmax)
 
 SmearMomentum (TVector3 &momentum)
 
void PoorManTracks ()
 

Variables

double fSigVx =0.005
 
double fSigVy =0.005
 
double fSigVz =0.005
 
double fSigPx =0.01
 
double fSigPy =0.01
 
double fSigPz =0.01
 
double fDtheta = 1.
 
int fPID =211
 
bool fAlternate =true
 
TParticlePDG * fPDG = 0
 
int fNumTrk =6
 
int niter =1
 
int fSeedPrefix =936650
 
int fVerbose =0
 
TClonesArray * fMcCands =0
 
TClonesArray * fCands =0
 

Function Documentation

TVector3 fVertex ( 0.  ,
0.  ,
0.   
)

Referenced by poormantracks(), and PoorManTracks().

int poormantracks ( int  nevt = 250,
int  laut = 0,
int  seed = -1 
)

Definition at line 29 of file poormantracks.C.

References Bool_t, RhoCandidate::Charge(), RhoCandidate::Combine(), RhoCandidate::Cov7(), RhoCandidate::Daughter(), fCands, RhoFitterBase::Fit(), RhoKalmanVtxFitter::FitVertexFast(), fMcCands, RhoCalculationTools::ForceConstantBz(), fSeedPrefix, fVerbose, fVertex(), RhoFitterBase::GetChi2(), RhoCandidate::GetFit(), RhoCandidate::GetMcTruth(), RhoFitterBase::GetNdf(), RhoVtxPoca::GetPocaVtx(), RhoCandidate::GetPosition(), RhoFitterBase::GetProb(), i, LoadPandaStyle(), max(), min(), mom, mommc, RhoCandidate::NDaughters(), RhoCandidate::P4(), RhoCandidate::P4Cov(), PoorManTracks(), RhoCandidate::Pos(), printf(), RollVertexBox(), seed, RhoKalmanVtxFitter::SetSilent(), SetVerbose(), sqrt(), and TString.

30 {
31  //gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
32  //gSystem->Load("libRho");
33  //gSystem->Load("libAnalysisTools");
34  fVerbose=laut;
35  gROOT->LoadMacro("$VMCWORKDIR/macro/run/Tools.C");
36  LoadPandaStyle();
37  int ranseed=fSeedPrefix+seed;
38  gRandom->SetSeed(ranseed);
39  RhoCalculationTools::ForceConstantBz(20.); // [kGs]=[0.1T]
40  PndAnalysisCalcTools::SetVerbose(0);
41  fMcCands = new TClonesArray("RhoCandidate");
42  fCands = new TClonesArray("RhoCandidate");
43 
44  //PArameters
45  double widx=1.5,
46  widy=1.5,
47  widz=1.5; //[cm] maximum vertex position off origin
48 
49  int bins=100;
50  double min=0,max=4.5,mdf=2.0; // GeV/c
51  double min2=0.13,max2=0.15;
52  double min3=1.,max3=1.5;
53  double max4=0.3;//cm
54  double maxchi = 5;
55  double vrange = 0.05; // cm
56  double pullrange =4.;
57  double momrange=0.1;//GeV/c
58  double angdiffran=0.05; //rad
59 
60  TList* histolist = new TList();
61 
62  // NO FIT
63  TH1F *hVtxDiffX=new TH1F("hvtxDiffx","Vertex Smearing;(x_{mc}-x)/cm",bins,-vrange,vrange);
64  TH1F *hVtxDiffY=new TH1F("hvtxDiffy","Vertex Smearing;(y_{mc}-y)/cm",bins,-vrange,vrange);
65  TH1F *hVtxDiffZ=new TH1F("hvtxDiffz","Vertex Smearing;(z_{mc}-z)/cm",bins,-vrange,vrange);
66 
67  TH1F *hVtxDiffPX=new TH1F("hvtxdiffpx","Momentum reco;(p_{x,mc}-p_{x} / GeV/c)",bins,-momrange,momrange);
68  TH1F *hVtxDiffPY=new TH1F("hvtxdiffpy","Momentum reco;(p_{y,mc}-p_{y} / GeV/c)",bins,-momrange,momrange);
69  TH1F *hVtxDiffPZ=new TH1F("hvtxdiffpz","Momentum reco;(p_{z,mc}-p_{z}) / GeV/c",bins,-momrange,momrange);
70  TH1F *hVtxDiffE=new TH1F("hvtxdiffe","Energy reco;(E_{mc}-E) / GeV",bins,-momrange,momrange);
71 
72  TH1F *hVtxPullPX=new TH1F("hvtxpullpx","Momentum pull distribution reco;(p_{x,mc}-p_{x})/#sigma_{p_{x}}",bins,-pullrange,pullrange);
73  TH1F *hVtxPullPY=new TH1F("hvtxpullpy","Momentum pull distribution reco;(p_{y,mc}-p_{y})/#sigma_{p_{y}}",bins,-pullrange,pullrange);
74  TH1F *hVtxPullPZ=new TH1F("hvtxpullpz","Momentum pull distribution reco;(p_{z,mc}-p_{z})/#sigma_{p_{z}}",bins,-pullrange,pullrange);
75  TH1F *hVtxPullE=new TH1F("hvtxpulle","Energy pull distribution reco;(E_{mc}-E)/#sigma_{E}",bins,-pullrange,pullrange);
76 
77  TH1F *hVtxDiffThe=new TH1F("hvtxdiffthe","Momentum reco;(#Theta_{mc}-#Theta}) / GeV/c",bins,-angdiffran,angdiffran);
78  TH1F *hVtxDiffPhi=new TH1F("hvtxdiffphi","Momentum reco;(#Phi_{mc}-#Phi}) / GeV/c",bins,-angdiffran,angdiffran);
79 
80  // POCA
81 
82  TH1F *hVtxPocaX=new TH1F("hvtxPocax","Vertex POCA;(x_{mc}-x)/cm",bins,-vrange,vrange);
83  TH1F *hVtxPocaY=new TH1F("hvtxPocay","Vertex POCA;(y_{mc}-y)/cm",bins,-vrange,vrange);
84  TH1F *hVtxPocaZ=new TH1F("hvtxPocaz","Vertex POCA;(z_{mc}-z)/cm",bins,-vrange,vrange);
85 
86  TH2F *hVtxPocaXY=new TH2F("hvtxPocaxy","Vertex POCA X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
87  TH2F *hVtxPocaRZ=new TH2F("hvtxPocarz","Vertex POCA R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
88 
89  TH1F *hVtxPullPocaX=new TH1F("hvtxpullpocax","Vertex pull distribution POCA;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
90  TH1F *hVtxPullPocaY=new TH1F("hvtxpullpocay","Vertex pull distribution POCA;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
91  TH1F *hVtxPullPocaZ=new TH1F("hvtxpullpocaz","Vertex pull distribution POCA;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);
92 
93  TH1F *hVtxPocas=new TH1F("hvtxpocas","Vertex POCA D value;D/cm",bins,0.,vrange);
94  TH1F *hVtxPocaEmpty=new TH1F("hvtxcpocaempty","",0,0.,1.);
95 
96  // FAST FIT
97 
98  TH1F *hVtxFastX=new TH1F("hvtxFastx","Vertex fast fit;(x_{mc}-x)/cm",bins,-vrange,vrange);
99  TH1F *hVtxFastY=new TH1F("hvtxTasty","Vertex fast fit;(y_{mc}-y)/cm",bins,-vrange,vrange);
100  TH1F *hVtxFastZ=new TH1F("hvtxFastz","Vertex fast fit;(z_{mc}-z)/cm",bins,-vrange,vrange);
101 
102  TH2F *hVtxFastXY=new TH2F("hvtxfastxy","Vertex fast fit X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
103  TH2F *hVtxFastRZ=new TH2F("hvtxfastrz","Vertex fast fit R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
104 
105  TH1F *hVtxErrFastX=new TH1F("hvtxerrfastx","Vertex error distribution fast;#sigma_{x}",bins,0.,vrange);
106  TH1F *hVtxErrFastY=new TH1F("hvtxerrfasty","Vertex error distribution fast;#sigma_{y}",bins,0.,vrange);
107  TH1F *hVtxErrFastZ=new TH1F("hvtxerrfastz","Vertex error distribution fast;#sigma_{z}",bins,0.,vrange);
108  TH1F *hVtxPullFastX=new TH1F("hvtxpullfastx","Vertex pull distribution fast;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
109  TH1F *hVtxPullFastY=new TH1F("hvtxpullfasty","Vertex pull distribution fast;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
110  TH1F *hVtxPullFastZ=new TH1F("hvtxpullfastz","Vertex pull distribution fast;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);
111 
112  TH1F *hVtxChi2Fast=new TH1F("hvtxchi2fast","Vertex Fast #Chi^{2};#Chi^{2}",bins,0,maxchi);
113  TH1F *hVtxChiProbFast=new TH1F("hvtxchipropfast","P(#Chi^{2}) Fast;p",bins,0.,1.);
114 
115  // FULL FIT
116  TH1F *hVtxFitX=new TH1F("hvtxFitx","Vertex full fit;(x_{mc}-x)/cm",bins,-vrange,vrange);
117  TH1F *hVtxFitY=new TH1F("hvtxFity","Vertex full fit;(y_{mc}-y)/cm",bins,-vrange,vrange);
118  TH1F *hVtxFitZ=new TH1F("hvtxFitz","Vertex full fit;(z_{mc}-z)/cm",bins,-vrange,vrange);
119 
120  TH2F *hVtxFitXY=new TH2F("hvtxfitxy","Vertex full fit X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
121  TH2F *hVtxFitRZ=new TH2F("hvtxfitrz","Vertex full fit R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
122 
123  TH1F *hVtxErrFitX=new TH1F("hvtxerrfitx","Vertex error distribution full fit;#sigma_{x}",bins,0.,vrange);
124  TH1F *hVtxErrFitY=new TH1F("hvtxerrfity","Vertex error distribution full fit;#sigma_{y}",bins,0.,vrange);
125  TH1F *hVtxErrFitZ=new TH1F("hvtxerrfitz","Vertex error distribution full fit;#sigma_{z}",bins,0.,vrange);
126  TH1F *hVtxPullFitX=new TH1F("hvtxpullfitx","Vertex pull distribution full fit;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
127  TH1F *hVtxPullFitY=new TH1F("hvtxpullfity","Vertex pull distribution full fit;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
128  TH1F *hVtxPullFitZ=new TH1F("hvtxpullfitz","Vertex pull distribution full fit;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);
129  TH1F *hVtxDiffFitPX=new TH1F("hvtxdifffitpx","Momentum full fit;(p_{x,mc}-p_{x} / GeV/c)",bins,-momrange,momrange);
130  TH1F *hVtxDiffFitPY=new TH1F("hvtxdifffitpy","Momentum full fit;(p_{y,mc}-p_{y} / GeV/c)",bins,-momrange,momrange);
131  TH1F *hVtxDiffFitPZ=new TH1F("hvtxdifffitpz","Momentum full fit;(p_{z,mc}-p_{z}) / GeV/c",bins,-momrange,momrange);
132  TH1F *hVtxDiffFitE=new TH1F("hvtxdifffite","Energy full fit;(E_{mc}-E) / GeV",bins,-momrange,momrange);
133  TH1F *hVtxDiffFitThe=new TH1F("hvtxdifffitthe","Momentum full fit;(#Theta_{mc}-#Theta}) / GeV/c",bins,-angdiffran,angdiffran);
134  TH1F *hVtxDiffFitPhi=new TH1F("hvtxdifffitphi","Momentum full fit;(#Phi_{mc}-#Phi}) / GeV/c",bins,-angdiffran,angdiffran);
135 
136  TH1F *hVtxPullFitPX=new TH1F("hvtxpullfitpx","Momentum pull distribution full fit;(p_{x,mc}-p_{x})/#sigma_{p_{x}}",bins,-pullrange,pullrange);
137  TH1F *hVtxPullFitPY=new TH1F("hvtxpullfitpy","Momentum pull distribution full fit;(p_{y,mc}-p_{y})/#sigma_{p_{y}}",bins,-pullrange,pullrange);
138  TH1F *hVtxPullFitPZ=new TH1F("hvtxpullfitpz","Momentum pull distribution full fit;(p_{z,mc}-p_{z})/#sigma_{p_{z}}",bins,-pullrange,pullrange);
139  TH1F *hVtxPullFitE=new TH1F("hvtxpullfite","Energy pull distribution full fit;(E_{mc}-E)/#sigma_{E}",bins,-pullrange,pullrange);
140 
141  TH1F *hVtxChi2Fit=new TH1F("hvtxchi2fit","Vertex Fit #Chi^{2};#Chi^{2}",bins,0,maxchi);
142  TH1F *hVtxChiProbFit=new TH1F("hvtxchipropfit","P(#Chi^{2}) Fit;p",bins,0.,1.);
143 
144 
145  // KIN Fit
146 
147  TH1F *hVtxKinX=new TH1F("hvtxKinx","Vertex full kin;(x_{mc}-x)/cm",bins,-vrange,vrange);
148  TH1F *hVtxKinY=new TH1F("hvtxKiny","Vertex full kin;(y_{mc}-y)/cm",bins,-vrange,vrange);
149  TH1F *hVtxKinZ=new TH1F("hvtxKinz","Vertex full kin;(z_{mc}-z)/cm",bins,-vrange,vrange);
150 
151  TH2F *hVtxKinXY=new TH2F("hvtxkinxy","Vertex full kin X-Y view;(x_{mc}-x)/#mum;(y_{mc}-y)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
152  TH2F *hVtxKinRZ=new TH2F("hvtxkinrz","Vertex full kin R-Z view;(z_{mc}-z)/#mum;(r_{mc}-r)/cm",bins,-vrange,vrange,bins,-vrange,vrange);
153 
154  TH1F *hVtxErrKinX=new TH1F("hvtxerrkinx","Vertex error distribution full kin;#sigma_{x}",bins,0.,vrange);
155  TH1F *hVtxErrKinY=new TH1F("hvtxerrkiny","Vertex error distribution full kin;#sigma_{y}",bins,0.,vrange);
156  TH1F *hVtxErrKinZ=new TH1F("hvtxerrkinz","Vertex error distribution full kin;#sigma_{z}",bins,0.,vrange);
157  TH1F *hVtxPullKinX=new TH1F("hvtxpullkinx","Vertex pull distribution full kin;(x_{mc}-x)/#sigma_{x}",bins,-pullrange,pullrange);
158  TH1F *hVtxPullKinY=new TH1F("hvtxpullkiny","Vertex pull distribution full kin;(y_{mc}-y)/#sigma_{y}",bins,-pullrange,pullrange);
159  TH1F *hVtxPullKinZ=new TH1F("hvtxpullkinz","Vertex pull distribution full kin;(z_{mc}-z)/#sigma_{z}",bins,-pullrange,pullrange);
160  TH1F *hVtxDiffKinPX=new TH1F("hvtxdiffkinpx","Momentum full kin;(p_{x,mc}-p_{x} / GeV/c)",bins,-momrange,momrange);
161  TH1F *hVtxDiffKinPY=new TH1F("hvtxdiffkinpy","Momentum full kin;(p_{y,mc}-p_{y} / GeV/c)",bins,-momrange,momrange);
162  TH1F *hVtxDiffKinPZ=new TH1F("hvtxdiffkinpz","Momentum full kin;(p_{z,mc}-p_{z}) / GeV/c",bins,-momrange,momrange);
163  TH1F *hVtxDiffKinE=new TH1F("hvtxdiffkine","Energy full kin;(E_{mc}-E) / GeV",bins,-momrange,momrange);
164  TH1F *hVtxDiffKinThe=new TH1F("hvtxdiffkinthe","Momentum full kin;(#Theta_{mc}-#Theta}) / GeV/c",bins,-angdiffran,angdiffran);
165  TH1F *hVtxDiffKinPhi=new TH1F("hvtxdiffkinphi","Momentum full kin;(#Phi_{mc}-#Phi}) / GeV/c",bins,-angdiffran,angdiffran);
166 
167  TH1F *hVtxPullKinPX=new TH1F("hvtxpullkinpx","Momentum pull distribution full kin;(p_{x,mc}-p_{x})/#sigma_{p_{x}}",bins,-pullrange,pullrange);
168  TH1F *hVtxPullKinPY=new TH1F("hvtxpullkinpy","Momentum pull distribution full kin;(p_{y,mc}-p_{y})/#sigma_{p_{y}}",bins,-pullrange,pullrange);
169  TH1F *hVtxPullKinPZ=new TH1F("hvtxpullkinpz","Momentum pull distribution full kin;(p_{z,mc}-p_{z})/#sigma_{p_{z}}",bins,-pullrange,pullrange);
170  TH1F *hVtxPullKinE=new TH1F("hvtxpullkine","Energy pull distribution full kin;(E_{mc}-E)/#sigma_{E}",bins,-pullrange,pullrange);
171 
172  TH1F *hVtxChi2Kin=new TH1F("hvtxchi2kin","Vertex Kin #Chi^{2};#Chi^{2}",bins,0,maxchi);
173  TH1F *hVtxChiProbKin=new TH1F("hvtxchipropkin","P(#Chi^{2}) Kin;p",bins,0.,1.);
174 
175  // Helix Params
176  TH1F *hPrgPull0 = new TH1F("hprgpull0","PRG Pull Distribution (0)",bins,-pullrange,pullrange);
177  TH1F *hPrgPull1 = new TH1F("hprgpull1","PRG Pull Distribution (1)",bins,-pullrange,pullrange);
178  TH1F *hPrgPull2 = new TH1F("hprgpull2","PRG Pull Distribution (2)",bins,-pullrange,pullrange);
179  TH1F *hPrgPull3 = new TH1F("hprgpull3","PRG Pull Distribution (3)",bins,-pullrange,pullrange);
180  TH1F *hPrgPull4 = new TH1F("hprgpull4","PRG Pull Distribution (4)",bins,-pullrange,pullrange);
181 
182  histolist->Add(hVtxDiffX);
183  histolist->Add(hVtxDiffY);
184  histolist->Add(hVtxDiffZ);
185 
186  histolist->Add(hVtxPocaX);
187  histolist->Add(hVtxPocaY);
188  histolist->Add(hVtxPocaZ);
189 
190  histolist->Add(hVtxFastX);
191  histolist->Add(hVtxFastY);
192  histolist->Add(hVtxFastZ);
193 
194  histolist->Add(hVtxFitX);
195  histolist->Add(hVtxFitY);
196  histolist->Add(hVtxFitZ);
197 
198  histolist->Add(hVtxKinX);
199  histolist->Add(hVtxKinY);
200  histolist->Add(hVtxKinZ);
201 
202  histolist->Add(hVtxErrFastX);
203  histolist->Add(hVtxErrFastY);
204  histolist->Add(hVtxErrFastZ);
205 
206  histolist->Add(hVtxErrFitX);
207  histolist->Add(hVtxErrFitY);
208  histolist->Add(hVtxErrFitZ);
209 
210  histolist->Add(hVtxErrKinX);
211  histolist->Add(hVtxErrKinY);
212  histolist->Add(hVtxErrKinZ);
213 
214  histolist->Add(hVtxPullPocaX);
215  histolist->Add(hVtxPullPocaY);
216  histolist->Add(hVtxPullPocaZ);
217 
218  histolist->Add(hVtxPullFastX);
219  histolist->Add(hVtxPullFastY);
220  histolist->Add(hVtxPullFastZ);
221 
222  histolist->Add(hVtxPullFitX);
223  histolist->Add(hVtxPullFitY);
224  histolist->Add(hVtxPullFitZ);
225 
226  histolist->Add(hVtxPullKinX);
227  histolist->Add(hVtxPullKinY);
228  histolist->Add(hVtxPullKinZ);
229 
230  histolist->Add(hVtxDiffPX);
231  histolist->Add(hVtxDiffPY);
232  histolist->Add(hVtxDiffPZ);
233 
234  histolist->Add(hVtxPullPX);
235  histolist->Add(hVtxPullPY);
236  histolist->Add(hVtxPullPZ);
237 
238  histolist->Add(hVtxDiffFitPX);
239  histolist->Add(hVtxDiffFitPY);
240  histolist->Add(hVtxDiffFitPZ);
241 
242  histolist->Add(hVtxPullFitPX);
243  histolist->Add(hVtxPullFitPY);
244  histolist->Add(hVtxPullFitPZ);
245 
246  histolist->Add(hVtxDiffKinPX);
247  histolist->Add(hVtxDiffKinPY);
248  histolist->Add(hVtxDiffKinPZ);
249 
250  histolist->Add(hVtxPullKinPX);
251  histolist->Add(hVtxPullKinPY);
252  histolist->Add(hVtxPullKinPZ);
253 
254  histolist->Add(hVtxDiffE);
255  histolist->Add(hVtxDiffFitE);
256  histolist->Add(hVtxDiffKinE);
257  histolist->Add(hVtxPullE);
258  histolist->Add(hVtxPullFitE);
259  histolist->Add(hVtxPullKinE);
260 
261  histolist->Add(hVtxDiffThe);
262  histolist->Add(hVtxDiffPhi);
263 
264  histolist->Add(hVtxDiffFitThe);
265  histolist->Add(hVtxDiffFitPhi);
266 
267  histolist->Add(hVtxDiffKinThe);
268  histolist->Add(hVtxDiffKinPhi);
269 
270  histolist->Add(hVtxPocas);
271  histolist->Add(hVtxChi2Fast);
272  histolist->Add(hVtxChi2Fit);
273  histolist->Add(hVtxChi2Kin);
274  histolist->Add(hVtxChiProbFast);
275  histolist->Add(hVtxChiProbFit);
276  histolist->Add(hVtxChiProbKin);
277 
278  histolist->Add(hVtxPocaXY);
279  histolist->Add(hVtxFastXY);
280  histolist->Add(hVtxFitXY);
281  histolist->Add(hVtxKinXY);
282  histolist->Add(hVtxPocaRZ);
283  histolist->Add(hVtxFastRZ);
284  histolist->Add(hVtxFitRZ);
285  histolist->Add(hVtxKinRZ);
286 
287  histolist->Add(hPrgPull0);
288  histolist->Add(hPrgPull1);
289  histolist->Add(hPrgPull2);
290  histolist->Add(hPrgPull3);
291  histolist->Add(hPrgPull4);
292 
293  RhoCandidate *aCand=0,*bCand=0,*cCand=0,*dCand=0;
294  RhoCandidate *tmpcand=0,*mccand=0;
295  TVector3 vertexMC,vertexPoc,vertexFast,vertexFit,vertexKin,vertexDiff,vertexDiffFast,vertexDiffFit,vertexDiffKin;
296  TLorentzVector mom,mommc,momdiff;
297  const TVector3 nullpunkt(0.,0.,0.);
298  double rnd1=0.,rnd2=0.,rnd3=0.;
299  RhoError amomcov(4);
300  TMatrixD aposcov(3,3);
301  TMatrixD helixCov(5,5);
302  TMatrixD jacobian(5,7);
303  Float_t mcparams[5];
304  Float_t helixparams[5];
305  TMatrixD covV(3,3);
306  TMatrixD covVfit(3,3);
307 
308  for(int i=0;i<nevt;i++)
309  {
310  printf("Event %i of %i. (seed=%i)\n",i,nevt,ranseed);
311  RhoCandidate *combiCand=0,*combiCand1=0,*combiCand2=0;
312  vertexMC = RollVertexBox(widx,widy,widz);
313  fVertex=vertexMC;
314 
315  PoorManTracks();
316  if (fCands->GetEntriesFast()<2) {
317  cout<<"Warning: too few particles created."<<endl;
318  continue;
319  }
320 
321  aCand=(RhoCandidate*)fCands->At(0);
322  bCand=(RhoCandidate*)fCands->At(1);
323  if(fCands->GetEntriesFast()==3)
324  {
325  cCand=(RhoCandidate*)fCands->At(2);
326  combiCand = aCand->Combine(bCand,cCand);
327  combiCand1 = aCand->Combine(bCand,cCand);
328  combiCand2 = aCand->Combine(bCand,cCand);
329  }else if(fCands->GetEntriesFast()==4)
330  {
331  cCand=(RhoCandidate*)fCands->At(2);
332  dCand=(RhoCandidate*)fCands->At(3);
333  combiCand = aCand->Combine(bCand,cCand,dCand);
334  combiCand1 = aCand->Combine(bCand,cCand,dCand);
335  combiCand2 = aCand->Combine(bCand,cCand,dCand);
336  }else{
337  combiCand = aCand->Combine(bCand);
338  combiCand1 = aCand->Combine(bCand);
339  combiCand2 = aCand->Combine(bCand);
340  }
341 
342  if(fVerbose>1) cout<<"aCand, bCand, comboCand, aPos, bPos, combipos"<<endl;
343  if(fVerbose>1) cout<<*aCand<<endl;
344  if(fVerbose>1) cout<<*bCand<<endl;
345  if(fVerbose>1) cout<<*combiCand<<endl;
346  if(fVerbose>1) aCand->GetPosition().Print();
347  if(fVerbose>1) bCand->GetPosition().Print();
348  if(fVerbose>1) combiCand->GetPosition().Print();
349 
350 
351  // CLEAN
352  for(int k=0;k<combiCand->NDaughters();k++)
353  {
354  tmpcand = combiCand->Daughter(k);
355  mccand = tmpcand->GetMcTruth();
356 
357  vertexDiff=mccand->Pos();vertexDiff-=tmpcand->Pos();
358  hVtxDiffX->Fill(vertexDiff.X());
359  hVtxDiffY->Fill(vertexDiff.Y());
360  hVtxDiffZ->Fill(vertexDiff.Z());
361 
362  mommc=mccand->P4();
363  mom=tmpcand->P4();
364  momdiff=mommc-mom;
365  amomcov=tmpcand->P4Cov();
366  hVtxDiffPX->Fill(momdiff.X());
367  hVtxDiffPY->Fill(momdiff.Y());
368  hVtxDiffPZ->Fill(momdiff.Z());
369  hVtxDiffE->Fill(momdiff.E());
370 
371  hVtxPullPX->Fill(momdiff.X()/sqrt(amomcov[0][0]));
372  hVtxPullPY->Fill(momdiff.Y()/sqrt(amomcov[1][1]));
373  hVtxPullPZ->Fill(momdiff.Z()/sqrt(amomcov[2][2]));
374  hVtxPullE->Fill(momdiff.E()/sqrt(amomcov[3][3]));
375 
376  hVtxDiffThe->Fill(mommc.Theta()-mom.Theta());
377  hVtxDiffPhi->Fill(mommc.Phi()-mom.Phi());
378 
379  // Now check the PRG description....
380  Bool_t prgcheckmc = PndAnalysisCalcTools::P7toPRG(mccand->Pos(),mccand->P4(),mccand->Charge(),mccand->Cov7(),nullpunkt, mcparams, helixCov, jacobian, kFALSE);
381  Bool_t prgchecktc = PndAnalysisCalcTools::P7toPRG(tmpcand->Pos(),tmpcand->P4(),tmpcand->Charge(),tmpcand->Cov7(),nullpunkt, helixparams, helixCov, jacobian, kFALSE);
382 // Bool_t prgcheckmc = PndAnalysisCalcTools::P7toHelix(mccand->Pos(),mccand->P4(),mccand->Charge(),mccand->Cov7(), mcparams, helixCov, kFALSE);
383 // Bool_t prgchecktc = PndAnalysisCalcTools::P7toHelix(tmpcand->Pos(),tmpcand->P4(),tmpcand->Charge(),tmpcand->Cov7(), helixparams, helixCov, kFALSE);
384  hPrgPull0->Fill((mcparams[0]-helixparams[0])/sqrt(helixCov[0][0]));
385  hPrgPull1->Fill((mcparams[1]-helixparams[1])/sqrt(helixCov[1][1]));
386  hPrgPull2->Fill((mcparams[2]-helixparams[2])/sqrt(helixCov[2][2]));
387  hPrgPull3->Fill((mcparams[3]-helixparams[3])/sqrt(helixCov[3][3]));
388  hPrgPull4->Fill((mcparams[4]-helixparams[4])/sqrt(helixCov[4][4]));
389  }
390 
391  // ********** POCA
392  RhoVtxPoca vPoca;
393  vertexPoc=vertexMC; //seed
394  if(fVerbose>0) cout<<"poca "<<flush;
395  double dist = vPoca.GetPocaVtx(vertexPoc,combiCand);
396  if(fVerbose>0) cout<<" - Dist = "<<dist<<" \t";if(fVerbose>1)cout<<endl;
397  vertexDiff=vertexMC - vertexPoc;
398  if(fVerbose>1) vertexMC.Print();
399  if(fVerbose>1) vertexPoc.Print();
400  if(fVerbose>1) vertexDiff.Print();
401 
402  hVtxPocaX->Fill(vertexDiff.X());
403  hVtxPocaY->Fill(vertexDiff.Y());
404  hVtxPocaZ->Fill(vertexDiff.Z());
405  hVtxPocaXY->Fill(vertexDiff.X(),vertexDiff.Y());
406  hVtxPocaRZ->Fill(vertexDiff.Z(),vertexMC.Perp()-vertexPoc.Perp());
407  hVtxPocas->Fill(dist);
408  hVtxPullPocaX->Fill(vertexDiff.X()/dist);
409  hVtxPullPocaY->Fill(vertexDiff.Y()/dist);
410  hVtxPullPocaZ->Fill(vertexDiff.Z()/dist);
411 
412  // ********** FAST FIT
413  PndKalmanVtxFitter vFastter(combiCand);
414  vFastter.SetSilent();
415  //vFastter.SetDebug();
416  if(fVerbose>0) cout<<"prg "<<flush;
417  vertexFast=nullpunkt; //seed
418  double checkfast = vFastter.FitVertexFast(vertexFast,covV);
419  double chiq=vFastter.GetChi2();
420  int ndf = vFastter.GetNdf();
421  double prob = vFastter.GetProb();
422  if(fVerbose>0) cout<<" - chiq = "<<chiq<<" \t";if(fVerbose>1)cout<<endl;
423  if(checkfast>0.)
424  {
425  hVtxChi2Fast->Fill(chiq);
426  hVtxChiProbFast->Fill(prob);
427 
428  vertexDiffFast=vertexMC - vertexFast;
429  if(fVerbose>1) vertexMC.Print();
430  if(fVerbose>1) vertexFast.Print();
431  if(fVerbose>0) vertexDiffFast.Print();
432 
433  hVtxFastX->Fill(vertexDiffFast.X());
434  hVtxFastY->Fill(vertexDiffFast.Y());
435  hVtxFastZ->Fill(vertexDiffFast.Z());
436  hVtxFastXY->Fill(vertexDiffFast.X(),vertexDiffFast.Y());
437  hVtxFastRZ->Fill(vertexDiffFast.Z(),vertexMC.Perp()-vertexFast.Perp());
438 
439  hVtxErrFastX->Fill(sqrt(covV[0][0]));
440  hVtxErrFastY->Fill(sqrt(covV[1][1]));
441  hVtxErrFastZ->Fill(sqrt(covV[2][2]));
442  hVtxPullFastX->Fill(vertexDiffFast.X()/sqrt(covV[0][0]));
443  hVtxPullFastY->Fill(vertexDiffFast.Y()/sqrt(covV[1][1]));
444  hVtxPullFastZ->Fill(vertexDiffFast.Z()/sqrt(covV[2][2]));
445  }
446  // ********** FULL FIT
447  PndKalmanVtxFitter vFitter(combiCand1);
448  vFitter.SetSilent();
449  //PndAnalysisCalcTools::SetVerbose(5);
450  //if(laut)vFitter.SetDebug(true);
451  vertexFit=nullpunkt; //seed
452  if(fVerbose>0) cout<<"prg "<<flush;
453  //if(laut>0) PndAnalysisCalcTools::SetVerbose(3);
454  bool check = vFitter.Fit();
455  double chiqfit = vFitter.GetChi2();
456  int ndffit = vFitter.GetNdf();
457  double probfit = vFitter.GetProb();
458  //if(laut>0) PndAnalysisCalcTools::SetVerbose(0);
459  if(fVerbose>0) cout<<" - chiq = "<<chiqfit<<" \t"<<flush;if(fVerbose>1)cout<<endl;
460  if (check)
461  {
462  //PndAnalysisCalcTools::SetVerbose(0);
463  hVtxChi2Fit->Fill(chiqfit);
464  hVtxChiProbFit->Fill(probfit);
465 
466  vertexFit=combiCand1->GetFit()->DecayVtx();
467  covVfit=combiCand1->GetFit()->DecayVtx().CovMatrix();
468  vertexDiffFit=vertexMC - vertexFit;
469  if(fVerbose>1) vertexMC.Print();
470  if(fVerbose>1) vertexFit.Print();
471  if(fVerbose>1) vertexDiffFit.Print();
472 
473  hVtxFitX->Fill(vertexDiffFit.X());
474  hVtxFitY->Fill(vertexDiffFit.Y());
475  hVtxFitZ->Fill(vertexDiffFit.Z());
476  hVtxFitXY->Fill(vertexDiffFit.X(),vertexDiffFit.Y());
477  hVtxFitRZ->Fill(vertexDiffFit.Z(),vertexMC.Perp()-vertexFit.Perp());
478 
479  hVtxErrFitX->Fill(sqrt(covVfit[0][0]));
480  hVtxErrFitY->Fill(sqrt(covVfit[1][1]));
481  hVtxErrFitZ->Fill(sqrt(covVfit[2][2]));
482  hVtxPullFitX->Fill(vertexDiffFit.X()/sqrt(covVfit[0][0]));
483  hVtxPullFitY->Fill(vertexDiffFit.Y()/sqrt(covVfit[1][1]));
484  hVtxPullFitZ->Fill(vertexDiffFit.Z()/sqrt(covVfit[2][2]));
485 
486  for(int k=0;k<combiCand1->NDaughters();k++)
487  {
488  tmpcand = combiCand1->Daughter(k);
489  mccand = tmpcand->GetMcTruth();
490  tmpcand = tmpcand->GetFit();
491  mommc=mccand->P4();
492  mom=tmpcand->P4();
493  momdiff=mommc-mom;
494  amomcov=tmpcand->P4Cov();
495  hVtxDiffFitPX->Fill(momdiff.X());
496  hVtxDiffFitPY->Fill(momdiff.Y());
497  hVtxDiffFitPZ->Fill(momdiff.Z());
498  hVtxDiffFitE->Fill(momdiff.E());
499 
500  hVtxPullFitPX->Fill(momdiff.X()/sqrt(amomcov[0][0]));
501  hVtxPullFitPY->Fill(momdiff.Y()/sqrt(amomcov[1][1]));
502  hVtxPullFitPZ->Fill(momdiff.Z()/sqrt(amomcov[2][2]));
503  hVtxPullFitE->Fill(momdiff.E()/sqrt(amomcov[3][3]));
504 
505  hVtxDiffFitThe->Fill(mommc.Theta()-mom.Theta());
506  hVtxDiffFitPhi->Fill(mommc.Phi()-mom.Phi());
507  }
508  }
509 
510  // ********** Kin FIT
511  RhoKinVtxFitter kFitter(combiCand2);
512  vertexFit=nullpunkt; //seed
513  if(fVerbose>0) cout<<"kin "<<flush;
514  bool test = kFitter.Fit();
515  double chiqkin = kFitter.GetChi2();
516  int ndfkin = kFitter.GetNdf();
517  double probkin = kFitter.GetProb();
518  //if(laut>0) PndAnalysisCalcTools::SetVerbose(0);
519  if(fVerbose>0) cout<<" - chiq = "<<chiqkin<<" \t"<<flush;if(fVerbose>1)cout<<endl;
520  if (test)
521  {
522  //PndAnalysisCalcTools::SetVerbose(0);
523  vertexKin=combiCand2->GetFit()->DecayVtx();
524  aposcov=combiCand2->GetFit()->DecayVtx().CovMatrix();
525  hVtxChi2Kin->Fill(chiqkin);
526  hVtxChiProbKin->Fill(probkin);
527 
528  vertexDiffKin=vertexMC - vertexKin;
529  if(fVerbose>1) vertexMC.Print();
530  if(fVerbose>1) vertexKin.Print();
531  if(fVerbose>0) vertexDiffKin.Print();
532 
533  hVtxKinX->Fill(vertexDiffKin.X());
534  hVtxKinY->Fill(vertexDiffKin.Y());
535  hVtxKinZ->Fill(vertexDiffKin.Z());
536  hVtxKinXY->Fill(vertexDiffKin.X(),vertexDiffKin.Y());
537  hVtxKinRZ->Fill(vertexDiffKin.Z(),vertexMC.Perp()-vertexKin.Perp());
538 
539  hVtxErrKinX->Fill(sqrt(aposcov[0][0]));
540  hVtxErrKinY->Fill(sqrt(aposcov[1][1]));
541  hVtxErrKinZ->Fill(sqrt(aposcov[2][2]));
542  hVtxPullKinX->Fill(vertexDiffKin.X()/sqrt(aposcov[0][0]));
543  hVtxPullKinY->Fill(vertexDiffKin.Y()/sqrt(aposcov[1][1]));
544  hVtxPullKinZ->Fill(vertexDiffKin.Z()/sqrt(aposcov[2][2]));
545 
546  for(int k=0;k<combiCand2->NDaughters();k++)
547  {
548  tmpcand = combiCand2->Daughter(k);
549  mccand = tmpcand->GetMcTruth();
550  tmpcand = tmpcand->GetFit();
551  mommc=mccand->P4();
552  mom=tmpcand->P4();
553  momdiff=mommc-mom;
554  amomcov=tmpcand->P4Cov();
555  hVtxDiffKinPX->Fill(momdiff.X());
556  hVtxDiffKinPY->Fill(momdiff.Y());
557  hVtxDiffKinPZ->Fill(momdiff.Z());
558  hVtxDiffKinE->Fill(momdiff.E());
559 
560  hVtxPullKinPX->Fill(momdiff.X()/sqrt(amomcov[0][0]));
561  hVtxPullKinPY->Fill(momdiff.Y()/sqrt(amomcov[1][1]));
562  hVtxPullKinPZ->Fill(momdiff.Z()/sqrt(amomcov[2][2]));
563  hVtxPullKinE->Fill(momdiff.E()/sqrt(amomcov[3][3]));
564 
565  hVtxDiffKinThe->Fill(mommc.Theta()-mom.Theta());
566  hVtxDiffKinPhi->Fill(mommc.Phi()-mom.Phi());
567  }
568  }
569  } // event loop
570 
571  cout<<" ---------- calculations done, writing histograms."<<endl;
572  //manipulate & fit histos now
573 
574  //Write and Plot all histos from List
575 
576  TString hfile="Data/HistoVertexing";
577  if(seed>=0){hfile+="-";hfile+=ranseed;}
578  hfile+=".root";
579  TFile* histofile = new TFile(hfile.Data(),"RECREATE");
580  histofile->cd();
581  TListIter iter(histolist);
582  TH1* tmph=NULL;
583  while( tmph=(TH1*)iter() ) tmph->Write();
584  histofile->Close();
585 
586  // plothistosfromfile(hfile.Data(),"pdf",3,3);
587 
588  //gROOT->Macro("drawhistos.C");//automatic draw 1D histos with gauss fit
589  return 0;
590 
591 }
fVerbose
int fVerbose
Definition: poormantracks.C:24
printf
printf("RealTime=%f seconds, CpuTime=%f seconds\n", rtime, ctime)
RhoCandidate::Combine
RhoCandidate * Combine(RhoCandidate *c)
Definition: RhoCandidate.cxx:1313
i
Int_t i
Definition: run_full.C:25
fCands
TClonesArray * fCands
Definition: poormantracks.C:27
sqrt
friend F32vec4 sqrt(const F32vec4 &a)
Definition: P4_F32vec4.h:29
RhoCandidate::Pos
TVector3 Pos() const
Definition: RhoCandidate.h:186
PoorManTracks
void PoorManTracks()
Definition: poormantracks.C:640
fMcCands
TClonesArray * fMcCands
Definition: poormantracks.C:26
RhoCandidate::Daughter
RhoCandidate * Daughter(Int_t n)
Definition: RhoCandidate.cxx:991
mommc
TVector3 mommc
Definition: anaLmdDigi.C:64
mom
Double_t mom
Definition: plot_dirc.C:14
max
friend F32vec4 max(const F32vec4 &a, const F32vec4 &b)
Definition: P4_F32vec4.h:26
RollVertexBox
TVector3 RollVertexBox(double widx, double widy, double widz)
Definition: poormantracks.C:596
fVertex
TVector3 fVertex(0., 0., 0.)
RhoCalculationTools::ForceConstantBz
static void ForceConstantBz(Double_t bz=0.)
Force a constant B field value for all positions.
Definition: RhoCalculationTools.h:37
fSeedPrefix
int fSeedPrefix
Definition: poormantracks.C:21
RhoCandidate::GetPosition
TVector3 GetPosition() const
Definition: RhoCandidate.h:185
RhoVtxPoca::GetPocaVtx
Double_t GetPocaVtx(TVector3 &vertex, RhoCandidate *composite)
Definition: RhoVtxPoca.cxx:27
seed
unsigned int seed
Definition: createPndLmdMisalignmentMatrices.C:5
RhoCandidate::P4
TLorentzVector P4() const
Definition: RhoCandidate.h:195
min
friend F32vec4 min(const F32vec4 &a, const F32vec4 &b)
Definition: P4_F32vec4.h:25
RhoCandidate::GetMcTruth
RhoCandidate * GetMcTruth() const
Definition: RhoCandidate.h:437
RhoCandidate::P4Cov
RhoError P4Cov() const
Definition: RhoCandidate.cxx:1579
RhoCandidate::Cov7
TMatrixD Cov7() const
Definition: RhoCandidate.cxx:1643
RhoCandidate::NDaughters
Int_t NDaughters() const
Definition: RhoCandidate.cxx:888
RhoCandidate::GetFit
RhoCandidate * GetFit() const
Definition: RhoCandidate.h:293
RhoCandidate::Charge
Double_t Charge() const
Definition: RhoCandidate.h:184
SetVerbose
drchit SetVerbose(iVerbose)
LoadPandaStyle
LoadPandaStyle()
TMatrixD
TMatrixT< double > TMatrixD
Definition: PndLmdDim.h:52
void PoorManTracks ( )

Definition at line 640 of file poormantracks.C.

References fAlternate, fCands, fDtheta, fMcCands, fNumTrk, fPDG, fPID, fSigPx, fSigPy, fSigPz, fSigVx, fSigVy, fSigVz, fVerbose, fVertex(), RhoCandidate::GetPosition(), p, RhoCandidate::P3Cov(), Pi, RollMomentumBox(), RhoCandidate::SetCov7(), RhoCandidate::SetMcTruth(), RhoCandidate::SetPos(), SmearMomentum(), and SmearVertex().

Referenced by poormantracks().

641 {
642  if(fVerbose>1) cout<<" --------------------------------------------------- "<<endl;
643  if(fVerbose>0) cout<<"Poor man tracks."<<endl;
644  fMcCands->Clear("C");
645  fCands->Clear("C");
646 
647  RhoError covPos(3), covP4(4);
648  covPos[0][0]=fSigVx*fSigVx;
649  covPos[1][1]=fSigVy*fSigVy;
650  covPos[2][2]=fSigVz*fSigVz;
651  covP4[0][0]=fSigPx*fSigPx;
652  covP4[1][1]=fSigPy*fSigPy;
653  covP4[2][2]=fSigPz*fSigPz;
654 
655  for (int iTr=0; iTr<fNumTrk; iTr++)
656  {
657  fPDG = TDatabasePDG::Instance()->GetParticle(fPID);
658  TVector3 momentum = RollMomentumBox(fVertex,fDtheta*TMath::Pi(),0.9,1.1);
659 
660  int sizeMc = fMcCands->GetEntriesFast();
661  RhoCandidate *pmc=new ((*fMcCands)[sizeMc]) RhoCandidate(momentum,fPDG);
662  //pmc->SetMcTruth(pmc);
663  pmc->SetPos(fVertex);
664  //pmc->SetType(fPID);
665  //pmc->SetMcMotherIdx(-1);
666  //if(fabs(fCharge)>0){
667  // bool rc = PndAnalysisCalcTools::FillHelixParams(pmc, true);
668  // if(!rc && fVerbose>0) {
669  // Warning("BuildMcCands()","Faild calculation helix parameters");
670  // std::cout<<*pmc<<std::endl;
671  // }
672  //}
673 
674  // Smeared "reco"
675  int size = fCands->GetEntriesFast();
676  TVector3 vertexMeasure = fVertex;
677  SmearVertex(vertexMeasure);
678  TVector3 momMeasure = momentum;
679  SmearMomentum(momMeasure);
680 
681 
682  RhoCandidate *p=new ((*fCands)[size]) RhoCandidate(momMeasure,fPDG);//TODO charge
683  p->SetMcTruth(pmc);
684  p->SetPos(vertexMeasure);
685  //p->SetType(fPID);
686  //p->SetMcMotherIdx(-1);
687  // setting cov's
688 // covP4[3][3]= p->Px()*p->Px()*covP4[0][0];
689 // covP4[3][3]+=p->Py()*p->Py()*covP4[1][1];
690 // covP4[3][3]+=p->Pz()*p->Pz()*covP4[2][2];
691 // covP4[3][3]/=(p->E()*p->E());
692  p->SetCov7(covPos,covP4);
693  //if(fabs(fCharge)>0){
694  // bool rc = PndAnalysisCalcTools::FillHelixParams(p, false);
695  // if(!rc && fVerbose>0) {
696  // Warning("BuildMcCands()","Faild calculation helix parameters");
697  // std::cout<<*p<<std::endl;
698  // }
699  //}
700  if(fVerbose>1) cout<<" ---------\n pmc, p, pmcpos, ppos, fVertex"<<endl;
701  if(fVerbose>1) std::cout<<*pmc<<std::endl;
702  if(fVerbose>1) std::cout<<*p<<std::endl;
703  if(fVerbose>1) p->P3Cov().Print();
704  if(fVerbose>1) pmc->GetPosition().Print();
705  if(fVerbose>1) p->GetPosition().Print();
706  if(fVerbose>1) fVertex.Print();
707  // print helix parameters and cov together with RhoCandidate V-P4 and cov7
708 
709  // alternate particle type
710  if(fAlternate) fPID*=-1;
711  }
712  // iterate again to set mother relations
713  // for (int i=0;i<fMcCands->GetEntriesFast();i++)
714  // {
715  // RhoCandidate* aMcCand=(RhoCandidate*)fMcCands->At(i);
716  // mcMotherID=aMcCand->GetMcMotherIdx();
717  // if (mcMotherID<0) continue;
718  // RhoCandidate* aMother=(RhoCandidate*)fMcCands->At(mcMotherID);
719  // if(0 == aMother) continue;
720  // aMcCand->SetMotherLink(aMother); // This should add the mother-daughter and daughter-mother relation
721  // }
722 
723 }
fVerbose
int fVerbose
Definition: poormantracks.C:24
SmearMomentum
SmearMomentum(TVector3 &momentum)
Definition: poormantracks.C:631
fNumTrk
int fNumTrk
Definition: poormantracks.C:19
RhoCandidate::P3Cov
RhoError P3Cov() const
Definition: RhoCandidate.cxx:1596
RhoCandidate::SetPos
void SetPos(const TVector3 &pos)
Definition: RhoCandidate.h:235
fCands
TClonesArray * fCands
Definition: poormantracks.C:27
SmearVertex
void SmearVertex(TVector3 &vertex)
Definition: poormantracks.C:606
fSigPy
double fSigPy
Definition: poormantracks.C:13
fSigVx
double fSigVx
Definition: poormantracks.C:9
fMcCands
TClonesArray * fMcCands
Definition: poormantracks.C:26
fSigVy
double fSigVy
Definition: poormantracks.C:10
fVertex
TVector3 fVertex(0., 0., 0.)
p
Double_t p
Definition: anasim.C:58
fSigVz
double fSigVz
Definition: poormantracks.C:11
RhoCandidate::GetPosition
TVector3 GetPosition() const
Definition: RhoCandidate.h:185
fPID
int fPID
Definition: poormantracks.C:16
fPDG
TParticlePDG * fPDG
Definition: poormantracks.C:18
fDtheta
double fDtheta
Definition: poormantracks.C:15
RhoCandidate::SetMcTruth
void SetMcTruth(RhoCandidate *mct)
Definition: RhoCandidate.h:436
fSigPx
double fSigPx
Definition: poormantracks.C:12
Pi
Double_t Pi
Definition: PndGeoHypGeCrystal.cxx:10
RollMomentumBox
TVector3 RollMomentumBox(const TVector3 &vtx, double dtheta, double ptmin, double ptmax)
Definition: poormantracks.C:617
fAlternate
bool fAlternate
Definition: poormantracks.C:17
RhoCandidate::SetCov7
void SetCov7(const TMatrixD &cov7)
Definition: RhoCandidate.cxx:413
fSigPz
double fSigPz
Definition: poormantracks.C:14
TVector3 RollMomentumBox ( const TVector3 &  vtx,
double  dtheta,
double  ptmin,
double  ptmax 
)

Definition at line 617 of file poormantracks.C.

References fVerbose, phi, pt(), theta, and CAMath::TwoPi().

Referenced by PoorManTracks().

618 {
619  // allow whole phi range
620  double phi=TMath::TwoPi()*gRandom->Rndm();
621  // allow theta at vertex pointing theta +- dtheta
622  double theta = vtx.Theta() + dtheta*2*(0.5-gRandom->Rndm());
623  // allow transverse momentum in [ptmin,ptmax]
624  double pt = ptmin + (ptmax-ptmin)*gRandom->Rndm();
625  TVector3 momentum;
626  momentum.SetPtThetaPhi(pt,theta,phi);
627  if(fVerbose>2) cout<<"RollMomentumBox: ("<<momentum.x()<<","<<momentum.y()<<","<<momentum.z()<<")"<<endl;
628  return momentum;
629 }
fVerbose
int fVerbose
Definition: poormantracks.C:24
theta
double theta
Definition: macro/detectors/emc/dedicated/reco_analys.C:24
pt
TString pt(TString pts, TString exts="px py pz")
Definition: invexp.C:133
dtheta
double dtheta
Definition: anaLmdCluster.C:54
phi
double phi
Definition: macro/detectors/emc/dedicated/reco_analys.C:24
CAMath::TwoPi
static float TwoPi()
Definition: PndCAMath.h:61
TVector3 RollVertexBox ( double  widx,
double  widy,
double  widz 
)

Definition at line 596 of file poormantracks.C.

References fVerbose.

Referenced by poormantracks().

597 {
598  double rnd1=widx*2*(0.5-gRandom->Rndm());
599  double rnd2=widy*2*(0.5-gRandom->Rndm());
600  double rnd3=widz*2*(0.5-gRandom->Rndm());
601  TVector3 vertex(rnd1,rnd2,rnd3);
602  if(fVerbose>2) cout<<"RollVertexBox: ("<<vertex.x()<<","<<vertex.y()<<","<<vertex.z()<<")"<<endl;
603  return vertex;
604 }
fVerbose
int fVerbose
Definition: poormantracks.C:24
SmearMomentum ( TVector3 &  momentum)

Definition at line 631 of file poormantracks.C.

References fSigPx, fSigPy, and fSigPz.

Referenced by PoorManTracks().

632 {
633  double rnd1=gRandom->Gaus(momentum.Px(),fSigPx);
634  double rnd2=gRandom->Gaus(momentum.Py(),fSigPy);
635  double rnd3=gRandom->Gaus(momentum.Pz(),fSigPz);
636  momentum.SetXYZ(rnd1,rnd2,rnd3);
637  return;
638 }
fSigPy
double fSigPy
Definition: poormantracks.C:13
fSigPx
double fSigPx
Definition: poormantracks.C:12
fSigPz
double fSigPz
Definition: poormantracks.C:14
void SmearVertex ( TVector3 &  vertex)

Definition at line 606 of file poormantracks.C.

References fSigVx, fSigVy, fSigVz, and fVerbose.

Referenced by PoorManTracks().

607 {
608  double rnd1=gRandom->Gaus(vertex.x(),fSigVx);
609  double rnd2=gRandom->Gaus(vertex.y(),fSigVy);
610  double rnd3=gRandom->Gaus(vertex.z(),fSigVz);
611  if(fVerbose>2) cout<<"SmearVertex (B): ("<<vertex.x()<<","<<vertex.y()<<","<<vertex.z()<<")"<<endl;
612  vertex.SetXYZ(rnd1,rnd2,rnd3);
613  if(fVerbose>2) cout<<"SmearVertex (A): ("<<vertex.x()<<","<<vertex.y()<<","<<vertex.z()<<")"<<endl;
614  return;
615 }
fVerbose
int fVerbose
Definition: poormantracks.C:24
fSigVx
double fSigVx
Definition: poormantracks.C:9
fSigVy
double fSigVy
Definition: poormantracks.C:10
fSigVz
double fSigVz
Definition: poormantracks.C:11

Variable Documentation

bool fAlternate =true

Definition at line 17 of file poormantracks.C.

Referenced by PoorManTracks().

TClonesArray* fCands =0

Definition at line 27 of file poormantracks.C.

Referenced by ana_dsdsj(), ana_pid(), poormantracks(), and PoorManTracks().

double fDtheta = 1.

Definition at line 15 of file poormantracks.C.

Referenced by PoorManTracks().

TClonesArray* fMcCands =0

Definition at line 26 of file poormantracks.C.

Referenced by ana_d0d0b_micro(), ana_dsdsj2(), ana_dsdsj_full(), ana_etac_3pi0(), poormantracks(), and PoorManTracks().

int fNumTrk =6

Definition at line 19 of file poormantracks.C.

Referenced by PoorManTracks().

TParticlePDG* fPDG = 0

Definition at line 18 of file poormantracks.C.

Referenced by PoorManTracks().

int fPID =211

Definition at line 16 of file poormantracks.C.

Referenced by PoorManTracks().

int fSeedPrefix =936650

Definition at line 21 of file poormantracks.C.

Referenced by poormantracks().

double fSigPx =0.01

Definition at line 12 of file poormantracks.C.

Referenced by PoorManTracks(), and SmearMomentum().

double fSigPy =0.01

Definition at line 13 of file poormantracks.C.

Referenced by PoorManTracks(), and SmearMomentum().

double fSigPz =0.01

Definition at line 14 of file poormantracks.C.

Referenced by PoorManTracks(), and SmearMomentum().

double fSigVx =0.005

Definition at line 9 of file poormantracks.C.

Referenced by PoorManTracks(), and SmearVertex().

double fSigVy =0.005

Definition at line 10 of file poormantracks.C.

Referenced by PoorManTracks(), and SmearVertex().

double fSigVz =0.005

Definition at line 11 of file poormantracks.C.

Referenced by PoorManTracks(), and SmearVertex().

int fVerbose =0

Definition at line 24 of file poormantracks.C.

Referenced by PndSdsNoiseProducer::AddDigiPixel(), PndSdsNoiseProducer::AddDigiStrip(), PndGemFindClusters::AddDigiToCluster(), PndSdsHybridHitProducer::AddHit(), PndSdsHybridHitProducer::AddHits(), PndSttCellTrackletGenerator::AddHitToBestCombi(), PndBarrelTrackFinder::AddHitToPreviousHits(), PndSttMvdGemTracking::AddHitToTrack(), PndLmdStripClusterTask::AddMSErr(), PndSdsDigiStripWriteoutBuffer::AddNewDataToTClonesArray(), PndEmcDigiWriteoutBuffer::AddNewDataToTClonesArray(), PndEmcWaveformWriteoutBuffer::AddNewDataToTClonesArray(), PndEmcWaveformSorterTask::AddNewDataToTClonesArray(), PndEmcDigiSorterTask::AddNewDataToTClonesArray(), PndMdtDigiSorterTask::AddNewDataToTClonesArray(), PndRichHitWriteoutBuffer::AddNewDataToTClonesArray(), PndSdsDigiPixelWriteoutBuffer::AddNewDataToTClonesArray(), PndGemDigiWriteoutBuffer::AddNewDataToTClonesArray(), PndFtsHitWriteoutBuffer::AddNewDataToTClonesArray(), PndSttHitWriteoutBuffer::AddNewDataToTClonesArray(), PndSdsDigiStripSorterTask::AddNewDataToTClonesArray(), PndFtsHitSorterTask::AddNewDataToTClonesArray(), PndSttHitSorterTask::AddNewDataToTClonesArray(), PndSdsDigiPixelSorterTask::AddNewDataToTClonesArray(), PndGemDigiSorterTask::AddNewDataToTClonesArray(), PndMvdTopix4RingSorterTask::AddNewDataToTClonesArray(), PndSciTHitSorterTask::AddNewDataToTClonesArray(), PndSciTHitWriteoutBuffer::AddNewDataToTClonesArray(), PndDrcDigiWriteoutBuffer::AddNewDataToTClonesArray(), PndEmcWaveformBuffer::AddNewDataToTClonesArray(), PndRichHitSorterTask::AddNewDataToTClonesArray(), PndDrcDigiSorterTask::AddNewDataToTClonesArray(), PndMdtWaveformWriteoutBuffer::AddNewDataToTClonesArray(), PndDiscWriteoutBuffer::AddNewDataToTClonesArray(), PndGemEventBuilderOnTracks::AddNewDataToTClonesArray(), PndSttMvdGemTracking::AddRemainingHits(), PndSttCellTrackletGenerator::AddRemainingHits(), PndEventBuilderManager::AnalyzeAndExtractEvents(), PndSttCellTrackletGenerator::AssignAmbiguousHits(), PndSttMvdGemTracking::AssignHits(), PndSttHitCorrector::CalcClassification(), PndHypIdealRecoTask::CalcDetPlane(), PndSdsRecoTask::CalcDetPlane(), PndMvdIdealRecoTask::CalcDetPlane(), PndSdsIdealRecoTask::CalcDetPlane(), PndSttCellTrackletGenerator::CalcDeviation(), PndSttCellTrackletGenerator::CalcDeviationOfRiemannTrack(), PndSttHitCorrector::CalcDifferencesBetweenAngles(), PndGemSmearingTask::CalcGFDetPlane(), PndSdsStripClusterTask::CalcMeanCharge(), PndLmdNoiseProducer::CalcReadoutCycles(), PndMvdNoiseProducer::CalcReadoutCycles(), PndHypGeGammaAna::CalculateCompton(), PndMvdQATask::CalculateDifference(), PndSttHitCorrector::CalculateTangentAngles(), PndSttMvdGemTracking::CheckCombinatorial(), PndSttCellTrackFinderAnalysisTask::CheckFirstTracklets(), PndSttCellTrackFinderAnalysisTask::CheckTrackletCombinations(), PndLmdSigCleanTask::CheckXY(), ClassImp(), PndSttCellTrackletGenerator::CombineTrackletsMultiStages(), PndSttCellTrackletGenerator::CombineTrackletsMultiStagesRecursive(), PndGemFindClustersTB::CompareDigiToClusters(), PndGemFindClusters::CompareDigiToClustersDigis(), PndGemEventBuilderOnTracks::CompareTrackToPreviousEvents(), PndGemFindHits::ConfirmHits(), PndGemFindHits::ConfirmHits2(), PndLmdTrackFinderCATask::CookAllCells(), PndLmdTrackFinderCATask::CookCells(), PndSttMvdGemTracking::Copy(), PndSttHitCorrector::CorrectHits(), PndSttHitCorrector::CorrectIsochrones(), PndSttMvdGemTracking::CountHitsInTrack(), PndSttCellTrackFinderAnalysisTask::CountMaxHitsCombi(), PndSttCellTrackFinderAnalysisTask::CountMaxHitsFirstStep(), PndGemFindClusters::CreateClusters(), PndGemDigiAna::CreateHistos(), PndGemFindHitsAna::CreateHistos(), PndSttCellTrackletGenerator::CreatePairCombis(), PndSttCellTrackletGenerator::CreatePndTrackCands(), PndSttCellTrackletGenerator::CreateRiemannTrack(), PndSttMvdGemTracking::DeleteHitFromTrack(), PndGemDigitize::DigitizeEvent(), PndGemDigitize::DigitizeRealisticEvent(), PndPidEmcBayesAssociatorTask::DoPidMatch(), PndSttCellTrackletGenerator::EvaluateMultiState(), PndSttMvdGemTracking::EvaluatePerformances(), PndSttCellTrackletGenerator::EvaluateState(), PndLmdTrackFinderCATask::Evolution(), PndEventCounterTask::Exec(), PndGapEventBuilderTask::Exec(), PndMvdReadInTBDataTask::Exec(), PndPidMvdAssociatorTask::Exec(), PndPidFtofAssociatorTask::Exec(), PndRichDbTask::Exec(), PndPidSciTAssociatorTask::Exec(), PndTrkAddTCA::Exec(), PndMvdSttGemRiemannTrackFinderTask::Exec(), PndRichRecoTask::Exec(), PndMcListConverter::Exec(), PndMvdRiemannTrackFinderTask::Exec(), PndLmdHitMergeTask::Exec(), PndLmdPixelClusterTask::Exec(), PndPidMlAssociatorTask::Exec(), TtCracowTask::Exec(), PndPidMdtHCAssociatorTask::Exec(), PndLmdNoiseProducer::Exec(), PndLmdStripClusterTask::Exec(), PndFtsHitProducerRealFast::Exec(), PndPidRichAssociatorTask::Exec(), PndSttCellTrackFinderAnalysisTask::Exec(), PndFtsDataAccessor::Exec(), PndFtsHitProducerMcPointCoordinates::Exec(), PndMdtHitProducer::Exec(), PndSttHitProducerRealFast::Exec(), PndPidDrcAssociatorTask::Exec(), PndPidDiscAssociatorTask::Exec(), PndHypGeGammaAna::Exec(), PndPidTestTask::Exec(), PndSttHelixHitProducer::Exec(), PndMvdEventMergerTask::Exec(), PndLmdNoiseTrkSuppressionTask::Exec(), PndSttMvdGemTracking::Exec(), PndPidIdealFwdOnlyAssociatorTask::Exec(), PndLmdTrackFinderTask::Exec(), PndLmdTrksFilterTask::Exec(), PndPidIdealAssociatorTask::Exec(), PndStraightLineTrackFinderTask::Exec(), PndFtsTrackerIdeal::Exec(), PndLmdLineTask::Exec(), PndLmdTrackFinderCATask::Exec(), PndEmcMakeClusterOnline::Exec(), PndLmdLinFitTask::Exec(), PndEmcFWEndcapTimebasedWaveforms::Exec(), PndTrackSmearTask::Exec(), PndMvdIdealTrackFinderTask::Exec(), PndLmdBPRungeKuttaTask::Exec(), PndGemFindHitsQA::Exec(), PndGemTrackFinderQA::Exec(), PndHypIdealTrackingTask::Exec(), PndSttHitProducerIdeal::Exec(), PndGemFindHitsAna::Exec(), PndCATracking::Exec(), PndGemDigiAna::Exec(), PndGemMatchHits::Exec(), PndLmdDigiProducer::Exec(), PndHypStripClusterTask::Exec(), PndLmdBPtestTask::Exec(), PndLmdGeaneTask::Exec(), PndRichHitFinder::Exec(), PndDrcHitFinder::Exec(), PndFtofHitProducerIdeal::Exec(), PndPidEmcBayesAssociatorTask::Exec(), PndLmdDigiQTask::Exec(), PndLmdSigCleanTask::Exec(), PndSdsRecoTask::Exec(), PndDrcHitProducerIdeal::Exec(), PndMvdConvertApvTask::Exec(), PndGemSmearingTask::Exec(), PndTrackDraw::Exec(), PndTrackingQualityBarrelTaskNewLinks::Exec(), PndTrkLegendreNew::Exec(), PndTrkLegendreSecTask::Exec(), PndTrkLegendreSecTask2::Exec(), PndTrkLegendreTask::Exec(), PndHypIdealRecoTask::Exec(), PndRecoMultiKalmanTask2::Exec(), PndLmdTrkQTask::Exec(), PndTrackingQATask::Exec(), PndLmdIdealFitTask::Exec(), PndSorterTask::Exec(), PndMvdIdealRecoTask::Exec(), PndTrkCombiLegendreTask::Exec(), PndSciTHitProducerIdeal::Exec(), PndSttIsochroneDraw::Exec(), PndMvdNoiseProducer::Exec(), PndEmcHitProducer::Exec(), PndTrkTrackFinder::Exec(), PndGemDigitize::Exec(), PndSciTDigiTask::Exec(), PndFtsHoughTrackerTaskQA::Exec(), PndPidMvaAssociatorTask::Exec(), PndLmdKalmanTask::Exec(), PndRecoKalmanTask2::Exec(), PndGemFindHits::Exec(), PndEmcExpClusterSplitter::Exec(), PndXYResidualTask::Exec(), PndGemFindClusters::Exec(), PndFtsHoughTrackerTask::Exec(), PndGemFindClustersTB::Exec(), PndRecoMultiKalmanTask::Exec(), PndSdsNoiseProducer::Exec(), PndPosCorrectorTask::Exec(), PndSdsHybridHitProducer::Exec(), PndSdsIdealRecoTask::Exec(), PndTSCorrectorTask::Exec(), PndBarrelTrackFinder::Exec(), PndMapSorterTask::Exec(), PndHypStripHitProducer::Exec(), PndRecoKalmanTask::Exec(), PndSdsStripHitProducerDif::Exec(), PndSdsStripHitProducer::Exec(), PndSdsPixelClusterTask::Exec(), PndSdsStripClusterTask::Exec(), PndTrkLegendreNew::ExtractLegendre(), PndTrkTrackFinder::ExtractLegendre(), PndTrkLegendreSecTask::ExtractLegendre(), PndTrkLegendreTask::ExtractLegendre(), PndTrkLegendreSecTask2::ExtractLegendre(), PndBarrelTrackFinder::ExtractMeanRPhiFromTrack(), PndBarrelTrackFinder::ExtractMeanZ_PFromTrack(), PndSdsStripClusterTask::FillClusterFinders(), PndLmdNoiseProducer::FillSensorLists(), PndMvdNoiseProducer::FillSensorLists(), PndEventBuilderOnMCEvents::FindEvents(), PndGemEventBuilderOnTracks::FindEvents(), PndGemFindHits::FindHits2(), PndLmdTrackFinderTask::FindHitsI(), PndStraightLineTrackFinderTask::FindHitsI(), PndLmdTrackFinderTask::FindHitsII(), PndStraightLineTrackFinderTask::FindHitsII(), PndLmdTrackFinderTask::FindHitsIII(), PndStraightLineTrackFinderTask::FindHitsIII(), PndTrkLegendreNew::FindMvdPixelReferenceHit(), PndTrkLegendreSecTask::FindMvdPixelReferenceHit(), PndTrkLegendreTask::FindMvdPixelReferenceHit(), PndTrkCombiLegendreTask::FindMvdPixelReferenceHit(), PndTrkLegendreSecTask2::FindMvdPixelReferenceHit(), PndTrkTrackFinder::FindMvdPixelReferenceHit(), PndTrkLegendreNew::FindMvdStripReferenceHit(), PndTrkLegendreSecTask::FindMvdStripReferenceHit(), PndTrkLegendreTask::FindMvdStripReferenceHit(), PndTrkCombiLegendreTask::FindMvdStripReferenceHit(), PndTrkLegendreSecTask2::FindMvdStripReferenceHit(), PndTrkTrackFinder::FindMvdStripReferenceHit(), PndHypGeGammaAna::FindPeaks(), PndTrkLegendreNew::FindReferenceHit(), PndTrkCombiLegendreTask::FindReferenceHit(), PndTrkTrackFinder::FindReferenceHit(), PndTrkLegendreSecTask2::FindReferenceHit(), PndTrkLegendreNew::FindSttReferenceHit(), PndTrkLegendreTask::FindSttReferenceHit(), PndTrkLegendreSecTask::FindSttReferenceHit(), PndTrkLegendreSecTask2::FindSttReferenceHit(), PndTrkTrackFinder::FindSttReferenceHit(), PndSdsDigiStripWriteoutBuffer::FindTimeForData(), PndSdsDigiPixelWriteoutBuffer::FindTimeForData(), PndSciTHitWriteoutBuffer::FindTimeForData(), PndSttCellTrackFinder::FindTracks(), PndSttCellTrackletGenerator::FindTracks(), PndHypGeGammaAna::Finish(), PndSttCellTrackFinderAnalysisTask::Finish(), PndLmdBPtestTask::Finish(), PndFtsHoughTrackerTaskQA::Finish(), PndFtsHoughTrackerTask::Finish(), PndEmcMakeClusterOnline::FinishClusters(), PndMdtDigitization::FinishTask(), PndMdtHitProducer::FinishTask(), PndEmcMakeClusterOnline::FinishTask(), PndSorterTask::FinishTask(), PndMapSorterTask::FinishTask(), PndHypGeGammaAna::FittingOutput(), PndSttMvdGemTracking::ForbidMultiAssignedHits(), PndLmdTrackFinderCATask::ForwardEvolution(), PndSttCellTrackletGenerator::GenerateTracklets(), PndSttCellTrackletGenerator::GenerateTrackletsGPU(), PndSttHitCorrector::GetBestCombinatedPhi(), PndSttCellTrackletGenerator::GetDeviationCount(), PndPidCorrelator::GetFMdtInfo(), PndPidCorrelator::GetFscInfo(), PndPidCorrelator::GetFtofInfo(), PndSttMvdGemTracking::GetHitsAssociatedToTrack(), PndSttMvdGemTracking::GetInitialParams(), PndSdsHybridHitProducer::GetLocalHitPoints(), PndSttHitCorrector::GetMiddleHitNeighbor(), PndPidCorrelator::GetMvdInfo(), PndSttCellTrackFinderAnalysisTask::GetNumLinksOfHits(), PndPidEmcBayesAssociatorTask::GetPdf(), PndHypGeGammaAna::GetPeakContent(), PndPidCorrelator::GetRichInfo(), PndSttHitCorrector::GetStraightNeighbors(), PndPidCorrelator::GetTrackInfo(), PndSttCellTrackletGenerator::GetTrackletInf(), PndSttMvdGemTracking::GetTracksAssociatedToHit(), PndSdsHybridHitProducer::GetTransformation(), PndFtsHoughTrackerTask::GetVerbose(), PndMvdReadInTBDataTask::Init(), PndLmdHybridHitProducer::Init(), PndMvdHybridHitProducer::Init(), PndLmdNoiseProducer::Init(), PndTrkGemCombiTask::Init(), PndSttCellTrackFinderAnalysisTask::Init(), PndSttMvdGemTracking::Init(), PndLmdPixelClusterTask::Init(), PndHypGeCOSYBackgroundAna::Init(), PndHypGeGammaAna::Init(), PndFtsDataAccessor::Init(), PndLmdLinFitTask::Init(), PndTrackCandDraw::Init(), PndFtsTrackerIdeal::Init(), PndMvdQATask::Init(), PndTrkLegendreTask::Init(), PndSdsRecoTask::Init(), PndTrkLegendreNew::Init(), PndTrkLegendreSecTask2::Init(), PndTrkLegendreSecTask::Init(), PndSorterTask::Init(), PndTrkCombiLegendreTask::Init(), PndMvdNoiseProducer::Init(), PndTrkTrackFinder::Init(), PndTrackDraw::Init(), PndFtsHoughTrackerTaskQA::Init(), PndSttIsochroneDraw::Init(), PndMvdConvertApvTask::Init(), PndRecoKalmanTask2::Init(), PndSdsNoiseProducer::Init(), PndFtsHoughTrackerTask::Init(), PndPosCorrectorTask::Init(), PndSdsHybridHitProducer::Init(), PndTSCorrectorTask::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndPidSttAssociatorTask::Init(), PndMapSorterTask::Init(), PndHypStripHitProducer::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndSdsStripHitProducer::Init(), PndSdsPixelClusterTask::Init(), PndSttCellTrackletGenerator::InitStartTracklets(), PndHypIdealRecoTask::InitTransMat(), PndGemSmearingTask::InitTransMat(), PndSdsRecoTask::InitTransMat(), PndMvdIdealRecoTask::InitTransMat(), PndSdsIdealRecoTask::InitTransMat(), PndSttCellTrackletGenerator::InsertCombination(), PndSttMvdGemTracking::IntersectionFinder(), PndEmc2DLocMaxFinder::isALocalMax(), PndSttMvdGemTracking::IsAssignable(), PndSttCellTrackletGenerator::IsEndTubeOfTracklet(), PndGemFindClusters::JoinTwoClusters(), PndSttHitCorrector::KeepBestAngle(), PndLmdLinFitTask::line3DfitMS(), PndBarrelTrackFinder::MatchHitToHit(), PndBarrelTrackFinder::MatchHitWithZInfoTT(), PndBarrelTrackFinder::MatchParallelSttHitTT(), PndGemTrackFinderQA::MatchRecoTracks(), PndBarrelTrackFinder::MatchSkewedSttHitTT(), PndPidCorrelator::MdtGeometry(), PndSdsDigiStripWriteoutBuffer::Modify(), PndSdsDigiPixelWriteoutBuffer::Modify(), PndGemDigiWriteoutBuffer::Modify(), PndEmcWaveformWriteoutBuffer::Modify(), PndDrcDigiWriteoutBuffer::Modify(), PndSciTHitWriteoutBuffer::Modify(), PndMdtWaveformWriteoutBuffer::Modify(), PndSttMvdGemTracking::OnlyOneHitToEachTrack(), PndSttMvdGemTracking::OrderGemHits(), plotHoughSpaces(), plotTrackCands(), PndDrcHitFinder::PndDrcHitFinder(), PndDrcHitProducerIdeal::PndDrcHitProducerIdeal(), PndEventBuilderOnMCEvents::PndEventBuilderOnMCEvents(), PndFtsDataAccessor::PndFtsDataAccessor(), PndFtsHoughTrackerTask::PndFtsHoughTrackerTask(), PndFtsHoughTrackerTaskQA::PndFtsHoughTrackerTaskQA(), PndFtsSignalOverlap::PndFtsSignalOverlap(), PndFtsTrackerIdeal::PndFtsTrackerIdeal(), PndGemEventBuilderOnTracks::PndGemEventBuilderOnTracks(), PndLmdDigiTask::PndLmdDigiTask(), PndLmdHybridHitProducer::PndLmdHybridHitProducer(), PndMdtHitProducer::PndMdtHitProducer(), PndMdtWaveformWriteoutBuffer::PndMdtWaveformWriteoutBuffer(), PndMixBackgroundEvents::PndMixBackgroundEvents(), PndMvdHybridHitProducer::PndMvdHybridHitProducer(), PndMvdPixelDigiSorterTask::PndMvdPixelDigiSorterTask(), PndRichDbTask::PndRichDbTask(), PndRichHitFinder::PndRichHitFinder(), PndRichRecoTask::PndRichRecoTask(), PndSdsHybridHitProducer::PndSdsHybridHitProducer(), PndSdsRecoTask::PndSdsRecoTask(), PndSttHelixHitProducer::PndSttHelixHitProducer(), PndSttMvdGemTracking::PndSttMvdGemTracking(), PndSttSignalOverlap::PndSttSignalOverlap(), PndSttTrackFitterQATask::PndSttTrackFitterQATask(), poormantracks(), PoorManTracks(), PndTrkLegendreTask::PrepareLegendre(), PndTrkLegendreSecTask::PrepareLegendre(), PndTrkLegendreSecTask2::PrepareLegendre(), PndGemFindClusters::PrintCluster(), PndGemFindClusters::PrintClusters(), PndSttMvdGemTracking::PropagateToGemPlane(), PndSttMvdGemTracking::PropagateToGemPlaneAsHelix(), PndTrackDraw::PropagateTrack(), PndTrackDraw::PropagateTrackHelix(), PndSttCellTrackletGenerator::RefitTracks(), PndTrkLegendreNew::Refresh(), PndTrkCombiLegendreTask::Refresh(), PndTrkTrackFinder::Refresh(), PndTrkLegendreSecTask::Refresh(), PndTrkLegendreTask::Refresh(), PndTrkLegendreSecTask2::Refresh(), PndFtsDataAccessor::Register(), PndFtsTrackerIdeal::Register(), PndTrkLegendreTask::RegisterTrack(), PndTrkLegendreSecTask::RegisterTrack(), PndTrkLegendreSecTask2::RegisterTrack(), PndFtsHoughTrackerTaskQA::ReInit(), PndFtsHoughTrackerTask::ReInit(), PndBarrelTrackFinder::RemoveCloneTracks(), PndSttMvdGemTracking::Reset(), PndTrkLegendreTask::Reset(), PndTrkLegendreNew::Reset(), PndTrkLegendreSecTask::Reset(), PndTrkLegendreSecTask2::Reset(), PndTrkCombiLegendreTask::Reset(), PndTrkTrackFinder::Reset(), PndSttMvdGemTracking::Retrack(), RollMomentumBox(), RollVertexBox(), PndSdsStripClusterTask::SelectSensorParams(), PndSdsStripHitProducerDif::SelectSensorParams(), PndSdsStripHitProducer::SelectSensorParams(), PndLmdPixelClusterTask::SetBackMapping(), PndMvdPixelClusterTask::SetBackMapping(), PndMvdHybridHitProducer::SetBranchNames(), PndLmdHybridHitProducer::SetBranchNames(), PndMvdPixelDigiSorterTask::SetBranchNames(), PndLmdStripClusterTask::SetCalculators(), PndLmdStripHitProducer::SetCalculators(), PndMvdStripClusterTask::SetCalculators(), PndMvdStripHitProducer::SetCalculators(), PndSdsStripClusterTask::SetCalculators(), PndSdsStripHitProducer::SetCalculators(), PndLmdPixelClusterTask::SetClusterFinder(), PndMvdPixelClusterTask::SetClusterFinder(), PndSttCellTrackletGenerator::SetCorrectedHits(), PndMdtTrkProducer::SetGeometry(), PndMvdStripClusterTask::SetParContainers(), PndMvdPixelClusterTask::SetParContainers(), PndRichHitFinder::SetParContainers(), PndDrcHitFinder::SetParContainers(), PndFtsHoughTrackerTask::SetParContainers(), PndMvdStripHitProducer::SetParContainers(), PndSdsHybridHitProducer::SetParContainers(), PndSdsStripHitProducer::SetParContainers(), PndSdsPixelClusterTask::SetParContainers(), PndSttMvdGemTracking::SetStartParameters(), PndSttMvdGemTracking::SetupGEMPlanes(), PndBlackBoxTask::SetVerbose(), PndEventCounterTask::SetVerbose(), PndMcListConverter::SetVerbose(), PndPidSciTAssociatorTask::SetVerbose(), PndUnassignedHitsTask::SetVerbose(), PndMvdRiemannTrackFinderTaskEff::SetVerbose(), PndPidFtofAssociatorTask::SetVerbose(), PndPidMvdAssociatorTask::SetVerbose(), PndMQStraightLineTrackFinder::SetVerbose(), PndMvdRiemannTrackFinderTaskCutPar::SetVerbose(), PndSttCellTrackFinderAnalysisTask::SetVerbose(), PndLmdTrackFinderTask::SetVerbose(), PndPidMlAssociatorTask::SetVerbose(), PndStraightLineTrackFinderTask::SetVerbose(), PndLmdTrackFinderCATask::SetVerbose(), PndMvdSttGemRiemannTrackFinderTask::SetVerbose(), PndPidMdtHCAssociatorTask::SetVerbose(), PndMvdRiemannTrackFinderTask::SetVerbose(), PndPidRichAssociatorTask::SetVerbose(), PndPidDrcAssociatorTask::SetVerbose(), PndPidDiscAssociatorTask::SetVerbose(), PndMvdIdealTrackFinderTask::SetVerbose(), PndGemTrackFinderQA::SetVerbose(), PndPidIdealFwdOnlyAssociatorTask::SetVerbose(), PndGemFindHitsQA::SetVerbose(), PndGemFindHitsAna::SetVerbose(), PndPidIdealAssociatorTask::SetVerbose(), PndGemDigiAna::SetVerbose(), PndTrackDraw::SetVerbose(), PndSttIsochroneDraw::SetVerbose(), PndPidEmcBayesAssociatorTask::SetVerbose(), PndPidMvaAssociatorTask::SetVerbose(), PndFtsHoughTrackerTask::SetVerbose(), PndPidSttAssociatorTask::SetVerbose(), PndHypIdealRecoTask::smear(), PndGemSmearingTask::smear(), PndMvdIdealRecoTask::smear(), PndSdsIdealRecoTask::smear(), PndHypIdealRecoTask::smearLocal(), PndGemSmearingTask::smearLocal(), PndMvdIdealRecoTask::smearLocal(), PndSdsIdealRecoTask::smearLocal(), SmearVertex(), PndGemFindClusters::SortCluster(), PndGemFindClusters::SortClusters(), PndGemFindClusters::SortDigis(), PndLmdTrackFinderCATask::SortHitsByDet(), PndLmdTrackFinderTask::SortHitsByDet(), PndStraightLineTrackFinderTask::SortHitsByDet(), PndLmdTrackFinderTask::SortHitsByDet2(), PndStraightLineTrackFinderTask::SortHitsByDet2(), PndLmdTrackFinderCATask::SortHitsByDetSimple(), PndLmdTrackFinderCATask::SortHitsByDetSimple2(), PndLmdTrackFinderCATask::SortHitsByZ(), PndLmdTrackFinderTask::SortHitsByZ(), PndStraightLineTrackFinderTask::SortHitsByZ(), PndSttCellTrackletGenerator::SplitData(), PndGemEventBuilderOnTracks::StoreEventData(), PndSttCellTrackFinderAnalysisTask::TestRecoQualityCombi(), PndSttCellTrackFinderAnalysisTask::TestRecoQualityFirstStep(), PndGemFindClusters::WriteClusters(), PndGemFindClustersTB::WriteClusters(), PndFtsHoughTrackerTask::~PndFtsHoughTrackerTask(), and PndFtsHoughTrackerTaskQA::~PndFtsHoughTrackerTaskQA().

int niter =1

Definition at line 20 of file poormantracks.C.

Referenced by RhoKalmanVtxFitter::FitVertexFast(), and PndLmdTrackFinderCATask::ForwardEvolution().

Generated on Sun Aug 18 2019 03:02:45 for FairRoot/PandaRoot by   doxygen 1.8.8