53 {
54 StatusCode result = StatusCode::SUCCESS;
55
56 ATH_MSG_INFO(
">>> CheckFlow_New_Minbias from Initialize");
57
58
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
71
72 const float eta_bin_max = 4.0;
73
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);
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);
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);
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);
96 hist_vec.push_back(m_hist_psi_corr_reco [ihar_i][ib_imp]);
97 }
98
99
100
101
102
103
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);
112 sprintf(name,"hist_v%d_b%d_ebe_ID2",ihar+1,ib_imp);
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);
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);
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);
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);
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);
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);
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);
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);
176 hist_vec.push_back(m_profile_b_dep_reco [ihar][ipt][ieta]);
177 }
178 }
179
180 sprintf(name,"profile_resolution_%d",ihar+1);
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
194 m_tesIO =
new TruthHelper::GenAccessIO();
196}
TProfile(*args, **kwargs)
1.16.1