Material distribution from GEANT geometry : HFT, run14

Use StarVMC/StarBASE to make plots (1D,2D) of radiation lengths per volume

preliminary check for the parameters with which the geantinos are thrown.
The eta and phi bin, as well as the # of geantinos per bin can be changed.

1) Default parameters.

Geantinos are thrown according :
const Float_t eta_min = -6.0;     // min pseudorapidity
const Float_t eta_max = +6.0;     // max pseudorapidity
const Float_t deta    = .1;       // size of etabin
const Float_t phi_min = -180.000; // min phi (in degrees)
const Float_t phi_max = +180.000; // max phi (in degrees)
const Float_t dphi    = 1;        // size of phibin (in degrees)
const Float_t nsample = 4.0;      // number of tracks / bin
Fig 1: rad length vs. phi (degree), z[cm], radius [cm], radius vs phi for PXL (PXMO) volume

Resolution is not good because it depends on the histograms binning too.



Explanation for the left plot of Fig1
Fig 2 : breakdown for PXMO volume
  • left : total radiation length per phi bin (with an eta bin = 0.1)
  • middle : # of geantinos per phi bin (with an eta bin = 0.1)
  • rigth : ratio left/middle

Explanation of the Y-axis
setup is now 
const Float_t eta_min = -6.0;     // min pseudorapidity
const Float_t eta_max = +6.0;     // max pseudorapidity
const Float_t deta    = 3;       // size of etabin
const Float_t phi_min = -180.000; // min phi (in degrees)
const Float_t phi_max = +180.000; // max phi (in degrees)
const Float_t dphi    = 1;        // size of phibin (in degrees)
const Float_t nsample = 4.0;      // number of tracks / bin
Fig 3 : deta=3, dphi=1

The number of entries makes sense : 4 (bin eta * range eta) * 360 (bin phi * range phi) *4 (# of geantinos) = 5760

Going back to the initial setup : 
Fig 4 : deta=0.1 , dphi =1