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  for ( int i = 0; i< 6; i++ ){
39  m_hist_Psi_n_true[i] = 0;
40  m_hist_Psi_n_reco[i] = 0;
41  m_hist_Psi_n_ebe[i] = 0;
42  m_hist_Psi_n_ebe_pt[i] = 0;
43  m_hist_vn_ebe[i] = 0;
44  }
45 
46  for ( int i = 0; i< 36; i++ ){
47  m_hist_psi_corr_true[i] = 0;
48  m_hist_psi_corr_reco[i] = 0;
49  }
50 }
51 
52 
53 StatusCode CheckFlow_New::initialize(){
54  ATH_CHECK(m_hijingKey.initialize());
55 
56  ATH_MSG_INFO(">>> CheckFlow_New from Initialize");
57 
58  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};
59  float eta_bin_max = 4.0;
60 
61 
62  SmartIF<ITHistSvc> rootHistSvc{Gaudi::svcLocator()->service("THistSvc")};
63  if (!rootHistSvc) {
64  ATH_MSG_ERROR( "Unable to locate THistSvc" );
65  return StatusCode::FAILURE;
66  }
67 
68  std::string StreamAndPath="/FlowOutPut/";
69  std::string histPath = StreamAndPath;
70  char name[100],name1[100];
71  for (int ihar=0;ihar<6;ihar++){
72 
73 // /*
74  sprintf(name,"hist_Psi_%d_true",ihar+1);
75  sprintf(name1,"Truth Psi_{%d} distribution;%dPsi_{%d} Truth;events",ihar+1,ihar+1,ihar+1);
76  m_hist_Psi_n_true [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
77  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_true[ihar]->GetName(),m_hist_Psi_n_true[ihar]).isFailure() ){
78  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_true[ihar]->GetName() << endmsg;
79  }
80  m_hist_Psi_n_true[ihar]->GetXaxis()->CenterTitle();
81  m_hist_Psi_n_true[ihar]->GetYaxis()->CenterTitle();
82 
83  sprintf(name,"hist_Psi_%d_reco",ihar+1);
84  sprintf(name1,"Reconstructed Psi_{%d} distribution;%dPsi_{%d} Reco;events",ihar+1,ihar+1,ihar+1);
85  m_hist_Psi_n_reco [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
86  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_reco[ihar]->GetName(),m_hist_Psi_n_reco[ihar]).isFailure() ){
87  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_reco[ihar]->GetName() << endmsg;
88  }
89  m_hist_Psi_n_reco[ihar]->GetXaxis()->CenterTitle();
90  m_hist_Psi_n_reco[ihar]->GetYaxis()->CenterTitle();
91 
92 
93  for (int ihar2=0;ihar2<6;ihar2++)
94  {
95  int ihar_i=ihar*6+ihar2;
96 
97  sprintf(name,"hist_Psi_corr_true_%d_%d",ihar+1,ihar2+1);
98  sprintf(name1,"true Psi_{%d} -Psi_{%d};%dPsi_{%d} -%dPsi_{%d} ;events",ihar+1,ihar2+1,ihar+1,ihar+1,ihar2+1,ihar2+1);
99  m_hist_psi_corr_true [ihar_i]=new TH1D (name,name1,1000,-2*M_PI,2*M_PI);
100  if ( rootHistSvc->regHist(histPath+m_hist_psi_corr_true[ihar_i]->GetName(),m_hist_psi_corr_true[ihar_i]).isFailure() ){
101  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_psi_corr_true[ihar_i]->GetName() << endmsg;
102  }
103  m_hist_psi_corr_true[ihar_i]->GetXaxis()->CenterTitle();
104  m_hist_psi_corr_true[ihar_i]->GetYaxis()->CenterTitle();
105 
106 
107  sprintf(name,"hist_Psi_corr_reco_%d_%d",ihar+1,ihar2+1);
108  sprintf(name1,"reco Psi_{%d} -Psi_{%d};%dPsi_{%d} -%dPsi_{%d} ;events",ihar+1,ihar2+1,ihar+1,ihar+1,ihar2+1,ihar2+1);
109  m_hist_psi_corr_reco [ihar_i]=new TH1D (name,name1,1000,-2*M_PI,2*M_PI);
110  if ( rootHistSvc->regHist(histPath+m_hist_psi_corr_reco[ihar_i]->GetName(),m_hist_psi_corr_reco[ihar_i]).isFailure() ){
111  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_psi_corr_reco[ihar_i]->GetName() << endmsg;
112  }
113  m_hist_psi_corr_reco[ihar_i]->GetXaxis()->CenterTitle();
114  m_hist_psi_corr_reco[ihar_i]->GetYaxis()->CenterTitle();
115  }
116 
117 
118 
119  //integrated vn event by event
120  sprintf(name,"hist_v%d_ebe",ihar+1);
121  sprintf(name1,"v%d;v%d;events",ihar+1,ihar+1);
122  m_hist_vn_ebe [ihar]=new TH1D (name,name1,1000,-0.5,0.5);
123  if ( rootHistSvc->regHist(histPath+m_hist_vn_ebe[ihar]->GetName(),m_hist_vn_ebe[ihar]).isFailure() ){
124  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_vn_ebe[ihar]->GetName() << endmsg;
125  }
126  m_hist_vn_ebe[ihar]->GetXaxis()->CenterTitle();
127  m_hist_vn_ebe[ihar]->GetYaxis()->CenterTitle();
128 
129  sprintf(name ,"hist_Psi%d_ebe",ihar+1);
130  sprintf(name1,"%d#Delta#Psi;%d(#Psi_{reco}-#Psi_{Truth});events",ihar+1,ihar+1);
131  m_hist_Psi_n_ebe [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
132  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_ebe[ihar]->GetName(),m_hist_Psi_n_ebe[ihar]).isFailure() ){
133  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_ebe[ihar]->GetName() << endmsg;
134  }
135  m_hist_Psi_n_ebe[ihar]->GetXaxis()->CenterTitle();
136  m_hist_Psi_n_ebe[ihar]->GetYaxis()->CenterTitle();
137 
138  sprintf(name ,"hist_Psi%d_ebe_pt",ihar+1);
139  sprintf(name1,"%d#Delta#Psi (pT weighted);%d(#Psi_{reco}-#Psi_{Truth});events",ihar+1,ihar+1);
140  m_hist_Psi_n_ebe_pt [ihar]=new TH1D (name,name1,1000,-M_PI,M_PI);
141  if ( rootHistSvc->regHist(histPath+m_hist_Psi_n_ebe_pt[ihar]->GetName(),m_hist_Psi_n_ebe_pt[ihar]).isFailure() ){
142  msg(MSG::WARNING) << "Can't book "<< histPath+m_hist_Psi_n_ebe_pt[ihar]->GetName() << endmsg;
143  }
144  m_hist_Psi_n_ebe_pt[ihar]->GetXaxis()->CenterTitle();
145  m_hist_Psi_n_ebe_pt[ihar]->GetYaxis()->CenterTitle();
146 
147 
148 
149 
150 
151 
152  for(int ieta=0;ieta<n_etabin;ieta++){
153  sprintf(name ,"profile_pt_dep_%d_eta%d" ,ihar+1,ieta);
154  sprintf(name1,"v%d vs pT (eta%d);pT;v%d",ihar+1,ieta,ihar+1);
155  m_profile_pt_dep [ihar][ieta]=new TProfile (name,name1,n_ptbin,pt_binvals);
156  if ( rootHistSvc->regHist(histPath+m_profile_pt_dep[ihar][ieta]->GetName(),m_profile_pt_dep[ihar][ieta]).isFailure() ){
157  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_pt_dep[ihar][ieta]->GetName() << endmsg;
158  }
159  m_profile_pt_dep [ihar][ieta]->GetXaxis()->CenterTitle();
160  m_profile_pt_dep [ihar][ieta]->GetYaxis()->CenterTitle();
161  }
162 
163  for(int ipt=0;ipt<n_ptbin;ipt++){
164  sprintf(name ,"profile_eta_dep_%d_pt%d",ihar+1,ipt);
165  sprintf(name1,"v%d vs #eta; (ipt%d)#eta;v%d",ihar+1,ipt,ihar+1);
166  m_profile_eta_dep [ihar][ipt]=new TProfile (name,name1,2*n_etabin, -eta_bin_max,eta_bin_max);
167  if ( rootHistSvc->regHist(histPath+m_profile_eta_dep[ihar][ipt]->GetName(),m_profile_eta_dep[ihar][ipt]).isFailure() ){
168  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_eta_dep[ihar][ipt]->GetName() << endmsg;
169  }
170  m_profile_eta_dep [ihar][ipt]->GetXaxis()->CenterTitle();
171  m_profile_eta_dep [ihar][ipt]->GetYaxis()->CenterTitle();
172  }
173 
174 
175  for(int ieta=0;ieta<n_etabin;ieta++){
176  sprintf(name ,"profile_pt_dep_reco_%d_eta%d",ihar+1,ieta);
177  sprintf(name1,"v%d vs pT (eta%d);pT;v%d",ihar+1,ieta,ihar+1);
178  m_profile_pt_dep_reco [ihar][ieta]=new TProfile (name,name1,n_ptbin,pt_binvals);
179  if ( rootHistSvc->regHist(histPath+m_profile_pt_dep_reco[ihar][ieta]->GetName(),m_profile_pt_dep_reco[ihar][ieta]).isFailure() ){
180  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_pt_dep_reco[ihar][ieta]->GetName() << endmsg;
181  }
182  m_profile_pt_dep_reco [ihar][ieta]->GetXaxis()->CenterTitle();
183  m_profile_pt_dep_reco [ihar][ieta]->GetYaxis()->CenterTitle();
184  }
185 
186 
187  for(int ipt=0;ipt<n_ptbin;ipt++){
188  sprintf(name ,"profile_eta_dep_reco_%d_pt%d",ihar+1,ipt);
189  sprintf(name1,"v%d vs #eta (pt%d);#eta;v%d",ihar+1,ipt,ihar+1);
190  m_profile_eta_dep_reco [ihar][ipt]=new TProfile (name,name1,2*n_etabin, -eta_bin_max,eta_bin_max);
191  if ( rootHistSvc->regHist(histPath+m_profile_eta_dep_reco[ihar][ipt]->GetName(),m_profile_eta_dep_reco[ihar][ipt]).isFailure() ){
192  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_eta_dep_reco[ihar][ipt]->GetName() << endmsg;
193  }
194  m_profile_eta_dep_reco [ihar][ipt]->GetXaxis()->CenterTitle();
195  m_profile_eta_dep_reco [ihar][ipt]->GetYaxis()->CenterTitle();
196  }
197 
198 
199  }
200  m_profile_resolution=new TProfile("profile_resolution","vn resolution;n;resolution",6, 0.5,6.5);
201  if(rootHistSvc->regHist(histPath+m_profile_resolution->GetName(),m_profile_resolution).isFailure() ){
202  msg(MSG::WARNING) << "Can't book "<< histPath+m_profile_resolution->GetName() << endmsg;
203  }
204 
205  msg(MSG::DEBUG) << "Histograms have been booked " << endmsg;
206  m_tesIO = new TruthHelper::GenAccessIO();
207  return StatusCode::SUCCESS;
208 }
209 
210 
211 
212 StatusCode CheckFlow_New::execute() {
213  msg(MSG::INFO) << ">>> CheckFlow_New from execute" << endmsg;
214 
215  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};
216  float eta_bin_max = 4.0;
217 
218  SG::ReadHandle<HijingEventParams> hijing_pars{m_hijingKey};
219  float b = hijing_pars->get_b();
220  float Psi_n[6],Psi_n_reco[6];
221  float Psi_n_reco_pos[6],Psi_n_reco_neg[6];
222  for(int ihar=0;ihar<6;ihar++){Psi_n[ihar]=hijing_pars->get_psi(ihar+1);}
223  msg(MSG::INFO)<<"SOUMYA "<<hijing_pars->get_psi(1)<<" "<<hijing_pars->get_psi(2)<<" "<<hijing_pars->get_psi(3)
224  <<hijing_pars->get_psi(4)<<" "<<hijing_pars->get_psi(5)<<" "<<hijing_pars->get_psi(6)<<" "<<b << endmsg;
225 
226 
227  // Check cut on impact parameter b
228  if(b<m_bcut_min || b>m_bcut_max) return StatusCode::SUCCESS;
229 
230 
231  double ngenerated_pos = 0,ngenerated_pt_pos=0;
232  double cos_n_pos[6],sin_n_pos[6],cos_n_pt_pos[6],sin_n_pt_pos[6];
233  double ngenerated_neg = 0,ngenerated_pt_neg=0;
234  double cos_n_neg[6],sin_n_neg[6],cos_n_pt_neg[6],sin_n_pt_neg[6];
235  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;
236  cos_n_neg[ihar]=0;sin_n_neg[ihar]=0; cos_n_pt_neg[ihar]=0;sin_n_pt_neg[ihar]=0;}
237 
238  // Iterate over all MC particles
239  std::vector<HepMC::ConstGenParticlePtr> particles;
240  StatusCode stat = m_tesIO->getMC(particles, false, m_key);
241  if (stat.isFailure()) {
242  msg(MSG::ERROR) << "Could not find " << m_key << endmsg;
243  return stat;
244  }
245  for (auto pitr: particles) {
246  int pid = pitr->pdg_id();
247  double pt = pitr->momentum().perp();
248  double rapid = pitr->momentum().pseudoRapidity();
249  double phi = pitr->momentum().phi();
250  msg(MSG::DEBUG)
251  << " PID = " << pid << " Status = " << pitr->status()
252  << " Eta = " << rapid << " Phi = " << phi<< endmsg;
253 
254  if( (std::abs(rapid) >= m_rapcut_min) && (std::abs(rapid) <= m_rapcut_max) &&
255  (std::abs(pt) >= m_ptcut_min) && (std::abs(pt) <= m_ptcut_max) ) {
256 
257  for(int ihar=0;ihar<6;ihar++){
258  float temp=(ihar+1)*(phi-Psi_n[ihar]);
259 
260  int ieta= (int)(std::abs(rapid)*n_etabin/eta_bin_max);
261  if(ieta>=0 && ieta<n_etabin) m_profile_pt_dep [ihar][ieta]->Fill(pt/1000,cos(temp));
262 
263 
264  float temp_pt=pt/1000;
265  for(int ipt=0;ipt<n_ptbin;ipt++){
266  if(temp_pt<pt_binvals[ipt+1]){
267  m_profile_eta_dep[ihar][ipt]->Fill(rapid ,cos(temp));
268  break;
269  }
270  }
271 
272  if( rapid >3.2 && rapid< 4.9){
273  cos_n_pos[ihar]+=std::cos( (ihar+1)*phi);
274  sin_n_pos[ihar]+=std::sin( (ihar+1)*phi);
275  ngenerated_pos++;
276 
277  cos_n_pt_pos[ihar]+=pt*std::cos( (ihar+1)*phi);
278  sin_n_pt_pos[ihar]+=pt*std::sin( (ihar+1)*phi);
279  ngenerated_pt_pos +=pt;
280  }
281  if( rapid <-3.2 && rapid >-4.9){
282  cos_n_neg[ihar]+=std::cos( (ihar+1)*phi);
283  sin_n_neg[ihar]+=std::sin( (ihar+1)*phi);
284  ngenerated_neg++;
285 
286  cos_n_pt_neg[ihar]+=pt*std::cos( (ihar+1)*phi);
287  sin_n_pt_neg[ihar]+=pt*std::sin( (ihar+1)*phi);
288  ngenerated_pt_neg +=pt;
289  }
290 
291  }
292  }
293  }
294 
295 
296 // Calculate the event by event vn and also the reconstructed Psi_n angles
297 // Also make correlation histos between Psi_n_truth and Psi_n_reco
298  float cos_n[6],sin_n[6],cos_n_pt[6],sin_n_pt[6];
299  for(int ihar=0;ihar<6;ihar++){
300  cos_n[ihar] = ( cos_n_pos[ihar]+ cos_n_neg[ihar] ) / (ngenerated_pos+ngenerated_neg);
301  sin_n[ihar] = ( sin_n_pos[ihar]+ sin_n_neg[ihar] ) / (ngenerated_pos+ngenerated_neg);
302 
303  float psi_reco=std::atan2(sin_n[ihar],cos_n[ihar])/(ihar+1);
304  m_hist_Psi_n_ebe[ihar]->Fill( (ihar+1)*(psi_reco-Psi_n[ihar]) );
305  m_hist_vn_ebe [ihar]->Fill(std::sqrt(cos_n[ihar]*cos_n[ihar] +sin_n[ihar]*sin_n[ihar] ));
306 
307  Psi_n_reco_pos[ihar]=std::atan2(sin_n_pos[ihar],cos_n_pos[ihar])/ (ihar+1);
308  Psi_n_reco_neg[ihar]=std::atan2(sin_n_neg[ihar],cos_n_neg[ihar])/ (ihar+1);
309  Psi_n_reco [ihar]=psi_reco;
310 
311 
312  cos_n_pt[ihar] = ( cos_n_pt_pos[ihar]+ cos_n_pt_neg[ihar] ) / (ngenerated_pt_pos+ngenerated_pt_neg);
313  sin_n_pt[ihar] = ( sin_n_pt_pos[ihar]+ sin_n_pt_neg[ihar] ) / (ngenerated_pt_pos+ngenerated_pt_neg);
314 
315  psi_reco=std::atan2(sin_n_pt[ihar],cos_n_pt[ihar])/(ihar+1);
316  m_hist_Psi_n_ebe_pt[ihar]->Fill( (ihar+1)*(psi_reco-Psi_n[ihar]) );
317  }
318 
319 
320 // Make the plots for the correlation between Psi_n truth (for different n) (same for Psi_n reco)
321  for(int ihar=0;ihar<6;ihar++){
322  m_hist_Psi_n_true[ihar]->Fill((ihar+1)*Psi_n[ihar]);
323  m_hist_Psi_n_reco[ihar]->Fill((ihar+1)*Psi_n_reco[ihar]);
324 
325  float psi1,psi2;
326  for(int ihar2=0;ihar2<6;ihar2++){
327  psi1=(ihar+1)*Psi_n[ihar];psi2=(ihar2+1)*Psi_n[ihar2];
328  m_hist_psi_corr_true[ihar*6+ihar2]->Fill( std::atan2( std::sin(psi1-psi2),std::cos(psi1-psi2) ) );
329 
330  psi1=(ihar+1)*Psi_n_reco[ihar];psi2=(ihar2+1)*Psi_n_reco[ihar2];
331  m_hist_psi_corr_reco[ihar*6+ihar2]->Fill( std::atan2( std::sin(psi1-psi2),std::cos(psi1-psi2) ) );
332  }
333  }
334 
335 
336 
337 
338 
339 // calculate the pt and eta dependence using the Psi_reco angles also fill the resolution TProfile
340  for(int ihar=0;ihar<6;ihar++){
341  m_profile_resolution->Fill( ihar+1, cos( (ihar+1) * (Psi_n_reco_pos[ihar] - Psi_n_reco_neg[ihar]) ) );
342  }
343  for (auto pitr: particles) {
344  double pt = pitr->momentum().perp();
345  double rapid = pitr->momentum().pseudoRapidity();
346  double phi = pitr->momentum().phi();
347  if( (std::abs(rapid) >= m_rapcut_min) && (std::abs(rapid) <= m_rapcut_max) &&
348  (std::abs(pt) >= m_ptcut_min) && (std::abs(pt) <= m_ptcut_max) ) {
349 
350  for(int ihar=0;ihar<6;ihar++){
351  float temp=(ihar+1)*(phi-Psi_n_reco_pos[ihar]);
352  if(rapid>0) temp=(ihar+1)*(phi-Psi_n_reco_neg[ihar]);
353 
354 
355  int ieta= (int)(std::abs(rapid)*n_etabin/eta_bin_max);
356  if(ieta>=0 && ieta<n_etabin) m_profile_pt_dep_reco [ihar][ieta]->Fill(pt/1000,std::cos(temp));
357 
358  float temp_pt=pt/1000;
359  for(int ipt=0;ipt<n_ptbin;ipt++){
360  if(temp_pt<pt_binvals[ipt+1]){
361  m_profile_eta_dep_reco[ihar][ipt]->Fill(rapid ,cos(temp));
362  break;
363  }
364  }
365  }
366  }
367  }
368 
369 
370 
371 
372 
373  return StatusCode::SUCCESS;
374 }
375 
376 StatusCode CheckFlow_New::finalize() {
377  msg(MSG::INFO) << ">>> CheckFlow_New from finalize" << endmsg;
378 
379  return StatusCode::SUCCESS;
380 }
381 
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
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
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:70
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
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
CheckFlow_New::n_etabin
@ n_etabin
Definition: CheckFlow_New.h:51
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
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
McEventCollection.h
lumiFormat.i
int i
Definition: lumiFormat.py:85
python.TrigEgammaMonitorHelper.TProfile
def TProfile(*args, **kwargs)
Definition: TrigEgammaMonitorHelper.py:81
CheckFlow_New::initialize
StatusCode initialize()
Definition: CheckFlow_New.cxx:53
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
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
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:266
AthAlgorithm
Definition: AthAlgorithm.h:47
CheckFlow_New::execute
StatusCode execute()
Definition: CheckFlow_New.cxx:212
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:221
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:77
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::n_ptbin
@ n_ptbin
Definition: CheckFlow_New.h:50
CheckFlow_New::m_bcut_max
DoubleProperty m_bcut_max
Definition: CheckFlow_New.h:43
CheckFlow_New::finalize
StatusCode finalize()
Definition: CheckFlow_New.cxx:376
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
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 Thu Nov 7 2024 21:11:54 for ATLAS Offline Software by   doxygen 1.8.18