00001 #include "StHbtMaker/Infrastructure/StHbtTrack.hh"
00002 #include "StHbtMaker/Infrastructure/StHbtTypes.hh"
00003 #include "StHbtMaker/Cut/trackCutMonitor_Y_vs_Pt.h"
00004 #include <cstdio>
00005
00006 #ifdef __ROOT__
00007 ClassImp(trackCutMonitor_Y_vs_Pt)
00008 #endif
00009
00010 trackCutMonitor_Y_vs_Pt::trackCutMonitor_Y_vs_Pt(){
00011 mHisto = new StHbt2DHisto("Y_vs_Pt","y vs Pt (GeV/c)",50,-2.,2,50,0.,2.);
00012 cout << " this " << this << endl;
00013 mHisto->SetDirectory(0);
00014 }
00015
00016 trackCutMonitor_Y_vs_Pt::trackCutMonitor_Y_vs_Pt(double aMass, const char* TitCutMoni, const char* title,
00017 int nbins1, double min1, double max1,
00018 int nbins2, double min2, double max2):mMass(aMass){
00019 mHisto = new StHbt2DHisto(TitCutMoni, title, nbins1 , min1, max1, nbins2 , min2, max2);
00020 mHisto->SetDirectory(0);
00021 }
00022
00023 trackCutMonitor_Y_vs_Pt::trackCutMonitor_Y_vs_Pt(const char* TitCutMoni, const char* title,
00024 int nbins1, double min1, double max1,
00025 int nbins2, double min2, double max2):mMass(0.139){
00026 mHisto = new StHbt2DHisto(TitCutMoni, title, nbins1 , min1, max1, nbins2 , min2, max2);
00027 mHisto->SetDirectory(0);
00028 }
00029
00030 trackCutMonitor_Y_vs_Pt::~trackCutMonitor_Y_vs_Pt(){
00031 delete mHisto;
00032 }
00033
00034
00035 void trackCutMonitor_Y_vs_Pt::Fill(const StHbtTrack* track){
00036 StHbtLorentzVector L;
00037 L.setVect(track->P());
00038 L.setE(track->P().massHypothesis(mMass));
00039 mHisto->Fill( L.rapidity(), L.vect().perp(), 1.);
00040 }
00041
00042
