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

#include <CheckFlow_New.h>

Inheritance diagram for CheckFlow_New:
Collaboration diagram for CheckFlow_New:

Public Member Functions

 CheckFlow_New (const std::string &name, ISvcLocator *pSvcLocator)
 
virtual StatusCode initialize () override
 
virtual StatusCode execute () override
 
virtual StatusCode finalize () override
 
virtual StatusCode sysInitialize () override
 Override sysInitialize. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. More...
 
ServiceHandle< StoreGateSvc > & evtStore ()
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & evtStore () const
 The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
 
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
 
virtual StatusCode sysStart () override
 Handle START transition. More...
 
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles. More...
 
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles. More...
 
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
 Declare a new Gaudi property. More...
 
Gaudi::Details::PropertyBase * declareProperty (const std::string &name, T &property, const std::string &doc="none")
 Declare a new Gaudi property. More...
 
void updateVHKA (Gaudi::Details::PropertyBase &)
 
MsgStream & msg () const
 
MsgStream & msg (const MSG::Level lvl) const
 
bool msgLvl (const MSG::Level lvl) const
 

Public Attributes

 InputCollections
 
 job = AlgSequence()
 
 EvtMax
 
 OutputLevel
 
 Output
 
 CheckFlow_New = job.CheckFlow_New
 
 McEventKey
 
 RapidityCutMax
 

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Private Types

enum  { n_ptbin =16, n_etabin =8 }
 
typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

StringProperty m_key {this, "McEventKey", "FLOW_EVENT"}
 
BooleanProperty m_produceHistogram {this, "HistogramFlag", true}
 
DoubleProperty m_bcut_min {this, "ImpactCutMin", 0.}
 
DoubleProperty m_bcut_max {this, "ImpactCutMax", 99.}
 
DoubleProperty m_ptcut_min {this, "PtCutMin", 0.}
 
DoubleProperty m_ptcut_max {this, "PtCutMax", 999999.}
 
DoubleProperty m_rapcut_min {this, "RapidityCutMin", 0.}
 
DoubleProperty m_rapcut_max {this, "RapidityCutMax", 5.5}
 
TH1D * m_hist_Psi_n_true [6] {}
 
TH1D * m_hist_Psi_n_reco [6] {}
 
TH1D * m_hist_psi_corr_true [36] {}
 
TH1D * m_hist_psi_corr_reco [36] {}
 
TH1D * m_hist_Psi_n_ebe [6] {}
 
TH1D * m_hist_Psi_n_ebe_pt [6] {}
 
TH1D * m_hist_vn_ebe [6] {}
 
TProfile * m_profile_pt_dep [6][n_etabin] {}
 
TProfile * m_profile_eta_dep [6][n_ptbin] {}
 
TProfile * m_profile_pt_dep_reco [6][n_etabin] {}
 
TProfile * m_profile_eta_dep_reco [6][n_ptbin] {}
 
TProfile * m_profile_resolution {}
 
SG::ReadHandleKey< HijingEventParams > m_hijingKey {this, "HijingEventParmsKey","Hijing_event_params"}
 
TruthHelper::GenAccessIO * m_tesIO {}
 
DataObjIDColl m_extendedExtraObjects
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 30 of file CheckFlow_New.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
private
Enumerator
n_ptbin 
n_etabin 

Definition at line 49 of file CheckFlow_New.h.

49  {
50  n_ptbin=16,
51  n_etabin=8
52  };

Constructor & Destructor Documentation

◆ CheckFlow_New()

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

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleKeyArrayType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleKeyType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  hndl,
const SG::VarHandleType &   
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > &  t,
const SG::NotHandleType &   
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleBase &  hndl,
const std::string &  doc,
const SG::VarHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKey &  hndl,
const std::string &  doc,
const SG::VarHandleKeyType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray &  hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType &   
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType &   
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
148  }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode CheckFlow_New::execute ( )
overridevirtual

Definition at line 200 of file CheckFlow_New.cxx.

200  {
201  msg(MSG::INFO) << ">>> CheckFlow_New from execute" << endmsg;
202 
203  float pt_binvals[n_ptbin+1]={0.0,0.25,0.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,6.0,8.0,12.0,16.0,20.0,30.0,40.0};
204  float eta_bin_max = 4.0;
205 
206  SG::ReadHandle<HijingEventParams> hijing_pars{m_hijingKey};
207  float b = hijing_pars->get_b();
208  float Psi_n[6],Psi_n_reco[6];
209  float Psi_n_reco_pos[6],Psi_n_reco_neg[6];
210  for(int ihar=0;ihar<6;ihar++){Psi_n[ihar]=hijing_pars->get_psi(ihar+1);}
211  msg(MSG::INFO)<<"SOUMYA "<<hijing_pars->get_psi(1)<<" "<<hijing_pars->get_psi(2)<<" "<<hijing_pars->get_psi(3)
212  <<hijing_pars->get_psi(4)<<" "<<hijing_pars->get_psi(5)<<" "<<hijing_pars->get_psi(6)<<" "<<b << endmsg;
213 
214 
215  // Check cut on impact parameter b
216  if(b<m_bcut_min || b>m_bcut_max) return StatusCode::SUCCESS;
217 
218 
219  double ngenerated_pos = 0,ngenerated_pt_pos=0;
220  double cos_n_pos[6],sin_n_pos[6],cos_n_pt_pos[6],sin_n_pt_pos[6];
221  double ngenerated_neg = 0,ngenerated_pt_neg=0;
222  double cos_n_neg[6],sin_n_neg[6],cos_n_pt_neg[6],sin_n_pt_neg[6];
223  for(int ihar=0;ihar<6;ihar++){cos_n_pos[ihar]=0;sin_n_pos[ihar]=0; cos_n_pt_pos[ihar]=0;sin_n_pt_pos[ihar]=0;
224  cos_n_neg[ihar]=0;sin_n_neg[ihar]=0; cos_n_pt_neg[ihar]=0;sin_n_pt_neg[ihar]=0;}
225 
226  // Iterate over all MC particles
227  std::vector<HepMC::ConstGenParticlePtr> particles;
228  StatusCode stat = m_tesIO->getMC(particles, false, m_key);
229  if (stat.isFailure()) {
230  msg(MSG::ERROR) << "Could not find " << m_key << endmsg;
231  return stat;
232  }
233  for (auto pitr: particles) {
234  int pid = pitr->pdg_id();
235  double pt = pitr->momentum().perp();
236  double rapid = pitr->momentum().pseudoRapidity();
237  double phi = pitr->momentum().phi();
238  msg(MSG::DEBUG)
239  << " PID = " << pid << " Status = " << pitr->status()
240  << " Eta = " << rapid << " Phi = " << phi<< endmsg;
241 
242  if( (std::abs(rapid) >= m_rapcut_min) && (std::abs(rapid) <= m_rapcut_max) &&
243  (std::abs(pt) >= m_ptcut_min) && (std::abs(pt) <= m_ptcut_max) ) {
244 
245  for(int ihar=0;ihar<6;ihar++){
246  float temp=(ihar+1)*(phi-Psi_n[ihar]);
247 
248  int ieta= (int)(std::abs(rapid)*n_etabin/eta_bin_max);
249  if(ieta>=0 && ieta<n_etabin) m_profile_pt_dep [ihar][ieta]->Fill(pt/1000,cos(temp));
250 
251 
252  float temp_pt=pt/1000;
253  for(int ipt=0;ipt<n_ptbin;ipt++){
254  if(temp_pt<pt_binvals[ipt+1]){
255  m_profile_eta_dep[ihar][ipt]->Fill(rapid ,cos(temp));
256  break;
257  }
258  }
259 
260  if( rapid >3.2 && rapid< 4.9){
261  cos_n_pos[ihar]+=std::cos( (ihar+1)*phi);
262  sin_n_pos[ihar]+=std::sin( (ihar+1)*phi);
263  ngenerated_pos++;
264 
265  cos_n_pt_pos[ihar]+=pt*std::cos( (ihar+1)*phi);
266  sin_n_pt_pos[ihar]+=pt*std::sin( (ihar+1)*phi);
267  ngenerated_pt_pos +=pt;
268  }
269  if( rapid <-3.2 && rapid >-4.9){
270  cos_n_neg[ihar]+=std::cos( (ihar+1)*phi);
271  sin_n_neg[ihar]+=std::sin( (ihar+1)*phi);
272  ngenerated_neg++;
273 
274  cos_n_pt_neg[ihar]+=pt*std::cos( (ihar+1)*phi);
275  sin_n_pt_neg[ihar]+=pt*std::sin( (ihar+1)*phi);
276  ngenerated_pt_neg +=pt;
277  }
278 
279  }
280  }
281  }
282 
283 
284 // Calculate the event by event vn and also the reconstructed Psi_n angles
285 // Also make correlation histos between Psi_n_truth and Psi_n_reco
286  float cos_n[6],sin_n[6],cos_n_pt[6],sin_n_pt[6];
287  for(int ihar=0;ihar<6;ihar++){
288  cos_n[ihar] = ( cos_n_pos[ihar]+ cos_n_neg[ihar] ) / (ngenerated_pos+ngenerated_neg);
289  sin_n[ihar] = ( sin_n_pos[ihar]+ sin_n_neg[ihar] ) / (ngenerated_pos+ngenerated_neg);
290 
291  float psi_reco=std::atan2(sin_n[ihar],cos_n[ihar])/(ihar+1);
292  m_hist_Psi_n_ebe[ihar]->Fill( (ihar+1)*(psi_reco-Psi_n[ihar]) );
293  m_hist_vn_ebe [ihar]->Fill(std::sqrt(cos_n[ihar]*cos_n[ihar] +sin_n[ihar]*sin_n[ihar] ));
294 
295  Psi_n_reco_pos[ihar]=std::atan2(sin_n_pos[ihar],cos_n_pos[ihar])/ (ihar+1);
296  Psi_n_reco_neg[ihar]=std::atan2(sin_n_neg[ihar],cos_n_neg[ihar])/ (ihar+1);
297  Psi_n_reco [ihar]=psi_reco;
298 
299 
300  cos_n_pt[ihar] = ( cos_n_pt_pos[ihar]+ cos_n_pt_neg[ihar] ) / (ngenerated_pt_pos+ngenerated_pt_neg);
301  sin_n_pt[ihar] = ( sin_n_pt_pos[ihar]+ sin_n_pt_neg[ihar] ) / (ngenerated_pt_pos+ngenerated_pt_neg);
302 
303  psi_reco=std::atan2(sin_n_pt[ihar],cos_n_pt[ihar])/(ihar+1);
304  m_hist_Psi_n_ebe_pt[ihar]->Fill( (ihar+1)*(psi_reco-Psi_n[ihar]) );
305  }
306 
307 
308 // Make the plots for the correlation between Psi_n truth (for different n) (same for Psi_n reco)
309  for(int ihar=0;ihar<6;ihar++){
310  m_hist_Psi_n_true[ihar]->Fill((ihar+1)*Psi_n[ihar]);
311  m_hist_Psi_n_reco[ihar]->Fill((ihar+1)*Psi_n_reco[ihar]);
312 
313  float psi1,psi2;
314  for(int ihar2=0;ihar2<6;ihar2++){
315  psi1=(ihar+1)*Psi_n[ihar];psi2=(ihar2+1)*Psi_n[ihar2];
316  m_hist_psi_corr_true[ihar*6+ihar2]->Fill( std::atan2( std::sin(psi1-psi2),std::cos(psi1-psi2) ) );
317 
318  psi1=(ihar+1)*Psi_n_reco[ihar];psi2=(ihar2+1)*Psi_n_reco[ihar2];
319  m_hist_psi_corr_reco[ihar*6+ihar2]->Fill( std::atan2( std::sin(psi1-psi2),std::cos(psi1-psi2) ) );
320  }
321  }
322 
323 
324 
325 
326 
327 // calculate the pt and eta dependence using the Psi_reco angles also fill the resolution TProfile
328  for(int ihar=0;ihar<6;ihar++){
329  m_profile_resolution->Fill( ihar+1, cos( (ihar+1) * (Psi_n_reco_pos[ihar] - Psi_n_reco_neg[ihar]) ) );
330  }
331  for (auto pitr: particles) {
332  double pt = pitr->momentum().perp();
333  double rapid = pitr->momentum().pseudoRapidity();
334  double phi = pitr->momentum().phi();
335  if( (std::abs(rapid) >= m_rapcut_min) && (std::abs(rapid) <= m_rapcut_max) &&
336  (std::abs(pt) >= m_ptcut_min) && (std::abs(pt) <= m_ptcut_max) ) {
337 
338  for(int ihar=0;ihar<6;ihar++){
339  float temp=(ihar+1)*(phi-Psi_n_reco_pos[ihar]);
340  if(rapid>0) temp=(ihar+1)*(phi-Psi_n_reco_neg[ihar]);
341 
342 
343  int ieta= (int)(std::abs(rapid)*n_etabin/eta_bin_max);
344  if(ieta>=0 && ieta<n_etabin) m_profile_pt_dep_reco [ihar][ieta]->Fill(pt/1000,std::cos(temp));
345 
346  float temp_pt=pt/1000;
347  for(int ipt=0;ipt<n_ptbin;ipt++){
348  if(temp_pt<pt_binvals[ipt+1]){
349  m_profile_eta_dep_reco[ihar][ipt]->Fill(rapid ,cos(temp));
350  break;
351  }
352  }
353  }
354  }
355  }
356 
357 
358 
359 
360 
361  return StatusCode::SUCCESS;
362 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

51 {
52  // If we didn't find any symlinks to add, just return the collection
53  // from the base class. Otherwise, return the extended collection.
54  if (!m_extendedExtraObjects.empty()) {
55  return m_extendedExtraObjects;
56  }
57  return Algorithm::extraOutputDeps();
58 }

◆ finalize()

StatusCode CheckFlow_New::finalize ( )
overridevirtual

Definition at line 364 of file CheckFlow_New.cxx.

364  {
365  msg(MSG::INFO) << ">>> CheckFlow_New from finalize" << endmsg;
366 
367  return StatusCode::SUCCESS;
368 }

◆ initialize()

StatusCode CheckFlow_New::initialize ( )
overridevirtual

Definition at line 41 of file CheckFlow_New.cxx.

41  {
42  ATH_CHECK(m_hijingKey.initialize());
43 
44  ATH_MSG_INFO(">>> CheckFlow_New from Initialize");
45 
46  float pt_binvals[n_ptbin+1]={0.0,0.25,0.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,6.0,8.0,12.0,16.0,20.0,30.0,40.0};
47  float eta_bin_max = 4.0;
48 
49 
50  SmartIF<ITHistSvc> rootHistSvc{Gaudi::svcLocator()->service("THistSvc")};
51  if (!rootHistSvc) {
52  ATH_MSG_ERROR( "Unable to locate THistSvc" );
53  return StatusCode::FAILURE;
54  }
55 
56  std::string StreamAndPath="/FlowOutPut/";
57  std::string histPath = StreamAndPath;
58  char name[100],name1[100];
59  for (int ihar=0;ihar<6;ihar++){
60 
61 // /*
62  sprintf(name,"hist_Psi_%d_true",ihar+1);
63  sprintf(name1,"Truth Psi_{%d} distribution;%dPsi_{%d} Truth;events",ihar+1,ihar+1,ihar+1);
64  m_hist_Psi_n_true [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
65  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_true[ihar]->GetName(),m_hist_Psi_n_true[ihar]).isFailure() ){
66  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_true[ihar]->GetName() << endmsg;
67  }
68  m_hist_Psi_n_true[ihar]->GetXaxis()->CenterTitle();
69  m_hist_Psi_n_true[ihar]->GetYaxis()->CenterTitle();
70 
71  sprintf(name,"hist_Psi_%d_reco",ihar+1);
72  sprintf(name1,"Reconstructed Psi_{%d} distribution;%dPsi_{%d} Reco;events",ihar+1,ihar+1,ihar+1);
73  m_hist_Psi_n_reco [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
74  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_reco[ihar]->GetName(),m_hist_Psi_n_reco[ihar]).isFailure() ){
75  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_reco[ihar]->GetName() << endmsg;
76  }
77  m_hist_Psi_n_reco[ihar]->GetXaxis()->CenterTitle();
78  m_hist_Psi_n_reco[ihar]->GetYaxis()->CenterTitle();
79 
80 
81  for (int ihar2=0;ihar2<6;ihar2++)
82  {
83  int ihar_i=ihar*6+ihar2;
84 
85  sprintf(name,"hist_Psi_corr_true_%d_%d",ihar+1,ihar2+1);
86  sprintf(name1,"true Psi_{%d} -Psi_{%d};%dPsi_{%d} -%dPsi_{%d} ;events",ihar+1,ihar2+1,ihar+1,ihar+1,ihar2+1,ihar2+1);
87  m_hist_psi_corr_true [ihar_i]=new TH1D (name,name1,1000,-2*M_PI,2*M_PI);
88  if ( rootHistSvc->regHist(histPath+m_hist_psi_corr_true[ihar_i]->GetName(),m_hist_psi_corr_true[ihar_i]).isFailure() ){
89  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_psi_corr_true[ihar_i]->GetName() << endmsg;
90  }
91  m_hist_psi_corr_true[ihar_i]->GetXaxis()->CenterTitle();
92  m_hist_psi_corr_true[ihar_i]->GetYaxis()->CenterTitle();
93 
94 
95  sprintf(name,"hist_Psi_corr_reco_%d_%d",ihar+1,ihar2+1);
96  sprintf(name1,"reco Psi_{%d} -Psi_{%d};%dPsi_{%d} -%dPsi_{%d} ;events",ihar+1,ihar2+1,ihar+1,ihar+1,ihar2+1,ihar2+1);
97  m_hist_psi_corr_reco [ihar_i]=new TH1D (name,name1,1000,-2*M_PI,2*M_PI);
98  if ( rootHistSvc->regHist(histPath+m_hist_psi_corr_reco[ihar_i]->GetName(),m_hist_psi_corr_reco[ihar_i]).isFailure() ){
99  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_psi_corr_reco[ihar_i]->GetName() << endmsg;
100  }
101  m_hist_psi_corr_reco[ihar_i]->GetXaxis()->CenterTitle();
102  m_hist_psi_corr_reco[ihar_i]->GetYaxis()->CenterTitle();
103  }
104 
105 
106 
107  //integrated vn event by event
108  sprintf(name,"hist_v%d_ebe",ihar+1);
109  sprintf(name1,"v%d;v%d;events",ihar+1,ihar+1);
110  m_hist_vn_ebe [ihar]=new TH1D (name,name1,1000,-0.5,0.5);
111  if ( rootHistSvc->regHist(histPath+m_hist_vn_ebe[ihar]->GetName(),m_hist_vn_ebe[ihar]).isFailure() ){
112  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_vn_ebe[ihar]->GetName() << endmsg;
113  }
114  m_hist_vn_ebe[ihar]->GetXaxis()->CenterTitle();
115  m_hist_vn_ebe[ihar]->GetYaxis()->CenterTitle();
116 
117  sprintf(name ,"hist_Psi%d_ebe",ihar+1);
118  sprintf(name1,"%d#Delta#Psi;%d(#Psi_{reco}-#Psi_{Truth});events",ihar+1,ihar+1);
119  m_hist_Psi_n_ebe [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
120  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_ebe[ihar]->GetName(),m_hist_Psi_n_ebe[ihar]).isFailure() ){
121  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_ebe[ihar]->GetName() << endmsg;
122  }
123  m_hist_Psi_n_ebe[ihar]->GetXaxis()->CenterTitle();
124  m_hist_Psi_n_ebe[ihar]->GetYaxis()->CenterTitle();
125 
126  sprintf(name ,"hist_Psi%d_ebe_pt",ihar+1);
127  sprintf(name1,"%d#Delta#Psi (pT weighted);%d(#Psi_{reco}-#Psi_{Truth});events",ihar+1,ihar+1);
128  m_hist_Psi_n_ebe_pt [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
129  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_ebe_pt[ihar]->GetName(),m_hist_Psi_n_ebe_pt[ihar]).isFailure() ){
130  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_ebe_pt[ihar]->GetName() << endmsg;
131  }
132  m_hist_Psi_n_ebe_pt[ihar]->GetXaxis()->CenterTitle();
133  m_hist_Psi_n_ebe_pt[ihar]->GetYaxis()->CenterTitle();
134 
135 
136 
137 
138 
139 
140  for(int ieta=0;ieta<n_etabin;ieta++){
141  sprintf(name ,"profile_pt_dep_%d_eta%d" ,ihar+1,ieta);
142  sprintf(name1,"v%d vs pT (eta%d);pT;v%d",ihar+1,ieta,ihar+1);
143  m_profile_pt_dep [ihar][ieta]=new TProfile (name,name1,n_ptbin,pt_binvals);
144  if ( rootHistSvc->regHist(histPath+m_profile_pt_dep[ihar][ieta]->GetName(),m_profile_pt_dep[ihar][ieta]).isFailure() ){
145  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_pt_dep[ihar][ieta]->GetName() << endmsg;
146  }
147  m_profile_pt_dep [ihar][ieta]->GetXaxis()->CenterTitle();
148  m_profile_pt_dep [ihar][ieta]->GetYaxis()->CenterTitle();
149  }
150 
151  for(int ipt=0;ipt<n_ptbin;ipt++){
152  sprintf(name ,"profile_eta_dep_%d_pt%d",ihar+1,ipt);
153  sprintf(name1,"v%d vs #eta; (ipt%d)#eta;v%d",ihar+1,ipt,ihar+1);
154  m_profile_eta_dep [ihar][ipt]=new TProfile (name,name1,2*n_etabin, -eta_bin_max,eta_bin_max);
155  if ( rootHistSvc->regHist(histPath+m_profile_eta_dep[ihar][ipt]->GetName(),m_profile_eta_dep[ihar][ipt]).isFailure() ){
156  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_eta_dep[ihar][ipt]->GetName() << endmsg;
157  }
158  m_profile_eta_dep [ihar][ipt]->GetXaxis()->CenterTitle();
159  m_profile_eta_dep [ihar][ipt]->GetYaxis()->CenterTitle();
160  }
161 
162 
163  for(int ieta=0;ieta<n_etabin;ieta++){
164  sprintf(name ,"profile_pt_dep_reco_%d_eta%d",ihar+1,ieta);
165  sprintf(name1,"v%d vs pT (eta%d);pT;v%d",ihar+1,ieta,ihar+1);
166  m_profile_pt_dep_reco [ihar][ieta]=new TProfile (name,name1,n_ptbin,pt_binvals);
167  if ( rootHistSvc->regHist(histPath+m_profile_pt_dep_reco[ihar][ieta]->GetName(),m_profile_pt_dep_reco[ihar][ieta]).isFailure() ){
168  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_pt_dep_reco[ihar][ieta]->GetName() << endmsg;
169  }
170  m_profile_pt_dep_reco [ihar][ieta]->GetXaxis()->CenterTitle();
171  m_profile_pt_dep_reco [ihar][ieta]->GetYaxis()->CenterTitle();
172  }
173 
174 
175  for(int ipt=0;ipt<n_ptbin;ipt++){
176  sprintf(name ,"profile_eta_dep_reco_%d_pt%d",ihar+1,ipt);
177  sprintf(name1,"v%d vs #eta (pt%d);#eta;v%d",ihar+1,ipt,ihar+1);
178  m_profile_eta_dep_reco [ihar][ipt]=new TProfile (name,name1,2*n_etabin, -eta_bin_max,eta_bin_max);
179  if ( rootHistSvc->regHist(histPath+m_profile_eta_dep_reco[ihar][ipt]->GetName(),m_profile_eta_dep_reco[ihar][ipt]).isFailure() ){
180  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_eta_dep_reco[ihar][ipt]->GetName() << endmsg;
181  }
182  m_profile_eta_dep_reco [ihar][ipt]->GetXaxis()->CenterTitle();
183  m_profile_eta_dep_reco [ihar][ipt]->GetYaxis()->CenterTitle();
184  }
185 
186 
187  }
188  m_profile_resolution=new TProfile("profile_resolution","vn resolution;n;resolution",6, 0.5,6.5);
189  if(rootHistSvc->regHist(histPath+m_profile_resolution->GetName(),m_profile_resolution).isFailure() ){
190  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_resolution->GetName() << endmsg;
191  }
192 
193  msg(MSG::DEBUG) << "Histograms have been booked " << endmsg;
194  m_tesIO = new TruthHelper::GenAccessIO();
195  return StatusCode::SUCCESS;
196 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ msg() [1/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level  lvl) const
inlineinherited

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

bool AthCommonMsg< Algorithm >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

66  {
67  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
68 
69  if (sc.isFailure()) {
70  return sc;
71  }
72  ServiceHandle<ICondSvc> cs("CondSvc",name());
73  for (auto h : outputHandles()) {
74  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
75  // do this inside the loop so we don't create the CondSvc until needed
76  if ( cs.retrieve().isFailure() ) {
77  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
78  return StatusCode::SUCCESS;
79  }
80  if (cs->regHandle(this,*h).isFailure()) {
81  sc = StatusCode::FAILURE;
82  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
83  << " with CondSvc");
84  }
85  }
86  }
87  return sc;
88 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ CheckFlow_New

CheckFlow_New.CheckFlow_New = job.CheckFlow_New
Todo:
Migrate to a GenBase class or at least AthAlgorithm

Definition at line 54 of file CheckFlow_New.py.

◆ EvtMax

CheckFlow_New.EvtMax

Definition at line 38 of file CheckFlow_New.py.

◆ InputCollections

CheckFlow_New.InputCollections

Definition at line 16 of file CheckFlow_New.py.

◆ job

CheckFlow_New.job = AlgSequence()

Definition at line 30 of file CheckFlow_New.py.

◆ m_bcut_max

DoubleProperty CheckFlow_New::m_bcut_max {this, "ImpactCutMax", 99.}
private

Definition at line 43 of file CheckFlow_New.h.

◆ m_bcut_min

DoubleProperty CheckFlow_New::m_bcut_min {this, "ImpactCutMin", 0.}
private

Definition at line 42 of file CheckFlow_New.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_hijingKey

SG::ReadHandleKey<HijingEventParams> CheckFlow_New::m_hijingKey {this, "HijingEventParmsKey","Hijing_event_params"}
private

Definition at line 71 of file CheckFlow_New.h.

◆ m_hist_psi_corr_reco

TH1D* CheckFlow_New::m_hist_psi_corr_reco[36] {}
private

Definition at line 58 of file CheckFlow_New.h.

◆ m_hist_psi_corr_true

TH1D* CheckFlow_New::m_hist_psi_corr_true[36] {}
private

Definition at line 57 of file CheckFlow_New.h.

◆ m_hist_Psi_n_ebe

TH1D* CheckFlow_New::m_hist_Psi_n_ebe[6] {}
private

Definition at line 60 of file CheckFlow_New.h.

◆ m_hist_Psi_n_ebe_pt

TH1D* CheckFlow_New::m_hist_Psi_n_ebe_pt[6] {}
private

Definition at line 61 of file CheckFlow_New.h.

◆ m_hist_Psi_n_reco

TH1D* CheckFlow_New::m_hist_Psi_n_reco[6] {}
private

Definition at line 56 of file CheckFlow_New.h.

◆ m_hist_Psi_n_true

TH1D* CheckFlow_New::m_hist_Psi_n_true[6] {}
private

Definition at line 55 of file CheckFlow_New.h.

◆ m_hist_vn_ebe

TH1D* CheckFlow_New::m_hist_vn_ebe[6] {}
private

Definition at line 62 of file CheckFlow_New.h.

◆ m_key

StringProperty CheckFlow_New::m_key {this, "McEventKey", "FLOW_EVENT"}
private

Definition at line 39 of file CheckFlow_New.h.

◆ m_produceHistogram

BooleanProperty CheckFlow_New::m_produceHistogram {this, "HistogramFlag", true}
private

Definition at line 40 of file CheckFlow_New.h.

◆ m_profile_eta_dep

TProfile* CheckFlow_New::m_profile_eta_dep[6][n_ptbin] {}
private

Definition at line 65 of file CheckFlow_New.h.

◆ m_profile_eta_dep_reco

TProfile* CheckFlow_New::m_profile_eta_dep_reco[6][n_ptbin] {}
private

Definition at line 67 of file CheckFlow_New.h.

◆ m_profile_pt_dep

TProfile* CheckFlow_New::m_profile_pt_dep[6][n_etabin] {}
private

Definition at line 64 of file CheckFlow_New.h.

◆ m_profile_pt_dep_reco

TProfile* CheckFlow_New::m_profile_pt_dep_reco[6][n_etabin] {}
private

Definition at line 66 of file CheckFlow_New.h.

◆ m_profile_resolution

TProfile* CheckFlow_New::m_profile_resolution {}
private

Definition at line 69 of file CheckFlow_New.h.

◆ m_ptcut_max

DoubleProperty CheckFlow_New::m_ptcut_max {this, "PtCutMax", 999999.}
private

Definition at line 45 of file CheckFlow_New.h.

◆ m_ptcut_min

DoubleProperty CheckFlow_New::m_ptcut_min {this, "PtCutMin", 0.}
private

Definition at line 44 of file CheckFlow_New.h.

◆ m_rapcut_max

DoubleProperty CheckFlow_New::m_rapcut_max {this, "RapidityCutMax", 5.5}
private

Definition at line 47 of file CheckFlow_New.h.

◆ m_rapcut_min

DoubleProperty CheckFlow_New::m_rapcut_min {this, "RapidityCutMin", 0.}
private

Definition at line 46 of file CheckFlow_New.h.

◆ m_tesIO

TruthHelper::GenAccessIO* CheckFlow_New::m_tesIO {}
private

Definition at line 72 of file CheckFlow_New.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ McEventKey

CheckFlow_New.McEventKey

Definition at line 55 of file CheckFlow_New.py.

◆ Output

CheckFlow_New.Output

Definition at line 47 of file CheckFlow_New.py.

◆ OutputLevel

CheckFlow_New.OutputLevel

Definition at line 39 of file CheckFlow_New.py.

◆ RapidityCutMax

CheckFlow_New.RapidityCutMax

Definition at line 56 of file CheckFlow_New.py.

The documentation for this class was generated from the following files:
CheckFlow_New::m_profile_eta_dep_reco
TProfile * m_profile_eta_dep_reco[6][n_ptbin]
Definition: CheckFlow_New.h:67
CheckFlow_New::m_tesIO
TruthHelper::GenAccessIO * m_tesIO
Definition: CheckFlow_New.h:72
ReadCellNoiseFromCool.name1
name1
Definition: ReadCellNoiseFromCool.py:233
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
CheckFlow_New::m_hist_psi_corr_reco
TH1D * m_hist_psi_corr_reco[36]
Definition: CheckFlow_New.h:58
CheckFlow_New::m_ptcut_max
DoubleProperty m_ptcut_max
Definition: CheckFlow_New.h:45
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:67
CheckFlow_New::m_ptcut_min
DoubleProperty m_ptcut_min
Definition: CheckFlow_New.h:44
test_pyathena.pt
pt
Definition: test_pyathena.py:11
M_PI
#define M_PI
Definition: ActiveFraction.h:11
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
CheckFlow_New::m_hijingKey
SG::ReadHandleKey< HijingEventParams > m_hijingKey
Definition: CheckFlow_New.h:71
CheckFlow_New::m_profile_resolution
TProfile * m_profile_resolution
Definition: CheckFlow_New.h:69
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
CheckFlow_New::m_hist_Psi_n_ebe
TH1D * m_hist_Psi_n_ebe[6]
Definition: CheckFlow_New.h:60
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
CheckFlow_New::m_hist_Psi_n_true
TH1D * m_hist_Psi_n_true[6]
Definition: CheckFlow_New.h:55
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
xAOD::phi
setEt phi
Definition: TrigEMCluster_v1.cxx:29
CheckFlow_New::m_profile_pt_dep_reco
TProfile * m_profile_pt_dep_reco[6][n_etabin]
Definition: CheckFlow_New.h:66
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
CheckFlow_New::m_hist_Psi_n_ebe_pt
TH1D * m_hist_Psi_n_ebe_pt[6]
Definition: CheckFlow_New.h:61
CheckFlow_New::m_profile_pt_dep
TProfile * m_profile_pt_dep[6][n_etabin]
Definition: CheckFlow_New.h:64
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
AthCommonDataStore
Definition: AthCommonDataStore.h:52
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthAlgorithm.cxx:66
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
python.TrigEgammaMonitorHelper.TProfile
def TProfile(*args, **kwargs)
Definition: TrigEgammaMonitorHelper.py:81
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
CheckFlow_New::n_ptbin
@ n_ptbin
Definition: CheckFlow_New.h:50
ParticleGun_EoverP_Config.pid
pid
Definition: ParticleGun_EoverP_Config.py:62
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
TruthHelper::GenAccessIO::getMC
StatusCode getMC(MCParticleCollection &mcParticles, const bool ifgen=false, const std::string &key="GEN_EVENT") const
Definition: Generators/FlowAfterburner/FlowAfterburner/GenAccessIO.h:29
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
beamspotman.stat
stat
Definition: beamspotman.py:264
CheckFlow_New::m_key
StringProperty m_key
Definition: CheckFlow_New.h:39
CheckFlow_New::m_profile_eta_dep
TProfile * m_profile_eta_dep[6][n_ptbin]
Definition: CheckFlow_New.h:65
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
CheckFlow_New::m_hist_vn_ebe
TH1D * m_hist_vn_ebe[6]
Definition: CheckFlow_New.h:62
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
CheckFlow_New::m_hist_psi_corr_true
TH1D * m_hist_psi_corr_true[36]
Definition: CheckFlow_New.h:57
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:76
TruthHelper::GenAccessIO
Definition: Generators/FlowAfterburner/FlowAfterburner/GenAccessIO.h:21
CheckFlow_New::m_hist_Psi_n_reco
TH1D * m_hist_Psi_n_reco[6]
Definition: CheckFlow_New.h:56
CheckFlow_New::m_bcut_max
DoubleProperty m_bcut_max
Definition: CheckFlow_New.h:43
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
a
TList * a
Definition: liststreamerinfos.cxx:10
h
python.CaloAddPedShiftConfig.int
int
Definition: CaloAddPedShiftConfig.py:45
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
python.Constants.INFO
int INFO
Definition: Control/AthenaCommon/python/Constants.py:15
DEBUG
#define DEBUG
Definition: page_access.h:11
AthCommonMsg< Algorithm >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
LArG4FSStartPointFilter.particles
list particles
Definition: LArG4FSStartPointFilter.py:84
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:629
CheckFlow_New::n_etabin
@ n_etabin
Definition: CheckFlow_New.h:51
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:801
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
CheckFlow_New::m_rapcut_min
DoubleProperty m_rapcut_min
Definition: CheckFlow_New.h:46
fitman.k
k
Definition: fitman.py:528
ServiceHandle< ICondSvc >
CheckFlow_New::m_rapcut_max
DoubleProperty m_rapcut_max
Definition: CheckFlow_New.h:47
Generated on Mon Sep 1 2025 21:30:59 for ATLAS Offline Software by   doxygen 1.8.18