Skip to content
Snippets Groups Projects
Commit 00f8a03a authored by Bertram Kopf's avatar Bertram Kopf
Browse files

JpsiGamKsKlKK: added some histograms

parent 73008909
No related branches found
No related tags found
No related merge requests found
...@@ -45,7 +45,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<const JpsiGamKsKlKKEventL ...@@ -45,7 +45,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<const JpsiGamKsKlKKEventL
plotKpKm( _KpKmMassDataHist, (*it), 1. ); plotKpKm( _KpKmMassDataHist, (*it), 1. );
plotCostPhiKs( _costKs_KsKlHeliDataHist, _phiKs_KsKlHeliDataHist,(*it), 1. ); plotCostPhiKs( _costKs_KsKlHeliDataHist, _phiKs_KsKlHeliDataHist,(*it), 1. );
plotCostPhiKs( _costKp_KpKmHeliDataHist, _phiKp_KpKmHeliDataHist,(*it), 1. ); plotCostPhiKs( _costKp_KpKmHeliDataHist, _phiKp_KpKmHeliDataHist,(*it), 1. );
plotCostGam( _costGamCmDataHist,(*it), 1. );
fillTuple(_dataTuple, (*it), 1.); fillTuple(_dataTuple, (*it), 1.);
++it; ++it;
...@@ -62,6 +62,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<const JpsiGamKsKlKKEventL ...@@ -62,6 +62,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<const JpsiGamKsKlKKEventL
plotKpKm( _KpKmMassMcHist, (*it), 1. ); plotKpKm( _KpKmMassMcHist, (*it), 1. );
plotCostPhiKs( _costKs_KsKlHeliMcHist, _phiKs_KsKlHeliMcHist,(*it), 1. ); plotCostPhiKs( _costKs_KsKlHeliMcHist, _phiKs_KsKlHeliMcHist,(*it), 1. );
plotCostPhiKp( _costKp_KpKmHeliMcHist, _phiKp_KpKmHeliMcHist, (*it), 1. ); plotCostPhiKp( _costKp_KpKmHeliMcHist, _phiKp_KpKmHeliMcHist, (*it), 1. );
plotCostGam( _costGamCmMcHist, (*it), 1. );
fillTuple(_mcTuple, (*it), 1.); fillTuple(_mcTuple, (*it), 1.);
...@@ -105,7 +106,9 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp ...@@ -105,7 +106,9 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp
plotKsKl( _KsKlMassDataHist, (*it), 1. ); plotKsKl( _KsKlMassDataHist, (*it), 1. );
plotKpKm( _KpKmMassDataHist, (*it), 1. ); plotKpKm( _KpKmMassDataHist, (*it), 1. );
plotCostPhiKs( _costKs_KsKlHeliDataHist, _phiKs_KsKlHeliDataHist, (*it), 1. ); plotCostPhiKs( _costKs_KsKlHeliDataHist, _phiKs_KsKlHeliDataHist, (*it), 1. );
plotCostPhiKp( _costKp_KpKmHeliDataHist, _phiKp_KpKmHeliDataHist, (*it), 1. ); plotCostPhiKp( _costKp_KpKmHeliDataHist, _phiKp_KpKmHeliDataHist, (*it), 1. );
plotCostGam( _costGamCmDataHist, (*it), 1. );
fillTuple(_dataTuple, (*it), 1.); fillTuple(_dataTuple, (*it), 1.);
++it; ++it;
...@@ -121,6 +124,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp ...@@ -121,6 +124,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp
plotKpKm( _KpKmMassMcHist, (*it), 1. ); plotKpKm( _KpKmMassMcHist, (*it), 1. );
plotCostPhiKs( _costKs_KsKlHeliMcHist, _phiKs_KsKlHeliMcHist, (*it), 1. ); plotCostPhiKs( _costKs_KsKlHeliMcHist, _phiKs_KsKlHeliMcHist, (*it), 1. );
plotCostPhiKp( _costKp_KpKmHeliMcHist, _phiKp_KpKmHeliMcHist,(*it), 1. ); plotCostPhiKp( _costKp_KpKmHeliMcHist, _phiKp_KpKmHeliMcHist,(*it), 1. );
plotCostGam( _costGamCmMcHist,(*it), 1. );
double evtWeight= theJpsiGamKsKlKKLh->calcEvtIntensity((*it), fitParam); double evtWeight= theJpsiGamKsKlKKLh->calcEvtIntensity((*it), fitParam);
plotDalitz(_dalitzFittedHist, (*it),evtWeight ); plotDalitz(_dalitzFittedHist, (*it),evtWeight );
...@@ -129,6 +133,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp ...@@ -129,6 +133,7 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp
plotKpKm( _KpKmMassFittedHist, (*it), evtWeight ); plotKpKm( _KpKmMassFittedHist, (*it), evtWeight );
plotCostPhiKs( _costKs_KsKlHeliFittedHist, _phiKs_KsKlHeliFittedHist,(*it), evtWeight ); plotCostPhiKs( _costKs_KsKlHeliFittedHist, _phiKs_KsKlHeliFittedHist,(*it), evtWeight );
plotCostPhiKp( _costKp_KpKmHeliFittedHist, _phiKp_KpKmHeliFittedHist,(*it), evtWeight ); plotCostPhiKp( _costKp_KpKmHeliFittedHist, _phiKp_KpKmHeliFittedHist,(*it), evtWeight );
plotCostGam( _costGamCmFittedHist,(*it), evtWeight );
fillTuple(_mcTuple, (*it), evtWeight); fillTuple(_mcTuple, (*it), evtWeight);
++it; ++it;
...@@ -152,7 +157,8 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp ...@@ -152,7 +157,8 @@ JpsiGamKsKlKKHist::JpsiGamKsKlKKHist(boost::shared_ptr<AbsJpsiGamKsKlKKLh> theJp
_costKs_KsKlHeliFittedHist->Scale(scaleFactor); _costKs_KsKlHeliFittedHist->Scale(scaleFactor);
_phiKs_KsKlHeliFittedHist->Scale(scaleFactor); _phiKs_KsKlHeliFittedHist->Scale(scaleFactor);
_costKp_KpKmHeliFittedHist->Scale(scaleFactor); _costKp_KpKmHeliFittedHist->Scale(scaleFactor);
_phiKp_KpKmHeliFittedHist->Scale(scaleFactor); _phiKp_KpKmHeliFittedHist->Scale(scaleFactor);
_costGamCmFittedHist->Scale(scaleFactor);
} }
JpsiGamKsKlKKHist::~JpsiGamKsKlKKHist() JpsiGamKsKlKKHist::~JpsiGamKsKlKKHist()
...@@ -208,6 +214,11 @@ void JpsiGamKsKlKKHist::initRootStuff() ...@@ -208,6 +214,11 @@ void JpsiGamKsKlKKHist::initRootStuff()
_phiKp_KpKmHeliDataHist= new TH1F("_phiKp_KpKmHeliDataHist", "#Phi_{K+} KsKlHeli data", 100, -1., 1.); _phiKp_KpKmHeliDataHist= new TH1F("_phiKp_KpKmHeliDataHist", "#Phi_{K+} KsKlHeli data", 100, -1., 1.);
_phiKp_KpKmHeliMcHist= new TH1F("_phiKp_KpKmHeliMcHist", "#Phi_{K+} KsKlHeli Mc", 100, -1., 1.); _phiKp_KpKmHeliMcHist= new TH1F("_phiKp_KpKmHeliMcHist", "#Phi_{K+} KsKlHeli Mc", 100, -1., 1.);
_phiKp_KpKmHeliFittedHist= new TH1F("_phiKp_KpKmHeliFittedHist", "#Phi_{K+} KsKlHeli Fit", 100, -1, 1); _phiKp_KpKmHeliFittedHist= new TH1F("_phiKp_KpKmHeliFittedHist", "#Phi_{K+} KsKlHeli Fit", 100, -1, 1);
_costGamCmDataHist= new TH1F("_costGamCmDataHist", "cos(#Theta_{#gamma}) CM data", 100, -1., 1.);
_costGamCmMcHist= new TH1F("_costGamCmMcHist", "cos(#Theta_{#gamma}) CM Mc", 100, -1., 1.);
_costGamCmFittedHist= new TH1F("_costGamCmFittedHist", "cos(#Theta_{#gamma}) CM Fit", 100, -1, 1);
std::string tupleVariables = "mKsKlKpKm:mKsKl:mKpKm:KsKlKpKmCostTheta:gamCosTheta:KsKlCosTheta:KpKmCosTheta:KsCosTheta:KpCosTheta:weight"; std::string tupleVariables = "mKsKlKpKm:mKsKl:mKpKm:KsKlKpKmCostTheta:gamCosTheta:KsKlCosTheta:KpKmCosTheta:KsCosTheta:KpCosTheta:weight";
...@@ -254,6 +265,12 @@ void JpsiGamKsKlKKHist::plotCostPhiKp(TH1F* theCostHisto, TH1F* thePhiHisto, co ...@@ -254,6 +265,12 @@ void JpsiGamKsKlKKHist::plotCostPhiKp(TH1F* theCostHisto, TH1F* thePhiHisto, co
thePhiHisto->Fill( v4.Phi(), weight ); thePhiHisto->Fill( v4.Phi(), weight );
} }
void JpsiGamKsKlKKHist::plotCostGam(TH1F* theCostHisto, const JpsiGamKsKlKKEvtData* theData, double weight){
Vector4<float> v4 = theData->V4_gamma_HeliPsi;
theCostHisto->Fill( v4.CosTheta(), weight );
}
void JpsiGamKsKlKKHist::fillTuple( TNtuple* theTuple, const JpsiGamKsKlKKEvtData* theData, double weight) void JpsiGamKsKlKKHist::fillTuple( TNtuple* theTuple, const JpsiGamKsKlKKEvtData* theData, double weight)
{ {
......
...@@ -73,6 +73,11 @@ private: ...@@ -73,6 +73,11 @@ private:
TH1F* _phiKp_KpKmHeliDataHist; TH1F* _phiKp_KpKmHeliDataHist;
TH1F* _phiKp_KpKmHeliMcHist; TH1F* _phiKp_KpKmHeliMcHist;
TH1F* _phiKp_KpKmHeliFittedHist; TH1F* _phiKp_KpKmHeliFittedHist;
TH1F* _costGamCmDataHist;
TH1F* _costGamCmMcHist;
TH1F* _costGamCmFittedHist;
TNtuple* _dataTuple; TNtuple* _dataTuple;
TNtuple* _mcTuple; TNtuple* _mcTuple;
...@@ -86,6 +91,8 @@ private: ...@@ -86,6 +91,8 @@ private:
void plotKpKm(TH1F* theHisto, const JpsiGamKsKlKKEvtData* theData, double weight); void plotKpKm(TH1F* theHisto, const JpsiGamKsKlKKEvtData* theData, double weight);
void plotCostPhiKs(TH1F* theCostHisto, TH1F* thePhiHisto, const JpsiGamKsKlKKEvtData* theData, double weight); void plotCostPhiKs(TH1F* theCostHisto, TH1F* thePhiHisto, const JpsiGamKsKlKKEvtData* theData, double weight);
void plotCostPhiKp(TH1F* theCostHisto, TH1F* thePhiHisto, const JpsiGamKsKlKKEvtData* theData, double weight); void plotCostPhiKp(TH1F* theCostHisto, TH1F* thePhiHisto, const JpsiGamKsKlKKEvtData* theData, double weight);
void plotCostGam(TH1F* theCostHisto, const JpsiGamKsKlKKEvtData* theData, double weight);
void fillTuple( TNtuple* theTuple, const JpsiGamKsKlKKEvtData* theData, double weight); void fillTuple( TNtuple* theTuple, const JpsiGamKsKlKKEvtData* theData, double weight);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment