T49HbtMomentumResolutionCalculator

class T49HbtMomentumResolutionCalculator : public TObject

    public:

                                          T49HbtMomentumResolutionCalculator(TH3D* MomentumDifferenceOfMCTracks, Int_t nKtBins, Double_t minKt, Double_t maxKt, Int_t nYBins, Double_t minY, Double_t maxY)
                                          T49HbtMomentumResolutionCalculator(const T49HbtMomentumResolutionCalculator&)
                                          ~T49HbtMomentumResolutionCalculator()
                           static TClass* Class()
                                     void FillDisplacedMomenta()
                                     void FillDisplacedMomenta(T49ParticleMCRoot* MC1, T49ParticleMCRoot* MC2)
                                    Int_t GetKtYIndex(Double_t lKtPair, Double_t lYPair)
                               TObjArray* GetResolution3dHisots()
                               TObjArray* GetResolutionDistributions(Int_t KtYIndex)
                          virtual TClass* IsA() const
      T49HbtMomentumResolutionCalculator& operator=(const T49HbtMomentumResolutionCalculator&)
                                     void PrintResolution3dHisots(Int_t KtYIndex)
                                     void SaveResolution(Char_t* Name)
                                     void SelectedAssocicatedTracks(TObjArray* MCTrackList, TObjArray* RETrackList)
                             virtual void ShowMembers(TMemberInspector& insp, char* parent)
                             virtual void Streamer(TBuffer& b)
                                     void StreamerNVirtual(TBuffer& b)
                                     void Use(TObjArray* MCTrackList, TObjArray* RETrackList)

Data Members

    private:

           TH3D* mMomentumDifferenceOfMCTracks  Master 3d histo
           Int_t mNumberOfCells                 Number of cells in master 3d histo
      TObjArray* mSelectedMCTracks              To hold MC tracks
      TObjArray* mSelectedRETracks              To hold RE tracks
      TObjArray* mListOfListOfDisplacedMomenta  List of list of displaced momenta
      TObjArray* mResolution                    Array holding the final result
           Int_t mNumberOfResolutionHistos      Number of final result histos
        Double_t mPionMass                      Constant
       Double_t* mKtBins                        
       Double_t* mYBins                         
           Int_t mNKtBins                       
           Int_t mNYBins                        
           Int_t mMaxKtY                        

Class Description

 The T49HbtMomentumResolutionCalculator class is intended to produce for each qSide,qLong,qOut bin in each kt-Y bin
 the momentum resolution in qSide,qOut,qLong and a possible shift in qSide,qOut,qLong. To obtain these values 3 vectors
 for each cell are filled holding the difference between a Monte Carlo pair (MC) and the reconstructed pair (RC)
 in qSide, qOut and qLong. A gaus fit to each vector gives the RMS (=resolution) and the shift, which will is needed
 by the smear procedure in T49HbtMomentumResolutionCorrection.

 The constructor wants to know the dimension of the qside,qout,qlong 3d object and the kt-Y binning.
 The mixing macro calls the Use function, enertering lists of MC and the according RC tracks.
 The result can be stored in a file which will be accessed by T49HbtMomentumResolutionCorrection
 via T49HbtMomentumResolutionProvider

 Structure of T49HbtMomentumResolutionCalculator:

             each kt-Y cell -> ListofDisplacedMomenta[#qside * #qout * #qlong] holding TVector3s of MomDiffs
                            -> TH3D* ResolutionMeanQside (out,long)  results
                            -> TH3D* ResolutionRMSQside  (out,long)  results

 Init master 3d histo

 Empty TObjArrays holding MC and Reconstructed tracks

 Select reconstructed tracks which are best associated with a given MC track -> fill lists with selected tracks

 Construct all possible pairs of MC tracks

 Calculate MC and RE momentum difference, add difference MC-RE to
 vector with index according to the MCqout,MCside,MClong position in master TH3D

 Kt

 Calculate and return resulting resolution histograms

 This member function plots for a given kty bin the mean and rms of each cell
 along with the cells position in the master TH3D.
 If commented out, also the vector is plotted from which rms and mean are derived

 This member function returns 6 TH1D's filled with the mean and the rms values for qSide, qOut & qLong of all cells
 of a given kty bin

 Store the calculated resolution TH3Ds in a root file,
 to have easy access using class T49HbtMomentumResolutionProvider

Inline Functions

                                    TClass* Class()
                                    TClass* IsA() const
                                       void ShowMembers(TMemberInspector& insp, char* parent)
                                       void Streamer(TBuffer& b)
                                       void StreamerNVirtual(TBuffer& b)
         T49HbtMomentumResolutionCalculator T49HbtMomentumResolutionCalculator(const T49HbtMomentumResolutionCalculator&)
        T49HbtMomentumResolutionCalculator& operator=(const T49HbtMomentumResolutionCalculator&)