Changeset 0b8551f in git for external/TrackCovariance/AcceptanceClx.cc
- Timestamp:
- Nov 29, 2021, 4:04:38 PM (3 years ago)
- Branches:
- master
- Children:
- 0c0c9af
- Parents:
- 9a7ea36 (diff), bd376e3 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the(diff)
links above to see all the changes relative to each parent. - git-author:
- Michele Selvaggi <michele.selvaggi@…> (11/29/21 16:04:38)
- git-committer:
- GitHub <noreply@…> (11/29/21 16:04:38)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
external/TrackCovariance/AcceptanceClx.cc
r9a7ea36 r0b8551f 1 1 #include <TVector3.h> 2 2 #include "AcceptanceClx.h" 3 #include <iostream>4 3 // 5 4 // Pt splitting routine … … 59 58 // Initializations 60 59 // 61 // std::cout << "Entered constructor of AccpeptanceClx" << std::endl;60 //cout << "Entered constructor of AccpeptanceClx" << endl; 62 61 // Setup grid 63 62 // Start grid parameters … … 74 73 TVectorF Tha(NpThInp, ThInit); 75 74 Int_t NpTh = Tha.GetNrows(); // Nr. of starting theta points 76 // std::cout << "AcceptanceClv:: Pta and Tha arrays defined" << std::endl;75 //cout << "AcceptanceClv:: Pta and Tha arrays defined" << endl; 77 76 // 78 77 // Grid splitting parameters … … 128 127 // Pt split 129 128 if (dPt > dPtMin && dAp > dAmin && Nsplits < MaxSplits) { 130 //std::cout << "Pt(" << ipt << ") = " << Pta(ipt) << ", dAp = " << dAp << std::endl;131 129 NsplitCnt++; // Total splits counter 132 130 Nsplits++; // Increase splits/cycle 133 NpPt++; // Increase #pt points 134 //std::cout << "New pt split: dAp = " << dAp << ", dPt = " << dPt << 135 // ", Nsplits = " << Nsplits << ", NpPt = " << NpPt << std::endl; 131 NpPt++; // Increase #pt points 136 132 Float_t newPt = 0.5 * (Pta(ipt + 1) + Pta(ipt)); 137 133 VecInsert(ipt, newPt, Pta); … … 156 152 // Theta split 157 153 if (dTh > dThMin && dAt > dAmin && Nsplits < MaxSplits) { 158 // std::cout << "Th(" << ith << ") = " << Tha(ith) << ", dAt = " << dAt << std::endl;154 //cout << "Th(" << ith << ") = " << Tha(ith) << ", dAt = " << dAt << endl; 159 155 NsplitCnt++; // Total splits counter 160 156 Nsplits++; // Increase splits 161 NpTh++; // Increase #pt points 162 //std::cout << "New th split: dAt = " << dAt << ", dTh = " << dTh << 163 // ", Nsplits = " << Nsplits << ", NpTh = " << NpTh << std::endl; 157 NpTh++; // Increase #pt points 164 158 Float_t newTh = 0.5 * (Tha(ith + 1) + Tha(ith)); 165 159 VecInsert(ith, newTh, Tha); … … 202 196 fThArray = Tha; // Array of Theta nodes 203 197 // 204 std::cout << "AcceptanceClx:: Acceptance encoding with " << fNPtNodes198 std::cout << "AcceptanceClx:: Acceptance encoding with " << fNPtNodes 205 199 <<" pt nodes and "<< fNThNodes <<" theta nodes"<< std::endl; 206 200 Int_t Nrows = fAcc.GetNrows(); 207 201 Int_t Ncols = fAcc.GetNcols(); 208 //std::cout<<"AcceptanceClx:: fAcc rows: "<<Nrows<<", cols: "<<Ncols<<std::endl;209 //std::cout<<"AcceptanceClx:: Pta size: "<<Pta.GetNrows()<<", Tha size: "<<Tha.GetNrows()<<std::endl;;210 //std::cout<<"AcceptanceClx:: pt nodes"<<std::endl; fPtArray.Print();211 //std::cout<<"AcceptanceClx:: th nodes"<<std::endl; fThArray.Print();212 202 } 213 203 … … 290 280 // 291 281 // Protect against values out of range 292 //std::cout << "AcceptanceClx::HitNumber: just in pt= " << pt << ", theta = " << theta << std::endl;293 //std::cout << "AcceptanceClx::HitNumber: ptArr(0)= " << fPtArray(0) << ", thArr(0) = " << fThArray(0) << std::endl;294 //std::cout << "AcceptanceClx::HitNumber: ptArr(E)= " << fPtArray(fNPtNodes - 1)295 // << ", thArr(E) = " << fThArray(fNThNodes - 1) << std::endl;296 282 Float_t eps = 1.0e-4; 297 283 Float_t pt0 = (Float_t)pt; … … 311 297 { 312 298 std::cout << "Search error: (ip, pt) = (" << ip << ", " << pt << "), pt0 = " << pt0 << std::endl; 313 std::cout << "Search error: pt nodes = " << fNPtNodes 299 std::cout << "Search error: pt nodes = " << fNPtNodes 314 300 << " , last value = " << fPtArray(fNPtNodes - 1) << std::endl; 315 301 } … … 343 329 } 344 330 // 331
Note:
See TracChangeset
for help on using the changeset viewer.