Fork me on GitHub

Changeset bd33fce in git for examples


Ignore:
Timestamp:
Mar 8, 2021, 9:52:15 AM (4 years ago)
Author:
GitHub <noreply@…>
Branches:
master
Children:
a10293e, b1cb322
Parents:
2eaa9fd (diff), 91ef0b8 (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@…> (03/08/21 09:52:15)
git-committer:
GitHub <noreply@…> (03/08/21 09:52:15)
Message:

Merge pull request #88 from fbedesch/master

Fixed conflicts and Reginv3

Location:
examples
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • examples/ExamplePVtxFit.C

    r2eaa9fd rbd33fce  
    4040        //
    4141        // Loop over all events
    42         Int_t Nev = TMath::Min(Nevent, (Int_t) numberOfEntries);
     42        Int_t Nev = TMath::Min(Nevent, (Int_t)numberOfEntries);
    4343        for (Int_t entry = 0; entry < Nev; ++entry)
    4444        {
     
    6161                                GenParticle* gp = (GenParticle*)trk->Particle.GetObject();
    6262                                //
    63                                 // Position of origin in meters
    64                                 Double_t x = 1.0e-3 * gp->X;
    65                                 Double_t y = 1.0e-3 * gp->Y;
    66                                 Double_t z = 1.0e-3 * gp->Z;
     63                                // Position of origin in mm
     64                                Double_t x = gp->X;
     65                                Double_t y = gp->Y;
     66                                Double_t z = gp->Z;
    6767                                //
    6868                                // group tracks originating from the primary vertex
     
    7878                                        Double_t oPar[5] = { obsD0, obsPhi, obsC, obsZ0, obsCtg };
    7979                                        TVectorD obsPar(5, oPar);       // Fill observed parameters
    80                                         TVector3 xv(x, y, z);
    8180                                        //
    8281                                        pr[Ntr] = new TVectorD(obsPar);
    83                                         cv[Ntr] = new TMatrixDSym(TrkUtil::CovToMm(trk->CovarianceMatrix()));
     82                                        //std::cout<<"Cov Matrix:"<<std::endl;
     83                                        //trk->CovarianceMatrix().Print();
     84                                        cv[Ntr] = new TMatrixDSym(trk->CovarianceMatrix());
    8485                                        Ntr++;
    8586                                }
    8687                        }               // End loop on tracks
    87                         //std::cout << "Total of " << Ntr << " primary tracks out of " << NtrG << " tracks" << std::endl;
    8888                }
    8989                //
     
    118118        // Show resulting histograms
    119119        //
    120         TCanvas* Cnv = new TCanvas("Cnv", "Delphes generated track plots", 50, 50, 900, 500);
     120        TCanvas* Cnv = new TCanvas("Cnv", "Delphes primary vertex pulls", 50, 50, 900, 500);
    121121        Cnv->Divide(2, 2);
    122122        Cnv->cd(1); gPad->SetLogy(1); gStyle->SetOptFit(1111);
Note: See TracChangeset for help on using the changeset viewer.