ATLAS Offline Software
Public Member Functions | Protected Attributes | Private Member Functions | Private Attributes | List of all members
SigAnalysis Class Reference

#include <SigAnalysis.h>

Inheritance diagram for SigAnalysis:
Collaboration diagram for SigAnalysis:

Public Member Functions

 SigAnalysis (const std::string &name)
 
 ~SigAnalysis ()
 
virtual void initialise ()
 standard operation interface More...
 
virtual void execute (const std::vector< TIDA::Track * > &reftracks, const std::vector< TIDA::Track * > &testtracks, TrackAssociator *matcher)
 
virtual void finalise ()
 
void setprint (bool p)
 
virtual void execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher)=0
 
virtual void execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher, TrigObjectMatcher *)
 
virtual void execute (const std::vector< TIDA::Track * > &, const std::vector< TIDA::Track * > &, TrackAssociator *, const TIDA::Event *)
 
virtual void execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher, const TIDA::Event *event, double *)
 
virtual void execute (const std::vector< TIDA::Track * > &, const std::vector< TIDA::Track * > &, TrackAssociator *, const TIDA::Event *, double *, TIDARoiDescriptor *)
 
std::string name () const
 return identifier More...
 
virtual void execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher, TrigObjectMatcher *)
 
virtual void execute (const std::vector< TIDA::Track * > &, const std::vector< TIDA::Track * > &, TrackAssociator *, const TIDA::Event *)
 
virtual void execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher, const TIDA::Event *event, double *)
 
virtual void execute (const std::vector< TIDA::Track * > &, const std::vector< TIDA::Track * > &, TrackAssociator *, const TIDA::Event *, double *, TIDARoiDescriptor *)
 
virtual void execute_vtx (const std::vector< TIDA::Vertex * > &, const std::vector< TIDA::Vertex * > &, const TIDA::Event *=0)
 
void addHistogram (TH1 *h)
 
const std::map< std::string, TH1 * > & getHistograms () const
 access the histograms More...
 
std::map< std::string, TH1 * >::const_iterator begin () const
 
std::map< std::string, TH1 * >::const_iterator end () const
 
void setBeamRef (double x, double y, double z=0)
 set the beamline positions More...
 
void setBeamRef (const std::vector< double > &v)
 
void setBeamTest (double x, double y, double z=0)
 
void setBeamTest (const std::vector< double > &v)
 
double beamRefx () const
 
double beamRefy () const
 
double beamRefz () const
 
double beamTestx () const
 
double beamTesty () const
 
double beamTestz () const
 
TIDA::FeatureStore & store ()
 
const TIDA::Event * event () const
 
void setevent (TIDA::Event *e)
 
const TIDARoiDescriptor * roi () const
 
void setroi (TIDARoiDescriptor *r)
 

Protected Attributes

std::string m_name
 identifier of the of the analysis - also used for the root directory into which the histograms are put More...
 
double m_xBeamReference
 beamline positions reference sample
More...
 
double m_yBeamReference
 
double m_zBeamReference
 
double m_xBeamTest
 test sample More...
 
double m_yBeamTest
 
double m_zBeamTest
 
TIDA::FeatureStore m_store
 
TIDA::Event * m_event
 
TIDARoiDescriptor * m_roi
 

Private Member Functions

void addHistogram (TH1F *h)
 
TH1F * find (const std::string &n)
 

Private Attributes

std::map< std::string, TH1F * > m_histos
 
Efficiency1D * m_eff_pt = nullptr
 
Efficiency1D * m_eff_ptp = nullptr
 
Efficiency1D * m_eff_ptm = nullptr
 
Efficiency1D * m_eff_eta = nullptr
 
Efficiency1D * m_eff_phi = nullptr
 
Efficiency1D * m_eff_z0 = nullptr
 
Efficiency1D * m_eff_d0 = nullptr
 
Efficiency1D * m_eff_a0 = nullptr
 
Efficiency1D * m_purity_pt = nullptr
 
Efficiency1D * m_purity_eta = nullptr
 
Efficiency1D * m_purity_phi = nullptr
 
Efficiency1D * m_purity_z0 = nullptr
 
Efficiency1D * m_purity_d0 = nullptr
 
Efficiency1D * m_purity_a0 = nullptr
 
TH2F * m_h2 = nullptr
 
TH2F * m_h2m = nullptr
 
TH2F * m_h2r = nullptr
 
int m_Nreco
 number of reconstructed tracks More...
 
bool m_print
 flag to print out the matched tracks etc More...
 
int m_icount
 Event counter. More...
 
TIDDirectory * m_dir
 

Detailed Description

Definition at line 28 of file SigAnalysis.h.

Constructor & Destructor Documentation

◆ SigAnalysis()

SigAnalysis::SigAnalysis ( const std::string &  name)
inline

Definition at line 36 of file SigAnalysis.h.

36  : TrackAnalysis( name ), m_Nreco(0), m_print(true), m_icount(0), m_dir(0) {
37  std::cout << "SigAnalysis::SigAnalysis() " << TrackAnalysis::name() << std::endl;
38  }

◆ ~SigAnalysis()

SigAnalysis::~SigAnalysis ( )
inline

Definition at line 40 of file SigAnalysis.h.

40  {
41  std::cout << "SigAnalysis::~SigAnalysis() " << name() << std::endl;
42  std::map<std::string, TH1F*>::iterator hitr=m_histos.begin();
43  std::map<std::string, TH1F*>::iterator hend=m_histos.end();
44  for ( ; hitr!=hend ; hitr++ ) delete hitr->second;
45 
46  Efficiency1D* heff[4] = { m_eff_pt, m_eff_eta, m_eff_phi, m_eff_z0 };
47  Efficiency1D* hpurity[4] = { m_purity_pt, m_purity_eta, m_purity_phi, m_purity_z0 };
48 
49  for ( int i=4 ; i-- ; ) { delete heff[i]; delete hpurity[i]; }
50 
51  delete m_eff_ptp;
52  delete m_eff_ptm;
53 
54  }

Member Function Documentation

◆ addHistogram() [1/2]

void TrackAnalysis::addHistogram ( TH1 *  h)
inlineinherited

Definition at line 97 of file TrackAnalysis.h.

97  {
98  std::string name = h->GetName();
99  m_histos.insert( std::map<std::string, TH1*>::value_type( name, h) );
100  }

◆ addHistogram() [2/2]

void SigAnalysis::addHistogram ( TH1F *  h)
inlineprivate

Definition at line 70 of file SigAnalysis.h.

70  {
71  std::string name = h->GetName();
72  m_histos.insert( std::map<std::string, TH1F*>::value_type( name, h) );
73  }

◆ beamRefx()

double TrackAnalysis::beamRefx ( ) const
inlineinherited

Definition at line 124 of file TrackAnalysis.h.

124 { return m_xBeamReference; }

◆ beamRefy()

double TrackAnalysis::beamRefy ( ) const
inlineinherited

Definition at line 125 of file TrackAnalysis.h.

125 { return m_yBeamReference; }

◆ beamRefz()

double TrackAnalysis::beamRefz ( ) const
inlineinherited

Definition at line 126 of file TrackAnalysis.h.

126 { return m_zBeamReference; }

◆ beamTestx()

double TrackAnalysis::beamTestx ( ) const
inlineinherited

Definition at line 128 of file TrackAnalysis.h.

128 { return m_xBeamTest; }

◆ beamTesty()

double TrackAnalysis::beamTesty ( ) const
inlineinherited

Definition at line 129 of file TrackAnalysis.h.

129 { return m_yBeamTest; }

◆ beamTestz()

double TrackAnalysis::beamTestz ( ) const
inlineinherited

Definition at line 130 of file TrackAnalysis.h.

130 { return m_zBeamTest; }

◆ begin()

std::map<std::string, TH1*>::const_iterator TrackAnalysis::begin ( ) const
inlineinherited

Definition at line 105 of file TrackAnalysis.h.

105 { return m_histos.begin(); }

◆ end()

std::map<std::string, TH1*>::const_iterator TrackAnalysis::end ( ) const
inlineinherited

Definition at line 106 of file TrackAnalysis.h.

106 { return m_histos.end(); }

◆ event()

const TIDA::Event* TrackAnalysis::event ( ) const
inlineinherited

Definition at line 135 of file TrackAnalysis.h.

135 { return m_event; }

◆ execute() [1/10]

virtual void TrackAnalysis::execute
inline

Definition at line 68 of file TrackAnalysis.h.

71  { }

◆ execute() [2/10]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > &  ,
const std::vector< TIDA::Track * > &  ,
TrackAssociator *  ,
const TIDA::Event *   
)
inlinevirtualinherited

Reimplemented in ConfAnalysis, and AnalysisR3_Tier0.

Definition at line 68 of file TrackAnalysis.h.

71  { }

◆ execute() [3/10]

virtual void TrackAnalysis::execute
inline

Definition at line 81 of file TrackAnalysis.h.

86  { }

◆ execute() [4/10]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > &  ,
const std::vector< TIDA::Track * > &  ,
TrackAssociator *  ,
const TIDA::Event *  ,
double *  ,
TIDARoiDescriptor *   
)
inlinevirtualinherited

Reimplemented in AnalysisR3_Tier0.

Definition at line 81 of file TrackAnalysis.h.

86  { }

◆ execute() [5/10]

void SigAnalysis::execute ( const std::vector< TIDA::Track * > &  reftracks,
const std::vector< TIDA::Track * > &  testtracks,
TrackAssociator *  matcher 
)
virtual

Implements TrackAnalysis.

Definition at line 205 of file SigAnalysis.cxx.

208 {
209  if ( m_print ) std::cout << "SigAnalysis::execute() \t " << name()
210  << "\tref " << reftracks.size()
211  << "\ttest " << testtracks.size() << std::endl;
212 
213  // std::cout << "\tx " << m_xBeamReference << "\ty " << m_yBeamReference
214  // << "\tx " << m_xBeamTest << "\ty " << m_yBeamTest << std::endl;
215 
216 
217  // std::cout << "SigAnalysis (resolutions really) filling " << std::endl;
218 
219  // should have these as a class variable
220  static const std::string varName[10] = { "pT", "eta", "phi", "z0", "d0", "a0", "nsct", "npix", "ntrt", "nstraw" };
221 
222  // std::cout << "SigAnalysis ref size " << reftracks.size() << "\ttest size " << testtracks.size() << std::endl;
223 
224  std::map<std::string, TH1F*>::iterator hmitr = m_histos.find("ntracks");
225  if ( hmitr!=m_histos.end() ) hmitr->second->Fill( reftracks.size() );
226 
227  hmitr = m_histos.find("ntracks_rec");
228  if ( hmitr!=m_histos.end() ) hmitr->second->Fill( testtracks.size() );
229 
230  bool dump = false;
231 
232  for ( int i=reftracks.size() ; i-- ; ) {
233 
234  double pTt = reftracks[i]->pT();
235  double z0t = reftracks[i]->z0();
236  double etat = reftracks[i]->eta();
237  double phit = reftracks[i]->phi();
238  double d0t = reftracks[i]->a0();
239  // this will be changed when we know the beam spot position
240  // double a0t = reftracks[i]->a0() + sin(phit)*m_xBeam - cos(phit)*m_yBeam;
241  double a0t = reftracks[i]->a0() + sin(phit)*m_xBeamReference - cos(phit)*m_yBeamReference;
242 
243  double nsctt = reftracks[i]->sctHits();
244  double npixt = reftracks[i]->pixelHits();
245 
246  double ntrtt = reftracks[i]->trHits();
247  double nstrawt = reftracks[i]->strawHits();
248 
249  // std::cout << "Fill h2 " << " " << m_h2m << " " << *reftracks[i] << std::endl;
250 
251  m_h2->Fill( phit, d0t );
252 
253  const TIDA::Track* matchedreco = matcher->matched(reftracks[i]);
254 
255  // std::cout << "\t\tSigAnalysis " << name() << "\t" << i << " " << *reftracks[i] << " -> ";
256 
257  // raw reference track distributions
258  double vpart[10] = { std::fabs(pTt), etat, phit, z0t, d0t, a0t, nsctt, npixt, ntrtt, nstrawt };
259  for ( int it=0 ; it<10 ; it++ ) {
260  // std::string hname = varName[it];
261  // std::map<std::string, TH1F*>::iterator hmitr = m_histos.find(hname);
262  // if ( hmitr!=m_histos.end() ) hmitr->second->Fill( vpart[it] );
263 
264  if ( TH1F* hptr = find( varName[it] ) ) hptr->Fill( vpart[it] );
265  else std::cerr << "hmmm histo " << varName[it] << " not found" << std::endl;
266 
267  }
268 
269 
270  if ( matchedreco ) {
271 
272  // efficiency histos
273  m_eff_pt->Fill(std::fabs(pTt));
274  m_eff_z0->Fill(z0t);
275  m_eff_eta->Fill(etat);
276  m_eff_phi->Fill(phit);
277  m_eff_d0->Fill(d0t);
278  m_eff_a0->Fill(a0t);
279 
280  // signed pT
281  if ( pTt<0 ) m_eff_ptm->Fill(std::fabs(pTt));
282  else m_eff_ptp->Fill(std::fabs(pTt));
283 
284  // residual histos
285  double pTr = matchedreco->pT();
286  double z0r = matchedreco->z0();
287  double etar = matchedreco->eta();
288  double phir = matchedreco->phi();
289  double d0r = matchedreco->a0();
290  double a0r = matchedreco->a0() + sin(phir)*m_xBeamTest - cos(phir)*m_yBeamTest; // this will be changed when we know the beam spot position
291 
292  if ( m_h2m ) m_h2m->Fill( phit, d0t );
293 
294  // if ( m_print ) std::cout << "SigAnalysis::execute() \t " << name() << "\t" << i << " "
295  if ( m_print ) std::cout << "SigAnalysis::execute() \t\t" << i << " "
296  << *reftracks[i] << " -> " << *matchedreco << "\t"
297  << pTr << " " << pTt << " " << d0r << " " << d0t << std::endl;
298 
299  double vres[6] = { 1/pTt-1/pTr, etat-etar, phit-phir, z0t-z0r, d0t-d0r, a0t-a0r };
300  for ( int it=0 ; it<6 ; it++ ) {
301  if ( TH1F* hptr = find(varName[it]+"_res") ) hptr->Fill( vres[it] );
302  else std::cerr << "hmmm histo " << varName[it]+"_res" << " not found" << std::endl;
303  }
304 
305  // in this loop over the reference tracks, could fill efficiency
306  // histograms
307 
308  }
309  else {
310  // fill efficiencies with unmatched histos
311  // std::cout << "NULL" << std::endl;
312  m_eff_pt->FillDenom(std::fabs(pTt));
313  m_eff_z0->FillDenom(z0t);
314  m_eff_eta->FillDenom(etat);
315  m_eff_phi->FillDenom(phit);
316  m_eff_d0->FillDenom(d0t);
317  m_eff_a0->FillDenom(a0t);
318 
319  // signed pT
320  if ( pTt<0 ) m_eff_ptm->FillDenom(std::fabs(pTt));
321  else m_eff_ptp->FillDenom(std::fabs(pTt));
322 
323  if ( std::fabs(pTt)>4000 ) dump = true;
324  }
325 
326  }
327 
328 
329  // for fake/purity histograms, loop over the test tracks
330  // and get the corresponding matched reference tracks from the
331  // reverse map in the TrackAscociator class - revmatched()
332 
333  if ( m_print ) std::cout << "SigAnalysis::execute() \t " << name() << "\t " << m_icount << " events\t " << testtracks.size() << " tracks (" << m_Nreco << ")" << "\n---------------" << std::endl;
334 
335  m_icount++;
336 
337  m_Nreco += testtracks.size();
338 
339  for ( int i=testtracks.size() ; i-- ; ) {
340 
341  // std::cout << "\t\tSigAnalysis purity " << name() << "\t" << i << " " << *testtracks[i] << " -> ";
342 
343  // double pTr = std::fabs(testtracks[i]->pT());
344  double pTr = testtracks[i]->pT();
345  double etar = testtracks[i]->eta();
346  double phir = testtracks[i]->phi();
347  double z0r = testtracks[i]->z0();
348  double d0r = testtracks[i]->a0();
349  double a0r = testtracks[i]->a0() + sin(phir)*m_xBeamTest - cos(phir)*m_yBeamTest; // this will be changed when we know the beam spot position
350  // double a0rp = testtracks[i]->a0() - sin(phir)*m_xBeam - cos(phir)*m_yBeam; // this will be changed when we know the beam spot position
351 
352  // std::cout << "d0 " << d0r << "\tphi " << phir << "\tx " << m_xBeamTest << "\ty " << m_yBeamTest << std::endl;
353 
354  double nsctr = testtracks[i]->sctHits();
355  double npixr = testtracks[i]->pixelHits();
356 
357  double ntrtr = testtracks[i]->trHits();
358  double nstrawr = testtracks[i]->strawHits();
359 
360  if ( m_h2r ) m_h2r->Fill( phir, d0r );
361 
362  const TIDA::Track* matchedref = matcher->revmatched(testtracks[i]);
363 
364  // if ( matchedref ) std::cout << *matchedref << std::endl;
365  // else std::cout << "NULL" << std::endl;
366 
367  // raw test track distributions
368  double vpart[10] = { std::fabs(pTr), etar, phir, z0r, d0r, a0r, nsctr, npixr, ntrtr, nstrawr };
369  for ( int it=0 ; it<10 ; it++ ) {
370  // std::string hname = name()+"_"+varName[it]+"_rec";
371  // std::string hname = varName[it]+"_rec";
372  // std::map<std::string, TH1F*>::iterator hmitr = m_histos.find(hname);
373  // if ( hmitr!=m_histos.end() ) hmitr->second->Fill( vpar[it] );
374  // else std::cerr << "hmmm histo " << hname << " not found" << std::endl;
375  if ( TH1F* hptr = find(varName[it]+"_rec") ) hptr->Fill( vpart[it] );
376  else std::cerr << "hmmm histo " << varName[it]+"_rec" << " not found" << std::endl;
377  }
378 
379  // purities
380  if ( matchedref ) {
381 
382  // std::cout << *matchedref << std::endl;
383 
384  m_purity_pt->Fill(std::fabs(pTr));
385  m_purity_z0->Fill(z0r);
386  m_purity_eta->Fill(etar);
387  m_purity_phi->Fill(phir);
388  m_purity_d0->Fill(d0r);
389  m_purity_a0->Fill(a0r);
390 
391  }
392  else {
393  // std::cout << "NULL" << std::endl;
394  m_purity_pt->FillDenom(std::fabs(pTr));
395  m_purity_z0->FillDenom(z0r);
396  m_purity_eta->FillDenom(etar);
397  m_purity_phi->FillDenom(phir);
398  m_purity_d0->FillDenom(d0r);
399  m_purity_a0->FillDenom(a0r);
400  }
401 
402  }
403 
404  if ( dump && m_print ) {
405 
406  std::cout << "SigAnalysis::execute() missed a high pT track - dumping tracks" << std::endl;
407 
408  for ( int i=reftracks.size() ; i-- ; ) {
409 
410  if ( std::fabs( reftracks[i]->pT() ) > 1000 ) {
411  std::cout << "\t dump " << *reftracks[i];
412  const TIDA::Track* matchedreco = matcher->matched(reftracks[i]);
413  if ( matchedreco ) std::cout << " <--> " << *matchedreco << std::endl;
414  else std::cout << std::endl;
415  }
416 
417  }
418 
419  for ( int i=testtracks.size() ; i-- ; ) {
420  const TIDA::Track* matchedref = matcher->revmatched(testtracks[i]);
421  if ( matchedref==0 ) std::cout << "\t\t\t\t\t " << *testtracks[i] << std::endl;
422  }
423 
424  }
425 
426  // std::cout << "SigAnalysis::execute() exiting" << std::endl;
427 
428 }

◆ execute() [6/10]

virtual void TrackAnalysis::execute

◆ execute() [7/10]

virtual void TrackAnalysis::execute
inline

Definition at line 73 of file TrackAnalysis.h.

77  {
78  execute( tracks1, tracks2, matcher, event );
79  }

◆ execute() [8/10]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > &  tracks1,
const std::vector< TIDA::Track * > &  tracks2,
TrackAssociator *  matcher,
const TIDA::Event *  event,
double *   
)
inlinevirtualinherited

Reimplemented in AnalysisR3_Tier0.

Definition at line 73 of file TrackAnalysis.h.

77  {
78  execute( tracks1, tracks2, matcher, event );
79  }

◆ execute() [9/10]

virtual void TrackAnalysis::execute
inline

Definition at line 61 of file TrackAnalysis.h.

64  {
65  execute( tracks1, tracks2, matcher );
66  }

◆ execute() [10/10]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > &  tracks1,
const std::vector< TIDA::Track * > &  tracks2,
TrackAssociator *  matcher,
TrigObjectMatcher *   
)
inlinevirtualinherited

Reimplemented in ConfAnalysis.

Definition at line 61 of file TrackAnalysis.h.

64  {
65  execute( tracks1, tracks2, matcher );
66  }

◆ execute_vtx()

virtual void TrackAnalysis::execute_vtx ( const std::vector< TIDA::Vertex * > &  ,
const std::vector< TIDA::Vertex * > &  ,
const TIDA::Event *  = 0 
)
inlinevirtualinherited

Reimplemented in AnalysisR3_Tier0, and Analysis_Tier0.

Definition at line 89 of file TrackAnalysis.h.

91  { }

◆ finalise()

void SigAnalysis::finalise ( )
virtual

Implements TrackAnalysis.

Definition at line 177 of file SigAnalysis.cxx.

177  {
178  std::cout << "SigAnalysis::finalise() " << name() << "\tNreco " << m_Nreco << " tracks" << std::endl;
179 
180  // if ( m_Nreco==0 ) return;
181 
182  m_dir->push();
183 
184  std::cout << "SigAnalysis::finalise() " << gDirectory->GetName() << std::endl;
185 
186  // std::map<std::string, TH1F*>::iterator hitr=m_histos.begin();
187  // std::map<std::string, TH1F*>::iterator hend=m_histos.end();
188  // for ( ; hitr!=hend ; hitr++ ) hitr->second->Write();
189 
190  // std::cout << "DBG >" << m_eff_pt->Hist()->GetName() << "< DBG" << std::endl;
191 
192  Efficiency1D* heff[8] = { m_eff_pt, m_eff_eta, m_eff_phi, m_eff_z0, m_eff_d0, m_eff_a0, m_eff_ptm, m_eff_ptp };
193  for ( int i=8 ; i-- ; ) { heff[i]->finalise(); } // heff[i]->Hist()->Write(); }
194 
195  // std::cout << "DBG >" << m_purity_pt->Hist()->GetName() << "< DBG" << std::endl;
196 
197  Efficiency1D* hpurity[6] = { m_purity_pt, m_purity_eta, m_purity_phi, m_purity_z0, m_purity_d0, m_purity_a0 };
198  for ( int i=6 ; i-- ; ) { hpurity[i]->finalise(); } // hpurity[i]->Hist()->Write(); }
199 
200  m_dir->pop();
201 
202 }

◆ find()

TH1F* SigAnalysis::find ( const std::string &  n)
inlineprivate

Definition at line 75 of file SigAnalysis.h.

75  {
76  std::map<std::string, TH1F*>::iterator hmitr = m_histos.find(n);
77  if ( hmitr!=m_histos.end() ) return hmitr->second;
78  else return 0;
79  }

◆ getHistograms()

const std::map<std::string, TH1*>& TrackAnalysis::getHistograms ( ) const
inlineinherited

access the histograms

Definition at line 103 of file TrackAnalysis.h.

103 { return m_histos; }

◆ initialise()

void SigAnalysis::initialise ( )
virtual

standard operation interface

Implements TrackAnalysis.

Definition at line 13 of file SigAnalysis.cxx.

13  {
14 
15  std::cout << "SigAnalysis::initialise() " << name() << std::endl;
16 
17  //+++ pT ranges
18  // double tmp_maxPt = 50000.;
19  double tmp_absResPt = 0.0005;
20 
21  const int pTResBins = 100;
22 
23  //+++ Eta ranges
24  double tmp_maxEta = 3.;
25  double tmp_absResEta = 0.04; // 0.0005;
26 
27  //+++ Phi ranges
28  double tmp_maxPhi = 3.142;
29  double tmp_absResPhi = 0.02; // 0.0001;
30 
31 
32  const int etaBins = 60;
33  const int etaResBins = 300;
34 
35  const int phiBins = 36;
36  const int phiResBins = 100;
37 
38  const int zBins = 50;
39  const double zMax = 400;
40 
41  const int zresBins = 100;
42  const double zresMax = 10;
43 
44  const int d0Bins = 100;
45  const double d0Max = 10;
46 
47  const int d0resBins = 100;
48  const double d0resMax = 10;
49 
50  // beamspot corrected position
51 
52  const int a0Bins = 100;
53  const double a0Max = 10;
54 
55  const int a0resBins = 100;
56  const double a0resMax = 5;
57 
58  //+++ Book histograms
59 
60  // calculate a logarithmic binning in pt
61 
62  const int ptnbins = 20;
63  // const int ptnbins = 30;
64  double ptbinlims[ptnbins+1];
65  for ( int i=0 ; i<=ptnbins ; i++ ) { ptbinlims[i] = std::pow(10, 2.0*i/ptnbins+2); }
66  // for ( int i=0 ; i<=ptnbins ; i++ ) { ptbinlims[i] = std::pow(10, 2.3*i/ptnbins+2); }
67 
68 
69 
70  TDirectory* dir = gDirectory;
71 
72  std::cout << "SigAnalysis::initialize() Directory " << gDirectory->GetName() << " " << name() << std::endl;
73 
74  m_dir = new TIDDirectory(name());
75  m_dir->push();
76 
77  std::cout << "SigAnalysis::initialize() Directory " << gDirectory->GetName() << " package directory, " << name() << std::endl;
78 
79  Efficiency1D* heff[8];
80  Efficiency1D* hpurity[6];
81 
82  // "reference" quantities
83  addHistogram( new TH1F( "pT", "pT", ptnbins, ptbinlims ) );
84  addHistogram( new TH1F( "eta", "eta", etaBins, -tmp_maxEta, tmp_maxEta ) );
85  addHistogram( new TH1F( "phi", "phi", phiBins, -tmp_maxPhi, tmp_maxPhi ) );
86  addHistogram( new TH1F( "z0", "z0", zBins, -zMax, zMax ) );
87  addHistogram( new TH1F( "d0", "d0", d0Bins, -d0Max, d0Max ) );
88  addHistogram( new TH1F( "a0", "a0", a0Bins, -a0Max, a0Max ) );
89 
90  // efficienies and purities
91  heff[0] = new Efficiency1D( find("pT"), "pT_eff" );
92  heff[1] = new Efficiency1D( find("eta"), "eta_eff" );
93  heff[2] = new Efficiency1D( find("phi"), "phi_eff" );
94  heff[3] = new Efficiency1D( find("z0"), "z0_eff" );
95  heff[4] = new Efficiency1D( find("d0"), "d0_eff" );
96  heff[5] = new Efficiency1D( find("a0"), "a0_eff" );
97 
98  heff[6] = new Efficiency1D( find("pT"), "pTm_eff" );
99  heff[7] = new Efficiency1D( find("pT"), "pTp_eff" );
100 
101  m_eff_pt = heff[0];
102  m_eff_eta = heff[1];
103  m_eff_phi = heff[2];
104  m_eff_z0 = heff[3];
105  m_eff_d0 = heff[4];
106  m_eff_a0 = heff[5];
107 
108  m_eff_ptm = heff[6];
109  m_eff_ptp = heff[7];
110 
111  hpurity[0] = new Efficiency1D( find("pT"), "pT_pur" );
112  hpurity[1] = new Efficiency1D( find("eta"), "eta_pur" );
113  hpurity[2] = new Efficiency1D( find("phi"), "phi_pur" );
114  hpurity[3] = new Efficiency1D( find("z0"), "z0_pur" );
115  hpurity[4] = new Efficiency1D( find("d0"), "d0_pur" );
116  hpurity[5] = new Efficiency1D( find("a0"), "a0_pur" );
117 
118  m_purity_pt = hpurity[0];
119  m_purity_eta = hpurity[1];
120  m_purity_phi = hpurity[2];
121  m_purity_z0 = hpurity[3];
122  m_purity_d0 = hpurity[4];
123  m_purity_a0 = hpurity[5];
124 
125  // "test" quantities
126  addHistogram( new TH1F( "pT_rec", "pT_rec", ptnbins, ptbinlims ) );
127  addHistogram( new TH1F( "eta_rec", "eta_rec", etaBins, -tmp_maxEta, tmp_maxEta ) );
128  addHistogram( new TH1F( "phi_rec", "phi_rec", phiBins, -tmp_maxPhi, tmp_maxPhi ) );
129  addHistogram( new TH1F( "z0_rec", "z0_rec", zBins, -zMax, zMax ) );
130  addHistogram( new TH1F( "d0_rec", "d0_rec", d0Bins, -d0Max, d0Max ) );
131  addHistogram( new TH1F( "a0_rec", "a0_rec", a0Bins, -a0Max, a0Max ) );
132 
133  // resolutions
134  addHistogram( new TH1F( "pT_res", "pT_res", pTResBins, -tmp_absResPt, tmp_absResPt ) );
135  addHistogram( new TH1F( "eta_res", "eta_res", etaResBins, -tmp_absResEta, tmp_absResEta ) );
136  addHistogram( new TH1F( "phi_res", "phi_res", phiResBins, -tmp_absResPhi, tmp_absResPhi ) );
137  addHistogram( new TH1F( "z0_res", "z0_res", zresBins, -zresMax, zresMax ) );
138  addHistogram( new TH1F( "d0_res", "d0_res", d0resBins, -0.5*d0resMax, 0.5*d0resMax ) );
139  addHistogram( new TH1F( "a0_res", "a0_res", a0resBins, -0.5*a0resMax, 0.5*a0resMax ) );
140 
141  // hit occupancies
142 
143  int NHits = 40;
144  int Ntracks = 50;
145 
146  addHistogram( new TH1F( "nsct", "nsct", NHits, -0.5, float(NHits-0.5) ) );
147  addHistogram( new TH1F( "nsct_rec", "nsct_rec", NHits, -0.5, float(NHits-0.5) ) );
148 
149  addHistogram( new TH1F( "npix", "npix", NHits, -0.5, float(NHits-0.5) ) );
150  addHistogram( new TH1F( "npix_rec", "npix_rec", NHits, -0.5, float(NHits-0.5) ) );
151 
152  addHistogram( new TH1F( "ntrt", "ntrt", NHits, -0.5, float(NHits-0.5) ) );
153  addHistogram( new TH1F( "ntrt_rec", "ntrt_rec", NHits, -0.5, float(NHits-0.5) ) );
154 
155  addHistogram( new TH1F( "nstraw", "nstraw", NHits*4, -0.5, float(4*NHits-0.5) ) );
156  addHistogram( new TH1F( "nstraw_rec", "nstraw_rec", NHits*4, -0.5, float(4*NHits-0.5) ) );
157 
158  addHistogram( new TH1F( "ntracks", "ntracks", Ntracks, -0.5, float(Ntracks+0.5) ) );
159  addHistogram( new TH1F( "ntracks_rec", "ntracks_rec", Ntracks, -0.5, float(Ntracks+0.5) ) );
160 
161 
162  // beam offset fitting histos
163  m_h2 = new TH2F( "d0vphi", "d0vphi", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
164  m_h2r = new TH2F( "d0vphi_rec", "d0vphi_rec", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
165  m_h2m = new TH2F( "d0vphi_match", "d0vphi_match", phiBins, -3.142, 3.142, d0Bins, -d0Max, d0Max );
166 
167  m_dir->pop();
168 
169  dir->cd();
170 
171  // std::cout << "initialize() Directory " << gDirectory->GetName() << " on leaving" << std::endl;
172 
173 }

◆ name()

std::string TrackAnalysis::name ( ) const
inlineinherited

return identifier

Definition at line 52 of file TrackAnalysis.h.

52 { return m_name; }

◆ roi()

const TIDARoiDescriptor* TrackAnalysis::roi ( ) const
inlineinherited

Definition at line 138 of file TrackAnalysis.h.

138 { return m_roi; }

◆ setBeamRef() [1/2]

void TrackAnalysis::setBeamRef ( const std::vector< double > &  v)
inlineinherited

Definition at line 112 of file TrackAnalysis.h.

112  {
113  if ( v.size()>0 ) m_xBeamReference = v[0];
114  if ( v.size()>1 ) m_yBeamReference = v[1];
115  if ( v.size()>2 ) m_zBeamReference = v[2];
116  }

◆ setBeamRef() [2/2]

void TrackAnalysis::setBeamRef ( double  x,
double  y,
double  z = 0 
)
inlineinherited

set the beamline positions

Definition at line 109 of file TrackAnalysis.h.

109 { m_xBeamReference = x; m_yBeamReference = y; m_zBeamReference = z; }

◆ setBeamTest() [1/2]

void TrackAnalysis::setBeamTest ( const std::vector< double > &  v)
inlineinherited

Definition at line 118 of file TrackAnalysis.h.

118  {
119  if ( v.size()>0 ) m_xBeamTest = v[0];
120  if ( v.size()>1 ) m_yBeamTest = v[1];
121  if ( v.size()>2 ) m_zBeamTest = v[2];
122  }

◆ setBeamTest() [2/2]

void TrackAnalysis::setBeamTest ( double  x,
double  y,
double  z = 0 
)
inlineinherited

Definition at line 110 of file TrackAnalysis.h.

110 { m_xBeamTest = x; m_yBeamTest = y; m_zBeamTest = z; }

◆ setevent()

void TrackAnalysis::setevent ( TIDA::Event *  e)
inlineinherited

Definition at line 136 of file TrackAnalysis.h.

136 { m_event=e; }

◆ setprint()

void SigAnalysis::setprint ( bool  p)
inline

Definition at line 66 of file SigAnalysis.h.

66 { m_print=p; }

◆ setroi()

void TrackAnalysis::setroi ( TIDARoiDescriptor *  r)
inlineinherited

Definition at line 139 of file TrackAnalysis.h.

139 { m_roi=r; }

◆ store()

TIDA::FeatureStore& TrackAnalysis::store ( )
inlineinherited

Definition at line 133 of file TrackAnalysis.h.

133 { return m_store; }

Member Data Documentation

◆ m_dir

TIDDirectory* SigAnalysis::m_dir
private

Definition at line 116 of file SigAnalysis.h.

◆ m_eff_a0

Efficiency1D* SigAnalysis::m_eff_a0 = nullptr
private

Definition at line 93 of file SigAnalysis.h.

◆ m_eff_d0

Efficiency1D* SigAnalysis::m_eff_d0 = nullptr
private

Definition at line 92 of file SigAnalysis.h.

◆ m_eff_eta

Efficiency1D* SigAnalysis::m_eff_eta = nullptr
private

Definition at line 89 of file SigAnalysis.h.

◆ m_eff_phi

Efficiency1D* SigAnalysis::m_eff_phi = nullptr
private

Definition at line 90 of file SigAnalysis.h.

◆ m_eff_pt

Efficiency1D* SigAnalysis::m_eff_pt = nullptr
private

Definition at line 85 of file SigAnalysis.h.

◆ m_eff_ptm

Efficiency1D* SigAnalysis::m_eff_ptm = nullptr
private

Definition at line 87 of file SigAnalysis.h.

◆ m_eff_ptp

Efficiency1D* SigAnalysis::m_eff_ptp = nullptr
private

Definition at line 86 of file SigAnalysis.h.

◆ m_eff_z0

Efficiency1D* SigAnalysis::m_eff_z0 = nullptr
private

Definition at line 91 of file SigAnalysis.h.

◆ m_event

TIDA::Event* TrackAnalysis::m_event
protectedinherited

Definition at line 164 of file TrackAnalysis.h.

◆ m_h2

TH2F* SigAnalysis::m_h2 = nullptr
private

Definition at line 103 of file SigAnalysis.h.

◆ m_h2m

TH2F* SigAnalysis::m_h2m = nullptr
private

Definition at line 104 of file SigAnalysis.h.

◆ m_h2r

TH2F* SigAnalysis::m_h2r = nullptr
private

Definition at line 105 of file SigAnalysis.h.

◆ m_histos

std::map<std::string, TH1F*> SigAnalysis::m_histos
private

Definition at line 83 of file SigAnalysis.h.

◆ m_icount

int SigAnalysis::m_icount
private

Event counter.

Definition at line 114 of file SigAnalysis.h.

◆ m_name

std::string TrackAnalysis::m_name
protectedinherited

identifier of the of the analysis - also used for the root directory into which the histograms are put

Definition at line 145 of file TrackAnalysis.h.

◆ m_Nreco

int SigAnalysis::m_Nreco
private

number of reconstructed tracks

Definition at line 108 of file SigAnalysis.h.

◆ m_print

bool SigAnalysis::m_print
private

flag to print out the matched tracks etc

Definition at line 111 of file SigAnalysis.h.

◆ m_purity_a0

Efficiency1D* SigAnalysis::m_purity_a0 = nullptr
private

Definition at line 100 of file SigAnalysis.h.

◆ m_purity_d0

Efficiency1D* SigAnalysis::m_purity_d0 = nullptr
private

Definition at line 99 of file SigAnalysis.h.

◆ m_purity_eta

Efficiency1D* SigAnalysis::m_purity_eta = nullptr
private

Definition at line 96 of file SigAnalysis.h.

◆ m_purity_phi

Efficiency1D* SigAnalysis::m_purity_phi = nullptr
private

Definition at line 97 of file SigAnalysis.h.

◆ m_purity_pt

Efficiency1D* SigAnalysis::m_purity_pt = nullptr
private

Definition at line 95 of file SigAnalysis.h.

◆ m_purity_z0

Efficiency1D* SigAnalysis::m_purity_z0 = nullptr
private

Definition at line 98 of file SigAnalysis.h.

◆ m_roi

TIDARoiDescriptor* TrackAnalysis::m_roi
protectedinherited

Definition at line 165 of file TrackAnalysis.h.

◆ m_store

TIDA::FeatureStore TrackAnalysis::m_store
protectedinherited

Definition at line 162 of file TrackAnalysis.h.

◆ m_xBeamReference

double TrackAnalysis::m_xBeamReference
protectedinherited

beamline positions reference sample

Definition at line 153 of file TrackAnalysis.h.

◆ m_xBeamTest

double TrackAnalysis::m_xBeamTest
protectedinherited

test sample

Definition at line 158 of file TrackAnalysis.h.

◆ m_yBeamReference

double TrackAnalysis::m_yBeamReference
protectedinherited

Definition at line 154 of file TrackAnalysis.h.

◆ m_yBeamTest

double TrackAnalysis::m_yBeamTest
protectedinherited

Definition at line 159 of file TrackAnalysis.h.

◆ m_zBeamReference

double TrackAnalysis::m_zBeamReference
protectedinherited

Definition at line 155 of file TrackAnalysis.h.

◆ m_zBeamTest

double TrackAnalysis::m_zBeamTest
protectedinherited

Definition at line 160 of file TrackAnalysis.h.

The documentation for this class was generated from the following files:
TrackAnalysis::m_roi
TIDARoiDescriptor * m_roi
Definition: TrackAnalysis.h:165
xAOD::iterator
JetConstituentVector::iterator iterator
Definition: JetConstituentVector.cxx:68
SigAnalysis::m_purity_d0
Efficiency1D * m_purity_d0
Definition: SigAnalysis.h:99
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
beamspotman.r
def r
Definition: beamspotman.py:676
CalculateHighPtTerm.pT
pT
Definition: ICHEP2016/CalculateHighPtTerm.py:57
TIDDirectory::pop
void pop()
Definition: TIDDirectory.h:79
TrackAnalysis::execute
virtual void execute(const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher)=0
TIDA::Track::a0
double a0() const
Definition: Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/Track.h:49
ConvertOldUJHistosToNewHistos.etaBins
list etaBins
Definition: ConvertOldUJHistosToNewHistos.py:145
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:20
SigAnalysis::m_eff_z0
Efficiency1D * m_eff_z0
Definition: SigAnalysis.h:91
skel.it
it
Definition: skel.GENtoEVGEN.py:396
T_Efficiency::finalise
void finalise(double scale=100)
actually calculate the efficiencies
Definition: T_Efficiency.h:80
TIDA::Track::pT
double pT() const
Definition: Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/Track.h:50
SigAnalysis::find
TH1F * find(const std::string &n)
Definition: SigAnalysis.h:75
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
SigAnalysis::m_purity_phi
Efficiency1D * m_purity_phi
Definition: SigAnalysis.h:97
TIDA::Track::z0
double z0() const
Definition: Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/Track.h:48
python.TrigEgammaMonitorHelper.TH2F
def TH2F(name, title, nxbins, bins_par2, bins_par3, bins_par4, bins_par5=None, bins_par6=None, path='', **kwargs)
Definition: TrigEgammaMonitorHelper.py:45
x
#define x
TIDA::Associator::matched
virtual const S * matched(T *t)
Definition: TIDAAssociator.h:45
run_Egamma1_LArStrip_Fex.dump
dump
Definition: run_Egamma1_LArStrip_Fex.py:88
TIDDirectory
Definition: TIDDirectory.h:25
Efficiency1D
Definition: Efficiency1D.h:19
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
ParseInputs.gDirectory
gDirectory
Definition: Final2012/ParseInputs.py:133
lumiFormat.i
int i
Definition: lumiFormat.py:85
z
#define z
SigAnalysis::m_purity_pt
Efficiency1D * m_purity_pt
Definition: SigAnalysis.h:95
beamspotman.n
n
Definition: beamspotman.py:731
PixelAthClusterMonAlgCfg.varName
string varName
end cluster ToT and charge
Definition: PixelAthClusterMonAlgCfg.py:125
Efficiency1D::FillDenom
virtual void FillDenom(double x, float w=1)
Definition: Efficiency1D.h:42
SigAnalysis::m_eff_phi
Efficiency1D * m_eff_phi
Definition: SigAnalysis.h:90
SigAnalysis::m_print
bool m_print
flag to print out the matched tracks etc
Definition: SigAnalysis.h:111
TrackAnalysis::name
std::string name() const
return identifier
Definition: TrackAnalysis.h:52
TIDA::Associator::revmatched
virtual const T * revmatched(S *t)
Definition: TIDAAssociator.h:52
TIDA::Track::phi
double phi() const
Definition: Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/Track.h:47
TrackAnalysis::event
const TIDA::Event * event() const
Definition: TrackAnalysis.h:135
Efficiency1D::Fill
virtual void Fill(double x, double w=1)
Definition: Efficiency1D.h:37
TrackAnalysis::m_histos
std::map< std::string, TH1 * > m_histos
lookup table for the histograms by name - does this need to be in the base class?
Definition: TrackAnalysis.h:149
SigAnalysis::execute
virtual void execute(const std::vector< TIDA::Track * > &reftracks, const std::vector< TIDA::Track * > &testtracks, TrackAssociator *matcher)
Definition: SigAnalysis.cxx:205
MakeTH3DFromTH2Ds.zBins
list zBins
Definition: MakeTH3DFromTH2Ds.py:86
TrackAnalysis::m_yBeamReference
double m_yBeamReference
Definition: TrackAnalysis.h:154
SigAnalysis::m_eff_eta
Efficiency1D * m_eff_eta
Definition: SigAnalysis.h:89
beamspotman.dir
string dir
Definition: beamspotman.py:623
TrackAnalysis::m_event
TIDA::Event * m_event
Definition: TrackAnalysis.h:164
SigAnalysis::m_eff_ptp
Efficiency1D * m_eff_ptp
Definition: SigAnalysis.h:86
SigAnalysis::m_purity_a0
Efficiency1D * m_purity_a0
Definition: SigAnalysis.h:100
TrackAnalysis::m_zBeamTest
double m_zBeamTest
Definition: TrackAnalysis.h:160
SigAnalysis::m_purity_eta
Efficiency1D * m_purity_eta
Definition: SigAnalysis.h:96
TrackAnalysis::m_store
TIDA::FeatureStore m_store
Definition: TrackAnalysis.h:162
TrackAnalysis::TrackAnalysis
TrackAnalysis(const std::string &name)
the beam test parts are not really usable in a multithreaded environment
Definition: TrackAnalysis.h:38
SigAnalysis::m_histos
std::map< std::string, TH1F * > m_histos
Definition: SigAnalysis.h:83
SigAnalysis::m_h2
TH2F * m_h2
Definition: SigAnalysis.h:103
TIDA::Track::eta
double eta() const
Definition: Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/Track.h:46
TIDDirectory::push
void push()
Definition: TIDDirectory.h:78
SigAnalysis::m_icount
int m_icount
Event counter.
Definition: SigAnalysis.h:114
SigAnalysis::m_eff_ptm
Efficiency1D * m_eff_ptm
Definition: SigAnalysis.h:87
python.PyAthena.v
v
Definition: PyAthena.py:154
y
#define y
h
SigAnalysis::m_purity_z0
Efficiency1D * m_purity_z0
Definition: SigAnalysis.h:98
SigAnalysis::m_h2r
TH2F * m_h2r
Definition: SigAnalysis.h:105
SigAnalysis::m_eff_a0
Efficiency1D * m_eff_a0
Definition: SigAnalysis.h:93
TrackAnalysis::m_xBeamTest
double m_xBeamTest
test sample
Definition: TrackAnalysis.h:158
SigAnalysis::addHistogram
void addHistogram(TH1F *h)
Definition: SigAnalysis.h:70
SigAnalysis::m_eff_d0
Efficiency1D * m_eff_d0
Definition: SigAnalysis.h:92
TIDA::Track
Definition: Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/Track.h:26
TrackAnalysis::m_yBeamTest
double m_yBeamTest
Definition: TrackAnalysis.h:159
python.TrigEgammaMonitorHelper.TH1F
def TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
Definition: TrigEgammaMonitorHelper.py:24
SigAnalysis::m_Nreco
int m_Nreco
number of reconstructed tracks
Definition: SigAnalysis.h:108
SigAnalysis::m_dir
TIDDirectory * m_dir
Definition: SigAnalysis.h:116
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
TrackAnalysis::m_zBeamReference
double m_zBeamReference
Definition: TrackAnalysis.h:155
SigAnalysis::m_h2m
TH2F * m_h2m
Definition: SigAnalysis.h:104
TrackAnalysis::m_name
std::string m_name
identifier of the of the analysis - also used for the root directory into which the histograms are pu...
Definition: TrackAnalysis.h:145
SigAnalysis::m_eff_pt
Efficiency1D * m_eff_pt
Definition: SigAnalysis.h:85
TrackAnalysis::m_xBeamReference
double m_xBeamReference
beamline positions reference sample
Definition: TrackAnalysis.h:153
Generated on Thu Nov 7 2024 21:56:26 for ATLAS Offline Software by   doxygen 1.8.18