ATLAS Offline Software
CheckFlow_New.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 // File: Generators/FlowAfterburnber/CheckFlow_New.h
6 // Description:
7 // This is a simple algorithm to histogram particle properties
8 // for diagnosing of flow generation
9 //
10 //************************** THIS PROGRAM ANALYZES FILES WITH FIXED b_imp*********************
11 //
12 // AuthorList:
13 // Andrzej Olszewski: Initial Code February 2006
14 // Andrzej Olszewski: Converted to ROOT histograms July 2007
15 // Soumya Mohapatra : Re-written to check the new Flow implementations (JUNE 2011)
16 
17 #include "FlowAfterburner/CheckFlow_New.h"
18 #include "GeneratorObjects/McEventCollection.h"
19 
20 #include "GaudiKernel/SmartDataPtr.h"
21 #include "GaudiKernel/DataSvc.h"
22 
23 #include "GaudiKernel/ITHistSvc.h"
24 
25 #include "TH1D.h"
26 #include "TProfile.h"
27 
28 #include "AtlasHepMC/GenEvent.h"
29 #include "AtlasHepMC/GenParticle.h"
30 #include "AtlasHepMC/GenVertex.h"
31 
32 #include "GeneratorObjects/HijingEventParams.h"
33 
35 CheckFlow_New::CheckFlow_New(const std::string& name, ISvcLocator* pSvcLocator) :
36  AthAlgorithm(name, pSvcLocator)
37 {
38 }
39 
40 
41 StatusCode CheckFlow_New::initialize(){
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 }
197 
198 
199 
200 StatusCode CheckFlow_New::execute() {
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 }
363 
364 StatusCode CheckFlow_New::finalize() {
365  msg(MSG::INFO) << ">>> CheckFlow_New from finalize" << endmsg;
366 
367  return StatusCode::SUCCESS;
368 }
369 
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
GenEvent.h
ReadCellNoiseFromCool.name1
name1
Definition: ReadCellNoiseFromCool.py:233
CheckFlow_New.h
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
HijingEventParams.h
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
GenVertex.h
test_pyathena.pt
pt
Definition: test_pyathena.py:11
M_PI
#define M_PI
Definition: ActiveFraction.h:11
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
CheckFlow_New::execute
virtual StatusCode execute() override
Definition: CheckFlow_New.cxx:200
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
GenParticle.h
CheckFlow_New::m_hist_Psi_n_ebe
TH1D * m_hist_Psi_n_ebe[6]
Definition: CheckFlow_New.h:60
CheckFlow_New::m_hist_Psi_n_true
TH1D * m_hist_Psi_n_true[6]
Definition: CheckFlow_New.h:55
CheckFlow_New::m_profile_pt_dep_reco
TProfile * m_profile_pt_dep_reco[6][n_etabin]
Definition: CheckFlow_New.h:66
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
CheckFlow_New::finalize
virtual StatusCode finalize() override
Definition: CheckFlow_New.cxx:364
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
McEventCollection.h
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
CheckFlow_New::n_ptbin
@ n_ptbin
Definition: CheckFlow_New.h:50
ParticleGun_EoverP_Config.pid
pid
Definition: ParticleGun_EoverP_Config.py:62
CheckFlow_New::initialize
virtual StatusCode initialize() override
Definition: CheckFlow_New.cxx:41
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
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
AthAlgorithm
Definition: AthAlgorithm.h:47
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
CheckFlow_New::m_hist_vn_ebe
TH1D * m_hist_vn_ebe[6]
Definition: CheckFlow_New.h:62
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
python.CaloAddPedShiftConfig.int
int
Definition: CaloAddPedShiftConfig.py:45
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
CheckFlow_New::n_etabin
@ n_etabin
Definition: CheckFlow_New.h:51
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
CheckFlow_New.CheckFlow_New
CheckFlow_New
Definition: CheckFlow_New.py:54
CheckFlow_New::m_rapcut_min
DoubleProperty m_rapcut_min
Definition: CheckFlow_New.h:46
CheckFlow_New::m_rapcut_max
DoubleProperty m_rapcut_max
Definition: CheckFlow_New.h:47
Generated on Mon Sep 1 2025 21:07:48 for ATLAS Offline Software by   doxygen 1.8.18