BOSS 7.0.8
BESIII Offline Software System
Loading...
Searching...
No Matches
Pi0EtaToGGRecAlg.h
Go to the documentation of this file.
1#ifndef Physics_Analysis_Pi0EtaToGGRec_Alg_H
2#define Physics_Analysis_Pi0EtaToGGRec_Alg_H
3
4#include "CLHEP/Vector/LorentzVector.h"
5#include "GaudiKernel/Algorithm.h"
6#include "GaudiKernel/SmartDataPtr.h"
7
8#include "EventModel/EventHeader.h"
9#include "EvtRecEvent/EvtRecEvent.h"
10#include "EvtRecEvent/EvtRecTrack.h"
11
12// predeclaration
13class RecEmcShower;
14class EvtRecTrack;
15
16class Pi0EtaToGGRecAlg : public Algorithm {
17
18 public:
19
20 Pi0EtaToGGRecAlg(const std::string& name, ISvcLocator* pSvcLocator);
21
22 StatusCode initialize();
23 StatusCode execute();
24 StatusCode finalize();
25
26 private:
27
28
29 /// Functions
30 HepLorentzVector getP4(RecEmcShower* gTrk);
31
32 bool validPhoton(EvtRecTrack* recTrk,
33 SmartDataPtr<EvtRecEvent> recEvt,
34 SmartDataPtr<EvtRecTrackCol> recTrkCol);
35
36 bool bothInEndcap(HepLorentzVector g1_P4, HepLorentzVector g2_P4);
37
38
39 /////////////////////////
40 //// declareProperty cuts
41
42 /// reject pi0/eta candidate if both photon children are in endcap
43 bool m_rejectEndEnd;
44
45 /// pi0 cuts
46 double m_pi0minmass_cut;
47 double m_pi0maxmass_cut;
48 double m_pi0chisq_cut;
49
50 /// eta cuts
51 double m_etaminmass_cut;
52 double m_etamaxmass_cut;
53 double m_etachisq_cut;
54
55 /// Individual photon cuts
56 double m_minEnergy;
57
58 bool m_useBarrelEndcap;
59 double m_maxCosThetaBarrel;
60 double m_minCosThetaEndcap;
61 double m_maxCosThetaEndcap;
62 double m_minEndcapEnergy;
63
64 bool m_applyTimeCut;
65 double m_minTime;
66 double m_maxTime;
67 double m_deltaTime;
68
69 bool m_applyDangCut;
70 double m_minDang;
71
72 int m_maxNGoodPhoton;
73 int m_maxNPi0;
74 //int m_maxNEta;
75
76 bool m_doClear;
77};
78#endif