Fork me on GitHub

Changeset 3a105e5 in git for external


Ignore:
Timestamp:
Feb 24, 2022, 4:15:47 PM (3 years ago)
Author:
michele <michele.selvaggi@…>
Branches:
master
Children:
2a2f310
Parents:
56fb0be
Message:

added first hit to track

Location:
external/TrackCovariance
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • external/TrackCovariance/ObsTrk.cc

    r56fb0be r3a105e5  
    120120        //if (pt > maxPt) std::cout << "Warning ObsTrk::GenToObsPar: pt " << pt << " is above grid range of " << maxPt << std::endl;
    121121        Double_t minAn = fGC->GetMinAng();
    122         //if (angd < minAn) std::cout << "Warning ObsTrk::GenToObsPar: angle " << angd 
     122        //if (angd < minAn) std::cout << "Warning ObsTrk::GenToObsPar: angle " << angd
    123123        //      << " is below grid range of " << minAn << std::endl;
    124124        Double_t maxAn = fGC->GetMaxAng();
     
    137137        Double_t Xfirst, Yfirst, Zfirst;
    138138        Int_t iLay = trk->FirstHit(Xfirst, Yfirst, Zfirst);
    139         TVector3 fXfirst(Xfirst, Yfirst, Zfirst);
     139        fXfirst = TVector3(Xfirst, Yfirst, Zfirst);
     140  //std::cout<<"obs trk: "<<Xfirst<<","<<Yfirst<<","<<Zfirst<<std::endl;
     141
    140142        if (inside)
    141143        {
  • external/TrackCovariance/SolTrack.cc

    r56fb0be r3a105e5  
    176176{
    177177        //
    178         // Return lists of hits associated to a track including all scattering layers. 
     178        // Return lists of hits associated to a track including all scattering layers.
    179179        // Return value is the total number of measurement hits
    180180        // kmh = total number of measurement layers hit for given type
     
    190190        {
    191191                Double_t R; Double_t phi; Double_t zz;
    192                 if (HitLayer(i, R, phi, zz)) 
     192                if (HitLayer(i, R, phi, zz))
    193193                {
    194194                        zhh[kh] = zz;
     
    207207{
    208208        //
    209         // Return lists of hits associated to a track for all measurement layers. 
     209        // Return lists of hits associated to a track for all measurement layers.
    210210        // Return value is the total number of measurement hits
    211211        // kmh = total number of measurement layers hit for given type
     
    236236Int_t SolTrack::FirstHit(Double_t &Xfirst, Double_t &Yfirst, Double_t &Zfirst)
    237237{
    238         Int_t iFirst = -1;     
     238        Int_t iFirst = -1;
    239239        Int_t iFirstLay = -1;   // Default return with no hits
    240240        Xfirst = 0.;
     
    249249                Double_t *dh = new Double_t[Nmh];
    250250                //
    251                 Int_t n = HitListXYZ(ih, Xh, Yh, Zh);   
     251                Int_t n = HitListXYZ(ih, Xh, Yh, Zh);
    252252                //
    253253                for(Int_t i=0; i<Nmh; i++){
    254                         Double_t rr = TMath::Sqrt(Xh[i]*Xh[i]+Yh[i]*Yh[i]);     // Hit radius           
     254                        Double_t rr = TMath::Sqrt(Xh[i]*Xh[i]+Yh[i]*Yh[i]);     // Hit radius
    255255                        dh[i] = TMath::ASin(C() * TMath::Sqrt((rr * rr - D() * D()) / (1. + 2 * C() * D()))) / C();     // Arc length traveled
    256256                }
     
    327327        //
    328328        //
    329         // Input flags: 
     329        // Input flags:
    330330        //                              Res = .TRUE. turn on resolution effects/Use standard resolutions
    331331        //                                        .FALSE. set all resolutions to 0
     
    410410                        nz = 1.0;
    411411                }
    412                 Double_t corr = TMath::Abs(pxi*nx + pyi * ny + pzi * nz) / p(); 
     412                Double_t corr = TMath::Abs(pxi*nx + pyi * ny + pzi * nz) / p();
    413413                Double_t Rlf = fG->lTh(i) / (corr*fG->lX0(i));                                  // Rad. length fraction
    414414                thms[ii] = 0.0136*TMath::Sqrt(Rlf)*(1.0 + 0.038*TMath::Log(Rlf)) / p();         // MS angle
     
    438438                Int_t ityp  = fG->lTyp(i);                      // Layer type Barrel or Z
    439439                Int_t nmeai = fG->lND(i);                       // # measurements in layer
    440                
     440
    441441                if (fG->isMeasure(i))
    442442                {
     
    470470                                                Rm(im, 2) = csa * dRphi(2) - ssa * dRz(2);      // C derivative
    471471                                                Rm(im, 3) = csa * dRphi(3) - ssa * dRz(3);      // z0 derivative
    472                                                 Rm(im, 4) = csa * dRphi(4) - ssa * dRz(4);      // cot(theta) derivative       
     472                                                Rm(im, 4) = csa * dRphi(4) - ssa * dRz(4);      // cot(theta) derivative
    473473                                        }
    474474                                        if (ityp == 2)          // Z type layer (Measure R-phi at const. Z)
     
    531531                                                {
    532532                                                        Double_t strk = 0;
    533                                                         if (nmk + 1 == 1) strk = fG->lStU(k);   // Stereo angle upper 
     533                                                        if (nmk + 1 == 1) strk = fG->lStU(k);   // Stereo angle upper
    534534                                                        if (nmk + 1 == 2) strk = fG->lStL(k);   // Stereo angle lower
    535535                                                        //if (im == km && Res) Sm(im, km) += sig*sig;   // Detector resolution on diagonal
     
    542542                                                        //
    543543                                                        // Loop on all layers below for MS contributions
    544                                                         for (Int_t jj = 0; jj < kk; jj++)               
     544                                                        for (Int_t jj = 0; jj < kk; jj++)
    545545                                                        {
    546546                                                                Double_t di = dik(ii, jj);
Note: See TracChangeset for help on using the changeset viewer.