1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | #include "DReaction_factory_p2pi_hists.h" |
9 | |
10 | |
11 | |
12 | |
13 | jerror_t DReaction_factory_p2pi_hists::brun(JEventLoop* locEventLoop, int32_t locRunNumber) |
14 | { |
15 | vector<double> locBeamPeriodVector; |
16 | locEventLoop->GetCalib("PHOTON_BEAM/RF/beam_period", locBeamPeriodVector); |
17 | dBeamBunchPeriod = locBeamPeriodVector[0]; |
18 | |
19 | return NOERROR; |
20 | } |
21 | |
22 | |
23 | |
24 | |
25 | jerror_t DReaction_factory_p2pi_hists::evnt(JEventLoop* locEventLoop, uint64_t locEventNumber) |
26 | { |
27 | |
28 | DReactionStep* locReactionStep = NULL__null; |
29 | DReaction* locReaction = new DReaction("p2pi_pmiss"); |
| Value stored to 'locReaction' during its initialization is never read |
30 | |
31 | |
32 | |
33 | |
34 | |
35 | |
36 | |
37 | bool unused = false; |
38 | locReaction = new DReaction("p2pi_preco"); |
39 | |
40 | |
41 | locReactionStep = new DReactionStep(); |
42 | locReactionStep->Set_InitialParticleID(Gamma); |
43 | locReactionStep->Set_TargetParticleID(Proton); |
44 | locReactionStep->Add_FinalParticleID(Proton); |
45 | locReactionStep->Add_FinalParticleID(PiPlus); |
46 | locReactionStep->Add_FinalParticleID(PiMinus); |
47 | locReaction->Add_ReactionStep(locReactionStep); |
48 | dReactionStepPool.push_back(locReactionStep); |
49 | |
50 | |
51 | |
52 | |
53 | locReaction->Set_EventStoreSkims("2q+,q-"); |
54 | |
55 | |
56 | |
57 | |
58 | |
59 | locReaction->Set_NumPlusMinusRFBunches(0); |
60 | |
61 | |
62 | |
63 | |
64 | |
65 | |
66 | |
67 | |
68 | locReaction->Add_AnalysisAction(new DHistogramAction_PID(locReaction, false)); |
69 | |
70 | |
71 | locReaction->Add_AnalysisAction(new DCustomAction_p2pi_hists(locReaction, false)); |
72 | |
73 | |
74 | locReaction->Add_AnalysisAction(new DHistogramAction_MissingMassSquared(locReaction, false, 1000, -0.1, 0.1)); |
75 | locReaction->Add_AnalysisAction(new DCutAction_MissingMassSquared(locReaction, false, -0.01, 0.005)); |
76 | |
77 | |
78 | deque<Particle_t> locRhoPIDs; locRhoPIDs.push_back(PiPlus); locRhoPIDs.push_back(PiMinus); |
79 | locReaction->Add_AnalysisAction(new DHistogramAction_InvariantMass(locReaction, 0, locRhoPIDs, false, 900, 0.3, 1.2, "Rho")); |
80 | |
81 | if(unused) |
82 | { |
83 | |
84 | locReaction->Add_AnalysisAction(new DCustomAction_p2pi_cuts(locReaction, false)); |
85 | |
86 | |
87 | locReaction->Add_AnalysisAction(new DCustomAction_p2pi_unusedHists(locReaction, false, "Unused")); |
88 | } |
89 | |
90 | |
91 | locReaction->Add_AnalysisAction(new DHistogramAction_ParticleComboKinematics(locReaction, false)); |
92 | |
93 | _data.push_back(locReaction); |
94 | |
95 | |
96 | |
97 | |
98 | |
99 | locReaction = new DReaction("p2pi_preco_kinfit"); |
100 | |
101 | |
102 | locReaction->Add_ReactionStep(locReactionStep); |
103 | |
104 | |
105 | |
106 | |
107 | locReaction->Set_EventStoreSkims("2q+,q-"); |
108 | |
109 | |
110 | locReaction->Set_KinFitType(d_P4AndVertexFit); |
111 | |
112 | |
113 | locReaction->Set_NumPlusMinusRFBunches(0); |
114 | |
115 | |
116 | |
117 | |
118 | |
119 | |
120 | |
121 | |
122 | locReaction->Add_AnalysisAction(new DHistogramAction_PID(locReaction, false)); |
123 | |
124 | |
125 | locReaction->Add_AnalysisAction(new DHistogramAction_MissingMassSquared(locReaction, false, 1000, -0.1, 0.1)); |
126 | locReaction->Add_AnalysisAction(new DHistogramAction_InvariantMass(locReaction, 0, locRhoPIDs, false, 900, 0.3, 1.2, "Rho")); |
127 | locReaction->Add_AnalysisAction(new DHistogramAction_InvariantMass(locReaction, 0, locRhoPIDs, true, 900, 0.3, 1.2, "Rho_KinFit")); |
128 | |
129 | |
130 | locReaction->Add_AnalysisAction(new DHistogramAction_KinFitResults(locReaction, 0.05, true)); |
131 | locReaction->Add_AnalysisAction(new DCutAction_KinFitFOM(locReaction, 5.73303E-7)); |
132 | |
133 | |
134 | locReaction->Add_AnalysisAction(new DHistogramAction_MissingMassSquared(locReaction, false, 1000, -0.1, 0.1, "PostKinFitCut")); |
135 | locReaction->Add_AnalysisAction(new DHistogramAction_InvariantMass(locReaction, 0, locRhoPIDs, false, 900, 0.3, 1.2, "Rho_PostKinFitCut")); |
136 | locReaction->Add_AnalysisAction(new DHistogramAction_InvariantMass(locReaction, 0, locRhoPIDs, true, 900, 0.3, 1.2, "Rho_KinFit_PostKinFitCut")); |
137 | |
138 | |
139 | locReaction->Add_AnalysisAction(new DHistogramAction_ParticleComboKinematics(locReaction, false)); |
140 | |
141 | _data.push_back(locReaction); |
142 | |
143 | return NOERROR; |
144 | } |
145 | |
146 | |
147 | |
148 | |
149 | jerror_t DReaction_factory_p2pi_hists::fini(void) |
150 | { |
151 | for(size_t loc_i = 0; loc_i < dReactionStepPool.size(); ++loc_i) |
152 | delete dReactionStepPool[loc_i]; |
153 | return NOERROR; |
154 | } |
155 | |