00001 #include "StHbtMaker/Infrastructure/StHbtTypes.hh"
00002 #include "StHbtMaker/Infrastructure/StHbtV0.hh"
00003 #include "StHbtMaker/Cut/v0CutMonitor_Minv.h"
00004 #include <cstdio>
00005
00006 #ifdef __ROOT__
00007 ClassImp(v0CutMonitor_Minv)
00008 #endif
00009
00010 v0CutMonitor_Minv::v0CutMonitor_Minv(){
00011 mHistoK0Short = new StHbt1DHisto("K0ShortMinv","invariant mass (GeV/c)",200,0.,1.);
00012 mHistoK0Short->SetDirectory(0);
00013 mHistoLambda = new StHbt1DHisto("LambdaMinv","invariant mass (GeV/c)",200,0.5,1.5);
00014 mHistoLambda->SetDirectory(0);
00015 mHistoAntiLambda = new StHbt1DHisto("AntiLambdaMinv","invariant mass (GeV/c)",200,0.5,1.5);
00016 mHistoAntiLambda->SetDirectory(0);
00017 }
00018
00019 v0CutMonitor_Minv::v0CutMonitor_Minv(const char* Titlek0,const char* Titlela,const char* Titlelab){
00020 mHistoK0Short = new StHbt1DHisto(Titlek0,"invariant mass (GeV/c)",200,0.,1.);
00021 mHistoK0Short->SetDirectory(0);
00022 mHistoLambda = new StHbt1DHisto(Titlela,"invariant mass (GeV/c)",200,0.5,1.5);
00023 mHistoLambda->SetDirectory(0);
00024 mHistoAntiLambda = new StHbt1DHisto(Titlelab,"invariant mass (GeV/c)",200,0.5,1.5);
00025 mHistoAntiLambda->SetDirectory(0);
00026 }
00027
00028 v0CutMonitor_Minv::v0CutMonitor_Minv(const char* TitCutMoni, const char* title,
00029 int nbins, double min, double max){
00030 mHistoK0Short = new StHbt1DHisto(TitCutMoni, title, nbins , min, max);
00031 mHistoK0Short->SetDirectory(0);
00032 mHistoLambda = new StHbt1DHisto(TitCutMoni, title, nbins , min, max);
00033 mHistoLambda->SetDirectory(0);
00034 mHistoAntiLambda = new StHbt1DHisto(TitCutMoni, title, nbins , min, max);
00035 mHistoAntiLambda->SetDirectory(0);
00036 }
00037
00038 v0CutMonitor_Minv::~v0CutMonitor_Minv(){
00039 delete mHistoK0Short;
00040 delete mHistoLambda;
00041 delete mHistoAntiLambda;
00042 }
00043
00044
00045 void v0CutMonitor_Minv::Fill(const StHbtV0* v0){
00046 mHistoK0Short->Fill( v0->massK0Short() , 1.);
00047 mHistoLambda->Fill( v0->massLambda() , 1.);
00048 mHistoAntiLambda->Fill( v0->massAntiLambda() , 1.);
00049 }
00050
00051
00052 void v0CutMonitor_Minv::Finish(){
00053 cout << " entries in K0Short histogram : " << mHistoK0Short->Integral() << endl;
00054 cout << " entries in Lambda histogram : " << mHistoLambda->Integral() << endl;
00055 cout << " entries in AntiLambda histogram : " << mHistoAntiLambda->Integral() << endl;
00056 }
00057
00058
00059 StHbtString v0CutMonitor_Minv::Report(){
00060 string Stemp;
00061 char Ctemp[100];
00062 sprintf(Ctemp," v0CutMonitor ");
00063 Stemp=Ctemp;
00064 StHbtString returnThis = Stemp;
00065 return returnThis;
00066 }
00067
