Changeset bd33fce in git for examples/ExamplePVtxFit.C

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

File:

• examples/ExamplePVtxFit.C (modified) (4 diffs)

examples/ExamplePVtxFit.C

@@ -40 +40 @@
         //
         // Loop over all events
-        Int_t Nev = TMath::Min(Nevent, (Int_t) numberOfEntries);
+        Int_t Nev = TMath::Min(Nevent, (Int_t)numberOfEntries);
         for (Int_t entry = 0; entry < Nev; ++entry)
         {
@@ -61 +61 @@
                                 GenParticle* gp = (GenParticle*)trk->Particle.GetObject();
                                 //
-                                // Position of origin in meters
-                                Double_t x = 1.0e-3 * gp->X;
-                                Double_t y = 1.0e-3 * gp->Y;
-                                Double_t z = 1.0e-3 * gp->Z;
+                                // Position of origin in mm
+                                Double_t x = gp->X;
+                                Double_t y = gp->Y;
+                                Double_t z = gp->Z;
                                 //
                                 // group tracks originating from the primary vertex
@@ -78 +78 @@
                                         Double_t oPar[5] = { obsD0, obsPhi, obsC, obsZ0, obsCtg };
                                         TVectorD obsPar(5, oPar);       // Fill observed parameters
-                                        TVector3 xv(x, y, z);
                                         //
                                         pr[Ntr] = new TVectorD(obsPar);
-                                        cv[Ntr] = new TMatrixDSym(TrkUtil::CovToMm(trk->CovarianceMatrix()));
+                                        //std::cout<<"Cov Matrix:"<<std::endl;
+                                        //trk->CovarianceMatrix().Print();
+                                        cv[Ntr] = new TMatrixDSym(trk->CovarianceMatrix());
                                         Ntr++;
                                 }
                         }               // End loop on tracks
-                        //std::cout << "Total of " << Ntr << " primary tracks out of " << NtrG << " tracks" << std::endl;
                 }
                 //
@@ -118 +118 @@
         // Show resulting histograms
         //
-        TCanvas* Cnv = new TCanvas("Cnv", "Delphes generated track plots", 50, 50, 900, 500);
+        TCanvas* Cnv = new TCanvas("Cnv", "Delphes primary vertex pulls", 50, 50, 900, 500);
         Cnv->Divide(2, 2);
         Cnv->cd(1); gPad->SetLogy(1); gStyle->SetOptFit(1111);