Skip to content
Snippets Groups Projects
Commit 821638d3 authored by Michael Leyhe's avatar Michael Leyhe
Browse files

Added new histograms to JpsiGamEtaPiPiHist

parent 0b520e78
No related branches found
No related tags found
No related merge requests found
......@@ -47,6 +47,18 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<const EvtDataBaseList>
_phiPip_PipPimHeliDataHist(0),
_phiPip_PipPimHeliMcHist(0),
_phiPip_PipPimHeliFittedHist(0),
_costPiPi_EtaPiPiHeliDataHist(0),
_costPiPi_EtaPiPiHeliMcHist(0),
_costPiPi_EtaPiPiHeliFittedHist(0),
_phiPiPi_EtaPiPiHeliDataHist(0),
_phiPiPi_EtaPiPiHeliMcHist(0),
_phiPiPi_EtaPiPiHeliFittedHist(0),
_costEtaPi_EtaPiPiHeliDataHist(0),
_costEtaPi_EtaPiPiHeliMcHist(0),
_costEtaPi_EtaPiPiHeliFittedHist(0),
_phiEtaPi_EtaPiPiHeliDataHist(0),
_phiEtaPi_EtaPiPiHeliMcHist(0),
_phiEtaPi_EtaPiPiHeliFittedHist(0),
_costGamCmDataHist(0),
_costGamCmMcHist(0),
_costGamCmFittedHist(0),
......@@ -76,7 +88,8 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<const EvtDataBaseList>
plotCostPhiEta(_costEta_EtaPipHeliDataHist , _phiEta_EtaPipHeliDataHist ,(*it), 1. );
plotCostPhiPip(_costPip_PipPimHeliDataHist, _phiPip_PipPimHeliDataHist,(*it), 1. );
plotCostGam(_costGamCmDataHist,(*it), 1. );
// plotCostPhi_PhiPhiHeli(_costPhi_KpKmDataHist, _phiPhi_KpKmDataHist, (*it)->FourVecs[enumJpsiGamEtaPiPiData::V4_KpKm_HeliKsKlKpKm], 1.);
plotCostPhiPip(_costPiPi_EtaPiPiHeliDataHist, _phiPiPi_EtaPiPiHeliDataHist,(*it), 1. );
plotCostPhiEta(_costEtaPi_EtaPiPiHeliDataHist , _phiEtaPi_EtaPiPiHeliDataHist ,(*it), 1. );
// fillTuple(_dataTuple, (*it), 1.);
++it;
}
......@@ -93,7 +106,8 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<const EvtDataBaseList>
plotCostPhiEta(_costEta_EtaPipHeliMcHist , _phiEta_EtaPipHeliMcHist ,(*it), 1. );
plotCostPhiPip(_costPip_PipPimHeliMcHist, _phiPip_PipPimHeliMcHist,(*it), 1. );
plotCostGam(_costGamCmMcHist,(*it), 1. );
// plotCostPhi_PhiPhiHeli(_costPhi_KpKmMcHist, _phiPhi_KpKmMcHist, (*it)->FourVecs[enumJpsiGamEtaPiPiMc::V4_KpKm_HeliKsKlKpKm], 1.);
plotCostPhiPip(_costPiPi_EtaPiPiHeliMcHist, _phiPiPi_EtaPiPiHeliMcHist,(*it), 1. );
plotCostPhiEta(_costEtaPi_EtaPiPiHeliMcHist , _phiEtaPi_EtaPiPiHeliMcHist ,(*it), 1. );
// fillTuple(_mcTuple, (*it), 1.);
++it;
}
......@@ -128,6 +142,18 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<AbsLh> theJpsiGamEtaPiP
_phiPip_PipPimHeliDataHist(0),
_phiPip_PipPimHeliMcHist(0),
_phiPip_PipPimHeliFittedHist(0),
_costPiPi_EtaPiPiHeliDataHist(0),
_costPiPi_EtaPiPiHeliMcHist(0),
_costPiPi_EtaPiPiHeliFittedHist(0),
_phiPiPi_EtaPiPiHeliDataHist(0),
_phiPiPi_EtaPiPiHeliMcHist(0),
_phiPiPi_EtaPiPiHeliFittedHist(0),
_costEtaPi_EtaPiPiHeliDataHist(0),
_costEtaPi_EtaPiPiHeliMcHist(0),
_costEtaPi_EtaPiPiHeliFittedHist(0),
_phiEtaPi_EtaPiPiHeliDataHist(0),
_phiEtaPi_EtaPiPiHeliMcHist(0),
_phiEtaPi_EtaPiPiHeliFittedHist(0),
_costGamCmDataHist(0),
_costGamCmMcHist(0),
_costGamCmFittedHist(0),
......@@ -160,7 +186,8 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<AbsLh> theJpsiGamEtaPiP
plotCostPhiEta(_costEta_EtaPipHeliDataHist , _phiEta_EtaPipHeliDataHist ,(*it), 1. );
plotCostPhiPip(_costPip_PipPimHeliDataHist, _phiPip_PipPimHeliDataHist,(*it), 1. );
plotCostGam(_costGamCmDataHist,(*it), 1. );
// plotCostPhi_PhiPhiHeli(_costPhi_KpKmDataHist, _phiPhi_KpKmDataHist, (*it)->FourVecs[enumJpsiGamEtaPiPiData::V4_KpKm_HeliKsKlKpKm], 1.);
plotCostPhiPip(_costPiPi_EtaPiPiHeliDataHist, _phiPiPi_EtaPiPiHeliDataHist,(*it), 1. );
plotCostPhiEta(_costEtaPi_EtaPiPiHeliDataHist , _phiEtaPi_EtaPiPiHeliDataHist ,(*it), 1. );
// fillTuple(_dataTuple, (*it), 1.);
++it;
}
......@@ -177,7 +204,8 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<AbsLh> theJpsiGamEtaPiP
plotCostPhiEta(_costEta_EtaPipHeliMcHist , _phiEta_EtaPipHeliMcHist ,(*it), 1. );
plotCostPhiPip(_costPip_PipPimHeliMcHist, _phiPip_PipPimHeliMcHist,(*it), 1. );
plotCostGam(_costGamCmMcHist,(*it), 1. );
plotCostPhiPip(_costPiPi_EtaPiPiHeliMcHist, _phiPiPi_EtaPiPiHeliMcHist,(*it), 1. );
plotCostPhiEta(_costEtaPi_EtaPiPiHeliMcHist , _phiEtaPi_EtaPiPiHeliMcHist ,(*it), 1. );
// fillTuple(_mcTuple, (*it), 1.);
......@@ -190,6 +218,8 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<AbsLh> theJpsiGamEtaPiP
plotCostPhiEta(_costEta_EtaPipHeliFittedHist , _phiEta_EtaPipHeliFittedHist ,(*it), evtWeight );
plotCostPhiPip(_costPip_PipPimHeliFittedHist, _phiPip_PipPimHeliFittedHist,(*it), evtWeight );
plotCostGam(_costGamCmFittedHist,(*it), evtWeight );
plotCostPhiPip(_costPiPi_EtaPiPiHeliFittedHist, _phiPiPi_EtaPiPiHeliFittedHist,(*it), evtWeight );
plotCostPhiEta(_costEtaPi_EtaPiPiHeliFittedHist , _phiEtaPi_EtaPiPiHeliFittedHist ,(*it), evtWeight );
++it;
}
......@@ -211,7 +241,10 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<AbsLh> theJpsiGamEtaPiP
_costPip_PipPimHeliFittedHist->Scale(scalingFactor);
_phiPip_PipPimHeliFittedHist->Scale(scalingFactor);
_costGamCmFittedHist->Scale(scalingFactor);
_costPiPi_EtaPiPiHeliFittedHist->Scale(scalingFactor);
_phiPiPi_EtaPiPiHeliFittedHist->Scale(scalingFactor);
_costEtaPi_EtaPiPiHeliFittedHist->Scale(scalingFactor);
_phiEtaPi_EtaPiPiHeliFittedHist->Scale(scalingFactor);
}
......@@ -312,6 +345,20 @@ void JpsiGamEtaPiPiHist::initRootStuff(){
_phiPip_PipPimHeliMcHist = new TH1F("_phiPip_PipPimHeliMcHist" , "#Phi_{#pi^{+}} #pi^{+}#pi^{-}Heli Mc" , 100, -TMath::Pi(), TMath::Pi());
_phiPip_PipPimHeliFittedHist= new TH1F("_phiPip_PipPimHeliFittedHist", "#Phi_{#pi^{+}} #pi^{+}#pi^{-}Heli Fit" , 100, -TMath::Pi(), TMath::Pi());
_costPiPi_EtaPiPiHeliDataHist = new TH1F("_costPiPi_EtaPiPiHeliDataHist" , "cos(#Theta_{#pi^{+}#pi^{-}}) #eta#pi^{+}#pi^{-}Heli data", 100, -1., 1.);
_costPiPi_EtaPiPiHeliMcHist = new TH1F("_costPiPi_EtaPiPiHeliMcHist" , "cos(#Theta_{#pi^{+}#pi^{-}}) #eta#pi^{+}#pi^{-}Heli Mc", 100, -1., 1.);
_costPiPi_EtaPiPiHeliFittedHist = new TH1F("_costPiPi_EtaPiPiHeliFittedHist" , "cos(#Theta_{#pi^{+}#pi^{-}}) #eta#pi^{+}#pi^{-}Heli Fit", 100, -1., 1.);
_phiPiPi_EtaPiPiHeliDataHist = new TH1F("_phiPiPi_EtaPiPiHeliDataHist" , "#Phi_{#pi^{+}#pi^{-}} #eta#pi^{+}#pi^{-}Heli data", 100, -TMath::Pi(), TMath::Pi());
_phiPiPi_EtaPiPiHeliMcHist = new TH1F("_phiPiPi_EtaPiPiHeliMcHist" , "#Phi_{#pi^{+}#pi^{-}} #eta#pi^{+}#pi^{-}Heli data", 100, -TMath::Pi(), TMath::Pi());
_phiPiPi_EtaPiPiHeliFittedHist = new TH1F("_phiPiPi_EtaPiPiHeliFittedHist" , "#Phi_{#pi^{+}#pi^{-}} #eta#pi^{+}#pi^{-}Heli data", 100, -TMath::Pi(), TMath::Pi());
_costEtaPi_EtaPiPiHeliDataHist = new TH1F("_costEtaPi_EtaPiPiHeliDataHist" , "cos(#Theta_{#eta#pi^{+-}}) #eta#pi^{+}#pi^{-}Heli data", 100, -1., 1.);
_costEtaPi_EtaPiPiHeliMcHist = new TH1F("_costEtaPi_EtaPiPiHeliDataHist" , "cos(#Theta_{#eta#pi^{+-}}) #eta#pi^{+}#pi^{-}Heli data", 100, -1., 1.);
_costEtaPi_EtaPiPiHeliFittedHist = new TH1F("_costEtaPi_EtaPiPiHeliDataHist" , "cos(#Theta_{#eta#pi^{+-}}) #eta#pi^{+}#pi^{-}Heli data", 100, -1., 1.);
_phiEtaPi_EtaPiPiHeliDataHist = new TH1F("_phiEtaPi_EtaPiPiHeliDataHist" , "#Phi_{#eta#pi^{+-}} #eta#pi^{+}#pi^{-}Heli data", 100, -TMath::Pi(), TMath::Pi());
_phiEtaPi_EtaPiPiHeliMcHist = new TH1F("_phiEtaPi_EtaPiPiHeliDataHist" , "#Phi_{#eta#pi^{+-}} #eta#pi^{+}#pi^{-}Heli data", 100, -TMath::Pi(), TMath::Pi());
_phiEtaPi_EtaPiPiHeliFittedHist = new TH1F("_phiEtaPi_EtaPiPiHeliDataHist" , "#Phi_{#eta#pi^{+-}} #eta#pi^{+}#pi^{-}Heli data", 100, -TMath::Pi(), TMath::Pi());
_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);
......@@ -364,12 +411,21 @@ void JpsiGamEtaPiPiHist::plotCostGam(TH1F* theCostHisto, EvtData* theData, doubl
Vector4<double> v4 = theData->FourVecs[enumJpsiGamEtaPiPiData::V4_gamma_HeliPsi];
theCostHisto->Fill( v4.CosTheta(), weight );
}
/*
void JpsiGamEtaPiPiHist::plotCostPhi_PhiPhiHeli(TH1F* theCostHisto, TH1F* thePhiHisto, const Vector4<double>& the4Vec, double weight){
theCostHisto->Fill( the4Vec.CosTheta(), weight);
thePhiHisto->Fill( the4Vec.Phi(), weight);
void JpsiGamEtaPiPiHist::plotCostPhi_PiPiHeli(TH1F* theCostHisto, TH1F* thePhiHisto, EvtData* theData, double weight){
Vector4<double> v4 = theData->FourVecs[enumJpsiGamEtaPiPiData::V4_PipPim_HeliEtaPipPim];
theCostHisto->Fill( v4.CosTheta(), weight);
thePhiHisto->Fill( v4.Phi(), weight);
}
void JpsiGamEtaPiPiHist::plotCostPhi_EtaPiHeli(TH1F* theCostHisto, TH1F* thePhiHisto, EvtData* theData, double weight){
Vector4<double> v4 = theData->FourVecs[enumJpsiGamEtaPiPiData::V4_EtaPip_HeliEtaPipPim];
theCostHisto->Fill( v4.CosTheta(), weight);
thePhiHisto->Fill( v4.Phi(), weight);
v4 = theData->FourVecs[enumJpsiGamEtaPiPiData::V4_EtaPim_HeliEtaPipPim];
theCostHisto->Fill( v4.CosTheta(), weight);
thePhiHisto->Fill( v4.Phi(), weight);
}
*/
......@@ -380,7 +436,7 @@ void JpsiGamEtaPiPiHist::fillTuple( TNtuple* theTuple, EvtData* theData, double
Vector4<double> V4_gamma_HeliPsi= theData->FourVecs[enumJpsiGamEtaPiPiData::V4_gamma_HeliPsi] ;
Vector4<double> V4_EtaPip_HeliEtaPipPim= theData->FourVecs[enumJpsiGamEtaPiPiData::V4_EtaPip_HeliEtaPipPim] ;
Vector4<double> V4_PipPim_HeliKsKlKpKm= theData->FourVecs[enumJpsiGamEtaPiPiData::V4_PipPim_HeliEtaPipPim] ;
Vector4<double> V4_PipPim_HeliEtaPipPim= theData->FourVecs[enumJpsiGamEtaPiPiData::V4_PipPim_HeliEtaPipPim] ;
Vector4<double> V4_Eta_HeliEtaPip= theData->FourVecs[enumJpsiGamEtaPiPiData::V4_Eta_HeliEtaPip] ;
Vector4<double> V4_Pip_HeliPipPim= theData->FourVecs[enumJpsiGamEtaPiPiData::V4_Pip_HeliPipPim] ;
......@@ -411,7 +467,7 @@ void JpsiGamEtaPiPiHist::fillTuple( TNtuple* theTuple, EvtData* theData, double
V4_EtaPipPim_HeliPsi.CosTheta(),
V4_gamma_HeliPsi.CosTheta(),
V4_EtaPip_HeliEtaPipPim.CosTheta(),
V4_PipPim_HeliKsKlKpKm.CosTheta(),
V4_PipPim_HeliEtaPipPim.CosTheta(),
V4_Eta_HeliEtaPip.CosTheta(),
V4_Pip_HeliPipPim.CosTheta(),
thePhiPhiDecayPlaneAngle,
......
......@@ -84,6 +84,20 @@ private:
TH1F* _phiPip_PipPimHeliMcHist;
TH1F* _phiPip_PipPimHeliFittedHist;
TH1F* _costPiPi_EtaPiPiHeliDataHist;
TH1F* _costPiPi_EtaPiPiHeliMcHist;
TH1F* _costPiPi_EtaPiPiHeliFittedHist;
TH1F* _phiPiPi_EtaPiPiHeliDataHist;
TH1F* _phiPiPi_EtaPiPiHeliMcHist;
TH1F* _phiPiPi_EtaPiPiHeliFittedHist;
TH1F* _costEtaPi_EtaPiPiHeliDataHist;
TH1F* _costEtaPi_EtaPiPiHeliMcHist;
TH1F* _costEtaPi_EtaPiPiHeliFittedHist;
TH1F* _phiEtaPi_EtaPiPiHeliDataHist;
TH1F* _phiEtaPi_EtaPiPiHeliMcHist;
TH1F* _phiEtaPi_EtaPiPiHeliFittedHist;
TH1F* _costGamCmDataHist;
TH1F* _costGamCmMcHist;
TH1F* _costGamCmFittedHist;
......@@ -104,7 +118,8 @@ private:
void plotCostPhiPip(TH1F* theCostHisto, TH1F* thePhiHisto, EvtData* theData, double weight);
// void plotCostPhi_PhiPhiHeli(TH1F* theCostHisto, TH1F* thePhiHisto, const Vector4<double>& the4Vec, double weight);
void plotCostGam(TH1F* theCostHisto, EvtData* theData, double weight);
void plotCostPhi_PiPiHeli(TH1F* theCostHisto, TH1F* thePhiHisto, EvtData* theData, double weight);
void plotCostPhi_EtaPiHeli(TH1F* theCostHisto, TH1F* thePhiHisto, EvtData* theData, double weight);
void fillTuple( TNtuple* theTuple, EvtData* 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