Fork me on GitHub

Changes in / [4e5ed02:c18dca6] in git


Ignore:
Files:
44 deleted
6 edited

Legend:

Unmodified
Added
Removed

  • classes/DelphesClasses.cc

    r4e5ed02 rc18dca6  
    139139  PositionError(0.0, 0.0, 0.0, 0.0),
    140140  Area(0.0, 0.0, 0.0, 0.0),
    141   TrackCovariance(5),
     141  TrackCovarianceACTS(6),
    142142  L(0),
    143143  D0(0), ErrorD0(0),
    144144  DZ(0), ErrorDZ(0),
    145145  P(0), ErrorP(0),
    146   C(0), ErrorC(0),
    147146  PT(0), ErrorPT(0),
    148147  CtgTheta(0), ErrorCtgTheta(0),
    … …  
    314313  object.P = P;
    315314  object.ErrorP = ErrorP;
    316   object.C = C;
    317   object.ErrorC = ErrorC;
    318315  object.PT = PT;
    319316  object.ErrorPT = ErrorPT;
    … …  
    380377  object.SoftDroppedSubJet1 = SoftDroppedSubJet1;
    381378  object.SoftDroppedSubJet2 = SoftDroppedSubJet2;
    382   object.TrackCovariance = TrackCovariance;
     379  object.TrackCovarianceACTS = TrackCovarianceACTS;
    383380  object.fFactory = fFactory;
    384381  object.fArray = 0;
    … …  
    437434  InitialPosition.SetXYZT(0.0, 0.0, 0.0, 0.0);
    438435  Area.SetXYZT(0.0, 0.0, 0.0, 0.0);
    439   TrackCovariance.Zero();
     436  TrackCovarianceACTS.Zero();
    440437  L = 0.0;
    441438  ErrorT = 0.0;
    … …  
    446443  P = 0.0;
    447444  ErrorP = 0.0;
    448   C = 0.0;
    449   ErrorC = 0.0;
    450445  PT = 0.0;
    451446  ErrorPT = 0.0;

  • classes/DelphesClasses.h

    r4e5ed02 rc18dca6  
    651651
    652652  Float_t L; // path length
    653   Float_t DZ;
    654   Float_t ErrorDZ;
    655653  Float_t ErrorT; // path length
    656654  Float_t D0;
    657655  Float_t ErrorD0;
    658   Float_t C;
    659   Float_t ErrorC;
     656  Float_t DZ;
     657  Float_t ErrorDZ;
    660658  Float_t P;
    661659  Float_t ErrorP;
    … …  
    702700  Float_t SumPt;
    703701
    704   // ACTS compliant 6x6 track covariance (D0, phi, Curvature, dz, ctg(theta))
    705 
    706   TMatrixDSym TrackCovariance;
     702  // ACTS compliant 6x6 track covariance (D, z0, phi0, theta, q/p, time)
     703
     704  TMatrixDSym TrackCovarianceACTS;
    707705
    708706  // vertex variables

  • external/TrackCovariance/ObsTrk.cc

    r4e5ed02 rc18dca6  
    2020        fGenPar.ResizeTo(5);
    2121        fGenParACTS.ResizeTo(6);
    22         fGenParILC.ResizeTo(5);
    2322        fObsPar.ResizeTo(5);
    2423        fObsParACTS.ResizeTo(6);
    25         fObsParILC.ResizeTo(5);
    2624        fCov.ResizeTo(5, 5);
    2725        fCovACTS.ResizeTo(6, 6);
    28         fCovILC.ResizeTo(5, 5);
    2926        fGenPar = XPtoPar(x,p,Q);
    3027        fGenParACTS = ParToACTS(fGenPar);
    31         fGenParILC = ParToILC(fGenPar);
    3228        /*
    3329        std::cout << "ObsTrk::ObsTrk: fGenPar";
    … …  
    3733        fObsPar = GenToObsPar(fGenPar, fGC);
    3834        fObsParACTS = ParToACTS(fObsPar);
    39         fObsParILC = ParToILC(fObsPar);
    4035        fObsX = ParToX(fObsPar);
    4136        fObsP = ParToP(fObsPar);
    4237        fObsQ = ParToQ(fObsPar);
    4338        fCovACTS = CovToACTS(fCov);
    44         fCovILC = CovToILC(fCov);
    4539}
    4640//
    … …  
    212206}
    213207
    214 // Parameter conversion to ILC format
    215 TVectorD ObsTrk::ParToILC(TVectorD Par)
    216 {
    217         TVectorD pILC(5);       // Return vector
    218         //
    219         pILC(0) = Par(0)*1.0e3;                 // d0 in mm
    220         pILC(1) = Par(1);                               // phi0 is unchanged
    221         pILC(2) = -2 * Par(2)*1.0e-3;   // w in mm^-1
    222         pILC(3) = Par(3)*1.0e3;                 // z0 in mm
    223         pILC(4) = Par(4);                               // tan(lambda) = cot(theta)
    224         //
    225         return pILC;
    226 }
    227 // Covariance conversion to ILC format
    228 TMatrixDSym ObsTrk::CovToILC(TMatrixDSym Cov)
    229 {
    230         TMatrixDSym cILC(5); cILC.Zero();
    231         //
    232         // Fill derivative matrix
    233         TMatrixD A(5, 5);       A.Zero();
    234         //
    235         A(0, 0) = 1.0e3;                // D-d0 in mm
    236         A(1, 1) = 1.0;          // phi0-phi0
    237         A(2, 2) = -2.0e-3;      // w-C
    238         A(3, 3) = 1.0e3;                // z0-z0 conversion to mm
    239         A(4, 4) = 1.0;          // tan(lambda) - cot(theta)
    240         //
    241         TMatrixDSym Cv = Cov;
    242         TMatrixD At(5, 5);
    243         At.Transpose(A);
    244         Cv.Similarity(At);
    245         cILC = Cv;
    246         //
    247         return cILC;
    248 }
    249 
    250208
    251209

  • external/TrackCovariance/ObsTrk.h

    r4e5ed02 rc18dca6  
    2929        TVector3 fObsP;                                 // Observed  track momentum @ track minimum approach
    3030        TVectorD fGenPar;                               // Generated helix track parameters (D, phi0, C, z0, cot(th))
    31         TVectorD fGenParACTS;                   // Generated helix track parameters (D, z0, phi0, th, q/p, time
    32         TVectorD fGenParILC;                            // Generated helix track parameters (w, phi0, d0, z0, tan(lambda))
     31        TVectorD fGenParACTS;                   // Generated helix track parameters (D, z0, phi0, th, q/p, time)
    3332        TVectorD fObsPar;                               // Observed  helix track parameters (D, phi0, C, z0, cot(th))
    34         TVectorD fObsParACTS;                   // Observed  helix track parameters (D, z0, phi0, th, q/p, time
    35         TVectorD fObsParILC;                            // Observed  helix track parameters (d0, phi0, w, z0, tan(lambda))
     33        TVectorD fObsParACTS;                   // Observed  helix track parameters (D, z0, phi0, th, q/p, time)
    3634        TMatrixDSym fCov;                               // INterpolated covariance of track parameters
    3735        TMatrixDSym fCovACTS;                   // Covariance of track parameters in ACTS format
    3836                                                                        // (D, z0, phi0, theta, q/p, time)
    39         TMatrixDSym fCovILC;                            // Covariance of track parameters in ILC format
    40                                                                         // (d0, phi0, w, z0, tan(lambda))
    4137        //
    4238        // Conversion to ACTS parametrization
    4339        //
    4440        TVectorD ParToACTS(TVectorD Par);               // Parameter conversion
    45         TMatrixDSym CovToACTS(TMatrixDSym Cov); // Covariance
    46         //
    47         // Conversion to ILC parametrization
    48         //
    49         TVectorD ParToILC(TVectorD Par);                // Parameter conversion
    50         TMatrixDSym CovToILC(TMatrixDSym Cov);  // Covariance conversion
     41        TMatrixDSym CovToACTS(TMatrixDSym Cov); // Covariance conversion
    5142        //
    5243public:
    … …  
    7768        // D, z0, phi0, theta, q/p, time
    7869        TVectorD GetGenParACTS()        { return fGenParACTS; }
    79         // d0, phi0, w, z0, tan(lambda)
    80         TVectorD GetGenParILC() { return fGenParILC; }
    8170        // Observed level X, P, Q
    8271        Double_t GetObsQ()      { return fObsQ; }
    … …  
    8776        // D, z0, phi0, theta, q/p, time
    8877        TVectorD GetObsParACTS()        { return fObsParACTS; }
    89         // d0, phi0, w, z0, tan(lambda)
    90         TVectorD GetObsParILC() { return fObsParILC; }
    9178        // Covariances
    9279        TMatrixDSym GetCov(){ return fCov; }
    93         TMatrixDSym GetCovACTS(){ return fCovACTS; }
    94         TMatrixDSym GetCovILC(){ return fCovILC; }
     80        TMatrixDSym GetCovACTS(){ return fCovACTS; };
    9581};
    9682

  • modules/FastJetFinder.cc

    r4e5ed02 rc18dca6  
    489489      fastjet::Filter trimmer(fastjet::JetDefinition(fastjet::kt_algorithm, fRTrim), fastjet::SelectorPtFractionMin(fPtFracTrim));
    490490      fastjet::PseudoJet trimmed_jet = trimmer(*itOutputList);
    491      
     491
     492      trimmed_jet = join(trimmed_jet.constituents());
     493
    492494      candidate->TrimmedP4[0].SetPtEtaPhiM(trimmed_jet.pt(), trimmed_jet.eta(), trimmed_jet.phi(), trimmed_jet.m());
    493495

  • modules/TrackCovariance.cc

    r4e5ed02 rc18dca6  
    9696  Candidate *candidate, *mother;
    9797  Double_t mass, p, pt, q, ct;
    98   Double_t dd0, ddz, dphi, dct, dp, dpt, dC;
     98  Double_t dd0, ddz, dphi, dct, dp, dpt;
    9999 
    100100
    … …  
    115115    candidate->InitialPosition.SetXYZT(track.GetObsX().X(),track.GetObsX().Y(),track.GetObsX().Z(),candidatePosition.T());
    116116   
    117     // save full covariance 5x5 matrix internally (D0, phi, Curvature, dz, ctg(theta))
    118     candidate->TrackCovariance = track.GetCov();
     117    // save full ACTS friendly covariance matrix internally (6x6: D, z0, phi0, theta, q/p, time)
     118    candidate->TrackCovarianceACTS = track.GetCovACTS();
    119119
    120120    pt = candidate->Momentum.Pt();
    … …  
    128128   
    129129    candidate->D0 = track.GetObsPar()[0];
     130    candidate->DZ = track.GetObsPar()[3];
     131    candidate->P  = track.GetObsP().Mag();
     132    candidate->CtgTheta = track.GetObsPar()[4];
    130133    candidate->Phi = track.GetObsPar()[1];
    131     candidate->C = track.GetObsPar()[2];
    132     candidate->DZ = track.GetObsPar()[3];
    133     candidate->CtgTheta = track.GetObsPar()[4];
    134     candidate->P  = track.GetObsP().Mag();
     134
    135135    candidate->PT = pt;
    136136    candidate->Charge = q;
    … …  
    142142    dpt       = 2 * TMath::Sqrt( track.GetCov()(2, 2))*pt*pt / (0.2998*fBz);
    143143    dp        = TMath::Sqrt((1.+ct*ct)*dpt*dpt + 4*pt*pt*ct*ct*dct*dct/(1.+ct*ct)/(1.+ct*ct));
    144     dC        = TMath::Sqrt(track.GetCov()(2, 2));
    145144
    146145    candidate->ErrorD0 = dd0;
    147146    candidate->ErrorDZ = ddz;
    148147    candidate->ErrorP = dp;
    149     candidate->ErrorC = dC;
    150148    candidate->ErrorCtgTheta = dct;
    151149    candidate->ErrorPhi = dphi;
Note: See TracChangeset for help on using the changeset viewer.