Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
CheckFlow_New_Minbias.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_Minbias.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 VARYING b_imp(Min Bias)*********************
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_Minbias.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 
29 #include "AtlasHepMC/GenEvent.h"
30 #include "AtlasHepMC/GenParticle.h"
31 #include "AtlasHepMC/GenVertex.h"
32 
33 #include "GeneratorObjects/HijingEventParams.h"
34 
35 
36 //typedef std::vector<HepMC::ConstGenParticlePtr> MCParticleCollection ;
37 
38 CheckFlow_New_Minbias::CheckFlow_New_Minbias(const std::string& name, ISvcLocator* pSvcLocator) :
39  AthAlgorithm(name, pSvcLocator)
40 {
41  //Declare the algorithm's properties
42  declareProperty("McEventKey" , m_key ="FLOW_EVENT");
43  declareProperty("HistogramFlag" , m_produceHistogram = true );
44  declareProperty("ImpactCutMin" , m_bcut_min = 0 );
45  declareProperty("ImpactCutMax" , m_bcut_max = 99 );
46  declareProperty("PtCutMin" , m_ptcut_min = 0 );
47  declareProperty("PtCutMax" , m_ptcut_max = 999999 );
48  declareProperty("RapidityCutMin", m_rapcut_min = 0 );
49  declareProperty("RapidityCutMax", m_rapcut_max = 5.5 );
50 }
51 
52 
53 StatusCode CheckFlow_New_Minbias::initialize(){
54  StatusCode result = StatusCode::SUCCESS;
55 
56  ATH_MSG_INFO(">>> CheckFlow_New_Minbias from Initialize");
57 
58 
59  ServiceHandle<ITHistSvc> rootHistSvc("THistSvc",name());
60  CHECK( rootHistSvc.retrieve() );
61  std::string histPath = "/FlowOutPut/";
62  std::vector<TH1*> hist_vec;
63  char name[100],name1[100];
64 
65 
66  for (int ihar=0;ihar<6;ihar++){
67  for(int ib_imp=0;ib_imp<n_b_bins;ib_imp++){
68 
69  const float pt_binvals[]={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};
70  //const float b_bin_vals[]={0.0,3.4955,4.9315,6.0375,6.9695,7.7895,8.5335,9.2135,9.8515,10.4485,11.0175,
71  // 11.554,12.070,12.560,13.033,13.492,13.944,14.409,14.929,15.6425};
72  const float eta_bin_max = 4.0;
73 //----------------------------------------The reconstructed and truth Psi and the correlations between them-----------------------------
74  sprintf(name,"hist_Psi_%d_true_b%d",ihar+1,ib_imp);
75  sprintf(name1,"Truth Psi_{%d} distribution;%dPsi_{%d} Truth;events",ihar+1,ihar+1,ihar+1);
76  m_hist_Psi_n_true [ihar][ib_imp]=new TH1D (name,name1,1000,-M_PI,M_PI);
77  hist_vec.push_back(m_hist_Psi_n_true [ihar][ib_imp]);
78 
79  sprintf(name,"hist_Psi_%d_reco_b%d",ihar+1,ib_imp);
80  sprintf(name1,"Reconstructed Psi_{%d} distribution;%dPsi_{%d} Reco;events",ihar+1,ihar+1,ihar+1);
81  m_hist_Psi_n_reco [ihar][ib_imp]=new TH1D (name,name1,1000,-M_PI,M_PI);
82  hist_vec.push_back(m_hist_Psi_n_reco [ihar][ib_imp]);
83 
84 
85  for (int ihar2=0;ihar2<6;ihar2++){
86  int ihar_i=ihar*6+ihar2;
87 
88  sprintf(name,"hist_Psi_corr_true_b%d_%d_%d",ib_imp,ihar+1,ihar2+1);
89  sprintf(name1,"true Psi_{%d} -Psi_{%d};%dPsi_{%d} -%dPsi_{%d} ;events",ihar+1,ihar2+1,ihar+1,ihar+1,ihar2+1,ihar2+1);
90  m_hist_psi_corr_true [ihar_i][ib_imp]=new TH1D (name,name1,1000,-2*M_PI,2*M_PI);
91  hist_vec.push_back(m_hist_psi_corr_true [ihar_i][ib_imp]);
92 
93  sprintf(name,"hist_Psi_corr_reco_%d_%d_%d",ib_imp,ihar+1,ihar2+1);
94  sprintf(name1,"reco Psi_{%d} -Psi_{%d};%dPsi_{%d} -%dPsi_{%d} ;events",ihar+1,ihar2+1,ihar+1,ihar+1,ihar2+1,ihar2+1);
95  m_hist_psi_corr_reco [ihar_i][ib_imp]=new TH1D (name,name1,1000,-2*M_PI,2*M_PI);
96  hist_vec.push_back(m_hist_psi_corr_reco [ihar_i][ib_imp]);
97  }
98 //------------------------------------------------------------------------------------------
99 
100 
101 
102 //-----------------------------------------------------------------------------------------
103  //integrated vn event by event
104  sprintf(name,"hist_v%d_b%d_ebe",ihar+1,ib_imp);
105  sprintf(name1,"v%d;v%d;events",ihar+1,ihar+1);
106  m_hist_vn_ebe [ihar][ib_imp]=new TH1D (name,name1,1000,-0.5,0.5);
107  hist_vec.push_back(m_hist_vn_ebe [ihar][ib_imp]);
108 
109  sprintf(name,"hist_v%d_b%d_ebe_ID1",ihar+1,ib_imp);
110  sprintf(name1,"v%d;v%d;events",ihar+1,ihar+1);
111  m_hist_vn_ebe_ID1 [ihar][ib_imp]=new TH1D (name,name1,800,0.0,0.4);
112  sprintf(name,"hist_v%d_b%d_ebe_ID2",ihar+1,ib_imp);
113  m_hist_vn_ebe_ID2 [ihar][ib_imp]=new TH1D (name,name1,800,0.0,0.4);
114  hist_vec.push_back(m_hist_vn_ebe_ID1 [ihar][ib_imp]);
115  hist_vec.push_back(m_hist_vn_ebe_ID2 [ihar][ib_imp]);
116 
117  sprintf(name ,"hist_Psi%d_b%d_ebe",ihar+1,ib_imp);
118  sprintf(name1,"%d#Delta#Psi;%d(#Psi_{reco}-#Psi_{Truth});events",ihar+1,ihar+1);
119  m_hist_Psi_n_ebe [ihar][ib_imp]=new TH1D (name,name1,1000,-M_PI,M_PI);
120  hist_vec.push_back(m_hist_Psi_n_ebe [ihar][ib_imp]);
121 
122  sprintf(name ,"hist_Psi%d_b%d_ebe_pt",ihar+1,ib_imp);
123  sprintf(name1,"%d#Delta#Psi (pT weighted);%d(#Psi_{reco}-#Psi_{Truth});events",ihar+1,ihar+1);
124  m_hist_Psi_n_ebe_pt [ihar][ib_imp]=new TH1D (name,name1,1000,-M_PI,M_PI);
125  hist_vec.push_back(m_hist_Psi_n_ebe_pt [ihar][ib_imp]);
126 //-----------------------------------------------------------------------------------------
127 
128 
129 
130 
131 
132 //-----------------------------------------------------------------------------------------
133  for(int ieta=0;ieta<n_etabin;ieta++){
134  sprintf(name ,"profile_pt_dep_%d_eta%d_b%d" ,ihar+1,ieta,ib_imp);
135  sprintf(name1,"v%d vs pT (eta%d);pT;v%d",ihar+1,ieta,ihar+1);
136  m_profile_pt_dep [ihar][ib_imp][ieta]=new TProfile (name,name1,n_ptbin,pt_binvals);
137  hist_vec.push_back(m_profile_pt_dep [ihar][ib_imp][ieta]);
138  }
139 
140  for(int ipt=0;ipt<n_ptbin;ipt++){
141  sprintf(name ,"profile_eta_dep_%d_pt%d_b%d",ihar+1,ipt,ib_imp);
142  sprintf(name1,"v%d vs #eta; (ipt%d)#eta;v%d",ihar+1,ipt,ihar+1);
143  m_profile_eta_dep [ihar][ib_imp][ipt]=new TProfile (name,name1,2*n_etabin, -eta_bin_max,eta_bin_max);
144  hist_vec.push_back(m_profile_eta_dep [ihar][ib_imp][ipt]);
145  }
146 
147 
148  for(int ieta=0;ieta<n_etabin;ieta++){
149  sprintf(name ,"profile_pt_dep_reco_%d_eta%d_b%d",ihar+1,ieta,ib_imp);
150  sprintf(name1,"v%d vs pT (eta%d);pT;v%d",ihar+1,ieta,ihar+1);
151  m_profile_pt_dep_reco [ihar][ib_imp][ieta]=new TProfile (name,name1,n_ptbin,pt_binvals);
152  hist_vec.push_back(m_profile_pt_dep_reco [ihar][ib_imp][ieta]);
153  }
154 
155 
156  for(int ipt=0;ipt<n_ptbin;ipt++){
157  sprintf(name ,"profile_eta_dep_reco_%d_pt%d_b%d",ihar+1,ipt,ib_imp);
158  sprintf(name1,"v%d vs #eta (pt%d);#eta;v%d",ihar+1,ipt,ihar+1);
159  m_profile_eta_dep_reco [ihar][ib_imp][ipt]=new TProfile (name,name1,2*n_etabin, -eta_bin_max,eta_bin_max);
160  hist_vec.push_back(m_profile_eta_dep_reco [ihar][ib_imp][ipt]);
161  }
162  }
163 
164 
165 
166  for(int ipt=0;ipt<n_ptbin;ipt++){
167  for(int ieta=0;ieta<n_etabin;ieta++){
168  sprintf(name ,"profile_b_dep_%d_pt%d_eta%d",ihar+1,ipt,ieta);
169  sprintf(name1,"v%d vs cent ;cent;v%d",ihar+1,ihar+1);
170  m_profile_b_dep [ihar][ipt][ieta]=new TProfile (name,name1,n_b_bins,-0.5,n_b_bins-0.5);
171  hist_vec.push_back(m_profile_b_dep [ihar][ipt][ieta]);
172 
173  sprintf(name ,"profile_b_dep_reco_%d_pt%d_eta%d",ihar+1,ipt,ieta);
174  sprintf(name1,"v%d vs cent ;cent;v%d",ihar+1,ihar+1);
175  m_profile_b_dep_reco [ihar][ipt][ieta]=new TProfile (name,name1,n_b_bins,-0.5,n_b_bins-0.5);
176  hist_vec.push_back(m_profile_b_dep_reco [ihar][ipt][ieta]);
177  }
178  }
179 
180  sprintf(name,"profile_resolution_%d",ihar+1);
181  m_profile_resolution[ihar]=new TProfile(name,"vn resolution;n;resolution",n_b_bins,-0.5,n_b_bins-0.5);
182  hist_vec.push_back(m_profile_resolution[ihar]);
183  }
184 //-----------------------------------------------------------------------------------------
185 
186  for(auto& hist: hist_vec){
187  CHECK(rootHistSvc->regHist(histPath+hist->GetName(),hist));
188  hist->GetXaxis()->CenterTitle();
189  hist->GetYaxis()->CenterTitle();
190  }
191 
192 
193  ATH_MSG_DEBUG("Histograms have been booked ");
194  m_tesIO = new TruthHelper::GenAccessIO();
195  return result;
196 }
197 
198 
199 
200 StatusCode CheckFlow_New_Minbias::execute() {
201 
202  ATH_MSG_INFO(">>> CheckFlow_New_Minbias from execute");
203  const float pt_binvals[]={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  const float b_bin_vals[]={0.0,3.4955,4.9315,6.0375,6.9695,7.7895,8.5335,9.2135,9.8515,10.4485,11.0175,11.554,12.070,12.560,13.033,13.492,13.944,14.409,14.929,15.6425};
205  const float eta_bin_max = 4.0;
206  const HijingEventParams *hijing_pars;
207 
208 
209  CHECK(evtStore()->retrieve(hijing_pars, "Hijing_event_params"));
210  float b = hijing_pars->get_b();
211  float Psi_n[6],Psi_n_reco[6];
212  for(int ihar=0;ihar<6;ihar++){Psi_n[ihar]=hijing_pars->get_psi(ihar+1);}
213  ATH_MSG_INFO("SOUMYA "<<hijing_pars->get_psi(1)<<" "<<hijing_pars->get_psi(2)<<" "\
214  <<hijing_pars->get_psi(3)<<" "<<hijing_pars->get_psi(4)<<" "\
215  <<hijing_pars->get_psi(5)<<" "<<hijing_pars->get_psi(6)<<" "\
216  <<b);
217 
218 
219 
220  int ib_imp=-1;
221  for(int ib=0;ib<n_b_bins;ib++){if(b<b_bin_vals[ib+1]) {ib_imp=ib;break;}}
222  if(ib_imp<0) return StatusCode::SUCCESS;
223  if(ib_imp==0) {std::cout<<"AAAAAAAAAAAAAAAAAAAA "<<b<<std::endl;}
224 
225  // Check cut on impact parameter b
226  if(b<m_bcut_min || b>m_bcut_max) return StatusCode::SUCCESS;
227 
228 
229  double ngenerated_pos = 0,ngenerated_pt_pos=0;
230  double ngenerated_neg = 0,ngenerated_pt_neg=0;
231  double cos_n_pos[6],sin_n_pos[6],cos_n_pt_pos[6],sin_n_pt_pos[6];
232  double cos_n_neg[6],sin_n_neg[6],cos_n_pt_neg[6],sin_n_pt_neg[6];
233  double cos_ID1[6],sin_ID1[6],tot_ID1=0.0;
234  double cos_ID2[6],sin_ID2[6],tot_ID2=0.0;
235  for(int ihar=0;ihar<6;ihar++){
236  cos_n_pos [ihar]=0;
237  sin_n_pos [ihar]=0;
238  cos_n_pt_pos[ihar]=0;
239  sin_n_pt_pos[ihar]=0;
240  cos_n_neg [ihar]=0;
241  sin_n_neg [ihar]=0;
242  cos_n_pt_neg[ihar]=0;
243  sin_n_pt_neg[ihar]=0;
244  cos_ID1 [ihar]=0;
245  sin_ID1 [ihar]=0;
246  cos_ID2 [ihar]=0;
247  sin_ID2 [ihar]=0;
248  }
249 
250 
251  // Iterate over MC particles
252  std::vector<HepMC::ConstGenParticlePtr> particles;
253  CHECK(m_tesIO->getMC(particles, false, m_key));
254 
255  for (auto pitr: particles) {
256  int pid = pitr->pdg_id();
257  int p_stat = pitr->status();
258  double pt = pitr->momentum().perp();
259  double rapid = pitr->momentum().pseudoRapidity();
260  double phi = pitr->momentum().phi();
261  ATH_MSG_DEBUG(" PID = " << pid << " Status = " << p_stat \
262  << " Eta = " << rapid << " Phi = " << phi);
263 
264  if( (std::abs(rapid) >= m_rapcut_min) &&
265  (std::abs(rapid) <= m_rapcut_max) &&
266  (std::abs(pt) >= m_ptcut_min ) &&
267  (std::abs(pt) <= m_ptcut_max)) {
268 
269  for(int ihar=0;ihar<6;ihar++){
270  float temp=(ihar+1)*(phi-Psi_n[ihar]);
271 
272  int ieta= (int)(std::abs(rapid)*n_etabin/eta_bin_max);
273  if(ieta>=0 && ieta<n_etabin) m_profile_pt_dep [ihar][ib_imp][ieta]->Fill(pt/1000,cos(temp));
274 
275 
276  float temp_pt=pt/1000;
277  for(int ipt=0;ipt<n_ptbin;ipt++){
278  if(temp_pt<pt_binvals[ipt+1]){
279  m_profile_eta_dep[ihar][ib_imp][ipt]->Fill(rapid ,cos(temp));
280  if(ieta>=0 && ieta<n_etabin) m_profile_b_dep [ihar][ipt][ieta]->Fill(ib_imp,cos(temp));
281  break;
282  }
283  }
284 
285  if( rapid >3.2 && rapid< 4.9){
286  cos_n_pos[ihar]+=std::cos( (ihar+1)*phi);
287  sin_n_pos[ihar]+=std::sin( (ihar+1)*phi);
288  ngenerated_pos++;
289 
290  cos_n_pt_pos[ihar]+=pt*cos( (ihar+1)*phi);
291  sin_n_pt_pos[ihar]+=pt*sin( (ihar+1)*phi);
292  ngenerated_pt_pos +=pt;
293  }
294  if( rapid <-3.2 && rapid >-4.9){
295  cos_n_neg[ihar]+=std::cos( (ihar+1)*phi);
296  sin_n_neg[ihar]+=std::sin( (ihar+1)*phi);
297  ngenerated_neg++;
298 
299  cos_n_pt_neg[ihar]+=pt*std::cos( (ihar+1)*phi);
300  sin_n_pt_neg[ihar]+=pt*std::sin( (ihar+1)*phi);
301  ngenerated_pt_neg +=pt;
302  }
303  }
304  }
305 
306  //EbE vn for ID acceptance end pt>0.5GeV
307  if(std::abs(pt)>=500){
308  tot_ID1++;
309  for(int ihar=0;ihar<6;ihar++){
310  cos_ID1[ihar]+=cos((ihar+1)*phi);
311  sin_ID1[ihar]+=sin((ihar+1)*phi);
312  }
313  }
314  if(std::abs(pt)>=0){
315  tot_ID2++;
316  for(int ihar=0;ihar<6;ihar++){
317  cos_ID2[ihar]+=cos((ihar+1)*phi);
318  sin_ID2[ihar]+=sin((ihar+1)*phi);
319  }
320  }
321  }
322 
323 
324 //EbE vn for ID acceptance end pt>0.5GeV
325  for(int ihar=0;ihar<6;ihar++){
326  if(tot_ID1>0.01){
327  double temp1= std::sqrt(cos_ID1[ihar]*cos_ID1[ihar] + sin_ID1[ihar]*sin_ID1[ihar])/tot_ID1;
328  m_hist_vn_ebe_ID1[ihar][ib_imp]->Fill(temp1);
329  }
330  if(tot_ID2>0.01){
331  double temp2= std::sqrt(cos_ID2[ihar]*cos_ID2[ihar] + sin_ID2[ihar]*sin_ID2[ihar])/tot_ID2;
332  m_hist_vn_ebe_ID2[ihar][ib_imp]->Fill(temp2);
333  }
334  }
335 
336 
337 // Calculate the event by event vn and also the reconstructed Psi_n angles
338 // Also make correlation histos between Psi_n_truth and Psi_n_reco
339  float cos_n[6],sin_n[6],cos_n_pt[6],sin_n_pt[6];
340  float Psi_n_reco_pos[6],Psi_n_reco_neg[6];
341  for(int ihar=0;ihar<6;ihar++){
342  cos_n[ihar] = ( cos_n_pos[ihar]+ cos_n_neg[ihar] ) / (ngenerated_pos+ngenerated_neg);
343  sin_n[ihar] = ( sin_n_pos[ihar]+ sin_n_neg[ihar] ) / (ngenerated_pos+ngenerated_neg);
344 
345  float psi_reco=std::atan2(sin_n[ihar],cos_n[ihar])/(ihar+1);
346  m_hist_Psi_n_ebe[ihar][ib_imp]->Fill( (ihar+1)*(psi_reco-Psi_n[ihar]) );
347  m_hist_vn_ebe [ihar][ib_imp]->Fill(std::sqrt(cos_n[ihar]*cos_n[ihar] +sin_n[ihar]*sin_n[ihar] ));
348 
349  Psi_n_reco_pos[ihar]=std::atan2(sin_n_pos[ihar],cos_n_pos[ihar])/ (ihar+1);
350  Psi_n_reco_neg[ihar]=std::atan2(sin_n_neg[ihar],cos_n_neg[ihar])/ (ihar+1);
351  Psi_n_reco [ihar]=psi_reco;
352 
353 
354  cos_n_pt[ihar] = ( cos_n_pt_pos[ihar]+ cos_n_pt_neg[ihar] ) / (ngenerated_pt_pos+ngenerated_pt_neg);
355  sin_n_pt[ihar] = ( sin_n_pt_pos[ihar]+ sin_n_pt_neg[ihar] ) / (ngenerated_pt_pos+ngenerated_pt_neg);
356 
357  psi_reco=std::atan2(sin_n_pt[ihar],cos_n_pt[ihar])/(ihar+1);
358  m_hist_Psi_n_ebe_pt[ihar][ib_imp]->Fill( (ihar+1)*(psi_reco-Psi_n[ihar]) );
359  }
360 
361 
362 // Make the plots for the correlation between Psi_n truth (for different n) (same for Psi_n reco)
363  for(int ihar=0;ihar<6;ihar++){
364  m_hist_Psi_n_true[ihar][ib_imp]->Fill((ihar+1)*Psi_n[ihar]);
365  m_hist_Psi_n_reco[ihar][ib_imp]->Fill((ihar+1)*Psi_n_reco[ihar]);
366 
367  float psi1,psi2;
368  for(int ihar2=0;ihar2<6;ihar2++){
369  psi1=(ihar+1)*Psi_n[ihar];psi2=(ihar2+1)*Psi_n[ihar2];
370  m_hist_psi_corr_true[ihar*6+ihar2][ib_imp]->Fill( std::atan2( std::sin(psi1-psi2),std::cos(psi1-psi2) ) );
371 
372  psi1=(ihar+1)*Psi_n_reco[ihar];psi2=(ihar2+1)*Psi_n_reco[ihar2];
373  m_hist_psi_corr_reco[ihar*6+ihar2][ib_imp]->Fill( std::atan2( std::sin(psi1-psi2),std::cos(psi1-psi2) ) );
374  }
375  }
376 
377 
378 
379 
380 
381 // calculate the pt and eta dependence using the Psi_reco angles also fill the resolution TProfile
382  for(int ihar=0;ihar<6;ihar++){
383  m_profile_resolution[ihar]->Fill( ib_imp, cos( (ihar+1) * (Psi_n_reco_pos[ihar] - Psi_n_reco_neg[ihar]) ) );
384  if(ib_imp==0) {std::cout<<"i11111111111111111111 "<<b<<std::endl;}
385  }
386  for (auto pitr: particles) {
387  double pt = pitr->momentum().perp();
388  double rapid = pitr->momentum().pseudoRapidity();
389  double phi = pitr->momentum().phi();
390  if( (std::abs(rapid) >= m_rapcut_min) && (std::abs(rapid) <= m_rapcut_max) &&
391  (std::abs(pt) >= m_ptcut_min) && (std::abs(pt) <= m_ptcut_max) ) {
392 
393  for(int ihar=0;ihar<6;ihar++){
394  float temp=(ihar+1)*(phi-Psi_n_reco_pos[ihar]);
395  if(rapid>0) temp=(ihar+1)*(phi-Psi_n_reco_neg[ihar]);
396 
397 
398  int ieta= (int)(std::abs(rapid)*n_etabin/eta_bin_max);
399  if(ieta>=0 && ieta<n_etabin) m_profile_pt_dep_reco [ihar][ib_imp][ieta]->Fill(pt/1000,cos(temp));
400 
401  float temp_pt=pt/1000;
402  for(int ipt=0;ipt<n_ptbin;ipt++){
403  if(temp_pt<pt_binvals[ipt+1]){
404  m_profile_eta_dep_reco[ihar][ib_imp][ipt]->Fill(rapid ,cos(temp));
405  if(ieta>=0 && ieta<n_etabin) m_profile_b_dep_reco [ihar][ipt][ieta]->Fill(ib_imp,cos(temp));
406  break;
407  }
408  }
409  }
410  }
411  }
412 
413  return StatusCode::SUCCESS;
414 }
415 
416 StatusCode CheckFlow_New_Minbias::finalize() {
417 
418  ATH_MSG_INFO(">>> CheckFlow_New from finalize");
419 
420  return StatusCode::SUCCESS;
421 }
422 
CheckFlow_New_Minbias::m_rapcut_min
double m_rapcut_min
Definition: CheckFlow_New_Minbias.h:48
CheckFlow_New_Minbias::m_hist_vn_ebe_ID1
TH1D * m_hist_vn_ebe_ID1[6][n_b_bins]
Definition: CheckFlow_New_Minbias.h:67
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
GenEvent.h
ReadCellNoiseFromCool.name1
name1
Definition: ReadCellNoiseFromCool.py:233
get_generator_info.result
result
Definition: get_generator_info.py:21
CheckFlow_New_Minbias::m_hist_Psi_n_reco
TH1D * m_hist_Psi_n_reco[6][n_b_bins]
Definition: CheckFlow_New_Minbias.h:60
CheckFlow_New_Minbias::n_b_bins
@ n_b_bins
Definition: CheckFlow_New_Minbias.h:55
CheckFlow_New_Minbias::m_hist_psi_corr_reco
TH1D * m_hist_psi_corr_reco[36][n_b_bins]
Definition: CheckFlow_New_Minbias.h:62
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
CheckFlow_New_Minbias.h
HijingEventParams.h
CheckFlow_New_Minbias::m_profile_resolution
TProfile * m_profile_resolution[6]
Definition: CheckFlow_New_Minbias.h:78
CheckFlow_New::m_ptcut_max
DoubleProperty m_ptcut_max
Definition: CheckFlow_New.h:45
CheckFlow_New_Minbias::m_bcut_max
double m_bcut_max
Definition: CheckFlow_New_Minbias.h:45
CheckFlow_New_Minbias::m_profile_eta_dep
TProfile * m_profile_eta_dep[6][n_b_bins][n_ptbin]
Definition: CheckFlow_New_Minbias.h:71
CheckFlow_New_Minbias::m_hist_Psi_n_ebe
TH1D * m_hist_Psi_n_ebe[6][n_b_bins]
Definition: CheckFlow_New_Minbias.h:64
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
plotmaker.hist
hist
Definition: plotmaker.py:148
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
PlotCalibFromCool.ib
ib
Definition: PlotCalibFromCool.py:419
HijingEventParams::get_psi
float get_psi(int n) const
Definition: HijingEventParams.h:95
drawFromPickle.cos
cos
Definition: drawFromPickle.py:36
CheckFlow_New::m_produceHistogram
BooleanProperty m_produceHistogram
Definition: CheckFlow_New.h:40
GenParticle.h
CheckFlow_New_Minbias::m_profile_b_dep
TProfile * m_profile_b_dep[6][n_ptbin][n_etabin]
Definition: CheckFlow_New_Minbias.h:75
AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
CheckFlow_New_Minbias::m_hist_psi_corr_true
TH1D * m_hist_psi_corr_true[36][n_b_bins]
Definition: CheckFlow_New_Minbias.h:61
CheckFlow_New_Minbias::finalize
StatusCode finalize()
Definition: CheckFlow_New_Minbias.cxx:416
McEventCollection.h
python.TrigEgammaMonitorHelper.TProfile
def TProfile(*args, **kwargs)
Definition: TrigEgammaMonitorHelper.py:81
CheckFlow_New_Minbias::n_etabin
@ n_etabin
Definition: CheckFlow_New_Minbias.h:54
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
HijingEventParams
Definition: HijingEventParams.h:23
python.LArMinBiasAlgConfig.int
int
Definition: LArMinBiasAlgConfig.py:59
ParticleGun_EoverP_Config.pid
pid
Definition: ParticleGun_EoverP_Config.py:62
CheckFlow_New_Minbias::m_ptcut_max
double m_ptcut_max
Definition: CheckFlow_New_Minbias.h:47
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition: Control/AthenaKernel/AthenaKernel/errorcheck.h:422
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
CheckFlow_New_Minbias::m_profile_b_dep_reco
TProfile * m_profile_b_dep_reco[6][n_ptbin][n_etabin]
Definition: CheckFlow_New_Minbias.h:76
CheckFlow_New_Minbias::execute
StatusCode execute()
Definition: CheckFlow_New_Minbias.cxx:200
CheckFlow_New_Minbias::m_profile_eta_dep_reco
TProfile * m_profile_eta_dep_reco[6][n_b_bins][n_ptbin]
Definition: CheckFlow_New_Minbias.h:73
CheckFlow_New_Minbias::initialize
StatusCode initialize()
Definition: CheckFlow_New_Minbias.cxx:53
AthAlgorithm
Definition: AthAlgorithm.h:47
CheckFlow_New::m_key
StringProperty m_key
Definition: CheckFlow_New.h:39
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
TruthHelper::GenAccessIO
Definition: Generators/FlowAfterburner/FlowAfterburner/GenAccessIO.h:21
CheckFlow_New_Minbias::m_ptcut_min
double m_ptcut_min
Definition: CheckFlow_New_Minbias.h:46
CheckFlow_New_Minbias::m_key
std::string m_key
Definition: CheckFlow_New_Minbias.h:41
CheckFlow_New::m_bcut_max
DoubleProperty m_bcut_max
Definition: CheckFlow_New.h:43
HijingEventParams::get_b
float get_b() const
Definition: HijingEventParams.h:91
CheckFlow_New_Minbias::m_hist_Psi_n_true
TH1D * m_hist_Psi_n_true[6][n_b_bins]
Definition: CheckFlow_New_Minbias.h:59
CheckFlow_New::m_bcut_min
DoubleProperty m_bcut_min
Definition: CheckFlow_New.h:42
CheckFlow_New_Minbias::m_profile_pt_dep
TProfile * m_profile_pt_dep[6][n_b_bins][n_etabin]
Definition: CheckFlow_New_Minbias.h:70
CheckFlow_New_Minbias::m_profile_pt_dep_reco
TProfile * m_profile_pt_dep_reco[6][n_b_bins][n_etabin]
Definition: CheckFlow_New_Minbias.h:72
LArG4FSStartPointFilter.particles
list particles
Definition: LArG4FSStartPointFilter.py:84
CheckFlow_New_Minbias::m_tesIO
TruthHelper::GenAccessIO * m_tesIO
Definition: CheckFlow_New_Minbias.h:80
CheckFlow_New_Minbias.CheckFlow_New_Minbias
CheckFlow_New_Minbias
Definition: CheckFlow_New_Minbias.py:51
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
CheckFlow_New_Minbias::m_rapcut_max
double m_rapcut_max
Definition: CheckFlow_New_Minbias.h:49
CheckFlow_New::m_rapcut_min
DoubleProperty m_rapcut_min
Definition: CheckFlow_New.h:46
CheckFlow_New_Minbias::n_ptbin
@ n_ptbin
Definition: CheckFlow_New_Minbias.h:53
CheckFlow_New_Minbias::m_hist_Psi_n_ebe_pt
TH1D * m_hist_Psi_n_ebe_pt[6][n_b_bins]
Definition: CheckFlow_New_Minbias.h:65
CheckFlow_New_Minbias::m_hist_vn_ebe_ID2
TH1D * m_hist_vn_ebe_ID2[6][n_b_bins]
Definition: CheckFlow_New_Minbias.h:68
ServiceHandle< ITHistSvc >
CheckFlow_New::m_rapcut_max
DoubleProperty m_rapcut_max
Definition: CheckFlow_New.h:47
CheckFlow_New_Minbias::m_hist_vn_ebe
TH1D * m_hist_vn_ebe[6][n_b_bins]
Definition: CheckFlow_New_Minbias.h:66
Generated on Mon Mar 31 2025 21:08:27 for ATLAS Offline Software by   doxygen 1.8.18