54 bool eventInCycleCut =
false;
62 eventInCycleCut =
true;
71 if(grif1->GetDetector() == grif2->GetDetector() && (
fAddback || grif1->GetCrystal() == grif2->GetCrystal())) {
77 if(angle < fAngles->Rounding())
continue;
78 if(
fFolding && angle > 90.) angle = 180. - angle;
82 if(angleIndex < 0)
continue;
84 double ggTime = TMath::Abs(grif2->GetTime() - grif1->GetTime());
86 fH2[slot].at(
"AngularCorrelation")->Fill(grif1->GetEnergy(), grif2->GetEnergy());
87 fH2[slot].at(
"AngularCorrelation")->Fill(grif2->GetEnergy(), grif1->GetEnergy());
88 fH2[slot].at(Form(
"AngularCorrelation%d", angleIndex))->Fill(grif1->GetEnergy(), grif2->GetEnergy());
89 fH2[slot].at(Form(
"AngularCorrelation%d", angleIndex))->Fill(grif2->GetEnergy(), grif1->GetEnergy());
91 fH2[slot].at(
"AngularCorrelationBG")->Fill(grif1->GetEnergy(), grif2->GetEnergy());
92 fH2[slot].at(
"AngularCorrelationBG")->Fill(grif2->GetEnergy(), grif1->GetEnergy());
93 fH2[slot].at(Form(
"AngularCorrelationBG%d", angleIndex))->Fill(grif1->GetEnergy(), grif2->GetEnergy());
94 fH2[slot].at(Form(
"AngularCorrelationBG%d", angleIndex))->Fill(grif2->GetEnergy(), grif1->GetEnergy());
102 for(
auto g2 = 0; g2 < (
fAddback ? fLastGriffin->GetSuppressedAddbackMultiplicity(fLastBgo) : fLastGriffin->GetSuppressedMultiplicity(fLastBgo)); ++g2) {
103 auto grif2 = (
fAddback ? fLastGriffin->GetSuppressedAddbackHit(g2) : fLastGriffin->GetSuppressedHit(g2));
105 if(
fSingleCrystal && fLastGriffin->GetNSuppressedAddbackFrags(g2) > 1)
continue;
109 if(grif1->GetDetector() == grif2->GetDetector() && (
fAddback || grif1->GetCrystal() == grif2->GetCrystal()))
continue;
112 if(angle < fAngles->Rounding())
continue;
113 if(
fFolding && angle > 90.) angle = 180. - angle;
117 if(angleIndex < 0)
continue;
119 fH2[slot].at(
"AngularCorrelationMixed")->Fill(grif1->GetEnergy(), grif2->GetEnergy());
120 fH2[slot].at(Form(
"AngularCorrelationMixed%d", angleIndex))->Fill(grif1->GetEnergy(), grif2->GetEnergy());
125 if(eventInCycleCut) {