00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056 #include "StHbtMaker/Cut/helensV0Cut.h"
00057 #include <cstdio>
00058
00059 #ifdef __ROOT__
00060 ClassImp(helensV0Cut)
00061 #endif
00062
00063
00064 helensV0Cut::helensV0Cut(){
00065 mNV0sPassed = mNV0sFailed = 0;
00066
00067 mV0MassRange[0] =0;
00068 mV0MassRange[1]=10000;
00069
00070 mdcaV0daughters[0]=0;
00071 mdcaV0daughters[1]=1000;
00072
00073 mdcaV0ToPrimVertex[0]=0;
00074 mdcaV0ToPrimVertex[1]=1000;
00075
00076 mdecayLengthV0[0]=0;
00077 mdecayLengthV0[1]=10000;
00078
00079 mtpcHitsPos[0]=0;
00080 mtpcHitsPos[1]=1000;
00081
00082 mtpcHitsNeg[0]=0;
00083 mtpcHitsNeg[1]=1000;
00084
00085
00086 mdcaPosToPrimVertex[0]=0;
00087 mdcaPosToPrimVertex[1]=1000;
00088
00089 mdcaNegToPrimVertex[0]=0;
00090 mdcaNegToPrimVertex[1]=10000;
00091
00092 mptArmV0[0]=0;
00093 mptArmV0[1]=100;
00094
00095 malphaV0[0]=-10;
00096 malphaV0[1]=10;
00097
00098 mChargedEdx=0;
00099 mdEdx[0]=0;
00100 mdEdx[1]=-10;
00101 mdEdx[2]=0;
00102 mdEdx[3]=-10;
00103
00104 mPt[0]=0;
00105 mPt[1]=100000;
00106
00107 mRapidity[0]=-100000;
00108 mRapidity[1]=100000;
00109
00110 V0Type = "K0Short";
00111 mMass = 0.498;
00112
00113
00114 }
00115
00116
00117
00118
00119
00120 bool helensV0Cut::Pass(const StHbtV0* V0){
00121 int inMassRange;
00122
00123 #ifdef STHBTDEBUG
00124 cout << endl;
00125 cout << " * dcaV0Daughters " << V0->dcaV0Daughters();
00126 cout << " * dcaV0ToPrimVertex " << V0->dcaV0ToPrimVertex();
00127 cout << " * decayLengthV0 " << V0->decayLengthV0();
00128 cout << " * tpcHitsPos " << V0->tpcHitsPos();
00129 cout << " * tpcHitsNeg " << V0->tpcHitsNeg();
00130 cout << " * dcaPosToPrimVertex " << V0->dcaPosToPrimVertex();
00131 cout << " * dcaNegToPrimVertex " << V0->dcaNegToPrimVertex();
00132 cout << " * ptArmV0 " << V0->ptArmV0();
00133 cout << " * alphaV0 " << V0->alphaV0();
00134 cout << " * dEdxPos " << V0->dedxPos();
00135 cout << " * dEdxNeg " << V0->dedxNeg();
00136 cout << endl;
00137 #endif
00138
00139
00140 inMassRange=0;
00141
00142
00143 if( strstr(V0Type,"k") || strstr(V0Type,"K")){
00144 if( V0->massK0Short() < (mV0MassRange[1]) &&
00145 V0->massK0Short() > (mV0MassRange[0]) ) inMassRange=1;
00146 }
00147 else if( (strstr(V0Type,"anti") || strstr(V0Type,"ANTI"))){
00148 if( V0->massAntiLambda() < (mV0MassRange[1]) &&
00149 V0->massAntiLambda() > (mV0MassRange[0]) ) inMassRange=1;
00150 }
00151 else if( (strstr(V0Type,"ambda") || strstr(V0Type,"AMBDA"))){
00152 if( V0->massLambda() < (mV0MassRange[1]) &&
00153 V0->massLambda() > (mV0MassRange[0]) ) inMassRange=1;
00154 }
00155
00156
00157
00158 bool goodPID = ( inMassRange &&
00159 (V0->dcaV0Daughters() > mdcaV0daughters[0]) &&
00160 (V0->dcaV0Daughters() < mdcaV0daughters[1]) &&
00161 (V0->dcaV0ToPrimVertex() > mdcaV0ToPrimVertex[0]) &&
00162 (V0->dcaV0ToPrimVertex() < mdcaV0ToPrimVertex[1]) &&
00163 (V0->decayLengthV0() > mdecayLengthV0[0]) &&
00164 (V0->decayLengthV0() < mdecayLengthV0[1]) &&
00165 (V0->tpcHitsPos() > mtpcHitsPos[0]) &&
00166 (V0->tpcHitsPos() < mtpcHitsPos[1]) &&
00167 (V0->tpcHitsNeg() > mtpcHitsNeg[0]) &&
00168 (V0->tpcHitsNeg() < mtpcHitsNeg[1]) &&
00169 (V0->dcaPosToPrimVertex() > mdcaPosToPrimVertex[0]) &&
00170 (V0->dcaPosToPrimVertex() < mdcaPosToPrimVertex[1]) &&
00171 (V0->dcaNegToPrimVertex() > mdcaNegToPrimVertex[0]) &&
00172 (V0->dcaNegToPrimVertex() < mdcaNegToPrimVertex[1]) &&
00173 (V0->ptArmV0() > mptArmV0[0]) &&
00174 (V0->ptArmV0() < mptArmV0[1]) &&
00175 (V0->alphaV0() > malphaV0[0]) &&
00176 (V0->alphaV0() < malphaV0[1]));
00177
00178 if(goodPID && mChargedEdx){
00179 if( mChargedEdx <0){
00180 goodPID = ( (V0->dedxNeg() > (mdEdx[0]*V0->ptNeg()+mdEdx[1])) &&
00181 (V0->dedxNeg() > (mdEdx[2]*V0->ptNeg()+mdEdx[3])));
00182 }
00183 if( mChargedEdx > 0){
00184 goodPID = ( (V0->dedxPos() > (mdEdx[0]*V0->ptPos()+mdEdx[1])) &&
00185 (V0->dedxPos() > (mdEdx[2]*V0->ptPos()+mdEdx[3])));
00186 }
00187 }
00188
00189 if (goodPID){
00190 float TEnergy = ::sqrt((V0->ptotV0())*(V0->ptotV0())+mMass*mMass);
00191 float TRapidity = 0.5*::log((TEnergy+V0->momV0().z())/
00192 (TEnergy-V0->momV0().z()));
00193
00194 float Pt = V0->ptV0();
00195
00196
00197
00198 #ifdef STHBTDEBUG
00199 cout << " * Pt " << Pt;
00200 cout << " * mPt[0] " << mPt[0];
00201 cout << " * mPt[1] " << mPt[1];
00202 cout << " * TRapidity " << TRapidity;
00203 cout << " * mRapidity[0] " << mRapidity[0];
00204 cout << " * mRapidity[1] " << mRapidity[1];
00205 cout << " * Pt " << (Pt > mPt[0]) && (Pt < mPt[1]);
00206 cout << " * y " << (TRapidity > mRapidity[0]) && (TRapidity < mRapidity[1]);
00207 cout << endl;
00208 #endif
00209
00210 bool goodV0=
00211 ((Pt > mPt[0]) &&
00212 (Pt < mPt[1]) &&
00213 (TRapidity > mRapidity[0]) &&
00214 (TRapidity < mRapidity[1]));
00215
00216 goodV0 ? mNV0sPassed++ : mNV0sFailed++;
00217 return (goodV0);
00218 }
00219 else{
00220 mNV0sFailed++;
00221 return (goodPID);
00222 }
00223 }
00224
00225 StHbtString helensV0Cut::Report(){
00226 string Stemp;
00227 char Ctemp[100];
00228 sprintf(Ctemp,"--helensV0Cut--\n Particle mass:\t%E\n",this->Mass());
00229 Stemp=Ctemp;
00230 sprintf(Ctemp,"V0 mass range:\t%E - %E\n",mV0MassRange[0],
00231 mV0MassRange[1]);
00232 Stemp+=Ctemp;
00233 sprintf(Ctemp,"dcaV0daughters:\t%E - %E\n",mdcaV0daughters[0],
00234 mdcaV0daughters[1]);
00235 Stemp+=Ctemp;
00236 sprintf(Ctemp,"dcaV0ToPrimVertex:\t%E - %E\n",mdcaV0ToPrimVertex[0],
00237 mdcaV0ToPrimVertex[1]);
00238 Stemp+=Ctemp;
00239 sprintf(Ctemp,"decayLengthV0:\t%E - %E\n",mdecayLengthV0[0],
00240 mdecayLengthV0[1]);
00241 Stemp+=Ctemp;
00242 sprintf(Ctemp,"tpcHitsPos:\t%d - %d\n",mtpcHitsPos[0],mtpcHitsPos[1]);
00243 Stemp+=Ctemp;
00244 sprintf(Ctemp,"tpcHitsNeg:\t%d - %d\n",mtpcHitsNeg[0],mtpcHitsNeg[1]);
00245 Stemp+=Ctemp;
00246 sprintf(Ctemp,"dcaPosToPrimVertex:\t%E - %E\n",mdcaPosToPrimVertex[0],
00247 mdcaPosToPrimVertex[1]);
00248 Stemp+=Ctemp;
00249 sprintf(Ctemp,"dcaNegToPrimVertex:\t%E - %E\n",mdcaNegToPrimVertex[0],
00250 mdcaNegToPrimVertex[1]);
00251 Stemp+=Ctemp;
00252 sprintf(Ctemp,"dedx>:\t%E pt+%E and \t%E pt+ %E for Charge %E\n ",mdEdx[0],mdEdx[1],mdEdx[2],mdEdx[3],mChargedEdx);
00253 Stemp+=Ctemp;
00254 sprintf(Ctemp,"ptArmV0:\t%E - %E\n",mptArmV0[0],mptArmV0[1]);
00255 Stemp+=Ctemp;
00256 sprintf(Ctemp,"alphaV0:\t%E - %E\n",malphaV0[0],malphaV0[1]);
00257 Stemp+=Ctemp;
00258 sprintf(Ctemp,"Particle pT:\t%E - %E\n",mPt[0],mPt[1]);
00259 Stemp+=Ctemp;
00260 sprintf(Ctemp,"Particle rapidity:\t%E - %E\n",mRapidity[0],mRapidity[1]);
00261 Stemp+=Ctemp;
00262 sprintf(Ctemp,"Number of V0s which passed:\t%ld Number which failed:\t%ld\n",mNV0sPassed,mNV0sFailed);
00263 Stemp += Ctemp;
00264 StHbtString returnThis = Stemp;
00265 return returnThis;
00266 }
