StRoot: StHbtMaker/ThCorrFctn/ThQinvCorrFctn.cxx Source File

00001 /***************************************************************************
00002  *
00003  * $Id: 
00004  *
00005  * Author: Laurent Conin, Fabrice Retiere, Subatech, France
00006  ***************************************************************************
00007  *
00008  * Description : Calculate the theoretical QInv correlation function 
00009  *
00010  ***************************************************************************
00011  *
00012  * $Log: 
00013  *
00014  ***************************************************************************/
00015 
00016 #include "StHbtMaker/ThCorrFctn/ThQinvCorrFctn.h"
00017 
00018 #ifdef __ROOT__
00019 ClassImp(ThQinvCorrFctn)
00020 #endif
00021 
00022 ThQinvCorrFctn::ThQinvCorrFctn(char* aTitle, int aNBins, 
00023                                double aHLo, double aHHi)
00024   : StHbtThCorrFctn() ,StHbtRoot1DCF(aTitle, aNBins, aHLo, aHHi) 
00025 { /* no-op */ };
00026 ThQinvCorrFctn::ThQinvCorrFctn(const ThQinvCorrFctn& ThCf) : StHbtThCorrFctn(ThCf),StHbtRoot1DCF( ThCf) 
00027 { /* no-op */ }
00028 
00029 ThQinvCorrFctn::~ThQinvCorrFctn() 
00030 { /* no-op */ };
00031 
00032 //____________________________
00033 void ThQinvCorrFctn::AddNum(StHbtThPair* aThPair){
00034 
00035 }
00036 void ThQinvCorrFctn::AddDen(StHbtThPair* aThPair){
00037 
00038   double tQInv = 2*fabs(aThPair->GetMeasPair()->KStar());
00039   //  cout << "Filling ThQinv CF with values: KStar:  " << tQInv << "  wDen:  " << aThPair->GetWeightDen() << " wNum:   " << aThPair->GetWeightNum() << endl;
00040   if(tQInv>=mHLo && tQInv<=mHHi){
00041     mDenominator->Fill(tQInv,aThPair->GetWeightDen());
00042   }
00043   if(tQInv>=mHLo && tQInv<=mHHi){
00044     mNumerator->Fill(tQInv,aThPair->GetWeightNum());
00045   }
00046 }
00047 
00048 void ThQinvCorrFctn::Finish(){
00049   StHbtRoot1DCF::Finish();
00050 }
00051 
00052 inline StHbt1DHisto* ThQinvCorrFctn::Numerator() const { cout << "Return Ratio" << endl; return mNumerator;};
00053 inline StHbt1DHisto* ThQinvCorrFctn::Denominator() const {cout << "Return Denimnatior" << endl; return mDenominator;};
00054 inline StHbt1DHisto* ThQinvCorrFctn::Ratio() const {cout << "Return Numerator" << endl; return mRatio;};
00055 inline void ThQinvCorrFctn::Write()  {mNumerator->Write();mDenominator->Write();mRatio->Write();};
00056 
00057 inline StHbtThCorrFctn* ThQinvCorrFctn::ThClone() const {return new ThQinvCorrFctn(*this);}