// $Id$ // // File: JEventProcessor_multiphoton.h // Created: Tue Oct 21 16:50:20 EDT 2014 // Creator: staylor (on Linux gluon05.jlab.org 2.6.32-358.18.1.el6.x86_64 x86_64) // #ifndef _JEventProcessor_multiphoton_ #define _JEventProcessor_multiphoton_ //#include //#include #include #include "TRACKING/DTrackFinder.h" #include #include #include #include #include #include #include #include #include #include #include #include class JEventProcessor_multiphoton:public jana::JEventProcessor{ public: JEventProcessor_multiphoton(); ~JEventProcessor_multiphoton(); const char* className(void){return "JEventProcessor_multiphoton";} private: enum particle_pair_t{ Pi0Pi0_, Pi0Eta_, EtaEta_, Pi0EtaPrime_, Pi0X_, EtaX_, }; jerror_t init(void); ///< Called once at program start. jerror_t brun(jana::JEventLoop *eventLoop, int32_t runnumber); ///< Called everytime a new run number is detected. jerror_t evnt(jana::JEventLoop *eventLoop, uint64_t eventnumber); ///< Called every event. jerror_t erun(void); ///< Called everytime run number changes, provided brun has been called. jerror_t fini(void); ///< Called after last event of last event source has been processed. void FillChargedParticleVectors(vector&tracks,map>&chargedParticles) const; void FillGammaParticleVector(double t0_rf,const DVector3 &vertex, vector&dirchits, vector&tofpoints, vector&neutrals, vector&gammaHyps, unsigned int &num_in_insert, bool &tof_veto,bool &got_bcal, bool &got_dirc,double &unused_energy, double &fcal_total_E, bool &close_shower_veto ) const; bool DoKinematicFit(double t0_rf,const DVector3 &vertex, const DBeamPhoton *beamphoton, map >&chargedParticles, vector&gammaParticles, DKinFitUtils_GlueX *dKinFitUtils, DKinFitter *dKinFitter, bool isMissing) const; void GetKF4vectors(DKinFitter *dKinFitter,DLorentzVector &beam_kf, DLorentzVector &missing_kf, map>&final_kf) const; void FillGJHisto(double weight,const DLorentzVector &beam, const DLorentzVector &meson,const DLorentzVector &analyzer, TH2F *histo) const; void MakePipPimGroups(vector&pips_kf, vector&pims_kf, vector>>&groups); void Fill5GammaHistos(double unused_energy,const DLorentzVector &beam_kf, const DLorentzVector &pair1, const DLorentzVector &pair2, const DLorentzVector &gam5, double weight); void Fill5GammaHistos_pippim(double unused_energy,const DLorentzVector &beam_kf, const DLorentzVector &pippim, const DLorentzVector &pair1, const DLorentzVector &pair2, const DLorentzVector &gam5,double weight); void Fill3GammaHistos(double unused_energy,const DLorentzVector &pair1, const DLorentzVector &gam3, const DLorentzVector &gam4, double weight); void TwoGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void ThreeGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void FourGammaAnalysis(double fcal_total_E,double unused_energy,bool got_segments,bool got_dirc,bool got_bcal,bool tof_veto,bool got_itof,int num_in_insert,int numPseudos,bool close_shower_veto,const DLorentzVector &beam_kf, map>&final_kf, double weight); void FiveGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void SixGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void SevenGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void EightGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void TenGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void TwelveGammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void PipPimAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void PipPim1GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void PipPim2GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void PipPim3GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void PipPim4GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight); void PipPim5GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void PipPim6GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void TwoPipTwoPimAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void TwoPipTwoPim1GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void TwoPipTwoPim2GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void TwoPipTwoPim3GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight) const; void TwoPipTwoPim4GammaAnalysis(double unused_energy,const DLorentzVector &beam_kf, map>&final_kf, double weight); void MissingProtonAnalysis(const DLorentzVector &beam_kf, const DTrackTimeBased *proton_track, const DLorentzVector &missing_kf, map>&final_kf, double weight) const; void MakeTriplets(); bool MakeMesonPairs(vector&gammas_kf,double weight, TH2F *histo); static bool SortParticleProbability(const paira, const pairb){ return a.first>b.first; } const DParticleID* pid_algorithm; TH2F *DircFcalCorrelation; TH1F *DircFcalDiff; TH1F *PipPimPi0_missing_proton,*PipPimEta_missing_proton; TH1F *PipPimEtaPrime_missing_proton,*Eta_3piq_ThetaLab; TH2F *Mandelstam_t_eta_3piq_missing_proton; TH2F *Mandelstam_t_etaprime_c_missing_proton; TH1F *MissingMass_off_2gamma,*TwoGamma_noproton,*TwoGamma_noproton_CLcut; TH1F *TwoGammaMass_kf_missing_proton,*Eta_ThetaLab,*TwoGammaTDiff; TH2F *TwoGamma_vs_t_noproton,*TwoGamma_PxDiff,* TwoGamma_PyDiff; TH1F *MissingMass_tdiffcut_CLcut; TH1F *MissingMass_tdiff_cut; TH1F *MissingMass_off_pippim2gamma,*MissingMass_off_pippim2gamma_CLcut; TH1F *TwoGamma_pippim_noproton_CLcut,*TwoGammaMass_pippim_missing_proton; TH1F *Eshower_BCAL_for_neutrons,*Eshower_FCAL_for_neutrons; TH2F *PipPim_vs_NPip_1G,*PipPi0_vs_t_N,*PipEta_vs_t_N,*PipEtaPrime_vs_t_N; TH2F *PipEta_vs_t_Ncut,*PipEta_vs_t_eta_Ncut,*PipEta_vs_t_pi_Ncut; TH2F *EtaCosThetaGF_npip2g_Ncut; TH1F *PipEta_Ncut; TH2F *Mandelstam_t_eta_missing_proton,*Mandelstam_t_pi0_missing_proton; TH2F *Mandelstam_t_etaprime_missing_proton,*EnergyBalance; TH1F *Pi02g_InsertCut; TH1F *TwoG_NPi3g,*ThreeGammaMass_NPi,*Pip3G_N,*Pi0G_NPi3g,*Pi0G_NPi3g_with_veto; TH1F *SixGammaMass_NPi,*ThreePi0Mass_NPi,*TwoPi0EtaMass_NPi; TH1F *Pip6G_N,*PipEta_N6g,*OmegaPip_N2g; TH1F *TwoEtaPi0Mass_NPi,*ThreeEtaMass_NPi; TH2F *TwoEtaPi0Dalitz_NPi,*Pi0gVsPi0g_NPi,*FourGamma2D_NPi6g; TH1F *MCMissingMass,*PipNCL,*Pip1G_N,*EtaPi0_N_eta_3pi; TH1F *OneCellEnergy,*Pi0N,*EtaN,*EtaPrimeN; TH2F *NeutrondTvsZ,*NeutrondXvsZ,*NeutrondYvsZ; TH2F *Pi0gVsPi0g,*FourGamma2D_6g; TH1F *TwoOmega_6g; TH2F *MCMissingMassVsEbeam,*OmegaDalitz1,*OmegaDalitz2; TH2F *Egamma_vs_t_a0,*Egamma_vs_t_f0,*Egamma_vs_t_s0,*Egamma_vs_t_a2; TH2F *Egamma_vs_t_f2,*ProtonThetaVsP_4g; TH2F *Egamma_vs_t_a0_acc,*Egamma_vs_t_f0_acc,*Egamma_vs_t_s0_acc,*Egamma_vs_t_a2_acc; TH2F *Egamma_vs_t_f2_acc,*ProtonThetaVsP_4g_acc; TH1F *TwoGamma_PKpKmPipPim,*PipPim_with_KsKs,*PPim_with_KsKs; TH2F *PipPim_vs_2Pi0,*PipPim_vs_2Eta,*PipPim_vs_EtaPi0; TH1F *PipPimMissingMass,*Pi0g_with_eta,*Pi0g_with_eta_bg; TH2F *RhoPimPVsTheta,*RhoPipPVsTheta,*RhoPipThetaVsPimTheta; TH1F *TwoGammaMassBCAL,*TwoGammaMassFCAL; TH2F *TwoEtaPi0Dalitz,*PipPim2D,*PipPimEtaMass_vs_Ebeam; TH1F *EtaCosTheta_a2_mc; TH1F *OmegaEta_pi0g,*OmegaEta_pi0g_bg; TH2F *EtaCosThetaH_a2_mc; TH1F *EtaCosTheta_a2; TH2F *EtaCosThetaH_a2; TH1F *PipPimDoca,*PipPim2gCL,*PipPim2gCL_acc; TH2F *PipPimFlightPath,*PipPimPoca,*PipPimMass_vs_p; TH1F *Pip2GammaMass; TH1F *LambdaLambdaBarMass; TH2F *PipPim2Gamma_PimVertexRvsZ,*PipPim2Gamma_PipVertexRvsZ; TH1F *TwoEplusTwoEminusCL,*TwoEplusTwoEminusMass_kf; TH1F *EplusEminusMass_1g,*OmegaEta_2pip2pim3g; TH1F *ThreePi0_10g,*NonLH2_missing_E; TH1F *ThreeEtaMass_kf,*TwoEtaPi0Mass_kf; TH1F *EightGammaKinFitCL; TH1F *TwoEta_2pip2pim,*TwoPi0_2pip2pim,*EtaPi0_2pip2pim; TH2F *PipPim_t; TH2F *MCProtonThetaVsP,*ProtonThetaVsP; TH1F *TwoGammaMass_6g_kf; TH1F *MCProtonAcceptedP,*ThreePi0Mass_kf,*TwoPi0EtaMass_kf; TH1F *MissingMass_1track_1gamma,*TwoGammaMass_kf,*ThreeGammaMass_kf; TH1F *TwoGammaMass_m; TH2F *Mandelstam_t_pi0_kf_NU,*Mandelstam_t_eta_kf_NU,*Mandelstam_t_etaprime_kf_NU; TH1F *TwoGammaMass_kf_NU; TH2F *TwoGammaMass_vs_Ebeam_NU,*TwoGammaMass_vs_Ebeam_theta_cut; TH2F *Ekin_vs_Emeas; TH1F *ThreeGammaKinFitCL,*TwoGammaMass_3g_kf,*OneGammaKinFitCL; TH1F *TwoGammaMass_ppippim6g,*MM2_p1gam,*FourGammaKinFitCL; TH1F *FourGamma_Pi02g,*FourGamma_Eta2g,*EtaPrimeEta_4g; TH2F *EtaCosThetaGF,*Pi0CosThetaGF,*EtaPrimeCosThetaGF_npip2g; TH2F *EtaCosThetaGF_npip2g,*Pi0CosThetaGF_npip2g; TH2F *TwoGamma_Pi02g,*TwoGamma_Pi02g_NoBCAL, *TwoGamma_Eta2g; TH2F *TwoGamma_Pi02g_NoBCAL_best; TH2F *Pi0g_Pi02g,*Pi0g_Pi02g_NoBCAL; TH1F *Pi0XMass_NoTOF, *Pi0XMass_NoBCAL_NoTOF,*Pi0XMass_NoBCAL_PseudoVeto; TH1F *Pi0XMass_NoITOF, *Pi0XMass_NoBCAL_NoITOF; TH1F *Pi0XMass_InsertCut_PseudoVeto,*Pi0XMass_PseudoVeto; TH1F *Pi0XMass_NoSegments,*Pi0XMass_NoBCAL_NoSegments; TH1F *Pi0XMass_NoBCAL,*Pi0XMass_InsertCut_NoBCAL,*Pi0XMass_InsertCut; TH1F *Pi0XMass_NoBCAL_NU,*Pi0XMass_InsertCut_NoBCAL_NU,*Pi0XMass_InsertCut_NU; TH1F *Pi0XMass_NU,*TwoPi0Mass_kf_NU,*EtaPi0Mass_kf_NU,*TwoEtaMass_kf_NU; TH1F *TwoPi0Mass_PseudoVeto; TH1F *EtaPrimePi0_4g_NU; TH2F *EtaPrimePi0Mass_t_NU,*TwoEtaMass_t_NU,*EtaPi0Mass_t_NU,*TwoPi0Mass_t_NU; TH2F *TwoPi0Mass_vs_Ebeam_NU,*EtaPi0Mass_vs_Ebeam_NU,*TwoEtaMass_vs_Ebeam_NU; TH2F *EtaCosThetaGF_NU,*Pi0CosThetaGF_NU; TH1F *Pi0XMass_NoBCAL_NoDIRC_NU,*Pi0XMass_InsertCut_NoBCAL_NoDIRC_NU; TH1F *Pi0XMass_InsertCut_NoDIRC_NU,*Pi0XMass_NoDIRC_NU; TH1F *FourGammaMass_NoBCAL,*FourGammaMass_InsertCut_NoBCAL; TH1F *FourGammaMass_kf,*ProtonGammaMass_kf,*FourGamma_NoMesonPairs; TH2F *FourGamma2d_kf,*FourGamma2d_bg,*FourGamma2d_NoMesonPairs; TH1F *EtaPi0Prob_kf,*TwoPi0Prob_kf, *TwoPi0Mass_kf,*EtaPi0Mass_kf; TH1F *EtaPi0Mass_bg,*TwoPi0Mass_bg; TH1F *FiveGammaKinFitCL,*FiveGammaMass_kf,*Pi0GammaMass_kf; TH1F *SixGammaKinFitCL,*SixGammaMass_kf; TH1F *EtaPrimePi0_8g,*EtaPrimePi0_4g,*Pi0Gamma_kpkm,*PhiOmegaMass_with_pi0; TH1F *TwoGammaMass_kpkm3g; TH2F *MissingMass_pippim6g_vs_Ebeam,*MissingMass_pippim5g_vs_Ebeam; TH1F *FiveGammaMass_ppippim,*SixGammaMass_ppippim,*FourGammaMass_ppippim; TH1F *Mandelstam_t_pi0,*Mandelstam_t_eta,*Sigma0GammaMass; TH2F *Mandelstam_t_etaprime_kf; TH1F *EplusEminusMissingGammaCL,*EtaPrime_sinphicosphi,*Eta_sinphicosphi; TH2F *SixGammaVertexRvsZ,*PipPim2GammaVertexRvsZ; TH1F *EplusEminusMass_with_2g,*EplusEminusMass_PipPim; TH1F *EplusEminusMass_with_pi0,*EplusEminusMass_with_eta; TH1F *EplusEminusMass_PipPim_eta_cut,*EplusEminusMass_PipPim_etaprime_cut; TH2F *Compton_CosThetaCM_vs_Ebeam_kf; TH1F *Pi0GammaMass_4g_kf; TH1F *BeamPhotonTimeDiff; TH2F *PiMinusdEdxFDCCut,*BeamPhotonProtonTimeDiff; TH1F *Pi0Pi0Gamma_7g,*FCAL2ShowerDistance; TH2F *ProtondEdxFDC,*PPi0Dalitz,*SCdIddT,*SCdIddTLowLevel; TH2F *ProtondEdxFDCwithTOFcut; TH2F *SCdIddTLowLevelHasFADC; TH2F *ProtondEdx,*NormProtondEdx,*DProtondEdx; TH1F *SCt; TH2F *dProtondEdx,*ProtondEdxCut,*ProtondEdxCutFDC; TH2F *ProbProtondEdx,*ProbKPlusdEdx,*ProbKMinusdEdx; TH1F *Pi0CosThetaH_f2_mc,*Pi0CosTheta_f2_mc; TH2F *ProtonTOFdtNodEdx,*ProtonSCtCorr,*ProtonSCt; TH1F *PPimMass,*PipPimMass,*PPimMass_TaggedPhoton,*PipPimMass_TaggedPhoton; TH1F *TimeDiffTaggerTracks,*TDiffGammaProton,*NumGammaPerEvent,*NumGammaPerEventMM2Cut; TH2F *ProtonTOFdt,*ProtonBCALdt,*ProtonBCALdtForward,*ProtonBCALdtWithdEdxCut; TH2F *ProtonBCALdtWithFDCdEdxCut,*ProtonFCALdt; TH2F *ProtonTOFdtWithdEdxCut,*ProtonFCALdtWithdEdxCut; TH1F *TwoGammaMassMissingE_with_CLcut; TH1F *PipPimPi0_no_proton,*PipPimEta_no_proton; TH1F *PipPimPi0_no_proton_cut2,*PipPimEta_no_proton_cut2; TH1F *TwoGammaMass_pippim_no_proton,*PiPlusPiMinusMass_epem; TH1F *TwoGammaMass_pippim_no_proton_cut2; TH1F *LambdaPip_3pip3pim_Rcut; TH2F *PipPim2D_3pip3pim_Rcut; TH2F *PiPlusdEdxFDC,*PiPlusTOFwithFDCdEdxCut,*PiPlusBCALTOFwithFDCdEdxCut; TH2F *PiPlusBCALdt,*PPimMassVsR; TH1F *TrigMask,*Protont_flagged; TH1F *FourGammaMassMissingE_with_CLcut; TH1F *PPimPi0Mass,*RhoGammaMass,*PipPim1GammaMass_cut,*RhoGammaMass_4pi; TH1F *TwoGammaMassWithPipPim,*TwoGammaMassWithPPim,*PipPimPi0Mass; TH2F *ProbProtonTOF,*NormProtonTOF; TH2F *PiPlusTOFwithdEdxCut,*ProtondEdxSC,*ProtondEdxSCNoCDCdEdx; TH2F *ProtondEdxAmp; TH1F *Protont; TH2F *ProbPiMinusTOF,*NormPiMinusTOF; TH2F *ProbPiPlusTOF,*NormPiPlusTOF,*PiPlusBCALTOFwithdEdxCut; TH2F *ProtondEdxFDCWithSC; TH1F *PipPim2Gamma_PimVertexZ_cut,*PipPim2Gamma_PipVertexZ_cut; TH2F *NormPiPlusdEdx,*NormPiMinusdEdx,*PiMinusdEdxFDC; TH2F *ProbPiPlusdEdx,*ProbPiMinusdEdx,*PiMinusBCALdt; TH1F *PipPim3GammaMass_cut_kf; TH2F *dPiPlusdEdx,*dPiMinusdEdx,*OmegaPiDalitz_kf; TH2F *OmegaPiDalitz_2pip2pim; TH2F *PiPlusdEdx,*PiMinusdEdx,*PiPlusdEdxCut,*PiMinusdEdxCut; TH2F *PiPlusdEdxFDCCut; TH1F *PipPimEtaMass_kf,*PipPimPi0Mass_kf,*KstarKstar; TH2F *Omega_t_mandelstam; TH2F *Mandelstam_t_omega_kf,*Mandelstam_t_eta_kf,*Mandelstam_t_pi0_kf; TH2F *MCt_vs_E,*MCt_vs_Ecoarse,*MCE_vs_t,*MCM_vs_t,*MCt2d,*MCt2d2; TH1F *PipPim1GammaMass_kf,*PipPim6GammaMass_kf,*PipPimEtaMass_6g_kf; TH1F *PipPimEtaPrime_6g,*B1Pip,*B1Pim; TH2F *Pi0GammaDalitz_kf; TH1F *PipPim5GammaKinFitCL,*FiveGammaMass_ppippim_kf,*PipPim5GammaMass_kf; TH1F *TwoGammaMass_no_tag,*OmegaGammaMass_4g_kf; TH1F *MCBeam,*MissingMass_2pi_cut1,*MissingMass_2pi_cut2,*MCBeamTagged; TH2F *Pi02GammaMass_vs_Pi0GammaMass_kf; TH1F *Pi02gMass_kf,*Pi02gMass_acc,*Eta2gMass_kf; TH2F *Pi02GammaMass_vs_Pi0GammaMass_acc; TH2F *MissingE_vs_6GammaMass; TH2F *MissingE_vs_5GammaMass; TH2F *MissingE_vs_4GammaMass; TH2F *MissingE_vs_3GammaMass; TH2F *MissingE_vs_2GammaMass; TH2F *MissingE_vs_8GammaMass; TH2F *PipPim2GCL_vs_Ptheta; TH2F *MissingE_vs_7GammaMass; TH1F *PipPimEtaPrimeMass_2pip2pim; TH1F *EplusEminusPiplusPiminusCL,*EplusEminusPiplusPiminusMass; TH1F *EplusEminusPiplusPiminusMass_cut; TH1F *TwoGammaMass_epem3g,* EplusEminus3GammaMass_kf,* EplusEminus3GammaCL; TH1F *EplusEminusMass_with_3g; TH1F *EplusEminusPi0GammaMass_kf; TH2F *MesonPairMultiplicity; TH1F *TwoEtaMass_fcal,*TwoPi0Mass_fcal,*EtaPi0Mass_fcal,*Pi02Gamma_fcal; TH1F *TwoEtaMass_fcal_acc,*TwoPi0Mass_fcal_acc,*EtaPi0Mass_fcal_acc,*Pi02Gamma_fcal_acc; TH2F *EOverPMinus,*EOverPPlus,*EOverPMinus_BCAL,*EOverPPlus_BCAL; TH2F *Pim_EOverP_with_timecut,*Pip_EOverP_with_timecut; TH1F *TwoGammaMass_eplus_eminus,*EplusEminusPi0Mass_kf; TH1F *EplusEminusPi0Mass_cut,*EplusEminusEtaMass_cut; TH2F *OneGammaE_vs_Ebeam_kf,*KsPi0_vs_KsPi0; TH1F *Mandelstam_t_1gamma_kf; TH1F *TwoPipTwoPim4GammaMass_with_2Ks; TH1F *PipPimEtaPi0Mass_bg,*OmegaEta_PipPim_bg; TH2F *EtaCosTheta_EtaOmega_bg,*Pi0CosTheta_EtaPrimePi0_bg; TH1F *Pi0PipPim_with_Eta_bg,*EtaPipPim_with_Pi0_bg,*F1Pi0Mass_bg; TH2F *PipPim_vs_EtaPi0_bg; TH1F *EtaPrimePi0_PipPim_bg,*TwoEta_PipPim_bg; TH2F *EtaCosTheta_EtaOmega_pi0g,*EtaCosTheta_EtaOmega_pi0g_bg; TH1F *EtaPi0Mass_p2pi_bg,*OmegaPi0Mass_bg; TH2F *PipPim_vs_2Pi0_bg,*PipPim_vs_2Eta_bg,*PipPi0Mass_vs_PimPi0Mass_bg; TH2F *KsPi0_vs_KsPi0_bg,*PipPimPi0_vs_PipPimPi0_bg; TH2F *PipPimEta_vs_PipPimEta_bg,*PipPimEta_vs_PipPimPi0_bg; TH1F *TwoPi0_2pip2pim_bg,*EtaPi0_2pip2pim_bg,*TwoEta_2pip2pim_bg; TH1F *OmegaOmega_bg,*EtaOmega_bg,*EtaPrimeOmegaMass_bg; TH2F *EtaCosTheta_EtaOmega_2pip2pim4g_bg; TH1F *PipPim2Pi0Mass_bg,*PipPim2EtaMass_bg,*Pi0PipPim_with_Pi0_bg; TH1F *EtaEtaPrimeMass_bg,*EtaPipPim_with_Eta_bg; TH1F *TwoPi0_NPi_bg,*EtaPi0_NPi_bg,*TwoEta_NPi_bg; TH1F *Pip2Pi0_N_bg,*Pip2Eta_N_bg,*PipEtaPi0_N_bg; TH2F *PipPi0_PipPi0_Dalitz_N_bg,*PipEta_PipPi0_Dalitz_N_bg; TH2F *PipEta_PipEta_Dalitz_N_bg; TH1F *TwoPi0_N3Pi_bg,*EtaPi0_N3Pi_bg,*PipPimPi0_N_with_Pi0_bg; TH1F *PipPimPi0_N_with_Eta_bg,*PipPimEta_N_with_Pi0_bg; TH2F *Omega2PI_N_Dalitz_bg; TH1F *EtaPi0_N_eta_3pi_bg; TH1F *TwoEta_NPi,*Pip2Eta_N; TH2F *PipEta_PipEta_Dalitz_N; TH2F *PipPim_vs_NPip_with_2pi0_bg,*PipPim_vs_NPip_with_etapi0_bg; TH2F *MissingE_vs_pippim6gMass,*FiveGammaMass_vs_Ebeam; TH2F *MissingE_vs_pippim5gMass; TH2F *MissingE_vs_pippim4gMass; TH2F *MissingE_vs_pippim3gMass; TH2F *MissingE_vs_pippim2gMass; TH2F *MissingE_vs_pippim1gMass; TH1F *PipPim2EtaMass; TH1F *FourGammaMass_ppippim_kf; TH1F *ThreeGammaMass_ppippim_kf,*TwoPi0Prob_pippim4g_kf; TH1F *Pi0Gamma_epem2g,*EtaGamma_epem2g; TH1F *Pi0Gamma_epem2g_cut,*EtaGamma_epem2g_cut; TH1F *PipPim3GammaMass_kf,*TwoGammaMass_pippim3g_kf; TH1F *PipPim3GammaKinFitCL; TH1F *TwoPi0Prob_5g_kf; TH2F *FourGamma2d_ppippim_kf,*FourGamma2d_ppippim_bg; TH2F *Pi0Pi0GammaMass_vs_Pi0GammaMass_kf,*FourGamma2d_5g_kf; TH1F *EplusEminus1GammaMass_with_2g; TH1F *PipPim4GammaKinFitCL,*Pi0GammaMass_ppippim_kf; TH1F *KpKmPi0Mass_baryon_cut; TH2F *PipPi0Mass_vs_PimPi0Mass_kf; TH1F *PipPim6GammaKinFitCL, *SixGammaMass_ppippim_kf; TH1F *EtaGammaMass_kf; TH1F *OmegaOmega_2pi5g; TH1F *PipPim2GammaMass_kf,*TwoGammaMass_p2pi_kf; TH2F *PipPim2g_vs_2g; TH1F *PipPim1GammaKinFitCL; TH3F *TwoGammaPair3d; TH1F *Eta2Pi0Prob,*PipPimKinFitCL,*PipPimMass_kf; TH1F *TwoPi0Prob,*EtaPi0Prob,*TwoPi0Prob_5g,*ThreePi0Prob; TH2F *MissingE_vs_MM2_ppippim6g; TH2F *MissingE_vs_MM2_ppippim5g; TH2F *MissingE_vs_MM2_ppippim4g; TH2F *MissingE_vs_MM2_ppippim3g; TH2F *MissingE_vs_MM2_ppippim2g; TH2F *MissingE_vs_MM2_ppippim1g; TH2F *MissingMass_pippim1g_vs_Ebeam,*PipPim1GammaMass_vs_Ebeam; TH2F *MissingMass_pippim2g_vs_Ebeam,*PipPim2GammaMass_vs_Ebeam; TH2F *MissingMass_pippim3g_vs_Ebeam,*PipPim3GammaMass_vs_Ebeam; TH2F *MissingMass_pippim4g_vs_Ebeam,*PipPim4GammaMass_vs_Ebeam; TH2F *MissingMassSq_for_2gam; TH1F *MissingMassSq,*MissingMassSqWithPipPim,*MissingMassSqOffPipPim; TH1F *MissingMassSqOffPipPimCut1,*MissingMassSqOffPipPimCut2; TH1F *MissingMassSqPCut; TH1F *MissingMassSqPCut1; TH1F *MissingMassSqPCut2; TH2F *TwoGammaVertex,*TwoGammaVertexCLcut; TH1F *PipPimEtaPi0Mass,*TwoGammaMass_off_pi0; TH2F *PipPi0Mass_vs_PimPi0Mass; TH1F *PipPim2GammaKinFitCL,*TwoGammaMass_p2pi_cut,*PipPim2GammaMass_cut; TH1F *TwoGammaKinFitCL,*TwoGammaMass_cut,*MissingMass_2gam_cut; TH1F *EtaGammaMass_7g; TH2F *TwoGammaMass_vs_ThreeGammaMass; TH1F *TwoGammaMass_omega_cut; TH1F *Pi0Prob_pippim,*EtaProb_pippim,*Pi0Prob_pippim3g,*EtaProb_pippim3g; TH1F *TwoPi0Prob_pippim4g,*EtaPi0Prob_pippim4g; TH2F *MissingMass_1track_gammas_vs_Ebeam; TH2F *MissingMass_1track_3gammas_vs_Ebeam; TH2F *MissingMass_1track_4gammas_vs_Ebeam; TH2F *MissingMass_1track_5gammas_vs_Ebeam; TH2F *MissingMass_1track_6gammas_vs_Ebeam; TH2F *MissingMass_1track_7gammas_vs_Ebeam; TH2F *MissingMass_1track_8gammas_vs_Ebeam; TH2F *TwoPi_2d,*ThreePiMassVsTwoPiMass; TH1F *MissingMass_3pi0_cut,*SixGammaMass_3pi0_cut; TH1F *MissingMass_eta_2pi0_cut,*SixGammaMass_eta_2pi0_cut; TH2F *TwoPi0Dalitz,*Pi0Pi0GammaMass_vs_Pi0GammaMass; TH1F *MissingMass_etapi0_cut,*MissingMass_2pi0_cut; TH1F *Pi0_cl_3g,*Eta_cl_3g; TH1F *MissingMass_1track_eta,*MissingMass_1track_eta_sb,*MissingMass_1track_pi0; TH1F *PhiDiff_p2gam,*PhiDiff_p3gam,*PhiDiff_p4gam,*PhiDiff_p5gam,*PhiDiff_p6gam; TH1F *ThreeGamma_2pip2pim,*TwoGamma_2pip2pim3g; TH1F *Pi0Gamma_2pip2pim,*EtaGamma_2pip2pim; TH1F *ThreePi0_against_eta,*Eta2Pi0_against_eta,*Eta2Pi0_against_pi0; TH1F *TwoPipTwoPim4GammaMass,*TwoPipTwoPim4GammaKinFitCL; TH1F *PhiDiff_p8gam,*PhiDiff_p7gam; TH2F *TwoGamma_vs_PipPimTwoGamma; TH1F *PhiDiff_pippim1gam,*PhiDiff_pippim2gam,*PhiDiff_pippim3gam; TH1F *PhiDiff_pippim4gam,*PhiDiff_pippim5gam,*PhiDiff_pippim6gam; TH2F *FourGamma2d,*FourGamma2d_ppippim,*ThreeGammaDalitz,*ProtonRvsZ; TH2F *FourGamma2d_5g,*FourGamma2d_2pip2pim; TH1F *OmegaEta,*PhiOmegaMass,*OmegaEtaSB; TH2F *PipRvsZ,*PimRvsZ; TH2F *FourGammaMassVsEbeam; TH1F *TwoPi0Mass,*EtaPi0Mass,*ProtonTwoGammaMissingMassSq; TH1F *TwoPi0Mass_acc,*EtaPi0Mass_acc,*TwoEtaMass_acc,*EtaPrimePi0_4g_acc; TH1F *OmegaGammaMass_4g_acc; TH2F *Pi0GammaDalitz_acc,*FourGamma2d_acc; TH2F *TwoPi0Mass_vs_Ebeam_acc,*TwoEtaMass_vs_Ebeam_acc,*EtaPi0Mass_vs_Ebeam_acc; TH1F *EMeasuredP4Gam,*EMeasuredPPipPimPi0; TH2F *Pi02GammaMass_vs_Ebeam_acc; TH1F *TwoGammaMass_PipPimNoBeam,*PipPimPi0_NoBeam; TH2F *ThreeGammaMass_vs_Ebeam,*SixGammaMass_vs_Ebeam; TH2F *EightGammaMass_vs_Ebeam,*SevenGammaMass_vs_Ebeam; TH2F *MissingE_vs_MM2_p8gam,*MissingE_vs_MM2_p7gam; TH1F *MissingMass_1track_3gammas_cut,*MissingMassSq_1gam; TH1F *MissingMass_1track_3gammas_no_missing_gamma; TH1F *FourGammaMass_with_pi0,*MissingMass_1track_4gammas_with_pi0; TH1F *MM2_p3gam,*MM2_p4gam,*ThreeGammaMass_with_pi0,*MissingMass_1track_gampi0; TH1F *ProtonFourGammaMissingMassSq,*Proton2PiMissingMassSq; TH1F *MissingMass_1track_gameta,*ThreeGammaMass_with_eta; TH1F *proton_pi0_mass,*proton_pim_mass,*TwoGammaMass_3gam,*TwoGammaMass_4gam; TH1F *TwoGammaMass_3gam_cut; TH1F *MM2_p4pi_0gam,*MM2_p4pi_2gam,*MissingMass_with_4pi_0gam; TH1F *MissingMassSq_p2pi_0gam,*MissingMass_2pi_0gam; TH1F *MissingMassSq_p2pi2g; TH1F *MissingMass_2pi_1gam,*MM2_p2pi1gam; TH1F *MissingMass_1track_3gammas,*MissingMass_1track_4gammas,*FourGammaMass; TH1F *MissingMass_1track,*MissingMass_1track_gammas; TH2F *TwoGammaMass_vs_Ebeam; TH1F *MissingMass_1track_gammas_cut,*ThreeGammaMass_4g; TH1F *MissingMass_2pi,*MissingMassSq_p2pi; TH1F *MissingMass_2pi_2gam,*TwoGammaMass_p2pi,*MissingMass_2pi_pi0; TH1F *TwoGammaMass,*MM2_p4pi,*TwoGammaMass_with_4pi,*FourGammaMass_no_kf; TH1F *PimPredMeasPhiDiff,*PipPredMeasPhiDiff; TH1F *PimPredMeasThetaDiff,*PipPredMeasThetaDiff; TH1F *MissingMass,*InvMass,*MM2_p2gam,*MissingMassCut,*InvMass_no_tag; TH1F *MissingMass_with_4pi,*SixGammaMass,*missing_E_p6gam,*TwoGammaMass_6g; TH1F *MM2_p6gam,*MissingMass_1track_6gammas; TH1F *MMppim,*MMppip,*FiveGammaMass,*FiveGammaMass_pi0_omega_cut; TH1F *PipPimMass_etaprime,*PipPimMass_etaprime_sb; TH1F *ThreeGammaMass,*MissingMass_with_4pi_2gam; TH1F *ThreeGammaMass_ppippim,*TwoGammaMass_ppippim3g,*ThreeGammaMass_vs_pi0; TH1F *MM2_p5gam,*MissingMass_with_pippim3g,*MissingMass_with_pippim4g; TH1F *pippim1g_mass_cut,*MissingMass_pippim1g_cut; TH2F *TwoGammaMassVsThreeGammaMass; TH1F *TwoGammaMass_5g,*MissingMass_1track_5gammas,*missing_E_p5gam; TH1F *MissingMass_1track_5gammas_no_missing_gamma; TH1F *PipPimEtaMass,*PipPimPi0Mass_3g; TH1F *PipPimMass_1g,*PipPimMass_1g_etaprime; TH1F *PPipPim3GammaMissingE,*PPipPim3GammaMissingMassSq; TH1F *PPipPim4GammaMissingE,*PPipPim4GammaMissingMassSq; TH1F *OmegaGammaMass,*OmegaSBGammaMass,*OmegaPi0Mass,*OmegaSBPi0Mass; TH1F *missing_E_p2gam,*missing_E_p3gam,*missing_E_p4gam; TH1F *missing_E_p2pi1gam,*missing_E_p2pi2gam; TH1F *Egamma; TH2F *MissingE_vs_MM2_p6gam,*MissingE_vs_MM2_p5gam,*MissingE_vs_MM2_p4gam; TH2F *MissingE_vs_MM2_p3gam,*MissingE_vs_MM2_p2gam; TH1F *PipPimGammaMass_3g_kf,* PipPimPi0Mass_3g_kf,*PipPimEtaMass_3g_kf; TH1F *PPimMass_2g,*Pip2GammaMas,*FourGammaMassMissingE; TH1F *FourGammaMassMissingPt,*TwoGammaMassMissingE,*TwoGammaMassMissingPt; TH1F *PipPimTimeDiff; TH2F *FourGammaMassT,*FourGammaMassT_kf,*TwoGammaMassT,*TwoPi0VsTwoPi0; TH2F *TwoGammaMassT_kf,*TwoGammaMassTheta_kf; TH1F *KShortKShort_4pi0; TH2F *PipEta_vs_t_eta_N,*PipEta_vs_t_pi_N; TH1F *TwoGammaMass_missingProton; TH2F *TwoGammaMass_vs_theta_missingProton,*TwoGammaMass_vs_t_missingProton; TH1F *PPim_3pip3pim; TH2F *PipPim_vs_PipPim_with_Lambda; TH1F *PipPim_with_KsKs_sb,*PPim_with_KsKs_sb; TH2F *KsPipDalitz_with_lambda_sb,*KsKs_vs_LambdaPip_sb,*KsKs_vs_pKs_sb; TH2F *KsKs_vs_LambdaPip,*KsKs_vs_pKs,*KsPipDalitz_with_lambda; TH2F *KsKs_vs_LambdaPip_1g,*KsKs_vs_pKs_1g,*KsPipDalitz_with_lambda_1g; TH2F *KsKs_vs_LambdaPip_2g,*KsKs_vs_pKs_2g,*KsPipDalitz_with_lambda_2g; TH1F *LambdaK0_3pip3pim,*K0Pip_3pip3pim,*LambdaPip_3pip3pim; TH1F *Sigma0K0_3pip3pim,*K0Pip_3pip3pimg,*Sigma0Pip_3pip3pim; TH1F *PipPim_with_LambdaK0Pip; TH1F *PipPim_with_Sigma0K0Pip; TH2F *PipPim_vs_PipPim_with_Sigma0; TH1F *LambdaGamma_3pip3pim,*PPim_with_KsKs_1g,*PipPim_with_KsKs_1g; TH1F *PipPim_3pip3pim_1g,*ThreePipThreePimMass_1g; TH2F *PipPim2D_1g,*ThreePi_vs_ThreePi_1g; TH1F *LambdaGamma_3pip3pim1g,*PPim_3pip3pim1g; TH2F *PipPim_vs_PipPim_with_Lambda1g; TH1F *TwoGamma_3pip3pim,*PPim_with_KsKs_2g,*PipPim_with_KsKs_2g; TH1F *PipPim_3pip3pim_2g,*ThreePipThreePimMass_2g; TH2F *PipPim2D_2g,*ThreePi_vs_ThreePi_2g; TH1F *OmegaGammaMass_pippim3g,*EtaGammaMass_pippim3g,*PhiGammaMass_pippim3g; TH1F *ProtonGammaTimeDiff; TH2F *FourGammaMass_all_fcal,*ThreePi_vs_ThreePi; TH2F *ThetaGamma_4g,*TwoPi0MassTheta,*TwoGammaMassTheta; TH1F *PipPimEtaPrimeMass_kf; TH1F *TwoEtaMass_kf,*TwoEtaMass_bg,*TwoEtaProb_kf; TH2F *PipEtaMassSq_vs_PimEtaMassSq_f1_low_sb; TH2F *PipEtaMassSq_vs_PimEtaMassSq_f1cut; TH2F *PipEtaMassSq_vs_PimEtaMassSq_f1_high_sb; TH2F *Pi0Pi0GammaMass_vs_Pi0GammaMass_pippim,*FourGamma2d_5g_pippim; TH1F *TwoPi0Prob_5g_pippim; TH1F *NumSCHits_for_3gam,*DTimeSC; TH2F *DeltaPhiVsDeltaTheta2; TH1F *OmegaStudy_TwoGammaMass; TH2F *MissingMassVsEbeam,*MissingMassVsEbeam_with_pi0,*MissingMassSqVsEbeam; TH1F *OmegaPiPiMass,*EtaGammaMass_ppippim_kf,*PipPimGammaMass_against_eta; TH2F *PipPimMassVsTwoPi0Mass,*ThreeGammaMass2d,*ThreeGammaMass2dCut; TH2F *PredPhiTheta2,*PipEtaPrimeMassSq_vs_PimEtaPrimeMassSq; TH2F *PipEtaMassSq_vs_PimEtaMassSq; TH2F *ReconPhiTheta; TH1F *TDiff2Gamma,*TDiff2GammaBeam,*EDiff2GammaBeam,*TwoGammaMissingMass; TH2F *TwoGammaMassVsTheta,*FourGammaMass_vs_Ebeam; TH1F *TwoPi0Mass_high,EtaPi0Mass_high,EtaEtaMass_high; TH1F *EightGammaMass_kf,*TwoGammaMass_8g; TH1F *FourPi0Mass,*OmegaOmega,*EtaOmega; TH1F *PPipPimPi0_MM2,*PPipPimPi0_missing_E,*PPipPimPi0_missing_E_cut; TH1F *SevenGammaMass_kf,*TwoGammaMass_7g,*ThreePi0Mass_7g,*SevenGammaKinFitCL; TH1F *PipPimMass_2g,*PipPimMass_with_pi0,*PipPimMass_with_eta; TH1F *FourNeutralTDiff,*TwoNeutralTDiff,*Beam2GammaTDiff,*Beam4GammaTDiff; double SPLIT_CUT,NUM_SIGMA_BG; TH1F *EtaDalitzZ; TH2F *EtaDalitzXY; TH2F *TwoGammaMass_vs_E2gamma; TH2F *Pi02GammaMass_vs_Ebeam, *EtaPi0Mass_vs_Ebeam,*TwoPi0Mass_vs_Ebeam; TH2F *TwoEtaMass_vs_Ebeam; TH2F *BCAL2GammaMassVsE_v2,*FCAL2GammaMassVsE_v2; TH1F *Pi0Gamma_with_epem; TH2F *PipPimMass_vs_Ebeam,*PipPimPi0Mass_vs_Ebeam; TH1F *TwoPipTwoPimKinFitCL,*TwoPipTwoPim1GammaKinFitCL; TH1F *TwoPipTwoPim2GammaKinFitCL; TH1F *TwoPipTwoPimMass,*TwoPipTwoPim1GammaMass,*TwoPipTwoPim2GammaMass; TH2F *PipPimMass_vs_PipPimMass,*PipPimMass_vs_PipPimMass_1g; TH2F *PipPimMass_vs_PipPimMass_with_omega; TH2F *PipPimMass_vs_PipPimMass_4g; TH2F *PipPimMass_vs_PipPimMass_2g; TH1F *PipPim1GammaMass_4pi,*PipPim2GammaMass_4pi,*TwoGammaMass_4pi; TH1F *PipPimPi0Mass_4pi,*PipPimEtaMass_4pi,*PipPimEtaMass_Eta_to_3pi; TH1F *PipPimPi0Mass_4pi3g,*PipPimEtaMass_4pi3g; TH2F *EtaPi0_2d,*EtaPi0_vs_Pi0Pi0; TH2F *Proton4GammaMissingE_vs_dphi_no_kf; TH2F *PipPimMass_NoKinFit,*PhiPi0Pi0Dalitz; TH1F *MissingMassSq_NoKinFit; TH2F *PipPimPi0_vs_t[48]; TH2F *EtaCosTheta_EtaOmega,*Pi0CosTheta_EtaPrimePi0; TH2F *EtaCosTheta_EtaOmega_5g,*EtaCosTheta_EtaOmega_2pip2pim4g; TH2F *EtaCosTheta_EtaOmega_2pi5g; TH2F *Egamma_vs_t_a0_bg,*Egamma_vs_t_f0_bg,*Egamma_vs_t_a2_bg; TH2F *Egamma_vs_t_f2_bg,*Egamma_vs_t_s0_bg,*EtaCosThetaGF_bg; TH2F *Pi0CosThetaGF_bg,*PPi0Dalitz_bg,*PPi0EtaDalitz_bg; TH1F *EtaPrimePi0_4g_bg,*EtaPrimeEta_4g_bg; TH2F *Proton6GammaMissingE_vs_dphi_no_kf; TH1F *SixGammaMass_no_kf,*OmegaEta_2pi5g; TH1F *MCProtonAcceptance,*MCProtonP,*BeamEnergy,*EtaCosTheta_phi_mc; TH1F *MCt,*MCtCut; TH2F *MCdeltat; TH1F *ComptonDPhi; TH1F *TwoGammaMass_missingPip; TH2F *DeltaPOverP_vs_P; TH2F *DeltaPhi_vs_P,*DeltaTheta_vs_P,*DeltaPOverP_vs_theta; TH2F *DeltaPx_vs_Phi,*DeltaPy_vs_Phi,*DeltaPz_vs_Phi; TH2F *PipPimEta_vs_PipPimEta, *PipPimEta_vs_PipPimPi0, *PipPimPi0_vs_PipPimPi0; TH1F *TwoPipTwoPim3GammaKinFitCL,*TwoPipTwoPim3GammaMass; TH2F *Pi0Gamma_vs_PipPimPi0,*PipPimEtaPrimeDalitz; TH1F *EtaEta_3pi0,*ThreePi0_4pi0,*Pi0Eta_eta_3pi0; TH1F *KShortKShort,*KShortKShortWithPi0,*KShortKShort_baryon_cut; TH2F *EplusEminusMomenta; TH2F *EplusEminusNoHadrons_dEcut; TH1F *EplusEminusMass,*EplusEminusNoHadrons,*EplusEminusBeamDeltaE; TH1F *EplusEminusGamma,*TwoGamma_eplus_eminus,*TwoGammaEplusEminus; TH1F *EplusEminusGammaWithProton,*EplusEminusMassWithProton; TH1F *TwoGamma_eplus_eminus_with_proton,*TwoGammaEplusEminusWithProton; TH1F *EplusEminusGammaCL,*EplusEminus2GammaCL,*EplusEminusCL; TH1F *EplusEminus1GammaMass_kf,*EplusEminus2GammaMass_kf,*EplusEminusMass_kf; TH1F *TwelveGammaKinFitCL,*TenGammaKinFitCL; TH1F *FivePi0Mass,*SixPi0Mass,*TwoPipTwoPimEta; TH1F *TenGammaMass_kf,*TwoGammaMass_10g; TH1F *TwelveGammaMass_kf,*TwoGammaMass_12g; TH1F *EplusEminusEtaMass_kf,*Eta4Pi0Mass; TH2F *FourGammaVertex; TH1F *EtaPrime2Pi0_10g; TH1F *EtaPrime2Pi0,*RhoPlusRhoMinus,*RhoPlusRhoMinus_bg,*RhoRho; TH1F *Pi0Pi0_against_omega,*Proton4GammaMMsq_no_kf; TH1F *EtaGamma_with_2pi0,*PipPim_3pip3pim,*Eta2Pi0_10g; TH1F *LambdaPimMass_2g; TH1F *Pi0CosTheta_f2,*LambdaKmMass_2kpkm,*LambdaPimPi0Mass; TH1F *FourGammaMass_acc,*Pi0GammaMass_4g_acc; TH2F *FourGammaMass_all_fcal_acc; TH1F *TwoGammaMass_acc,*TwoPipTwoPim2Pi0; TH1F *Pi0XMass,*Pi0XMass_bg,*Pi0XMass_HarderVeto,*Pi0XMass_HarderVeto_NoBCAL; TH1F *Pi0XMass_PsuedoVeto,*Pi0XMass_InsertCut_PsuedoVeto; TH1F *Pi0XMass_InsertCut_NoBCAL_PseudoVeto; TH1F *Pi0XMass_FCALEcut,*Pi0XMass_NoBCAL_FCALEcut; TH1F *EtaXMass,*EtaXMass_bg; TH2F *TwoGammaMass_vs_Ebeam_acc; TH1F *ThreeGammaMass_acc,*TwoGammaMass_6g_acc; TH2F *ThreeGammaMass_vs_Ebeam_acc; TH1F *FiveGammaMass_acc,*SixGammaMass_acc,*ThreePi0Mass_acc; TH2F *FiveGammaMass_vs_Ebeam_acc,*SixGammaMass_vs_Ebeam_acc; TH1F *TwoPi0EtaMass_acc,*TwoEtaPi0Mass_acc,*ThreeEtaMass_acc; TH2F *FourGamma2d_5g_acc,*TwoGammaMass_vs_ThreeGammaMass_acc; TH1F *OmegaEta_acc,*Pi0GammaMass_7g,*TwoPhiMass; TH2F *Pi0Pi0GammaMass_vs_Pi0GammaMass_acc; TH1F *PipPimPi0Mass_no_kf,*TwoGammaMass_pippim2g_no_kf; TH2F *PipPimPi0_NoBeam_vs_E,*FourGamma2d_Pi02g; TH1F *TwoGammaMass_against_pi0,*PhiPipPimMass; TH1F *EtaCosTheta_phi; TH2F *FourGamma2d_no_kf,*PipPimPi0_vs_PipPim1Gamma; TH1F *TwoGammaMass_4g_no_kf,*Pi02gMass_no_kf; TH1F *EtaPrimeOmegaMass; TH1F *TwoGammaMass_PipPim6g_kf,*ThreePi0Mass_PipPim6g_kf,*TwoPi0EtaMass_PipPim6g_kf; TH1F *TwoEtaPi0Mass_PipPim6g_kf,*ThreeEtaMass_PipPim6g_kf; TH1F *PipPimMass_2pip2pim_no_kf; TH2F *PipPimMass_vs_PipPimMass_no_kf,*PipPimMass_vs_PipPimMass_no_kf_cut; TH1F *MissingMSqPPipPim,*MissingEMeasuredPPipPim; TH2F *MissingPzMeasuredPPipPim,*MissingPyMeasuredPPipPim,*MissingPxMeasuredPPipPim; TH2F *PPi0EtaDalitz,*TwoPipTwoPimDalitz; TH1F *OmegaPi0Mass_epem,*PPi0Mass,*PipPimMass_2pip2pim; TH2F *PipPimPi0Mass_vs_Ebeam_acc,*Phi3PiCosThetaVsPhiH,*Omega3PiCosThetaVsPhiH; TH1F *F1Pi0Mass_8g,*F1Pi0Mass_8g_cut,*EtaPrimePi0_8g_cut; TH1F *Eta2Pi0_against_pi0_cut; TH1F *Omega2Pi0_pippim6g; TH1F *PipPimPi0_6g,*Eta2Pi0_pippim6g,*PipPimEta_with_2pi0; TH2F *OmegaPi0Pi0Dalitz,*PipPimMass_vs_p_v2,*KShortKShort_vs_E,*KShortKShort_vs_t; TH1F *Pi0PipPim_with_Eta,*Pi0PipPim_with_Pi0; TH1F *PipPimPi0_with_2eta,*PipPimEta_with_etapi0; TH2F *PipPimMass_vs_PipPimMass_pi0,*PipPimMass_vs_PipPimMass_eta; TH1F *KShortKShortPi0,*ThreePipThreePimMass; TH1F *PipPimEtaMass_4pig,*PipPimEtaPrimeMass_4pig,*EShower_2pip2pim; TH1F *KShortKShort_1g,*UnusedEnergy; TH2F *PipPimMass_vs_Ebeam_tcut; TH1F *PipPim_CL,*PipPim_chisq; TH1F *NeutralBeta,*KinFitCL,*KinFitChiSq,*KinFitChiSqMissingProton; TH1F *KinFitCLMissingProton,*MissingMass_off_2gamma_CLcut; TH1F *TwoGamma_no_proton_CLcut,*TwoGammaMass_kf_no_proton; // RooRealVar *RooMass,*RooE,*RooT,*RooWeight; // RooDataSet *RooData; TH1F *OmegaEtaSBv2; TH1F *EtaPi0Mass_p2pi_kf; TH1F *PipPim4GammaMass_kf,*PipPim2Pi0Mass,*PipPim2Pi0Mass_RhoMinusCut; TH1F *PipPim2Pi0Mass_RhoPlusCut,*OmegaEta_PipPim; TH2F *OmegaEtaSB_PipPim; TH2F *OmegaEtaSBv2_PipPim,*Pi0CosTheta_EtaPrimePi0_NoOverlaps; TH2F *FourGamma2d_ppippim_NoOverlaps,*PipPi0Mass_vs_PimPi0Mass_NoOverlaps; TH1F *OmegaPi0Mass_kf; TH2F *PipPimMassVsTwoPi0Mass_NoOverlaps; TH1F *EtaPipPim_with_Eta; TH1F *EtaEtaPrimeMass,*EtaPipPim_with_Pi0; TH1F *TwoEta_PipPim; TH1F *EtaPrimePi0_PipPim,*F1Pi0Mass; TH2F *EtaPrimePi0SB_PipPim,*TwoGamma_vs_PipPimTwoGamma_NoOverlaps; TH2F *SigmaPlusPi0_SigmaPlusK0_Dalitz,*SigmaPlusEta_SigmaPlusK0_Dalitz; TH2F *PipPim_vs_NPip_with_2pi0,*PipPim_vs_NPip_with_etapi0; TH1F *FourGamma_NoMesonPairs_InsertCut,*TwoPi0_1combo,*TwoEta_1combo,*Pi0Eta_1combo; TH1F *FourGammaMass_InsertCut; TH1F *TOFdE,*FCALTOFDiff; TH2F *FCALTOFtDiff; TH2F *FCALTOFDiff2D,*FCALTOFDiff2D_dt2,*TOFdE2D; TH1F *FourGammaMass_NoVeto; TH1F *PPimDOCA_3pip3pim; TH2F *PipPimPi0_vs_PPi0; TH2F *TwoPi0Mass_t,*TwoEtaMass_t,*EtaPi0Mass_t,*EtaPrimePi0Mass_t; TH1F *NumFDCPseudos,*NumPseudos4gamma; TH1F *Pi0XMass_NoCloseShowers,*Pi0XMass_NoBCAL_NoCloseShowers; TH2F *TwoGamma_Pi02g_NoBCAL_NoCloseShowers,*Pi0g_Pi02g_NoCloseShowers; TH2F *TwoGamma_Pi02g_NoCloseShowers,*Pi0g_Pi02g_NoBCAL_NoCloseShowers; TH2F *TwoGamma_Pi02g_NoBCAL_best_NoCloseShowers; // For deuteron target TH2F *PPim_vs_PipPim_2p; TTree *threePiQTree; float treeval[100]; bool EMULATE_TRIGGER; bool FILL_ROO_DATASET; double FCAL_POS_CUT,BETA_CUT; double BCAL_R_CUT; double BCAL_THRESHOLD,FCAL_THRESHOLD,BCAL_Z_CUT,FCAL_RADIAL_CUT; double PI0_VETO_CUT; int MIN_CLUSTERS_IN_INSERT; double BCAL_DTRACK_CUT; double TRACK_CL_CUT,CL_CUT,MIN_BEAM_E,MAX_BEAM_E; bool INSERT_ONLY,FCAL_ONLY,DROP_ONE_CELL_BCAL_SHOWERS; int DEBUG_LEVEL; double UNUSED_ENERGY_CUT,PROTON_GAMMA_DT_CUT; double PROTON_CL_CUT,PIM_CL_CUT,PIP_CL_CUT,ANTIPROTON_CL_CUT; double ETA_CUT_VALUE,PI0_CUT_VALUE,ETAPRIME_CUT_VALUE; double DIRC_FCAL_TCUT,MIN_MASS_CUT; bool OUTPUT_EVENT_NUMBER; double FCAL_SUM_THRESHOLD,GAMMA_THETA_CUT; double MIN_VERTEX_Z,MAX_VERTEX_Z,MAX_VERTEX_R; const DMagneticFieldMap *bfield; DGeometry *dgeom; DVector3 mFCALCenter; int NUM_PSEUDO_CUT=8,NUM_IN_INSERT; bool VETO_GAP_EVENTS,FILL_TREES; bool TEST_PARTICLES; map>>particle_pairs; map>>bg_particle_pairs; map>bg_weights; map>>bachelor_gammas; vector> >doublets; vector> > triplets; }; #endif // _JEventProcessor_multiphoton_