Changeset 3b3071a in git for external/TrackCovariance/TrkUtil.h

Timestamp:

Apr 17, 2021, 10:28:12 AM (4 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

master

Children:

27197df

Parents:

9e08f27 (diff), 64294db (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.

Message:

Merge branch 'master' of ​https://github.com/delphes/delphes

File:

• external/TrackCovariance/TrkUtil.h (modified) (4 diffs)

external/TrackCovariance/TrkUtil.h

@@ -6 +6 @@
 #include <TVectorD.h>
 #include <TMatrixDSym.h>
+#include <TRandom.h>
 //
 //
@@ -15 +16 @@
 protected:
         Double_t fBz;                                                   // Solenoid magnetic field
+        //
+        Int_t fGasSel;                                                  // Gas selection: 0: He-Iso, 1: He, 2:Ar-Eth, 3: Ar
+        Double_t fRmin;                                                 // Lower                DCH radius
+        Double_t fRmax;                                                 // Higher       DCH radius
+        Double_t fZmin;                                                 // Lower                DCH z
+        Double_t fZmax;                                                 // Higher       DCH z
         //
         // Service routines
@@ -31 +38 @@
         TVectorD ParToILC(TVectorD Par);                // Parameter conversion
         TMatrixDSym CovToILC(TMatrixDSym Cov);  // Covariance conversion
+        //
 
 public:
@@ -60 +68 @@
         static TVectorD ParToMm(TVectorD Par);                  // Parameter conversion
         static TMatrixDSym CovToMm(TMatrixDSym Cov);    // Covariance conversion
-
+        //
+        // Cluster counting in gas
+        //
+        // Define gas volume (units = meters) 
+        void SetDchBoundaries(Double_t Rmin, Double_t Rmax, Double_t Zmin, Double_t Zmax);
+        // Gas mixture selection
+        void SetGasMix(Int_t Opt);
+        // Get number of ionization clusters
+        Bool_t IonClusters(Double_t &Ncl, Double_t mass, TVectorD Par);
+        Double_t Nclusters(Double_t bgam);      // mean clusters/meter vs beta*gamma
+        static Double_t Nclusters(Double_t bgam, Int_t Opt);    // mean clusters/meter vs beta*gamma
+        Double_t funcNcl(Double_t *xp, Double_t *par);
+        Double_t TrkLen(TVectorD Par);                                  // Track length inside chamber
 };