JpsiGamEtaPiPiHist modified

Examples/JpsiGamEtaPiPi/JpsiGamEtaPiPi.cfg

@@ -13,7 +13,9 @@ startHypo  = production
 
 enableHyp = etaToPiPiEtaHyp
 enableHyp = etaToa980PiHyp
-enableHyp = etaToa2_1320PiHyp
+enableHyp = etaTof0_980EtaHyp
+enableHyp = etaToSigmaEtaHyp
+#enableHyp = etaToa2_1320PiHyp
 enableHyp = f1ToPiPiEtaHyp
 enableHyp = f1Toa980PiHyp
 enableHyp = usePhasespace
@@ -21,8 +23,8 @@ enableHyp = usePhasespace
 
 mnParFix = J1P-1C-1L1S1PsiToEtaGammaPhi
 
-#massRangeMin = 1.26
-#massRangeMax = 1.30
+massRangeMin = 1.30
+massRangeMax = 1.35
 
 #massRangeMin = 0.
 #massRangeMax = 5.
@@ -33,9 +35,9 @@ mnParFix = J1P-1C-1L1S1PsiToEtaGammaPhi
 #massRangeMin = 1.39
 #massRangeMax = 1.46
 
-massRangeMin = 1.75
-massRangeMax = 1.85
+#massRangeMin = 1.75
+#massRangeMax = 1.85
 
-mode = qaMode
+mode = pwa
 # Determines whether information about configuration options should be emitted
 verbose = 1

Examples/JpsiGamEtaPiPi/JpsiGamEtaPiPiHist.cc

@@ -2,6 +2,7 @@
 #include <fstream>
 #include <sstream>
 #include <string>
+#include <TStyle.h>
 
 #include "Examples/JpsiGamEtaPiPi/JpsiGamEtaPiPiHist.hh"
 #include "Examples/JpsiGamEtaPiPi/JpsiGamEtaPiPiEventList.hh"
@@ -15,6 +16,7 @@
 #include "TH2F.h"
 #include "TMath.h"
 #include "TNtuple.h"
+#include "TCanvas.h"
 #include "TLorentzVector.h"
 #include "ErrLogger/ErrLogger.hh"
 
@@ -214,6 +216,46 @@ JpsiGamEtaPiPiHist::JpsiGamEtaPiPiHist(boost::shared_ptr<AbsLh> theJpsiGamEtaPiP
 
 
 JpsiGamEtaPiPiHist::~JpsiGamEtaPiPiHist(){
+	gStyle->SetPalette(1);
+	std::stringstream ss;
+	ss << _massRange.first << "_" << _massRange.second;
+
+	TCanvas* c_output=new TCanvas( "c_output","c_output",1000,600);
+	c_output->Divide(2,2);
+	c_output->cd(1);
+	_EtaPiMassDataHist->SetLineColor(kRed);
+	_EtaPiMassDataHist->DrawCopy("e");
+	_EtaPiMassFittedHist->DrawCopy("same");
+	c_output->cd(2);
+	_PipPimMassDataHist->SetLineColor(kRed);
+	_PipPimMassDataHist->DrawCopy("e");
+	_PipPimMassFittedHist->DrawCopy("same");
+	c_output->cd(3);
+	_dalitzDataHist->DrawCopy("colz");
+	c_output->cd(4);
+	_dalitzFittedHist->DrawCopy("colz");
+/*	c_output->cd(3);
+	_costEta_EtaPipHeliDataHist->SetLineColor(kRed);
+	_costEta_EtaPipHeliDataHist->DrawCopy("e");
+	_costEta_EtaPipHeliFittedHist->DrawCopy("same");
+	c_output->cd(4);
+	_phiEta_EtaPipHeliDataHist->SetLineColor(kRed);
+	_phiEta_EtaPipHeliDataHist->DrawCopy("e");
+	_phiEta_EtaPipHeliFittedHist->DrawCopy("same");
+	c_output->cd(5);
+	_costPip_PipPimHeliDataHist->SetLineColor(kRed);
+	_costPip_PipPimHeliDataHist->DrawCopy("e");
+	_costPip_PipPimHeliFittedHist->DrawCopy("same");
+	c_output->cd(6);
+	_phiPip_PipPimHeliDataHist->SetLineColor(kRed);
+	_phiPip_PipPimHeliDataHist->DrawCopy("e");
+	_phiPip_PipPimHeliFittedHist->DrawCopy("same");
+*/
+
+	TString outputfilename;
+	outputfilename = "JpsiGamEtaPiPi_" + ss.str() + ".pdf";
+	c_output->Print(outputfilename);
+	c_output->Close();
 	_theTFile->Write();
 	_theTFile->Close();
 }

Examples/JpsiGamEtaPiPi/MJpsiGamEtaPiPiApp.cc

@@ -313,7 +313,7 @@ int main(int __argc,char *__argv[]){
 		//std::ofstream theCompStream ( "componentIntensity.dat");
 		//theProdLh->dumpComponentIntensity( theCompStream, finalFitParams, theErrMatrix );
  		JpsiGamEtaPiPiHist theHist(theLhPtr, finalFitParams);
-// 		theHist.setMassRange(theAppParams.massRange() );
+ 		theHist.setMassRange(theAppParams.massRange() );
 // 		theHist.fill();
 
 		return 0;

Examples/JpsiGamEtaPiPi/PsiToXGamStates.cc

@@ -15,10 +15,10 @@ JpsiToXGamStates::JpsiToXGamStates() :
   _eta2JPC(new jpcRes(2,-1,1) )
 {
   //LS combinations for Psi decays to X gamma
-  fillJPCLS(_psiJPC, _etaJPC, _gammaJPC, _JPCLS_PsiToEtaGamma);
-  fillJPCLS(_psiJPC, _f0JPC, _gammaJPC, _JPCLS_PsiToF0Gamma);
-  fillJPCLS(_psiJPC, _f1JPC, _gammaJPC, _JPCLS_PsiToF1Gamma);
-  fillJPCLS(_psiJPC, _f2JPC, _gammaJPC, _JPCLS_PsiToF2Gamma);
+  fillJPCLS(_psiJPC, _etaJPC , _gammaJPC, _JPCLS_PsiToEtaGamma );
+  fillJPCLS(_psiJPC, _f0JPC  , _gammaJPC, _JPCLS_PsiToF0Gamma  );
+  fillJPCLS(_psiJPC, _f1JPC  , _gammaJPC, _JPCLS_PsiToF1Gamma  );
+  fillJPCLS(_psiJPC, _f2JPC  , _gammaJPC, _JPCLS_PsiToF2Gamma  );
   fillJPCLS(_psiJPC, _eta2JPC, _gammaJPC, _JPCLS_PsiToEta2Gamma);
   
 }