Difference between revisions of "Constants to Review"
(→Reconstruction constants) |
(→Reconstruction constants) |
||
Line 42: | Line 42: | ||
gas_gain=1e5; | gas_gain=1e5; | ||
hit->dist = hit->tdrift*55.0E-4; // Use number hardwired in simulation for now | hit->dist = hit->tdrift*55.0E-4; // Use number hardwired in simulation for now | ||
+ | |||
+ | '''FDC''' | ||
+ | |||
+ | FDC/DFDCGeometry.h has a bunch of constants that need to be reconciled with HDGeant/hitFDC.c | ||
'''FCAL''' | '''FCAL''' |
Revision as of 17:03, 19 May 2014
These are constants whose values need to be revisited
Digitization-level constants
CDC
a_scale = 1.0E6/1.3E5
FDC
a_scale = 2.4E4/1.3E5; // for cathodes, what about wires?
FCAL
a_scale = 4.0E1/2.5E5
BCAL
a_scale = 0.0001; // 100 keV/count (?)
TOF
a_scale = 0.2/5.2E5
ADC: hit->sigma_t = 4.0; // ns (what is the fADC time resolution?)
RDC: hit->sigma_t = 0.160; // ns (what is the TOF TDC time resolution?)
t0 offset for fADC?
Start Counter
a_scale = 2.0E-2/5.2E-5
ADC: hit->sigma_t = 4.0; // ns (what is the fADC time resolution?)
RDC: hit->sigma_t = 0.160; // ns (what is the TOF TDC time resolution?)
t0 offset for fADC?
Reconstruction constants
CDC/DCDCTrackHit_factory.h
w_eff=29.5e-9; gas_gain=1e5; hit->dist = hit->tdrift*55.0E-4; // Use number hardwired in simulation for now
FDC
FDC/DFDCGeometry.h has a bunch of constants that need to be reconciled with HDGeant/hitFDC.c
FCAL
Various constants in DFCALCluster::shower_profile()
BCAL
DBCALCluster.cc: // The method below for determining sig_theta works better than the one // above. sigma_z is determined using errors when reconstructing MC data. double sigma_z = sqrt(1.394*1.394/m_E + 0.859*0.859); m_sig_theta = sigma_z*sin(m_theta)*sin(m_theta)/DBCALGeometry::BCALINNERRAD; DBCALCluster_factory.cc: In DBCALCluster_factory::evnt(), energy threshold is 5*k_MeV In DBCALCluster_factory::clusterize() float seedThresh = 1*k_GeV; float minSeed = 10*k_MeV; In DBCALCluster_factory::overlap( const DBCALCluster& highEClust, const DBCALCluster& lowEClust ) //deltaPhi_force_merge and delta_z_force_merge were determined by looking //at the separation of decay photons from pi0's from a pythia sample. //There are no events where the decay photons have separation //(delta_phi < 0.2 && delta_z < 25 cm), so in most cases it should be safe //to merge clusters together if they are so close. const double deltaPhi_force_merge = 0.1; //radians const double delta_z_force_merge = 15.0*k_cm; //A major cause of extra clusters are lower energy hits, which have poor //z-resolution and so are not properly merged. Treat low energy //clusters (< 40 MeV) as a special case. Again, hopefully this is only //a temporary fix until we have a more comprehensive solution. const double delta_z_force_merge_low_E = 40.0*k_cm; const double low_E = .04*k_GeV;