StRoot: StHbtMaker/ThCorrFctn/StHbtThPairGaussFit.cxx Source File

00001 /***************************************************************************
00002  *
00003  * $Id: 
00004  *
00005  * Author: Laurent Conin, Fabrice Retiere, Subatech, France
00006  ***************************************************************************
00007  *
00008  * Description : Implementaton of StHbtThPairGaussFit
00009  *
00010  ***************************************************************************
00011  *
00012  * $Log: 
00013  *
00014  ***************************************************************************/
00015 
00016 #include "StHbtMaker/ThCorrFctn/StHbtThPairGaussFit.h"
00017 
00018 StHbtThPairGaussFit::StHbtThPairGaussFit(): StHbtThPairGauss() 
00019 { /* no-op */};  
00020 
00021 StHbtThPairGaussFit::~StHbtThPairGaussFit()
00022 { /* no-op */};  
00023 
00024 void StHbtThPairGaussFit::Set(const StHbtPair* aPair){
00025   SetMomentum_PID (aPair);
00026   SetPosition();
00027   if (mRef!=RCMS) {
00028     mSourceDist1=*mEmPoint1;
00029     mSourceDist2=*mEmPoint2;
00030     BoostPosition();
00031   }
00032   mMeasPair=aPair;
00033   mWeightOk=false;
00034 }
00035 
00036 double StHbtThPairGaussFit::GetRejectionProb2Size (double aX, double aY, double aZ, double aT){
00037   return 1.0;
00038   /*
00039   double tCorSizeX,tCorSizeY,tCorSizeZ,tCorSizeT;
00040   if (mSizeX==aX) { tCorSizeX=0.;}
00041   else {tCorSizeX=1./(aX*aX)-1./(mSizeX*mSizeX);}
00042   if (mSizeY==aY) { tCorSizeY=0.;}
00043   else {tCorSizeY=1./(aY*aY)-1./(mSizeY*mSizeY);}
00044   if (mSizeZ==aZ) { tCorSizeZ=0.;}
00045   else {tCorSizeZ=1./(aZ*aZ)-1./(mSizeZ*mSizeZ);}
00046   if (mTime==aT) { tCorSizeT=0.;}
00047   else {tCorSizeT=1./(aT*aT)-1./(mTime*mTime);}
00048   const StHbtLorentzVector  *tEM1;
00049   const StHbtLorentzVector *tEM2;
00050   if(mRef==RCMS){
00051     tEM1=mEmPoint1;
00052     tEM2=mEmPoint2;
00053   } else {
00054     tEM1= &mSourceDist1;
00055     tEM2= &mSourceDist2;
00056   }
00057   return exp(-0.5*((::pow(tEM1->x(),2)+::pow(tEM2->x(),2))*tCorSizeX
00058                    +(::pow(tEM1->y(),2)+::pow(tEM2->y(),2))*tCorSizeY
00059                    +(::pow(tEM1->z(),2)+::pow(tEM2->z(),2))*tCorSizeZ
00060                    +(::pow(tEM1->t(),2)+::pow(tEM2->t(),2))*tCorSizeT));
00061   */
00062 }
00063 
00064 void StHbtThPairGaussFit::setVariables(const StHbtPair*){
00065   
00066 }