StRoot: StHighptPool/macros/common/Bin.h Source File

00001 /*
00002   standard histogram bins and ranges
00003 */
00004 
00005 #ifndef Bin_H
00006 #define Bin_H
00007 
00008 #include "TObject.h"
00009 #include "TArrayD.h"
00010 #include <cmath>
00011 #include <cstdlib>
00012 
00013 namespace Bin
00014 {
00015   //
00016   // event stuff
00017   //
00018 
00019   // 0-9 9 is overflow(most central)
00020   const Double_t flowCentMin = -0.5;
00021   const Double_t flowCentMax = 9.5;
00022   const Int_t    nFlowCentBin = 10;
00023   
00024   const Double_t vertexZEvtMin = -200;
00025   const Double_t vertexZEvtMax = 200;
00026   const Int_t    nVertexZEvtBin = 
00027     (Int_t) ceil((vertexZEvtMax-vertexZEvtMin)/5);
00028   const Int_t    nVertexZEvtThinBin = (Int_t)ceil((vertexZEvtMax-vertexZEvtMin)/1);
00029 
00030   // 0-9 0 is most central (kFive)
00031    const Double_t zdcCentMin = -0.5; 
00032    const Double_t zdcCentMax = 9.5;
00033    const Int_t    nZdcCentBin = 10;
00034 
00035    const Double_t vertexXMin = -1;
00036    const Double_t vertexXMax = 1;
00037    const Int_t nVertexXBin = (Int_t)ceil((vertexXMax-vertexXMin)/.05);
00038 
00039    const Double_t vertexREvtMin = 0;
00040    const Double_t vertexREvtMax = 5;
00041    const Int_t nVertexREvtBin = (Int_t)ceil((vertexREvtMax-vertexREvtMin)/.1);
00042 
00043    const Double_t zdcMin=-0.5;
00044    const Double_t zdcMax=200.5;
00045    const Int_t nZdcBin=(Int_t)ceil((zdcMax-zdcMin));
00046    
00047    const Double_t hMinusMin=-0.5;
00048    const Double_t hMinusMax=400.5;
00049    const Int_t nHMinusBin=(Int_t)ceil((hMinusMax-hMinusMin));
00050 
00051    const Double_t nChMin=-0.5;
00052    const Double_t nChMax=800.5;
00053    const Int_t nNchBin=(Int_t)ceil((nChMax-nChMin));
00054 
00055 
00056    const Double_t ctbMin = -0.5;
00057    const Double_t ctbMax = 30000.5;
00058    const Int_t nCtbBin=(Int_t)ceil((ctbMax-ctbMin)/100);
00059 
00060   //
00061   //****** track stuff
00062   //
00063   void initPtAry(TArrayD* a, int type);
00064 
00065    const Double_t ptMin = 0.5;
00066    const Double_t ptMax = 8.0;
00067    const Int_t    nPtBin = (Int_t)ceil((ptMax-ptMin)/.5);
00068   
00069    const Double_t ptThinMin = 1.2;
00070    const Double_t ptThinMax = 8.0;
00071    const Int_t    nPtThinBin = (Int_t)ceil((ptThinMax-ptThinMin)/.2);
00072 
00073    const Double_t ptTinyMin = 0;
00074    const Double_t ptTinyMax = 8;
00075    const Int_t nPtTinyBin = (Int_t)ceil((ptTinyMax-ptTinyMin)/.2);
00076 
00077    // will rebin with these
00078    const double ptRebinMin = 1;
00079    const double ptRebinMax = 6;
00080    const double nPtRebinBin = (int)ceil((ptRebinMax-ptRebinMin)/.05);
00081    
00082 
00083    const Double_t resPtMin = -1.5;
00084    const Double_t resPtMax = 1.5;
00085    const Int_t nResPtBin = (Int_t)ceil((resPtMax-resPtMin)/.005);
00086 
00087    const Double_t resEtaMin = -.5;
00088    const Double_t resEtaMax = .5;
00089    const Int_t nResEtaBin = (Int_t)ceil((resEtaMax-resEtaMin)/.02);
00090 
00091    const Double_t resCurvMin = -2;
00092    const Double_t resCurvMax = 1;
00093    const Int_t nResCurvBin = (Int_t)ceil((resCurvMax-resCurvMin)/.02);
00094                           
00095    const Double_t etaMin = -1.0;
00096    const Double_t etaMax = 1.0;
00097    const Int_t nEtaBin = (Int_t)ceil((etaMax-etaMin)/.1);
00098 
00099    const Double_t etaBigMin = -1.0;
00100    const Double_t etaBigMax = 1.0;
00101    const Int_t nEtaBigBin = (Int_t)ceil((etaBigMax-etaBigMin)/.25);
00102 
00103    const Double_t etaThinMin = -1.;
00104    const Double_t etaThinMax = 1;
00105    const Int_t    nEtaThinBin = (Int_t)ceil((etaThinMax-etaThinMin)/.025);
00106 
00107    const Double_t etaSmallMin = -1.0;
00108    const Double_t etaSmallMax = 1.0;
00109    const Int_t nEtaSmallBin = (Int_t)ceil((etaSmallMax-etaSmallMin)/.05); 
00110 
00111    const Double_t sDcaMin = -3;
00112    const Double_t sDcaMax = 3;
00113    const Int_t nSDcaBin = (Int_t)ceil((sDcaMax-sDcaMin)/.1);
00114 
00115    const Double_t dcaXYGlMin = -3;
00116    const Double_t dcaXYGlMax = 3;
00117    const Double_t nDcaXYGlBin =(Int_t)ceil((dcaXYGlMax-dcaXYGlMin)/.05);
00118 
00119    const Double_t dcaXYGlWideMin = -3;
00120    const Double_t dcaXYGlWideMax = 3;
00121    const Double_t nDcaXYGlWideBin =(Int_t)ceil((dcaXYGlMax-dcaXYGlMin)/.1);
00122 
00123 
00124   
00125    const Double_t dcaMin = 0;
00126    const Double_t dcaMax = 3;
00127    const Int_t nDcaBin = (Int_t)ceil((dcaMax-dcaMin)/.05);
00128 
00129    const Double_t dcaThinMin = 0;
00130    const Double_t dcaThinMax = 3;
00131    const Int_t nDcaThinBin = (Int_t)ceil((dcaThinMax-dcaThinMin)/.05);
00132 
00133    const Double_t commonFracMin = 0;
00134    const Double_t commonFracMax = 1;
00135    const Int_t nCommonFracBin = (Int_t)ceil((commonFracMax-commonFracMin)/.02);
00136 
00137   const Double_t fitPtsMin = 0.5;
00138   const Double_t fitPtsMax = 45.5;
00139   const Int_t nFitPtsBin = 45;
00140 
00141   const Double_t fitPtsWideMin=25.5; // 26-29,30-33,34-37,38-41,42-45
00142   const Double_t fitPtsWideMax=45.5;
00143   const Int_t nFitPtsWideBin=5;
00144 
00145   const Double_t fitPtsCentMin = 0.5;
00146   const Double_t fitPtsCentMax = 45.5;
00147   const Int_t nFitPtsCentBin = 45;
00148 
00149    const Double_t fitHitMin = 0.5;
00150    const Double_t fitHitMax = 45.5;
00151    const Int_t nFitHitBin = 1;
00152 
00153    const Double_t fracHitMin = 0;
00154    const Double_t fracHitMax = 1;
00155    const Int_t nFracHitBin = (Int_t)ceil((fracHitMax-fracHitMin)/.02);
00156 
00157    const Double_t fitHitThinMin = 0.5;
00158    const Double_t fitHitThinMax = 45.5;
00159    const Int_t nFitHitThinBin = 45;
00160   
00161    const Double_t vertexZMin = -200;
00162    const Double_t vertexZMax = 200;
00163    const Int_t nVertexZBin = (Int_t)ceil((vertexZMax-vertexZMin)/10);
00164    const Int_t nEnterZBin  = (Int_t)ceil((vertexZMax-vertexZMin)/10);
00165 
00166   const Double_t vertexZWideMin = -200;
00167   const Double_t vertexZWideMax = 200;
00168   const Int_t nVertexZWideBin = (Int_t)ceil((vertexZWideMax-vertexZWideMin)/25);
00169 
00170   const Double_t vertexZBigMin = -200;
00171   const Double_t vertexZBigMax = 200;
00172   const Int_t nVertexZBigBin = (Int_t) ceil((vertexZBigMax-vertexZBigMin)/50);
00173     
00174    const Double_t midZBigMin = -200;
00175    const Double_t midZBigMax = 200;
00176    const Int_t nMidZBigBin = (Int_t) ceil((midZBigMax-midZBigMin)/50);
00177 
00178    const Double_t absZBigMin = 0;
00179    const Double_t absZBigMax = 200;
00180    const Double_t nAbsZBigBin = 4;
00181   
00182    const Double_t sectorWestMin = 0.5;
00183    const Double_t sectorWestMax = 12.5;
00184 
00185    const Double_t sectorEastMin = 12.5;
00186    const Double_t sectorEastMax = 24.5;
00187   
00188    const Int_t nSectorBin = 12;
00189 
00190 
00191    const Double_t dPhiMin = -.02;
00192    const Double_t dPhiMax = .02;
00193    const Int_t nDPhiBin = (Int_t)ceil((dPhiMax-dPhiMin)/.001);
00194 
00195    const Double_t dEtaMin = -.1;
00196    const Double_t dEtaMax = .1;
00197    const Int_t nDEtaBin = (Int_t)ceil((dEtaMax-dEtaMin)/.001);
00198 
00199   //   const Double_t dPtMin = -6;
00200   //  const Double_t dPtMax = 6;
00201   //  const Int_t nDPtBin = (Int_t)ceil((dPtMax-dPtMin)/.2);
00202 
00203    const Double_t phiMin = -3.14159;
00204    const Double_t phiMax = 3.14159;
00205    const Double_t phiDegMin = -165;
00206    const Double_t phiDegMax = 195;
00207    const Int_t nPhiBin = 12;
00208    const Int_t nPhiSmallBin = 24;
00209    const Int_t nPhiThinBin = 24*5;
00210 
00211    const Double_t dipMin = -.5;
00212    const Double_t dipMax = .5;
00213    const Int_t nDipBin = (Int_t)ceil((dipMax-dipMin)/.02);
00214 
00215    const Double_t dipWideMin = -.5;
00216    const Double_t dipWideMax = .5;
00217    const Int_t nDipWideBin = (Int_t)ceil((dipMax-dipMin)/.1);
00218 
00219    const Double_t pzMin = 5;
00220    const Double_t pzMax = -5;
00221    const Int_t nPzBin = (Int_t)ceil((pzMax-pzMin)/1);
00222 
00223    const Double_t dcaZGlMin = -.5;
00224    const Double_t dcaZGlMax = .5;
00225    const Int_t nDcaZGlBin = (Int_t)ceil((dcaZGlMax-dcaZGlMin)/.02);
00226 
00227    const Double_t dcaZGlWideMin = -.5;
00228    const Double_t dcaZGlWideMax = .5;
00229    const Int_t nDcaZGlWideBin = (Int_t)ceil((dcaZGlWideMax-dcaZGlWideMin)/.1);
00230 
00231    const Double_t dDipMin = -.003;
00232    const Double_t dDipMax = .003;
00233    const Int_t nDDipBin = (Int_t)ceil((dDipMax-dDipMin)/.0001);
00234 
00235    const Double_t dcaXYPrMin = -.2;
00236    const Double_t dcaXYPrMax = .2;
00237    const Int_t    nDcaXYPrBin = (Int_t)ceil((dcaXYPrMax-dcaXYPrMin)/.01);
00238 
00239    const Double_t dedxMin = 0;
00240    const Double_t dedxMax = 1e-5;
00241    const Int_t nDedxBin = (Int_t)ceil((dedxMax-dedxMin)/.2e-6);
00242 
00243 }
00244 
00245 
00246 
00247 
00248 
00249 
00250 #endif