Fork me on GitHub

Ignore:
Timestamp:
Nov 29, 2021, 4:04:38 PM (3 years ago)
Author:
GitHub <noreply@…>
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)
Message:

Merge pull request #102 from fbedesch/master

Major update to handle highly displaced tracks

File:
1 edited

Legend:

Unmodified
Added
Removed
  • external/TrackCovariance/AcceptanceClx.cc

    r9a7ea36 r0b8551f  
    11#include <TVector3.h>
    22#include "AcceptanceClx.h"
    3 #include <iostream>
    43//
    54// Pt splitting routine
     
    5958        // Initializations
    6059        //
    61         //std::cout << "Entered constructor of AccpeptanceClx" << std::endl;
     60        //cout << "Entered constructor of AccpeptanceClx" << endl;
    6261        // Setup grid
    6362        // Start grid parameters
     
    7473        TVectorF Tha(NpThInp, ThInit);
    7574        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;
    7776        //
    7877        // Grid splitting parameters
     
    128127                                // Pt split
    129128                                if (dPt > dPtMin && dAp > dAmin && Nsplits < MaxSplits) {
    130                                         //std::cout << "Pt(" << ipt << ") = " << Pta(ipt) << ", dAp = " << dAp << std::endl;
    131129                                        NsplitCnt++;    // Total splits counter
    132130                                        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
    136132                                        Float_t newPt = 0.5 * (Pta(ipt + 1) + Pta(ipt));
    137133                                        VecInsert(ipt, newPt, Pta);
     
    156152                                // Theta split
    157153                                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;
    159155                                        NsplitCnt++;    // Total splits counter
    160156                                        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
    164158                                        Float_t newTh = 0.5 * (Tha(ith + 1) + Tha(ith));
    165159                                        VecInsert(ith, newTh, Tha);
     
    202196        fThArray = Tha; // Array of Theta nodes
    203197        //
    204         std::cout << "AcceptanceClx:: Acceptance encoding with " << fNPtNodes
     198                std::cout << "AcceptanceClx:: Acceptance encoding with " << fNPtNodes
    205199                <<" pt nodes and "<< fNThNodes <<" theta nodes"<< std::endl;
    206200        Int_t Nrows = fAcc.GetNrows();
    207201        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();
    212202}
    213203
     
    290280        //
    291281        // 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;
    296282        Float_t eps = 1.0e-4;
    297283        Float_t pt0 = (Float_t)pt;
     
    311297        {
    312298                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 
    314300                        << " , last value = " << fPtArray(fNPtNodes - 1) << std::endl;
    315301        }
     
    343329}
    344330//
     331
Note: See TracChangeset for help on using the changeset viewer.