StRoot: StAssociationMaker/StMcParameterDB.cxx Source File

00001 /*****************************************
00002  *
00003  * $Id: StMcParameterDB.cxx,v 1.9 2005/11/22 21:44:16 fisyak Exp $
00004  *
00005  * Changed extension to .cxx so that it
00006  * would be accessible from Root macro
00007  *
00008  * $Log: StMcParameterDB.cxx,v $
00009  * Revision 1.9  2005/11/22 21:44:16  fisyak
00010  * Add Ssd to Associator, add IdTruth options for Svt and Ssd
00011  *
00012  * Revision 1.8  2004/03/30 03:02:02  calderon
00013  * -Changed default value of the z cut to 5 mm.
00014  *
00015  * Revision 1.7  2003/09/02 17:55:28  perev
00016  * gcc 3.2 updates + WarnOff
00017  *
00018  * Revision 1.6  2003/06/27 03:01:19  calderon
00019  * The z cut now depends on z_mc.
00020  * The parameterization is done in the parameter DB
00021  * with a linearly increasing rms, symmetric in +/- z.
00022  *
00023  * Revision 1.5  1999/12/14 07:07:41  calderon
00024  * Added Ratio Number of Common Hits / Number of Reconstructed Hits for
00025  * each detector.
00026  * Numbering scheme from StEvent & StMcEvent as per SVT request
00027  * Added Kink, V0 and Xi vertex associations.
00028  *
00029  * Revision 1.4  1999/12/08 00:00:25  calderon
00030  * New version of StAssociationMaker.
00031  * -Uses new StEvent / StMcEvent
00032  * -Includes maps using reconstructed and monte carlo objects as keys for:
00033  *   TPC Hits
00034  *   SVT Hits
00035  *   SSD Hits
00036  *   FTPC Hits
00037  *   Tracks (using all 3 hit multimaps)
00038  *
00039  * Revision 1.3  1999/10/01 14:08:58  calderon
00040  * Added Local Hit resolution Histogram. It is made by default
00041  * without any requirement of association, to serve
00042  * as a diagnostic.
00043  * Before building track multimap, check the size of the
00044  * tpc hit map.  If it is too small, print out a warning
00045  * and exit.
00046  *
00047  * Revision 1.2  1999/09/23 21:25:21  calderon
00048  * Added Log & Id
00049  * Modified includes according to Yuri
00050  *
00051  *****************************************/
00052 #include <Stiostream.h>
00053 #include "StMcParameterDB.h"
00054 #include "PhysicalConstants.h"
00055 #include "SystemOfUnits.h"
00056 
00057 ClassImp(StMcParameterDB)
00058     
00059 StMcParameterDB* StMcParameterDB::mParamDB = 0;
00060 
00061 StMcParameterDB* StMcParameterDB::instance()
00062 {
00063     if (!mParamDB) mParamDB = new StMcParameterDB;
00064     return mParamDB;
00065 }
00066 
00067 StMcParameterDB::StMcParameterDB()
00068     :mXCutTpc(5.*millimeter),mYCutTpc(5.*millimeter),mZCutTpc(5.*millimeter),mReqCommonHitsTpc(3),
00069      mXCutSvt(1.*millimeter),mYCutSvt(1.*millimeter),mZCutSvt(1.*millimeter),mReqCommonHitsSvt(1),
00070      mXCutSsd(1.*millimeter),mYCutSsd(1.*millimeter),mZCutSsd(1.*millimeter),mReqCommonHitsSsd(1),
00071      mRCutFtpc(3.*millimeter),mPhiCutFtpc(5.*degree),mReqCommonHitsFtpc(2)
00072 { /*noop*/ }
00073 
00074 void StMcParameterDB::setXCutTpc(float val) { mXCutTpc = val ;}
00075 void StMcParameterDB::setYCutTpc(float val) { mYCutTpc = val ;}
00076 void StMcParameterDB::setZCutTpc(float val) { mZCutTpc = val ;}
00077 void StMcParameterDB::setReqCommonHitsTpc(unsigned int val) { mReqCommonHitsTpc = val;}
00078 
00079 void StMcParameterDB::setXCutSvt(float val) { mXCutSvt = val ;}
00080 void StMcParameterDB::setYCutSvt(float val) { mYCutSvt = val ;}
00081 void StMcParameterDB::setZCutSvt(float val) { mZCutSvt = val ;}
00082 void StMcParameterDB::setReqCommonHitsSvt(unsigned int val) { mReqCommonHitsSvt = val;}
00083 
00084 void StMcParameterDB::setXCutSsd(float val) { mXCutSsd = val ;}
00085 void StMcParameterDB::setYCutSsd(float val) { mYCutSsd = val ;}
00086 void StMcParameterDB::setZCutSsd(float val) { mZCutSsd = val ;}
00087 void StMcParameterDB::setReqCommonHitsSsd(unsigned int val) { mReqCommonHitsSsd = val;}
00088 
00089 void StMcParameterDB::setRCutFtpc(float val)   { mRCutFtpc   = val ;}
00090 void StMcParameterDB::setPhiCutFtpc(float val) { mPhiCutFtpc = val ;}
00091 void StMcParameterDB::setReqCommonHitsFtpc(unsigned int val) { mReqCommonHitsFtpc = val;}
00092 
00093 
00094 float StMcParameterDB::zCutTpc(float z) const {
00095     return (mZCutTpc/0.166657) * (0.166657+2.07482e-4*fabs(z));
00096 }
00097 
00098 
00099 ostream& operator<<(ostream &os, const StMcParameterDB& mcDb)
00100 {
00101     os << " TPC Cuts " << endl;
00102     os << " X Cut   : " << mcDb.xCutTpc()/millimeter << " mm" << endl; 
00103     os << " Y Cut   : " << mcDb.yCutTpc()/millimeter << " mm" << endl; 
00104     os << " Z Cut   : " << mcDb.zCutTpc()/millimeter << " mm" << endl; 
00105     os << " Required Hits for Associating Tracks: " << mcDb.reqCommonHitsTpc() << endl;
00106     os << " SVT Cuts " << endl;
00107     os << " X Cut   : " << mcDb.xCutSvt()/millimeter << " mm" << endl; 
00108     os << " Y Cut   : " << mcDb.yCutSvt()/millimeter << " mm" << endl; 
00109     os << " Z Cut   : " << mcDb.zCutSvt()/millimeter << " mm" << endl; 
00110     os << " Required Hits for Associating Tracks: " << mcDb.reqCommonHitsSvt() << endl;
00111     os << " SSD Cuts " << endl;
00112     os << " X Cut   : " << mcDb.xCutSsd()/millimeter << " mm" << endl; 
00113     os << " Y Cut   : " << mcDb.yCutSsd()/millimeter << " mm" << endl; 
00114     os << " Z Cut   : " << mcDb.zCutSsd()/millimeter << " mm" << endl; 
00115     os << " Required Hits for Associating Tracks: " << mcDb.reqCommonHitsSsd() << endl;
00116     os << " FTPC Cuts " << endl;
00117     os << " R Cut   : " << mcDb.rCutFtpc()/millimeter << " mm" << endl; 
00118     os << " Phi Cut : " << mcDb.phiCutFtpc()/degree << " degrees" << endl; 
00119     os << " Required Hits for Associating Tracks: " << mcDb.reqCommonHitsFtpc() << endl;
00120     
00121     return os;
00122 }