#ifndef DSelector_kpi0lambda_h #define DSelector_kpi0lambda_h #include #include "DSelector/DSelector.h" #include "DSelector/DHistogramActions.h" #include "DSelector/DCutActions.h" #include "TH1I.h" #include "TH2I.h" class DSelector_kpi0lambda : public DSelector { public: DSelector_kpi0lambda(TTree* locTree = NULL) : DSelector(locTree){} virtual ~DSelector_kpi0lambda(){} void Init(TTree *tree); Bool_t Process(Long64_t entry); private: void Get_ComboWrappers(void); void Finalize(void); // BEAM POLARIZATION INFORMATION UInt_t dPreviousRunNumber; bool dIsPolarizedFlag; //else is AMO bool dIsPARAFlag; //else is PERP or AMO //CREATE REACTION-SPECIFIC PARTICLE ARRAYS //Step 0 DParticleComboStep* dStep0Wrapper; DBeamParticle* dComboBeamWrapper; DChargedTrackHypothesis* dKPlusWrapper; //Step 1 DParticleComboStep* dStep1Wrapper; DKinematicData* dDecayingLambdaWrapper; DChargedTrackHypothesis* dProtonWrapper; DChargedTrackHypothesis* dPiMinusWrapper; //Step 2 DParticleComboStep* dStep2Wrapper; DKinematicData* dDecayingPi0Wrapper; DNeutralParticleHypothesis* dPhoton1Wrapper; DNeutralParticleHypothesis* dPhoton2Wrapper; ClassDef(DSelector_kpi0lambda, 0); }; void DSelector_kpi0lambda::Get_ComboWrappers(void) { //Step 0 dStep0Wrapper = dComboWrapper->Get_ParticleComboStep(0); dComboBeamWrapper = static_cast(dStep0Wrapper->Get_InitialParticle()); dKPlusWrapper = static_cast(dStep0Wrapper->Get_FinalParticle(0)); //Step 1 dStep1Wrapper = dComboWrapper->Get_ParticleComboStep(1); dDecayingLambdaWrapper = dStep1Wrapper->Get_InitialParticle(); dProtonWrapper = static_cast(dStep1Wrapper->Get_FinalParticle(0)); dPiMinusWrapper = static_cast(dStep1Wrapper->Get_FinalParticle(1)); //Step 2 dStep2Wrapper = dComboWrapper->Get_ParticleComboStep(2); dDecayingPi0Wrapper = dStep2Wrapper->Get_InitialParticle(); dPhoton1Wrapper = static_cast(dStep2Wrapper->Get_FinalParticle(0)); dPhoton2Wrapper = static_cast(dStep2Wrapper->Get_FinalParticle(1)); } #endif // DSelector_kpi0lambda_h