13 m_2occupancy(nullptr),
32 type(std::move(recObj))
37 std::string Zsig(
"Z");
38 std::string Jsig(
"Jpsi");
40 bool doZ =
false;
bool doJ =
false;
42 if (!
type.compare(Zsig)) {
name=
"Z"; doZ =
true; }
43 else if (!
type.compare(Jsig)) {
name=
"J/#psi"; doJ =
true; }
44 else {
name=
"other"; }
48 double maxBin = 1000.0;
96 void RecoPhysPlots::fill(
const std::vector<std::pair<const xAOD::Muon*, const xAOD::Muon*> >& mumucandidates){
97 for (
auto mumu: mumucandidates){
98 const float invariant_mass = (mumu.first->p4() + mumu.second->p4()).M();
99 fill(mumu.first->eta(), mumu.second->eta(), invariant_mass);
110 if((etaRegionCoarsePlus == 3.5 && etaRegionCoarseMinus == 3.5)) {
m_M_EA_EA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(16);}
111 if((etaRegionCoarsePlus == 3.5 && etaRegionCoarseMinus == 2.5)) {
m_M_EA_BA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(15);}
112 if((etaRegionCoarsePlus == 3.5 && etaRegionCoarseMinus == 1.5)) {
m_M_EA_BC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(14);}
113 if((etaRegionCoarsePlus == 3.5 && etaRegionCoarseMinus == 0.5)) {
m_M_EA_EC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(13);}
115 if((etaRegionCoarsePlus == 2.5 && etaRegionCoarseMinus == 3.5)) {
m_M_BA_EA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(12);}
116 if((etaRegionCoarsePlus == 2.5 && etaRegionCoarseMinus == 2.5)) {
m_M_BA_BA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(11);}
117 if((etaRegionCoarsePlus == 2.5 && etaRegionCoarseMinus == 1.5)) {
m_M_BA_BC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(10);}
118 if((etaRegionCoarsePlus == 2.5 && etaRegionCoarseMinus == 0.5)) {
m_M_BA_EC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(9);}
120 if((etaRegionCoarsePlus == 1.5 && etaRegionCoarseMinus == 3.5)) {
m_M_BC_EA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(8);}
121 if((etaRegionCoarsePlus == 1.5 && etaRegionCoarseMinus == 2.5)) {
m_M_BC_BA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(7);}
122 if((etaRegionCoarsePlus == 1.5 && etaRegionCoarseMinus == 1.5)) {
m_M_BC_BC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(6);}
123 if((etaRegionCoarsePlus == 1.5 && etaRegionCoarseMinus == 0.5)) {
m_M_BC_EC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(5);}
125 if((etaRegionCoarsePlus == 0.5 && etaRegionCoarseMinus == 3.5)) {
m_M_EC_EA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(4);}
126 if((etaRegionCoarsePlus == 0.5 && etaRegionCoarseMinus == 2.5)) {
m_M_EC_BA->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(3);}
127 if((etaRegionCoarsePlus == 0.5 && etaRegionCoarseMinus == 1.5)) {
m_M_EC_BC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(2);}
128 if((etaRegionCoarsePlus == 0.5 && etaRegionCoarseMinus == 0.5)) {
m_M_EC_EC->Fill(invariant_mass/1000.0);
m_2occupancy->Fill(1);}
130 m_Mass->Fill(invariant_mass/1000.0);
157 if (
eta < -2.50) {
return -1.;}
158 if (
eta < -2.00) {
return 0.5;}
159 if (
eta < -1.70) {
return 1.5;}
160 if (
eta < -1.05) {
return 2.5;}
161 if (
eta < 0.00) {
return 3.5;}
162 if (
eta > 2.50) {
return 8.5;}
163 if (
eta > 2.00) {
return 7.5;}
164 if (
eta > 1.70) {
return 6.5;}
165 if (
eta > 1.05) {
return 5.5;}
166 if (
eta > 0.00) {
return 4.5;}
173 if(
eta < -2.50) {
return -1.;}
174 if(
eta < -1.05) {
return 0.5;}
175 if(
eta < 0.00) {
return 1.5;}
176 if(
eta < 1.05) {
return 2.5;}
177 if(
eta < 2.50) {
return 3.5;}
183 axis->SetBinLabel(1,
"CSC C");
184 axis->SetBinLabel(2,
"EC C");
185 axis->SetBinLabel(3,
"TRANS C");
186 axis->SetBinLabel(4,
"BARR C");
187 axis->SetBinLabel(5,
"BARR A");
188 axis->SetBinLabel(6,
"TRANS A");
189 axis->SetBinLabel(7,
"EC A");
190 axis->SetBinLabel(8,
" CSC A");
191 axis->SetTitle(
"Detector #eta Regions");
198 axis->SetBinLabel( 1,
"EC-EC");
199 axis->SetBinLabel( 2,
"EC-BC");
200 axis->SetBinLabel( 3,
"EC-BA");
201 axis->SetBinLabel( 4,
"EC-EA");
202 axis->SetBinLabel( 5,
"BC-EC");
203 axis->SetBinLabel( 6,
"BC-BC");
204 axis->SetBinLabel( 7,
"BC-BA");
205 axis->SetBinLabel( 8,
"BC-EA");
206 axis->SetBinLabel( 9,
"BA-EC");
207 axis->SetBinLabel(10,
"BA-BC");
208 axis->SetBinLabel(11,
"BA-BA");
209 axis->SetBinLabel(12,
"BA-EA");
210 axis->SetBinLabel(13,
"EA-EC");
211 axis->SetBinLabel(14,
"EA-BC");
212 axis->SetBinLabel(15,
"EA-BA");
213 axis->SetBinLabel(16,
"EA-EA");
215 axis->SetTitle(
"#eta Region Permutations [+#mu, -#mu]");
222 axis->SetTitle(
"#phi Sectors 1-16");
223 if(
axis->GetNbins() == 16)
axis->CenterLabels(kTRUE);
226 if(
axis->GetNbins()%16){
230 int nSectorDevisions =
axis->GetNbins()/16;
231 int skip_counter = 1;
233 axis->SetNdivisions(16,nSectorDevisions,0,kFALSE);
236 for (
int i = 1;
i <=
axis->GetNbins();
i++) {
238 if (skip_counter == nSectorDevisions) {
239 axis->SetBinLabel(
i,Form(
"%d",
i/nSectorDevisions));
242 axis->SetBinLabel(
i,
" ");
253 m_M_Mean-> SetBinContent(iBin, pInputHist->GetMean(1));
254 m_M_Mean-> SetBinError( iBin, pInputHist->GetMeanError(1));
255 m_M_Sigma->SetBinContent(iBin, pInputHist->GetRMS(1));