BOSS 7.0.8
BESIII Offline Software System
Loading...
Searching...
No Matches
DimuPreSelect Class Reference

#include <DimuPreSelect.h>

+ Inheritance diagram for DimuPreSelect:

Public Member Functions

 DimuPreSelect (const std::string &name, ISvcLocator *pSvcLocator)
 
StatusCode initialize ()
 
StatusCode execute ()
 
StatusCode finalize ()
 
 DimuPreSelect ()
 
 ~DimuPreSelect ()
 
void BookNtuple (NTuple::Tuple *&tuple)
 
int IsDimu ()
 
void Print ()
 

Detailed Description

Constructor & Destructor Documentation

◆ DimuPreSelect() [1/2]

DimuPreSelect::DimuPreSelect ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Definition at line 49 of file DimuPreSelect/DimuPreSelect-00-00-09/src/DimuPreSelect.cxx.

49 :
50 Algorithm(name, pSvcLocator)
51{
52 declareProperty("CmsEnergy", m_ecms = 3.686);
53 declareProperty("Vr0Cut", m_vr0cut = 1.0);
54 declareProperty("Vz0Cut", m_vz0cut = 5.0);
55 declareProperty("PUpCut", m_pcut_up = 2.0);
56 declareProperty("PDownCut", m_pcut_down = 0.5);
57 declareProperty("PSymCut", m_psymcut = 0.5);
58 declareProperty("TCut", m_tcut = 4);
59 declareProperty("EUpCut", m_ecut_up = 1.0);
60 declareProperty("EDownCut", m_ecut_down = 0.1);
61 declareProperty("DThetaCut", m_dthetacut = 0.05);
62 declareProperty("DPhiCut", m_dphicut = 0.4);
63 declareProperty("PartSelect", m_partselect= 0);
64 declareProperty("MuDigiCut", m_mudigicut = 6);
65 declareProperty("MuTrkCut", m_mutrkcut = 1);
66 declareProperty("UseMDC", m_useFlag[0]= 1);
67 declareProperty("UseTOF", m_useFlag[1]= 1);
68 declareProperty("UseEMC", m_useFlag[2]= 1);
69 declareProperty("UseMUC", m_useFlag[3]= 1);
70 declareProperty ("Output", m_output = false);
71
72 declareProperty ("SelectDimu",m_selectFlag= true);
73 //Declare the properties
74}

◆ DimuPreSelect() [2/2]

DimuPreSelect::DimuPreSelect ( )

Definition at line 24 of file EventPreSelect/EventPreSelect-00-00-21/src/DimuPreSelect.cxx.

25{
26 m_output=false;
27 m_totevent = m_currun = m_curevent = 0;
28 for(int i = 0; i < 20; i++) m_cutpass[i] = 0;
29 m_subpass[0] = m_subpass[1] = m_subpass[2] = m_subpass[3] = 0;
30 m_totalpass = 0;
31
32 m_propMgr.declareProperty("CmsEnergy", m_ecms = 3.686);
33 m_propMgr.declareProperty("Vr0Cut", m_vr0cut = 1.0);
34 m_propMgr.declareProperty("Vz0Cut", m_vz0cut = 5.0);
35 m_propMgr.declareProperty("PUpCut", m_pcut_up = 2.0);
36 m_propMgr.declareProperty("PDownCut", m_pcut_down = 0.5);
37 m_propMgr.declareProperty("PSymCut", m_psymcut = 0.5);
38 m_propMgr.declareProperty("TCut", m_tcut = 4);
39 m_propMgr.declareProperty("EUpCut", m_ecut_up = 1.0);
40 m_propMgr.declareProperty("EDownCut", m_ecut_down = 0.1);
41 m_propMgr.declareProperty("DThetaCut", m_dthetacut = 0.05);
42 m_propMgr.declareProperty("DPhiCut", m_dphicut = 0.4);
43 m_propMgr.declareProperty("PartSelect", m_partselect= 0);
44 m_propMgr.declareProperty("MuDigiCut", m_mudigicut = 6);
45 m_propMgr.declareProperty("MuTrkCut", m_mutrkcut = 1);
46
47 IJobOptionsSvc* jobSvc;
48 Gaudi::svcLocator()->service("JobOptionsSvc", jobSvc);
49 jobSvc->setMyProperties("DimuSelect", &m_propMgr);
50}

◆ ~DimuPreSelect()

DimuPreSelect::~DimuPreSelect ( )
inline

Member Function Documentation

◆ BookNtuple()

void DimuPreSelect::BookNtuple ( NTuple::Tuple *&  tuple)

Definition at line 52 of file EventPreSelect/EventPreSelect-00-00-21/src/DimuPreSelect.cxx.

53{
54 m_output=true;
55 m_passtuple=tuple;
56 if(!m_passtuple) {
57 std::cerr << "invalid ntuple in DimuPreSelect!" << std::endl;
58 } else {
59 m_passtuple->addItem ("run", m_run);
60 m_passtuple->addItem ("event", m_event);
61 m_passtuple->addItem ("part", m_part);
62 m_passtuple->addItem ("c1", m_c1);
63 m_passtuple->addItem ("c2", m_c2);
64 m_passtuple->addItem ("r1", m_r1);
65 m_passtuple->addItem ("r2", m_r2);
66 m_passtuple->addItem ("z1", m_z1);
67 m_passtuple->addItem ("z2", m_z2);
68 m_passtuple->addItem ("p1", m_p1);
69 m_passtuple->addItem ("p2", m_p2);
70 m_passtuple->addItem ("t1", m_t1);
71 m_passtuple->addItem ("t2", m_t2);
72 m_passtuple->addItem ("e1", m_e1);
73 m_passtuple->addItem ("e2", m_e2);
74 m_passtuple->addItem ("theta1",m_theta1);
75 m_passtuple->addItem ("theta2",m_theta2);
76 m_passtuple->addItem ("phi1", m_phi1);
77 m_passtuple->addItem ("phi2", m_phi2);
78 m_passtuple->addItem ("mudigi",m_mudigi);
79 m_passtuple->addItem ("mutrk", m_mutrk);
80
81 m_passtuple->addItem ("zeroC", m_zeroCFlag);
82 m_passtuple->addItem ("vtRZ", m_vtRZFlag);
83 m_passtuple->addItem ("pLim", m_pLimFlag);
84 m_passtuple->addItem ("pSym", m_pSymFlag);
85 m_passtuple->addItem ("tLim", m_tLimFlag);
86 m_passtuple->addItem ("eLim", m_eLimFlag);
87 m_passtuple->addItem ("eBB", m_eBBFlag);
88 m_passtuple->addItem ("partSlt",m_partFlag);
89 m_passtuple->addItem ("muDigiN",m_mudigiFlag);
90 m_passtuple->addItem ("muTrkN", m_mutrkFlag);
91
92 m_passtuple->addItem ("mdc", m_mdcFlag);
93 m_passtuple->addItem ("tof", m_tofFlag);
94 m_passtuple->addItem ("emc", m_emcFlag);
95 m_passtuple->addItem ("muc", m_mucFlag);
96 }
97}

Referenced by EventPreSelect::initialize().

◆ execute()

StatusCode DimuPreSelect::execute ( )

Definition at line 149 of file DimuPreSelect/DimuPreSelect-00-00-09/src/DimuPreSelect.cxx.

150{
151 if(m_selectFlag==false) return( StatusCode::SUCCESS );
152
153 MsgStream log(msgSvc(),name());
154 log<<MSG::INFO<<"in execute()"<<endreq;
155
156 m_totevent ++;
157 if(m_totevent%1000==0) std::cout << m_totevent << "\tdone!" <<std::endl;
158
159 setFilterPassed(false);
160 m_mdcPass = m_tofPass = m_emcPass = m_mucPass = false;
161 for(int i=0; i<4; i++) m_passFlag[i] = false;
162
163 SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
164 m_currun = eventHeader->runNumber();
165 m_curevent = eventHeader->eventNumber();
166 if( m_output ) { m_run = m_currun; m_event = m_curevent; }
167
168 // Select by MDC Info
169 SmartDataPtr<RecMdcTrackCol> mdcTrackCol(eventSvc(),"/Event/Recon/RecMdcTrackCol");
170 if(!mdcTrackCol) {
171 log << MSG::FATAL << "Could not find RecMdcTrackCol!" << endreq;
172 return( StatusCode::FAILURE);
173 }
174 log << MSG::INFO << "MDC tracks:\t" << mdcTrackCol->size() << endreq;
175
176 if( mdcTrackCol->size() != 2 ) return ( StatusCode::SUCCESS );
177
178 m_cutpass[0] += 1;
179
180 double c1, c2, r1, r2, z1, z2, p1, p2;
181 c1 = c2 = r1 = r2 = z1 = z2 = p1 = p2 = 0.;
182
183 c1 = (*mdcTrackCol)[0]->charge(); c2 = (*mdcTrackCol)[1]->charge();
184 r1 = (*mdcTrackCol)[0]->r(); r2 = (*mdcTrackCol)[1]->r();
185 z1 = (*mdcTrackCol)[0]->z(); z2 = (*mdcTrackCol)[1]->z();
186 p1 = (*mdcTrackCol)[0]->p(); p2 = (*mdcTrackCol)[1]->p();
187
188 if( m_output ) { m_c1 = c1; m_c2 = c2; m_r1 = r1; m_r2 = r2; m_z1 = z1; m_z2 = z2; m_p1 = p1; m_p2 = p2; }
189
190 bool mdcflag1 = c1 + c2 == 0;
191 bool mdcflag2 = fabs(r1)<=m_vr0cut && fabs(z1)<m_vz0cut;
192 bool mdcflag3 = fabs(r2)<=m_vr0cut && fabs(z2)<m_vz0cut;
193 bool mdcflag4 = p1<m_pcut_up && p2<m_pcut_up;
194 //bool mdcflag4 = p1<2 && p1>1;
195 //bool mdcflag5 = p2<2 && p2>1;
196 bool mdcflag5 = fabs( p1-p2 )/( p1+p2 ) < m_psymcut;
197
198 log << MSG::INFO << "r1:\t" << r1 << "\tz1:" << z1 << endreq;
199 log << MSG::INFO << "r2:\t" << r2 << "\tz2:" << z2 << endreq;
200 log << MSG::INFO << "p1:\t" << p1 << "\tp2:" << p2 << endreq;
201
202 if( mdcflag1 ) { m_cutpass[1] += 1; if(m_output) m_zeroCFlag = 1; }
203 if( mdcflag2 && mdcflag3 ) { m_cutpass[2] += 1; if(m_output) m_vtRZFlag = 1; }
204 if( mdcflag4 ) { m_cutpass[3] += 1; if(m_output) m_pLimFlag = 1; }
205 if( mdcflag5 ) { m_cutpass[4] += 1; if(m_output) m_pSymFlag = 1; }
206 if( mdcflag1 && mdcflag2 && mdcflag3 && mdcflag4 && mdcflag5 ) { m_mdcPass = true; m_subpass[0] += 1; }
207 if( !m_useFlag[0] ) m_passFlag[0] = true;
208 else m_passFlag[0] = m_mdcPass;
209 log << MSG::INFO << "MDC selection done!" << endreq;
210
211 // Select by TOF Info
212 SmartDataPtr<RecTofTrackCol> tofTrackCol(eventSvc(),"/Event/Recon/RecTofTrackCol");
213 if(!tofTrackCol) {
214 log << MSG::FATAL << "Could not find RecTofTrackCol!" << endreq;
215 return( StatusCode::FAILURE);
216 }
217 log << MSG::INFO << "TOF tracks:\t" << tofTrackCol->size() << endreq;
218
219 double t1, t2;
220 t1 = 0., t2 = 1000;
221 // if(tofTrackCol->size() < 8 && tofTrackCol->size() > 20)
222 if(tofTrackCol->size() > 7 && tofTrackCol->size() < 21)
223 {
224 int goodtof = 0;
225 for(int itof = 0; itof < tofTrackCol->size(); itof++)
226 {
227 TofHitStatus *status = new TofHitStatus;
228 status->setStatus((*tofTrackCol)[itof]->status());
229
230 if( !(status->is_cluster()) ) { delete status; continue; }
231 if(goodtof==0) t1 = (*tofTrackCol)[itof]->tof();
232 if(goodtof==1) t2 = (*tofTrackCol)[itof]->tof();
233
234 goodtof++;
235 delete status;
236 }
237 }
238
239 if( m_output ) { m_t1 = t1; m_t2 = t2; }
240
241 bool tofflag1 = fabs( t1-t2 ) < m_tcut; // ns
242 log << MSG::INFO << "t1:\t" << t1 << "\tt2:" << t2 << "dt:\t" << fabs(t1-t2) << endreq;
243 if( tofflag1 ) { m_cutpass[5] += 1; if(m_output) m_tLimFlag = 1; m_tofPass = true; m_subpass[1] += 1; }
244 if( !m_useFlag[1] ) m_passFlag[1] = true;
245 else m_passFlag[1] = m_tofPass;
246 log << MSG::INFO << "TOF selection done!" << endreq;
247
248 // Select by EMC Info
249 SmartDataPtr<RecEmcShowerCol> emcShowerCol(eventSvc(),"/Event/Recon/RecEmcShowerCol");
250 if (!emcShowerCol) {
251 log << MSG::FATAL << "Could not find RecEmcShowerCol!" << endreq;
252 return( StatusCode::FAILURE);
253 }
254 log << MSG::INFO << "EMC showers:\t" << emcShowerCol->size() << endreq;
255
256 if( emcShowerCol->size() < 2 ) return ( StatusCode::SUCCESS );
257
258 double e1, e2, theta1, theta2, phi1, phi2;
259 int part;
260
261 e1 = (*emcShowerCol)[0]->energy(); e2 = (*emcShowerCol)[1]->energy();
262 theta1 = (*emcShowerCol)[0]->theta(); theta2 = (*emcShowerCol)[1]->theta();
263 phi1 = (*emcShowerCol)[0]->phi(); phi2 = (*emcShowerCol)[1]->phi();
264 part = (*emcShowerCol)[0]->module();
265
266 if( m_output ) { m_e1 = e1; m_e2 = e2; m_theta1 = theta1; m_theta2 = theta2; m_phi1 = phi1; m_phi2 = phi2;
267 m_part = part;
268 }
269
270 bool emcFlag1 = e1 < m_ecut_up && e1 > m_ecut_down;
271 bool emcFlag2 = e2 < m_ecut_up && e2 > m_ecut_down;
272 bool emcFlag3 = fabs(theta1 + theta2 - PI) < m_dthetacut;
273 bool emcFlag4 = fabs(phi1 - phi2) - PI < m_dphicut;
274 bool emcFlag5 = !m_partselect || (m_partselect==1&&part==1) || (m_partselect==2&&part!=1);
275
276 log << MSG::INFO << "e1:\t" << e1 << "\te2:\t" << e2 << endreq;
277 log << MSG::INFO << "theta1:\t" << theta1 << "\ttheta2:\t" << theta2 << endreq;
278 log << MSG::INFO << "phi1:\t" << phi1 << "\tphi2:\t" << phi2 << endreq;
279 log << MSG::INFO << "part:\t" << part << "\tpartFlag:\t" << emcFlag5 << endreq;
280
281 if( emcFlag1 && emcFlag2 ) { m_cutpass[6] += 1; if(m_output) m_eLimFlag = 1; }
282 if( emcFlag3 && emcFlag4 ) { m_cutpass[7] += 1; if(m_output) m_eBBFlag = 1; }
283 if( emcFlag5 ) { m_cutpass[8] += 1; if(m_output) m_partFlag = 1; }
284 if( emcFlag1 && emcFlag2 && emcFlag3 && emcFlag4 && emcFlag5 ) { m_emcPass = true; m_subpass[2] += 1; }
285 if( !m_useFlag[2] ) m_passFlag[2] = true;
286 else m_passFlag[2] = m_emcPass;
287 log << MSG::INFO << "EMC selection done!" << endreq;
288
289 // Select by MUC Info
290 SmartDataPtr<MucDigiCol> mucDigiCol(eventSvc(),"/Event/Digi/MucDigiCol");
291 if(!mucDigiCol) {
292 log << MSG::FATAL << "Could not find MucDigiCol!" << endreq;
293 return( StatusCode::FAILURE);
294 }
295 SmartDataPtr<RecMucTrackCol> mucTrackCol(eventSvc(),"/Event/Recon/RecMucTrackCol");
296 if (!mucTrackCol) {
297 log << MSG::FATAL << "Could not find RecMucTrackCol" << endreq;
298 return( StatusCode::FAILURE);
299 }
300
301 int mudigiNum, mutrkNum;
302 mudigiNum = mutrkNum = 0;
303 mudigiNum = mucDigiCol->size(); mutrkNum = mucTrackCol->size();
304 if(m_output) { m_mudigi = mudigiNum; m_mutrk = mutrkNum; }
305
306 bool mucflag1 = mudigiNum >= m_mudigicut;
307 bool mucflag2 = mutrkNum >= m_mutrkcut;
308
309 log << MSG::INFO << "MUC digi:\t" << mudigiNum << "\tMUC track:\t" << mutrkNum << endreq;
310
311 if( mucflag1 ) { m_cutpass[9] += 1; if(m_output) m_mudigiFlag = 1; }
312 if( mucflag2 ) { m_cutpass[10] += 1; if(m_output) m_mutrkFlag = 1; }
313 if( mucflag1 && mucflag2 ) { m_mucPass = true; m_subpass[3] += 1; }
314 if( !m_useFlag[3] ) m_passFlag[3] = true;
315 else m_passFlag[3] = m_mucPass;
316 log << MSG::INFO << "MUC selection done!" << endreq;
317
318
319 // All selections
320 if( m_passFlag[0] && m_passFlag[1] && m_passFlag[2] && m_passFlag[3] )
321 {
322 m_totalpass += 1;
323 setFilterPassed(true);
324 }
325 log << MSG::INFO << "Set filter passed!" << endreq;
326
327 if(m_output) { m_mdcFlag = m_mdcPass; m_tofFlag = m_tofPass; m_emcFlag = m_emcPass; m_mucFlag = m_mucPass; }
328
329 if( m_output ) m_passtuple->write();
330
331 return StatusCode::SUCCESS;
332}
Double_t phi2
Double_t phi1
Double_t e1
Double_t e2
const double PI
Definition: PipiJpsi.cxx:55
IMessageSvc * msgSvc()
@ theta2
Definition: TrkKalDeriv.h:24
@ theta1
Definition: TrkKalDeriv.h:24
bool is_cluster() const
Definition: TofHitStatus.h:25
void setStatus(unsigned int status)

◆ finalize()

StatusCode DimuPreSelect::finalize ( )

Definition at line 334 of file DimuPreSelect/DimuPreSelect-00-00-09/src/DimuPreSelect.cxx.

335{
336 if(m_selectFlag==false) return StatusCode::SUCCESS;
337
338 MsgStream log(msgSvc(), name());
339 log << MSG::INFO << "in finalize()" << endmsg;
340 const string str[3] = {"all","barrel","endcap"};
341
342 cout << "pass 0 - 2 MDC tracks : " << m_cutpass[0] << endl;
343 cout << "pass 1 - 0 charge : " << m_cutpass[1] << endl;
344 cout << "pass 2 - |r|<"<<m_vr0cut<<" && |z|<"<<m_vz0cut<<" : " << m_cutpass[2] << endl;
345 cout << "pass 3 - p < "<<m_pcut_up<<" GeV/c : " << m_cutpass[3] << endl;
346 cout << "pass 4 - p_sym < "<<m_psymcut<<" : " << m_cutpass[4] << endl;
347 cout << "pass 5 - |t1-t1| < "<<m_tcut<<" ns : " << m_cutpass[5] << endl;
348 cout << "pass 6 - "<<m_ecut_down<<" < e < "<<m_ecut_up<<" : " << m_cutpass[6] << endl;
349 cout << "pass 7 - |dth|<"<<m_dthetacut<<" && |dphi|<"<<m_dphicut<<": " << m_cutpass[7] << endl;
350 cout << "pass 8 - "<< str[(int)m_partselect] <<" part is selected : " << m_cutpass[8] << endl;
351 cout << "pass 9 - mudigi >= "<<m_mudigicut<<" : " << m_cutpass[9] << endl;
352 cout << "pass 10- mutrk >= "<<m_mutrkcut<<" : " << m_cutpass[10]<< endl;
353
354 cout << "pass MDC : " << m_subpass[0] << "\tUsed: " << (m_useFlag[0]?"Yes":"No") << endl;
355 cout << "pass TOF : " << m_subpass[1] << "\tUsed: " << (m_useFlag[1]?"Yes":"No") << endl;
356 cout << "pass EMC : " << m_subpass[2] << "\tUsed: " << (m_useFlag[2]?"Yes":"No") << endl;
357 cout << "pass MUC : " << m_subpass[3] << "\tUsed: " << (m_useFlag[3]?"Yes":"No") << endl;
358
359
360 cout<<"Total event: " << m_totevent <<endl;
361 cout<<"Dimu event: " << m_totalpass <<endl;
362
363 return StatusCode::SUCCESS;
364}

◆ initialize()

StatusCode DimuPreSelect::initialize ( )

Definition at line 77 of file DimuPreSelect/DimuPreSelect-00-00-09/src/DimuPreSelect.cxx.

78{
79 MsgStream log(msgSvc(), name());
80 log << MSG::INFO << "in initialize()" << endmsg;
81
82 m_totevent = m_currun = m_curevent = 0;
83 for(int i = 0; i < 20; i++) m_cutpass[i] = 0;
84 m_subpass[0] = m_subpass[1] = m_subpass[2] = m_subpass[3] = 0;
85 m_totalpass = 0;
86
87 if(m_output)
88 {
89 StatusCode status;
90 NTuplePtr nt(ntupleSvc(), "FILE1/dimu");
91 if ( nt ) m_passtuple = nt;
92 else {
93 m_passtuple = ntupleSvc()->book ("FILE1/dimu", CLID_ColumnWiseTuple, "DimuPreSelect N-Tuple example");
94 if ( m_passtuple )
95 {
96 status =m_passtuple->addItem ("run", m_run);
97 status =m_passtuple->addItem ("event", m_event);
98 status =m_passtuple->addItem ("part", m_part);
99 status =m_passtuple->addItem ("c1", m_c1);
100 status =m_passtuple->addItem ("c2", m_c2);
101 status =m_passtuple->addItem ("r1", m_r1);
102 status =m_passtuple->addItem ("r2", m_r2);
103 status =m_passtuple->addItem ("z1", m_z1);
104 status =m_passtuple->addItem ("z2", m_z2);
105 status =m_passtuple->addItem ("p1", m_p1);
106 status =m_passtuple->addItem ("p2", m_p2);
107 status =m_passtuple->addItem ("t1", m_t1);
108 status =m_passtuple->addItem ("t2", m_t2);
109 status =m_passtuple->addItem ("e1", m_e1);
110 status =m_passtuple->addItem ("e2", m_e2);
111 status =m_passtuple->addItem ("theta1",m_theta1);
112 status =m_passtuple->addItem ("theta2",m_theta2);
113 status =m_passtuple->addItem ("phi1", m_phi1);
114 status =m_passtuple->addItem ("phi2", m_phi2);
115 status =m_passtuple->addItem ("mudigi",m_mudigi);
116 status =m_passtuple->addItem ("mutrk", m_mutrk);
117
118 status =m_passtuple->addItem ("zeroC", m_zeroCFlag);
119 status =m_passtuple->addItem ("vtRZ", m_vtRZFlag);
120 status =m_passtuple->addItem ("pLim", m_pLimFlag);
121 status =m_passtuple->addItem ("pSym", m_pSymFlag);
122 status =m_passtuple->addItem ("tLim", m_tLimFlag);
123 status =m_passtuple->addItem ("eLim", m_eLimFlag);
124 status =m_passtuple->addItem ("eBB", m_eBBFlag);
125 status =m_passtuple->addItem ("partSlt",m_partFlag);
126 status =m_passtuple->addItem ("muDigiN",m_mudigiFlag);
127 status =m_passtuple->addItem ("muTrkN", m_mutrkFlag);
128
129 status =m_passtuple->addItem ("mdc", m_mdcFlag);
130 status =m_passtuple->addItem ("tof", m_tofFlag);
131 status =m_passtuple->addItem ("emc", m_emcFlag);
132 status =m_passtuple->addItem ("muc", m_mucFlag);
133
134 }
135 else
136 {
137 log << MSG::ERROR << "Cannot book N-tuple:"<<long(m_passtuple)<<endmsg;
138 return StatusCode::FAILURE;
139 }
140 }
141 }
142
143 log << MSG::INFO << "Initialize done!" << endmsg;
144
145 return StatusCode::SUCCESS;
146}
INTupleSvc * ntupleSvc()

◆ IsDimu()

int DimuPreSelect::IsDimu ( )

Definition at line 99 of file EventPreSelect/EventPreSelect-00-00-21/src/DimuPreSelect.cxx.

100{
101 IMessageSvc* msgSvc;
102 Gaudi::svcLocator()->service("MessageSvc", msgSvc);
103 MsgStream log(msgSvc, "DimuPreSelect");
104
105 IDataProviderSvc* eventSvc;
106 Gaudi::svcLocator()->service("EventDataSvc", eventSvc);
107
108 m_mdcPass = false;
109 m_tofPass = false;
110 m_emcPass = false;
111 m_mucPass = false;
112
113 // Select by MDC Info
114 SmartDataPtr<RecMdcTrackCol> mdcTrackCol(eventSvc,"/Event/Recon/RecMdcTrackCol");
115 if(!mdcTrackCol) {
116 log << MSG::ERROR << "Could not find RecMdcTrackCol!" << endreq;
117 return false;
118 }
119 log << MSG::INFO << "MDC tracks:\t" << mdcTrackCol->size() << endreq;
120
121 if( mdcTrackCol->size() != 2 ) return 0;
122
123 m_cutpass[0] += 1;
124
125 double c1, c2, r1, r2, z1, z2, p1, p2;
126 c1 = c2 = r1 = r2 = z1 = z2 = p1 = p2 = 0.;
127
128 c1 = (*mdcTrackCol)[0]->charge(); c2 = (*mdcTrackCol)[1]->charge();
129 r1 = (*mdcTrackCol)[0]->r(); r2 = (*mdcTrackCol)[1]->r();
130 z1 = (*mdcTrackCol)[0]->z(); z2 = (*mdcTrackCol)[1]->z();
131 p1 = (*mdcTrackCol)[0]->p(); p2 = (*mdcTrackCol)[1]->p();
132
133
134 bool mdcflag1 = c1 + c2 == 0;
135 bool mdcflag2 = fabs(r1)<=m_vr0cut && fabs(z1)<m_vz0cut;
136 bool mdcflag3 = fabs(r2)<=m_vr0cut && fabs(z2)<m_vz0cut;
137 bool mdcflag4 = p1<m_pcut_up && p2<m_pcut_up;
138 //bool mdcflag4 = p1<2 && p1>1;
139 //bool mdcflag5 = p2<2 && p2>1;
140 bool mdcflag5 = fabs( p1-p2 )/( p1+p2 ) < m_psymcut;
141
142 log << MSG::INFO << "r1:\t" << r1 << "\tz1:" << z1 << endreq;
143 log << MSG::INFO << "r2:\t" << r2 << "\tz2:" << z2 << endreq;
144 log << MSG::INFO << "p1:\t" << p1 << "\tp2:" << p2 << endreq;
145
146 double zeroCFlag,vtRZFlag,pLimFlag,pSymFlag;
147 if( mdcflag1 ) { m_cutpass[1] += 1; zeroCFlag = 1; }
148 if( mdcflag2 && mdcflag3 ) { m_cutpass[2] += 1; vtRZFlag = 1; }
149 if( mdcflag4 ) { m_cutpass[3] += 1; pLimFlag = 1; }
150 if( mdcflag5 ) { m_cutpass[4] += 1; pSymFlag = 1; }
151
152 if( mdcflag1 && mdcflag2 && mdcflag3 && mdcflag4 && mdcflag5 ) { m_mdcPass = true; m_subpass[0] += 1; }
153 log << MSG::INFO << "MDC selection done!" << endreq;
154
155 // Select by TOF Info
156 SmartDataPtr<RecTofTrackCol> tofTrackCol(eventSvc,"/Event/Recon/RecTofTrackCol");
157 if(!tofTrackCol) {
158 log << MSG::ERROR << "Could not find RecTofTrackCol!" << endreq;
159 return 0;
160 }
161 log << MSG::INFO << "TOF tracks:\t" << tofTrackCol->size() << endreq;
162
163 double t1, t2;
164 t1 = 0., t2 = 1000;
165 // if(tofTrackCol->size() < 8 && tofTrackCol->size() > 20)
166 if(tofTrackCol->size() > 7 && tofTrackCol->size() < 21)
167 {
168 int goodtof = 0;
169 for(int itof = 0; itof < tofTrackCol->size(); itof++)
170 {
171 TofHitStatus *status = new TofHitStatus;
172 status->setStatus((*tofTrackCol)[itof]->status());
173
174 if( !(status->is_cluster()) ) { delete status; continue; }
175 if(goodtof==0) t1 = (*tofTrackCol)[itof]->tof();
176 if(goodtof==1) t2 = (*tofTrackCol)[itof]->tof();
177
178 goodtof++;
179 delete status;
180 }
181 }
182
183 double tLimFlag;
184 bool tofflag1 = fabs( t1-t2 ) < m_tcut; // ns
185 log << MSG::INFO << "t1:\t" << t1 << "\tt2:" << t2 << "dt:\t" << fabs(t1-t2) << endreq;
186 if( tofflag1 ) { m_cutpass[5] += 1; tLimFlag = 1; m_tofPass = true; m_subpass[1] += 1; }
187 log << MSG::INFO << "TOF selection done!" << endreq;
188
189 // Select by EMC Info
190 SmartDataPtr<RecEmcShowerCol> emcShowerCol(eventSvc,"/Event/Recon/RecEmcShowerCol");
191 if (!emcShowerCol) {
192 log << MSG::ERROR << "Could not find RecEmcShowerCol!" << endreq;
193 return 0;
194 }
195 log << MSG::INFO << "EMC showers:\t" << emcShowerCol->size() << endreq;
196
197 if( emcShowerCol->size() < 2 ) return 0;
198
199 double e1, e2, theta1, theta2, phi1, phi2;
200 int part;
201
202 e1 = (*emcShowerCol)[0]->energy(); e2 = (*emcShowerCol)[1]->energy();
203 theta1 = (*emcShowerCol)[0]->theta(); theta2 = (*emcShowerCol)[1]->theta();
204 phi1 = (*emcShowerCol)[0]->phi(); phi2 = (*emcShowerCol)[1]->phi();
205 part = (*emcShowerCol)[0]->module();
206
207 bool emcFlag1 = e1 < m_ecut_up && e1 > m_ecut_down;
208 bool emcFlag2 = e2 < m_ecut_up && e2 > m_ecut_down;
209 bool emcFlag3 = fabs(theta1 + theta2 - CLHEP::pi) < m_dthetacut;
210 bool emcFlag4 = fabs(phi1 - phi2) - CLHEP::pi < m_dphicut;
211 bool emcFlag5 = !m_partselect || (m_partselect==1&&part==1) || (m_partselect==2&&part!=1);
212
213 log << MSG::INFO << "e1:\t" << e1 << "\te2:\t" << e2 << endreq;
214 log << MSG::INFO << "theta1:\t" << theta1 << "\ttheta2:\t" << theta2 << endreq;
215 log << MSG::INFO << "phi1:\t" << phi1 << "\tphi2:\t" << phi2 << endreq;
216 log << MSG::INFO << "part:\t" << part << "\tpartFlag:\t" << emcFlag5 << endreq;
217
218 double eLimFlag,eBBFlag,partFlag;
219 if( emcFlag1 && emcFlag2 ) { m_cutpass[6] += 1; eLimFlag = 1; }
220 if( emcFlag3 && emcFlag4 ) { m_cutpass[7] += 1; eBBFlag = 1; }
221 if( emcFlag5 ) { m_cutpass[8] += 1; partFlag = 1; }
222 if( emcFlag1 && emcFlag2 && emcFlag3 && emcFlag4 && emcFlag5 ) { m_emcPass = true; m_subpass[2] += 1; }
223 log << MSG::INFO << "EMC selection done!" << endreq;
224
225 // Select by MUC Info
226 SmartDataPtr<MucDigiCol> mucDigiCol(eventSvc,"/Event/Digi/MucDigiCol");
227 if(!mucDigiCol) {
228 log << MSG::ERROR << "Could not find MucDigiCol!" << endreq;
229 return 0;
230 }
231 SmartDataPtr<RecMucTrackCol> mucTrackCol(eventSvc,"/Event/Recon/RecMucTrackCol");
232 if (!mucTrackCol) {
233 log << MSG::ERROR << "Could not find RecMucTrackCol" << endreq;
234 return 0;
235 }
236
237 int mudigiNum, mutrkNum;
238 mudigiNum = mutrkNum = 0;
239 mudigiNum = mucDigiCol->size(); mutrkNum = mucTrackCol->size();
240
241 bool mucflag1 = mudigiNum >= m_mudigicut;
242 bool mucflag2 = mutrkNum >= m_mutrkcut;
243
244 log << MSG::INFO << "MUC digi:\t" << mudigiNum << "\tMUC track:\t" << mutrkNum << endreq;
245
246 double mudigiFlag,mutrkFlag;
247 if( mucflag1 ) { m_cutpass[9] += 1; mudigiFlag = 1; }
248 if( mucflag2 ) { m_cutpass[10] += 1; mutrkFlag = 1; }
249 if( mucflag1 && mucflag2 ) { m_mucPass = true; m_subpass[3] += 1; }
250 log << MSG::INFO << "MUC selection done!" << endreq;
251
252 if(m_output) {
253 m_c1 = c1; m_c2 = c2; m_r1 = r1; m_r2 = r2; m_z1 = z1; m_z2 = z2; m_p1 = p1; m_p2 = p2;
254 m_t1 = t1; m_t2 = t2;
255 m_e1 = e1; m_e2 = e2; m_theta1 = theta1; m_theta2 = theta2; m_phi1 = phi1; m_phi2 = phi2;
256 m_part = part;
257 m_mudigi = mudigiNum; m_mutrk = mutrkNum;
258 m_zeroCFlag = zeroCFlag; m_vtRZFlag = vtRZFlag; m_pLimFlag = pLimFlag; m_pSymFlag = pSymFlag;
259 m_tLimFlag = tLimFlag; m_eLimFlag = eLimFlag; m_eBBFlag = eBBFlag; m_partFlag = partFlag;
260 m_mudigiFlag = mudigiFlag; m_mutrkFlag = mutrkFlag;
261
262 m_mdcFlag = m_mdcPass; m_tofFlag = m_tofPass; m_emcFlag = m_emcPass; m_mucFlag = m_mucPass;
263
264 m_passtuple->write();
265 }
266
267 // All selections
268 if( m_mdcPass && m_tofPass && m_emcPass && m_mucPass )
269 {
270 m_totalpass += 1;
271 if( part==1 ) return 1; //barrel
272 else return 2; //endcap
273 //setFilterPassed(true);
274 }
275
276 return 0;
277
278}

Referenced by EventPreSelect::execute().

◆ Print()

void DimuPreSelect::Print ( )

Definition at line 280 of file EventPreSelect/EventPreSelect-00-00-21/src/DimuPreSelect.cxx.

281{
282 const string str[3] = {"all","barrel","endcap"};
283
284 cout << "pass 0 - 2 MDC tracks : " << m_cutpass[0] << endl;
285 cout << "pass 1 - 0 charge : " << m_cutpass[1] << endl;
286 cout << "pass 2 - |r|<"<<m_vr0cut<<" && |z|<"<<m_vz0cut<<" : " << m_cutpass[2] << endl;
287 cout << "pass 3 - p < "<<m_pcut_up<<" GeV/c : " << m_cutpass[3] << endl;
288 cout << "pass 4 - p_sym < "<<m_psymcut<<" : " << m_cutpass[4] << endl;
289 cout << "pass 5 - |t1-t1| < "<<m_tcut<<" ns : " << m_cutpass[5] << endl;
290 cout << "pass 6 - "<<m_ecut_down<<" < e < "<<m_ecut_up<<" : " << m_cutpass[6] << endl;
291 cout << "pass 7 - |dth|<"<<m_dthetacut<<" && |dphi|<"<<m_dphicut<<": " << m_cutpass[7] << endl;
292 cout << "pass 8 - "<< str[(int)m_partselect] <<" part is selected : " << m_cutpass[8] << endl;
293 cout << "pass 9 - mudigi >= "<<m_mudigicut<<" : " << m_cutpass[9] << endl;
294 cout << "pass 10- mutrk >= "<<m_mutrkcut<<" : " << m_cutpass[10]<< endl;
295
296 cout << "pass MDC : " << m_subpass[0] << endl;
297 cout << "pass TOF : " << m_subpass[1] << endl;
298 cout << "pass EMC : " << m_subpass[2] << endl;
299 cout << "pass MUC : " << m_subpass[3] << endl;
300
301 cout<<"Total event: " << m_totevent <<endl;
302 cout<<"Dimu event: " << m_totalpass <<endl;
303}

Referenced by EventPreSelect::finalize().


The documentation for this class was generated from the following files: