36#include "FairRootManager.h"
37#include "KFParticle/KFParticle.h"
38#include "TDirectory.h"
512 FairRootManager* ioman = FairRootManager::Instance();
513 if (
nullptr == ioman) { Fatal(
"CbmKresEta::Init",
"RootManager not instantised!"); }
515 fMcTracks = (TClonesArray*) ioman->GetObject(
"MCTrack");
516 if (
nullptr ==
fMcTracks) { Fatal(
"CbmKresEta::Init",
"No MCTrack array!"); }
520 if (
nullptr ==
fPrimVertex) { LOG(fatal) <<
"CbmKresEta::Init No PrimaryVertex array!"; }
522 fGlobalTracks = (TClonesArray*) ioman->GetObject(
"GlobalTrack");
523 if (
nullptr ==
fGlobalTracks) { Fatal(
"CbmKresEta::Init",
"No GlobalTrack array!"); }
525 fStsTracks = (TClonesArray*) ioman->GetObject(
"StsTrack");
526 if (
nullptr ==
fStsTracks) { Fatal(
"CbmKresEta::Init",
"No StsTrack array!"); }
529 if (
nullptr ==
fStsTrackMatches) { Fatal(
"CbmKresEta::Init",
"No StsTrackMatch array!"); }
532 if (
nullptr ==
fRichProjections) { Fatal(
"CbmKresEta::Init",
"No RichProjection array!"); }
534 fRichRings = (TClonesArray*) ioman->GetObject(
"RichRing");
535 if (
nullptr ==
fRichRings) { Fatal(
"CbmKresEta::Init",
"No RichRing array!"); }
538 if (
nullptr ==
fRichRingMatches) { Fatal(
"CbmKresEta::Init",
"No RichRingMatch array!"); }
540 fRichHits = (TClonesArray*) ioman->GetObject(
"RichHit");
541 if (
nullptr ==
fRichHits) { Fatal(
"CbmKresEta::Init",
"No RichHit array!"); }
543 fArrayMvdHit = (TClonesArray*) ioman->GetObject(
"MvdHit");
544 if (
nullptr ==
fArrayMvdHit) { Fatal(
"CbmKresEta::Init",
"No MvdHit array!"); }
546 fArrayStsHit = (TClonesArray*) ioman->GetObject(
"StsHit");
547 if (
nullptr ==
fArrayStsHit) { Fatal(
"CbmKresEta::Init",
"No StsHit array!"); }
554void CbmKresEta::Exec(
int fEventNumEta,
double OpeningAngleCut,
double GammaInvMassCut,
int RealPID)
560 Fatal(
"CbmKresConversionManual::Exec",
"No PrimaryVertex array!");
629 for (
Int_t i = 0; i < ngTracks; i++) {
631 if (
nullptr == gTrack)
continue;
635 if (stsInd < 0)
continue;
637 if (stsTrack ==
nullptr)
continue;
639 if (stsMatch ==
nullptr)
continue;
642 if (stsMcTrackId < 0)
continue;
644 if (mcTrack ==
nullptr)
continue;
647 if (richInd < 0 && proj->GetX() > -115 && proj->GetX() < 115
648 && ((proj->GetY() < -120 && proj->GetY() > -200) || (proj->GetY() > 120 && proj->GetY() < 200)))
658 double charge = track_par->GetQp();
660 if (chi2 != chi2)
continue;
661 if (chi2 == 0)
continue;
663 if (chi2 > 3) {
SaveOutsideTracks(mcTrack, stsTrack, charge, stsInd, richInd, stsMcTrackId, Ring); }
664 if (chi2 > 3)
continue;
666 SaveTargetTracks(mcTrack, stsTrack, Momentum, charge, stsInd, richInd, stsMcTrackId, Ring);
704 if (fEventNumEta % numformix == 0) {
723 if (fEventNumEta % numformix == 0) {
730 if (fEventNumEta % numformix == 0) {
742 int InRich =
FindInRich(richInd, stsMcTrackId);
763 int stsInd,
int richInd,
int stsMcTrackId,
CbmRichRing* RING)
765 int InRich =
FindInRich(richInd, stsMcTrackId);
791 if (richMatch !=
nullptr && richMatch->
GetNofLinks() > 0) {
793 if (richMcTrackId > 0) {
794 if (stsMcTrackId == richMcTrackId) {
796 if (mcTrack2 !=
nullptr) {
798 if (TMath::Abs(pdgRICH) == 11) RingsInRich++;
811 if (
nullptr != ring) {
814 for (
int i = 0; i < nofHits; i++) {
815 int hitInd = ring->
GetHit(i);
817 if (
nullptr == hit)
continue;
823 && momentum > 0.2 && momentum < 4.)
832 double FinalAngle = 400;
847 Xpart1 = stsHit1->
GetX();
848 Ypart1 = stsHit1->
GetY();
849 Zpart1 = stsHit1->
GetZ();
854 Xpart1 = mvdHit1->
GetX();
855 Ypart1 = mvdHit1->
GetY();
856 Zpart1 = mvdHit1->
GetZ();
862 Xpart2 = stsHit2->
GetX();
863 Ypart2 = stsHit2->
GetY();
864 Zpart2 = stsHit2->
GetZ();
869 Xpart2 = mvdHit2->
GetX();
870 Ypart2 = mvdHit2->
GetY();
871 Zpart2 = mvdHit2->
GetZ();
875 if (TMath::Abs(Zpart1 - Zpart2) > 2
876 && Zpart1 > Zpart2) {
877 for (
int i = hits1sts - 2; i > -1; i--) {
878 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
881 Xpart1 = stsHit->
GetX();
882 Ypart1 = stsHit->
GetY();
883 Zpart1 = stsHit->
GetZ();
885 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
886 for (
int i = hits1mvd - 2; i > -1; i--) {
887 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
890 Xpart1 = mvdHit->
GetX();
891 Ypart1 = mvdHit->
GetY();
892 Zpart1 = mvdHit->
GetZ();
897 if (TMath::Abs(Zpart1 - Zpart2) > 2
898 && Zpart1 < Zpart2) {
899 for (
int i = hits2sts - 2; i > -1; i--) {
900 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
903 Xpart2 = stsHit->
GetX();
904 Ypart2 = stsHit->
GetY();
905 Zpart2 = stsHit->
GetZ();
907 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
908 for (
int i = hits2mvd - 2; i > -1; i--) {
909 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
912 Xpart2 = mvdHit->
GetX();
913 Ypart2 = mvdHit->
GetY();
914 Zpart2 = mvdHit->
GetZ();
920 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
921 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
924 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
929 double FinalAngle = 400;
944 Xpart1 = mvdHit1->
GetX();
945 Ypart1 = mvdHit1->
GetY();
946 Zpart1 = mvdHit1->
GetZ();
951 Xpart1 = stsHit1->
GetX();
952 Ypart1 = stsHit1->
GetY();
953 Zpart1 = stsHit1->
GetZ();
959 Xpart2 = mvdHit2->
GetX();
960 Ypart2 = mvdHit2->
GetY();
961 Zpart2 = mvdHit2->
GetZ();
966 Xpart2 = stsHit2->
GetX();
967 Ypart2 = stsHit2->
GetY();
968 Zpart2 = stsHit2->
GetZ();
972 if (TMath::Abs(Zpart1 - Zpart2) > 2
973 && Zpart1 < Zpart2) {
974 for (
int i = 1; i < hits1mvd; i++) {
975 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
978 Xpart1 = mvdHit->
GetX();
979 Ypart1 = mvdHit->
GetY();
980 Zpart1 = mvdHit->
GetZ();
982 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
983 for (
int i = 0; i < hits1sts; i++) {
984 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
987 Xpart1 = stsHit->
GetX();
988 Ypart1 = stsHit->
GetY();
989 Zpart1 = stsHit->
GetZ();
994 if (TMath::Abs(Zpart1 - Zpart2) > 2
995 && Zpart1 > Zpart2) {
996 for (
int i = 1; i < hits2mvd; i++) {
997 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1000 Xpart2 = mvdHit->
GetX();
1001 Ypart2 = mvdHit->
GetY();
1002 Zpart2 = mvdHit->
GetZ();
1004 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1005 for (
int i = 0; i < hits2sts; i++) {
1006 if (TMath::Abs(Zpart1 - Zpart2) < 2)
continue;
1009 Xpart2 = stsHit->
GetX();
1010 Ypart2 = stsHit->
GetY();
1011 Zpart2 = stsHit->
GetZ();
1017 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1018 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
1021 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1026 vector<CbmMCTrack*> MCtracks_minus, vector<CbmMCTrack*> MCtracks_plus,
1027 vector<CbmStsTrack*> StsTrack_minus, vector<CbmStsTrack*> StsTrack_plus,
1028 vector<TVector3> Momenta_minus, vector<TVector3> Momenta_plus,
1029 std::vector<int> Rings_minus, std::vector<int> Rings_plus,
1030 std::vector<int> stsIndex_minus, std::vector<int> stsIndex_plus,
1031 vector<CbmRichRing*> richRing_minus, vector<CbmRichRing*> richRing_plus,
1032 vector<Int_t> MCIndex_minus, vector<Int_t> MCIndex_plus)
1034 for (
size_t i = 0; i < Momenta_minus.size(); i++) {
1035 for (
size_t j = 0; j < Momenta_plus.size(); j++) {
1040 TVector3 part1 = Momenta_minus[i];
1041 TVector3 part2 = Momenta_plus[j];
1045 int richcheck_0 = 0;
1046 int richcheck_1 = 0;
1049 richcheck_0 =
CheckIfElectron(richRing_minus[i], Momenta_minus[i].Mag());
1050 richcheck_1 =
CheckIfElectron(richRing_plus[j], Momenta_plus[j].Mag());
1054 richcheck_0 = Rings_minus[i];
1055 richcheck_1 = Rings_plus[j];
1057 int richcheck = richcheck_0 + richcheck_1;
1071 if (mcTrackgrmama->
GetPdgCode() == 221) { TruePair = 1; }
1075 if (TruePair == 1) {
1076 if (richcheck == 0) {
1091 if (richcheck == 1) {
1106 if (richcheck == 2) {
1129 if (richcheck == 0) {
1144 if (richcheck == 1) {
1159 if (richcheck == 2) {
1178 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
1179 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
1181 if (richcheck == 1 && TruePair == 1) {
1185 if (richcheck == 1 && TruePair == 0) {
1194 frefId.push_back(stsIndex_minus[i]);
1195 frefId.push_back(stsIndex_plus[j]);
1197 fMCId.push_back(MCIndex_minus[i]);
1198 fMCId.push_back(MCIndex_plus[j]);
1210 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1218 if (richcheck == 2) {
1226 if (richcheck == 1 || richcheck == 2) {
1238 vector<CbmMCTrack*> MCtracks_minus_Outside,
1239 vector<CbmMCTrack*> MCtracks_plus_Outside,
1240 vector<CbmStsTrack*> StsTrack_minus_Outside,
1241 vector<CbmStsTrack*> StsTrack_plus_Outside, std::vector<int> Rings_minus_Outside,
1242 std::vector<int> Rings_plus_Outside, std::vector<int> stsIndex_minus_Outside,
1243 std::vector<int> stsIndex_plus_Outside, vector<CbmRichRing*> richRing_minus_Outside,
1244 vector<CbmRichRing*> richRing_plus_Outside, vector<Int_t> MCIndex_minus_Outside,
1245 vector<Int_t> MCIndex_plus_Outside)
1247 for (
size_t i = 0; i < StsTrack_minus_Outside.size(); i++) {
1248 for (
size_t j = 0; j < StsTrack_plus_Outside.size(); j++) {
1253 CbmMCTrack* part1MC = MCtracks_minus_Outside[i];
1254 CbmMCTrack* part2MC = MCtracks_plus_Outside[j];
1256 KFParticle electron;
1258 KFParticle positron;
1260 const KFParticle* daughters[2] = {&electron, &positron};
1261 KFParticle intersection;
1262 intersection.Construct(daughters, 2);
1264 if (intersection.GetZ() > 75 || intersection.GetZ() < -5)
continue;
1270 double chi_electron = 0;
1271 double chi_positron = 0;
1273 intersection.GetZ(), chi_electron);
1275 intersection.GetZ(), chi_positron);
1276 if (chi_electron != chi_electron)
continue;
1277 if (chi_electron == 0)
continue;
1278 if (chi_electron > 3)
continue;
1279 if (chi_positron != chi_positron)
continue;
1280 if (chi_positron == 0)
continue;
1281 if (chi_positron > 3)
continue;
1284 int richcheck_0 = 0;
1285 int richcheck_1 = 0;
1288 richcheck_0 =
CheckIfElectron(richRing_minus_Outside[i], part1.Mag());
1290 richcheck = richcheck_0 + richcheck_1;
1294 richcheck = Rings_minus_Outside[i] + Rings_plus_Outside[j];
1309 if (mcTrackgrmama->
GetPdgCode() == 221) { TruePair = 1; }
1313 if (TruePair == 1) {
1314 if (richcheck == 0) {
1329 if (richcheck == 1) {
1344 if (richcheck == 2) {
1373 if (richcheck == 0) {
1388 if (richcheck == 1) {
1403 if (richcheck == 2) {
1422 if (TMath::Abs(OpeningAngle) > AngleCut)
continue;
1423 if (TMath::Abs(InvmassReco) > InvMassCut)
continue;
1426 if (richcheck == 1 && TruePair == 1) {
1430 if (richcheck == 1 && TruePair == 0) {
1440 frefId.push_back(stsIndex_minus_Outside[i]);
1441 frefId.push_back(stsIndex_plus_Outside[j]);
1443 fMCId.push_back(MCIndex_minus_Outside[i]);
1444 fMCId.push_back(MCIndex_plus_Outside[j]);
1456 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1464 if (richcheck == 2) {
1472 if (richcheck == 1 || richcheck == 2) {
1529 vector<vector<int>> StsIndex, vector<vector<int>> MCIndex,
1530 vector<vector<CbmMCTrack*>> GammasMC, vector<TH1*> gg, vector<TH1*> rap_pt_separation)
1532 if (Gammas.size() < 2)
return;
1533 for (
size_t gamma1 = 0; gamma1 < Gammas.size() - 1; gamma1++) {
1534 for (
size_t gamma2 = gamma1 + 1; gamma2 < Gammas.size(); gamma2++) {
1536 TVector3 e11 = Gammas[gamma1][0];
1537 TVector3 e12 = Gammas[gamma1][1];
1538 TVector3 e21 = Gammas[gamma2][0];
1539 TVector3 e22 = Gammas[gamma2][1];
1547 if (StsIndex[gamma1][0] == StsIndex[gamma2][0] || StsIndex[gamma1][0] == StsIndex[gamma2][1]
1548 || StsIndex[gamma1][1] == StsIndex[gamma2][0] || StsIndex[gamma1][1] == StsIndex[gamma2][1])
1552 double openingAngleBetweenGammasReco =
1556 int STSmcId1 = MCIndex[gamma1][0];
1557 int STSmcId2 = MCIndex[gamma1][1];
1558 int STSmcId3 = MCIndex[gamma2][0];
1559 int STSmcId4 = MCIndex[gamma2][1];
1560 if (STSmcId1 != STSmcId2 && STSmcId1 != STSmcId3 && STSmcId1 != STSmcId4 && STSmcId2 != STSmcId3
1561 && STSmcId2 != STSmcId4 && STSmcId3 != STSmcId4) {
1597 if (TrueEta == 1) gg[0]->Fill(openingAngleBetweenGammasReco);
1598 if (TrueEta == 0) gg[1]->Fill(openingAngleBetweenGammasReco);
1600 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40)
continue;
1602 gg[2]->Fill(params.
fMinv);
1604 gg[3]->Fill(TMath::Abs(mcTrack1->
GetPdgCode()));
1605 gg[3]->Fill(TMath::Abs(mcTrack2->
GetPdgCode()));
1606 gg[3]->Fill(TMath::Abs(mcTrack3->
GetPdgCode()));
1607 gg[3]->Fill(TMath::Abs(mcTrack4->
GetPdgCode()));
1610 if (position ==
"Both") {
1612 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation[1]->Fill(params.
fMinv);
1613 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation[2]->Fill(params.
fMinv);
1614 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation[3]->Fill(params.
fMinv);
1615 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation[4]->Fill(params.
fMinv);
1618 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation[5]->Fill(params.
fMinv);
1619 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation[6]->Fill(params.
fMinv);
1620 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation[7]->Fill(params.
fMinv);
1621 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation[8]->Fill(params.
fMinv);
1624 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation[9]->Fill(params.
fMinv);
1625 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation[10]->Fill(params.
fMinv);
1626 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation[11]->Fill(params.
fMinv);
1627 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation[12]->Fill(params.
fMinv);
1630 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation[13]->Fill(params.
fMinv);
1631 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation[14]->Fill(params.
fMinv);
1632 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation[15]->Fill(params.
fMinv);
1635 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation[16]->Fill(params.
fMinv);
1636 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation[17]->Fill(params.
fMinv);
1643 gg[8]->Fill(params.
fMinv);
1644 gg[10]->Fill(openingAngleBetweenGammasReco);
1650 gg[9]->Fill(params.
fMinv);
1651 gg[11]->Fill(openingAngleBetweenGammasReco);
1661 cout <<
"Mixing for Eta(Target) - nof entries " << nof_Target << endl;
1662 for (
Int_t a = 0; a < nof_Target - 1; a++) {
1663 for (
Int_t b = a + 1; b < nof_Target; b++) {
1671 double openingAngleBetweenGammasReco =
1673 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40)
continue;
1690 cout <<
"Mixing for Eta(Outside) - nof entries " << nof_Outside << endl;
1691 for (
Int_t a = 0; a < nof_Outside - 1; a++) {
1692 for (
Int_t b = a + 1; b < nof_Outside; b++) {
1700 double openingAngleBetweenGammasReco =
1702 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40)
continue;
1715 vector<TH1*> rap_pt_separation_two)
1719 cout <<
"Mixing for Eta(Both) - nof entries " << nof_Both << endl;
1720 for (
Int_t a = 0; a < nof_Both - 1; a++) {
1721 for (
Int_t b = a + 1; b < nof_Both; b++) {
1729 double openingAngleBetweenGammasReco =
1731 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40)
continue;
1742 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_all[18]->Fill(params.
fMinv);
1743 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_all[19]->Fill(params.
fMinv);
1744 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_all[20]->Fill(params.
fMinv);
1745 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_all[21]->Fill(params.
fMinv);
1748 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_all[22]->Fill(params.
fMinv);
1749 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_all[23]->Fill(params.
fMinv);
1750 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_all[24]->Fill(params.
fMinv);
1751 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_all[25]->Fill(params.
fMinv);
1754 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_all[26]->Fill(params.
fMinv);
1755 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_all[27]->Fill(params.
fMinv);
1756 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_all[28]->Fill(params.
fMinv);
1757 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_all[29]->Fill(params.
fMinv);
1760 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_all[30]->Fill(params.
fMinv);
1761 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_all[31]->Fill(params.
fMinv);
1762 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_all[32]->Fill(params.
fMinv);
1765 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_all[33]->Fill(params.
fMinv);
1766 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_all[34]->Fill(params.
fMinv);
1773 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_onetwo[18]->Fill(params.
fMinv);
1774 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_onetwo[19]->Fill(params.
fMinv);
1775 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_onetwo[20]->Fill(params.
fMinv);
1776 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_onetwo[21]->Fill(params.
fMinv);
1779 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_onetwo[22]->Fill(params.
fMinv);
1780 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_onetwo[23]->Fill(params.
fMinv);
1781 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_onetwo[24]->Fill(params.
fMinv);
1782 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_onetwo[25]->Fill(params.
fMinv);
1785 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_onetwo[26]->Fill(params.
fMinv);
1786 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_onetwo[27]->Fill(params.
fMinv);
1787 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_onetwo[28]->Fill(params.
fMinv);
1788 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_onetwo[29]->Fill(params.
fMinv);
1791 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_onetwo[30]->Fill(params.
fMinv);
1792 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_onetwo[31]->Fill(params.
fMinv);
1793 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_onetwo[32]->Fill(params.
fMinv);
1796 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_onetwo[33]->Fill(params.
fMinv);
1797 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_onetwo[34]->Fill(params.
fMinv);
1804 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_two[18]->Fill(params.
fMinv);
1805 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_two[19]->Fill(params.
fMinv);
1806 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_two[20]->Fill(params.
fMinv);
1807 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_two[21]->Fill(params.
fMinv);
1810 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_two[22]->Fill(params.
fMinv);
1811 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_two[23]->Fill(params.
fMinv);
1812 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_two[24]->Fill(params.
fMinv);
1813 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_two[25]->Fill(params.
fMinv);
1816 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_two[26]->Fill(params.
fMinv);
1817 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_two[27]->Fill(params.
fMinv);
1818 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_two[28]->Fill(params.
fMinv);
1819 if (params.
fPt > 1.2 && params.
fPt <= 1.6) rap_pt_separation_two[29]->Fill(params.
fMinv);
1822 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_two[30]->Fill(params.
fMinv);
1823 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_two[31]->Fill(params.
fMinv);
1824 if (params.
fPt > 0.8 && params.
fPt <= 1.2) rap_pt_separation_two[32]->Fill(params.
fMinv);
1827 if (params.
fPt > 0.0 && params.
fPt <= 0.4) rap_pt_separation_two[33]->Fill(params.
fMinv);
1828 if (params.
fPt > 0.4 && params.
fPt <= 0.8) rap_pt_separation_two[34]->Fill(params.
fMinv);
1838 gDirectory->mkdir(
"Eta");
1839 gDirectory->cd(
"Eta");
1841 gDirectory->mkdir(
"Cuts_for_Eta_gg");
1842 gDirectory->cd(
"Cuts_for_Eta_gg");
1844 gDirectory->mkdir(
"Cuts_for_Eta_gg_Target");
1845 gDirectory->cd(
"Cuts_for_Eta_gg_Target");
1846 gDirectory->mkdir(
"Cuts_for_Eta_gg_zeroInRich_Target");
1847 gDirectory->cd(
"Cuts_for_Eta_gg_zeroInRich_Target");
1851 gDirectory->cd(
"..");
1852 gDirectory->mkdir(
"Cuts_for_Eta_gg_oneInRich_Target");
1853 gDirectory->cd(
"Cuts_for_Eta_gg_oneInRich_Target");
1857 gDirectory->cd(
"..");
1858 gDirectory->mkdir(
"Cuts_for_Eta_gg_twoInRich_Target");
1859 gDirectory->cd(
"Cuts_for_Eta_gg_twoInRich_Target");
1863 gDirectory->cd(
"..");
1867 gDirectory->cd(
"..");
1869 gDirectory->mkdir(
"Cuts_for_Eta_gg_Outside");
1870 gDirectory->cd(
"Cuts_for_Eta_gg_Outside");
1871 gDirectory->mkdir(
"Cuts_for_Eta_gg_zeroInRich_Outside");
1872 gDirectory->cd(
"Cuts_for_Eta_gg_zeroInRich_Outside");
1876 gDirectory->cd(
"..");
1877 gDirectory->mkdir(
"Cuts_for_Eta_gg_oneInRich_Outside");
1878 gDirectory->cd(
"Cuts_for_Eta_gg_oneInRich_Outside");
1882 gDirectory->cd(
"..");
1883 gDirectory->mkdir(
"Cuts_for_Eta_gg_twoInRich_Outside");
1884 gDirectory->cd(
"Cuts_for_Eta_gg_twoInRich_Outside");
1888 gDirectory->cd(
"..");
1892 gDirectory->cd(
"..");
1894 gDirectory->mkdir(
"Cuts_for_Eta_gg_Both");
1895 gDirectory->cd(
"Cuts_for_Eta_gg_Both");
1896 gDirectory->mkdir(
"Cuts_for_Eta_gg_zeroInRich_Both");
1897 gDirectory->cd(
"Cuts_for_Eta_gg_zeroInRich_Both");
1901 gDirectory->cd(
"..");
1902 gDirectory->mkdir(
"Cuts_for_Eta_gg_oneInRich_Both");
1903 gDirectory->cd(
"Cuts_for_Eta_gg_oneInRich_Both");
1907 gDirectory->cd(
"..");
1908 gDirectory->mkdir(
"Cuts_for_Eta_gg_twoInRich_Both");
1909 gDirectory->cd(
"Cuts_for_Eta_gg_twoInRich_Both");
1913 gDirectory->cd(
"..");
1917 gDirectory->cd(
"..");
1919 gDirectory->cd(
"..");
1922 gDirectory->mkdir(
"Eta_gg_Target");
1923 gDirectory->cd(
"Eta_gg_Target");
1925 gDirectory->mkdir(
"all");
1926 gDirectory->cd(
"all");
1930 gDirectory->cd(
"..");
1932 gDirectory->mkdir(
"onetwo");
1933 gDirectory->cd(
"onetwo");
1937 gDirectory->cd(
"..");
1939 gDirectory->mkdir(
"two");
1940 gDirectory->cd(
"two");
1944 gDirectory->cd(
"..");
1946 gDirectory->cd(
"..");
1949 gDirectory->mkdir(
"Eta_gg_Outside");
1950 gDirectory->cd(
"Eta_gg_Outside");
1952 gDirectory->mkdir(
"all");
1953 gDirectory->cd(
"all");
1957 gDirectory->cd(
"..");
1959 gDirectory->mkdir(
"onetwo");
1960 gDirectory->cd(
"onetwo");
1964 gDirectory->cd(
"..");
1966 gDirectory->mkdir(
"two");
1967 gDirectory->cd(
"two");
1971 gDirectory->cd(
"..");
1973 gDirectory->cd(
"..");
1976 gDirectory->mkdir(
"Eta_gg_Both");
1977 gDirectory->cd(
"Eta_gg_Both");
1979 gDirectory->mkdir(
"all");
1980 gDirectory->cd(
"all");
1981 gDirectory->mkdir(
"multi_all_Both");
1982 gDirectory->cd(
"multi_all_Both");
1986 gDirectory->cd(
"..");
1990 gDirectory->cd(
"..");
1992 gDirectory->mkdir(
"onetwo");
1993 gDirectory->cd(
"onetwo");
1994 gDirectory->mkdir(
"multi_onetwo_Both");
1995 gDirectory->cd(
"multi_onetwo_Both");
1999 gDirectory->cd(
"..");
2003 gDirectory->cd(
"..");
2005 gDirectory->mkdir(
"two");
2006 gDirectory->cd(
"two");
2007 gDirectory->mkdir(
"multi_two_Both");
2008 gDirectory->cd(
"multi_two_Both");
2012 gDirectory->cd(
"..");
2016 gDirectory->cd(
"..");
2018 gDirectory->cd(
"..");
2026 gDirectory->cd(
"..");
2032 new TH1D(
"Particle_pull_px_Target",
"Target track, pull Px: mc - reco; P_{x}(mc) - P_{x}(reco)", 100, -0.5, 0.5);
2035 new TH1D(
"Particle_pull_py_Target",
"Target track, pull Px: mc - reco; P_{y}(mc) - P_{y}(reco)", 100, -0.5, 0.5);
2038 new TH1D(
"Particle_pull_pz_Target",
"Target track, pull Px: mc - reco; P_{z}(mc) - P_{z}(reco)", 100, -0.5, 0.5);
2041 new TH1D(
"Particle_pull_px_Outside",
"Outside track, pull Px: mc - reco; P_{x}(mc) - P_{x}(reco)", 100, -0.5, 0.5);
2044 new TH1D(
"Particle_pull_py_Outside",
"Outside track, pull Px: mc - reco; P_{y}(mc) - P_{y}(reco)", 100, -0.5, 0.5);
2047 new TH1D(
"Particle_pull_pz_Outside",
"Outside track, pull Px: mc - reco; P_{z}(mc) - P_{z}(reco)", 100, -0.5, 0.5);
2050 new TH1D(
"Particle_pull_X_Outside",
"Outside track, pull Px: mc - reco; X(mc) - X(reco)", 500, -5, 5);
2053 new TH1D(
"Particle_pull_Y_Outside",
"Outside track, pull Px: mc - reco; Y(mc) - Y(reco)", 500, -5, 5);
2056 new TH1D(
"Particle_pull_Z_Outside",
"Outside track, pull Px: mc - reco; Z(mc) - Z(reco)", 500, -5, 5);
2061 new TH1D(
"GammasInvMass_fromEta_beforeCuts_zeroInRich_Target",
2062 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Target; invariant "
2063 "mass in GeV/c^{2}",
2067 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target",
2068 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target; opening "
2073 new TH1D(
"GammasMomentum_fromEta_beforeCuts_zeroInRich_Target",
2074 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Target; p in GeV/c", 400, 0, 4);
2077 new TH1D(
"GammasRapidity_fromEta_beforeCuts_zeroInRich_Target",
2078 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Target; rapidity", 400, 0., 4.);
2081 new TH1D(
"PlaneAngles_last_fromEta_zeroInRich_Target",
2082 "PlaneAngles_last_fromEta_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2085 new TH1D(
"PlaneAngles_first_fromEta_zeroInRich_Target",
2086 "PlaneAngles_first_fromEta_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2089 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target",
2090 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target; invariant "
2091 "mass in GeV/c^{2}",
2095 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target",
2096 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target; opening "
2101 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target",
2102 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target; p in GeV/c", 400, 0, 4);
2105 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target",
2106 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target; rapidity", 400, 0., 4.);
2109 new TH1D(
"PlaneAngles_last_wrongpairs_zeroInRich_Target",
2110 "PlaneAngles_last_wrongpairs_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2113 new TH1D(
"PlaneAngles_first_wrongpairs_zeroInRich_Target",
2114 "PlaneAngles_first_wrongpairs_zeroInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2118 new TH1D(
"GammasInvMass_fromEta_beforeCuts_oneInRich_Target",
2119 "GammasInvMass_fromEta_beforeCuts_oneInRich_Target; invariant "
2120 "mass in GeV/c^{2}",
2124 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target",
2125 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target; opening "
2130 new TH1D(
"GammasMomentum_fromEta_beforeCuts_oneInRich_Target",
2131 "GammasMomentum_fromEta_beforeCuts_oneInRich_Target; p in GeV/c", 400, 0, 4);
2134 new TH1D(
"GammasRapidity_fromEta_beforeCuts_oneInRich_Target",
2135 "GammasRapidity_fromEta_beforeCuts_oneInRich_Target; rapidity", 400, 0., 4.);
2138 new TH1D(
"PlaneAngles_last_fromEta_oneInRich_Target",
2139 "PlaneAngles_last_fromEta_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2142 new TH1D(
"PlaneAngles_first_fromEta_oneInRich_Target",
2143 "PlaneAngles_first_fromEta_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2146 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target",
2147 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target; invariant "
2148 "mass in GeV/c^{2}",
2152 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target",
2153 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target; opening "
2158 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target",
2159 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target; p in GeV/c", 400, 0, 4);
2162 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target",
2163 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target; rapidity", 400, 0., 4.);
2166 new TH1D(
"PlaneAngles_last_wrongpairs_oneInRich_Target",
2167 "PlaneAngles_last_wrongpairs_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2170 new TH1D(
"PlaneAngles_first_wrongpairs_oneInRich_Target",
2171 "PlaneAngles_first_wrongpairs_oneInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2175 new TH1D(
"GammasInvMass_fromEta_beforeCuts_twoInRich_Target",
2176 "GammasInvMass_fromEta_beforeCuts_twoInRich_Target; invariant "
2177 "mass in GeV/c^{2}",
2181 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target",
2182 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target; opening "
2187 new TH1D(
"GammasMomentum_fromEta_beforeCuts_twoInRich_Target",
2188 "GammasMomentum_fromEta_beforeCuts_twoInRich_Target; p in GeV/c", 400, 0, 4);
2191 new TH1D(
"GammasRapidity_fromEta_beforeCuts_twoInRich_Target",
2192 "GammasRapidity_fromEta_beforeCuts_twoInRich_Target; rapidity", 400, 0., 4.);
2195 new TH1D(
"PlaneAngles_last_fromEta_twoInRich_Target",
2196 "PlaneAngles_last_fromEta_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2199 new TH1D(
"PlaneAngles_first_fromEta_twoInRich_Target",
2200 "PlaneAngles_first_fromEta_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2203 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target",
2204 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target; invariant "
2205 "mass in GeV/c^{2}",
2209 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target",
2210 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target; opening "
2215 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target",
2216 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target; p in GeV/c", 400, 0, 4);
2219 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target",
2220 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target; rapidity", 400, 0., 4.);
2223 new TH1D(
"PlaneAngles_last_wrongpairs_twoInRich_Target",
2224 "PlaneAngles_last_wrongpairs_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2227 new TH1D(
"PlaneAngles_first_wrongpairs_twoInRich_Target",
2228 "PlaneAngles_first_wrongpairs_twoInRich_Target; #theta angle in degree;#", 720, -1., 179.);
2234 new TH1D(
"GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside",
2235 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside; invariant "
2236 "mass in GeV/c^{2}",
2240 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside",
2241 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside; opening "
2246 new TH1D(
"GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside",
2247 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside; p in GeV/c", 400, 0, 4);
2250 new TH1D(
"GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside",
2251 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside; rapidity", 400, 0., 4.);
2254 new TH1D(
"PlaneAngles_last_fromEta_zeroInRich_Outside",
2255 "PlaneAngles_last_fromEta_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2258 new TH1D(
"PlaneAngles_first_fromEta_zeroInRich_Outside",
2259 "PlaneAngles_first_fromEta_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2262 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside",
2263 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside; "
2264 "invariant mass in GeV/c^{2}",
2268 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside",
2269 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside; "
2270 "opening angle in degree",
2274 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside",
2275 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside; p in GeV/c", 400, 0, 4);
2278 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside",
2279 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside; rapidity", 400, 0., 4.);
2282 new TH1D(
"PlaneAngles_last_wrongpairs_zeroInRich_Outside",
2283 "PlaneAngles_last_wrongpairs_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2286 new TH1D(
"PlaneAngles_first_wrongpairs_zeroInRich_Outside",
2287 "PlaneAngles_first_wrongpairs_zeroInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2291 new TH1D(
"GammasInvMass_fromEta_beforeCuts_oneInRich_Outside",
2292 "GammasInvMass_fromEta_beforeCuts_oneInRich_Outside; invariant "
2293 "mass in GeV/c^{2}",
2297 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside",
2298 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside; opening "
2303 new TH1D(
"GammasMomentum_fromEta_beforeCuts_oneInRich_Outside",
2304 "GammasMomentum_fromEta_beforeCuts_oneInRich_Outside; p in GeV/c", 400, 0, 4);
2307 new TH1D(
"GammasRapidity_fromEta_beforeCuts_oneInRich_Outside",
2308 "GammasRapidity_fromEta_beforeCuts_oneInRich_Outside; rapidity", 400, 0., 4.);
2311 new TH1D(
"PlaneAngles_last_fromEta_oneInRich_Outside",
2312 "PlaneAngles_last_fromEta_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2315 new TH1D(
"PlaneAngles_first_fromEta_oneInRich_Outside",
2316 "PlaneAngles_first_fromEta_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2319 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside",
2320 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside; invariant "
2321 "mass in GeV/c^{2}",
2325 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside",
2326 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside; opening "
2331 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside",
2332 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside; p in GeV/c", 400, 0, 4);
2335 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside",
2336 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside; rapidity", 400, 0., 4.);
2339 new TH1D(
"PlaneAngles_last_wrongpairs_oneInRich_Outside",
2340 "PlaneAngles_last_wrongpairs_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2343 new TH1D(
"PlaneAngles_first_wrongpairs_oneInRich_Outside",
2344 "PlaneAngles_first_wrongpairs_oneInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2348 new TH1D(
"GammasInvMass_fromEta_beforeCuts_twoInRich_Outside",
2349 "GammasInvMass_fromEta_beforeCuts_twoInRich_Outside; invariant "
2350 "mass in GeV/c^{2}",
2354 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside",
2355 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside; opening "
2360 new TH1D(
"GammasMomentum_fromEta_beforeCuts_twoInRich_Outside",
2361 "GammasMomentum_fromEta_beforeCuts_twoInRich_Outside; p in GeV/c", 400, 0, 4);
2364 new TH1D(
"GammasRapidity_fromEta_beforeCuts_twoInRich_Outside",
2365 "GammasRapidity_fromEta_beforeCuts_twoInRich_Outside; rapidity", 400, 0., 4.);
2368 new TH1D(
"PlaneAngles_last_fromEta_twoInRich_Outside",
2369 "PlaneAngles_last_fromEta_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2372 new TH1D(
"PlaneAngles_first_fromEta_twoInRich_Outside",
2373 "PlaneAngles_first_fromEta_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2376 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside",
2377 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside; invariant "
2378 "mass in GeV/c^{2}",
2382 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside",
2383 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside; opening "
2388 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside",
2389 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside; p in GeV/c", 400, 0, 4);
2392 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside",
2393 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside; rapidity", 400, 0., 4.);
2396 new TH1D(
"PlaneAngles_last_wrongpairs_twoInRich_Outside",
2397 "PlaneAngles_last_wrongpairs_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2400 new TH1D(
"PlaneAngles_first_wrongpairs_twoInRich_Outside",
2401 "PlaneAngles_first_wrongpairs_twoInRich_Outside; #theta angle in degree;#", 720, -1., 179.);
2407 new TH1D(
"GammasInvMass_fromEta_beforeCuts_zeroInRich_Both",
2408 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Both; invariant mass "
2413 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both",
2414 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both; opening "
2419 new TH1D(
"GammasMomentum_fromEta_beforeCuts_zeroInRich_Both",
2420 "GammasMomentum_fromEta_beforeCuts_zeroInRich_Both; p in GeV/c", 400, 0, 4);
2423 new TH1D(
"GammasRapidity_fromEta_beforeCuts_zeroInRich_Both",
2424 "GammasRapidity_fromEta_beforeCuts_zeroInRich_Both; rapidity", 400, 0., 4.);
2427 new TH1D(
"PlaneAngles_last_fromEta_zeroInRich_Both",
2428 "PlaneAngles_last_fromEta_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2431 new TH1D(
"PlaneAngles_first_fromEta_zeroInRich_Both",
2432 "PlaneAngles_first_fromEta_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2435 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both",
2436 "GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both; invariant "
2437 "mass in GeV/c^{2}",
2441 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both",
2442 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both; opening "
2447 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both",
2448 "GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both; p in GeV/c", 400, 0, 4);
2451 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both",
2452 "GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both; rapidity", 400, 0., 4.);
2455 new TH1D(
"PlaneAngles_last_wrongpairs_zeroInRich_Both",
2456 "PlaneAngles_last_wrongpairs_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2459 new TH1D(
"PlaneAngles_first_wrongpairs_zeroInRich_Both",
2460 "PlaneAngles_first_wrongpairs_zeroInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2464 new TH1D(
"GammasInvMass_fromEta_beforeCuts_oneInRich_Both",
2465 "GammasInvMass_fromEta_beforeCuts_oneInRich_Both; invariant mass "
2470 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both",
2471 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both; opening angle "
2476 new TH1D(
"GammasMomentum_fromEta_beforeCuts_oneInRich_Both",
2477 "GammasMomentum_fromEta_beforeCuts_oneInRich_Both; p in GeV/c", 400, 0, 4);
2480 new TH1D(
"GammasRapidity_fromEta_beforeCuts_oneInRich_Both",
2481 "GammasRapidity_fromEta_beforeCuts_oneInRich_Both; rapidity", 400, 0., 4.);
2484 new TH1D(
"PlaneAngles_last_fromEta_oneInRich_Both",
2485 "PlaneAngles_last_fromEta_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2488 new TH1D(
"PlaneAngles_first_fromEta_oneInRich_Both",
2489 "PlaneAngles_first_fromEta_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2492 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both",
2493 "GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both; invariant "
2494 "mass in GeV/c^{2}",
2498 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both",
2499 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both; opening "
2504 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both",
2505 "GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both; p in GeV/c", 400, 0, 4);
2508 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both",
2509 "GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both; rapidity", 400, 0., 4.);
2512 new TH1D(
"PlaneAngles_last_wrongpairs_oneInRich_Both",
2513 "PlaneAngles_last_wrongpairs_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2516 new TH1D(
"PlaneAngles_first_wrongpairs_oneInRich_Both",
2517 "PlaneAngles_first_wrongpairs_oneInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2520 new TH1D(
"PlaneAngles_last_fromEta_afterCuts_oneInRich_Both",
2521 "PlaneAngles_last_fromEta_afterCuts_oneInRich_Both; #theta angle "
2526 new TH1D(
"PlaneAngles_first_fromEta_afterCuts_oneInRich_Both",
2527 "PlaneAngles_first_fromEta_afterCuts_oneInRich_Both; #theta angle "
2532 new TH1D(
"PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both",
2533 "PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both; #theta "
2534 "angle in degree;#",
2538 new TH1D(
"PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both",
2539 "PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both; #theta "
2540 "angle in degree;#",
2545 new TH1D(
"GammasInvMass_fromEta_beforeCuts_twoInRich_Both",
2546 "GammasInvMass_fromEta_beforeCuts_twoInRich_Both; invariant mass "
2551 new TH1D(
"GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both",
2552 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both; opening angle "
2557 new TH1D(
"GammasMomentum_fromEta_beforeCuts_twoInRich_Both",
2558 "GammasMomentum_fromEta_beforeCuts_twoInRich_Both; p in GeV/c", 400, 0, 4);
2561 new TH1D(
"GammasRapidity_fromEta_beforeCuts_twoInRich_Both",
2562 "GammasRapidity_fromEta_beforeCuts_twoInRich_Both; rapidity", 400, 0., 4.);
2565 new TH1D(
"PlaneAngles_last_fromEta_twoInRich_Both",
2566 "PlaneAngles_last_fromEta_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2569 new TH1D(
"PlaneAngles_first_fromEta_twoInRich_Both",
2570 "PlaneAngles_first_fromEta_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2573 new TH1D(
"GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both",
2574 "GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both; invariant "
2575 "mass in GeV/c^{2}",
2579 new TH1D(
"GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both",
2580 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both; opening "
2585 new TH1D(
"GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both",
2586 "GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both; p in GeV/c", 400, 0, 4);
2589 new TH1D(
"GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both",
2590 "GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both; rapidity", 400, 0., 4.);
2593 new TH1D(
"PlaneAngles_last_wrongpairs_twoInRich_Both",
2594 "PlaneAngles_last_wrongpairs_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2597 new TH1D(
"PlaneAngles_first_wrongpairs_twoInRich_Both",
2598 "PlaneAngles_first_wrongpairs_twoInRich_Both; #theta angle in degree;#", 720, -1., 179.);
2605 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_all_Target",
2606 "AngleBetweenGammas_trueEta_before_cuts_all_Target (between "
2607 "#gamma#gamma from #eta); angle [deg]",
2611 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_all_Target",
2612 "AngleBetweenGammas_wrongEta_before_cuts_all_Target (between "
2613 "#gamma#gamma); angle [deg]",
2617 new TH1D(
"InvMass_GammaCombinations_all_Target",
2618 "InvMass_GammaCombinations_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2623 new TH2D(
"TrueEta_pt_vs_rap_all_Target",
"TrueEta_pt_vs_rap_all_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2627 new TH2D(
"WrongEta_pt_vs_rap_all_Target",
"WrongEta_pt_vs_rap_all_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2631 new TH2D(
"TrueEta_pt_vs_rap_est_all_Target",
"TrueEta_pt_vs_rap_est_all_Target; rapidity y; p_{t} in GeV/c ", 10,
2635 new TH2D(
"WrongEta_pt_vs_rap_est_all_Target",
"WrongEta_pt_vs_rap_est_all_Target; rapidity y; p_{t} in GeV/c ", 10,
2639 new TH1D(
"TrueEta_InvMass_after_cuts_all_Target",
2640 "TrueEta_InvMass_after_cuts_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2643 new TH1D(
"WrongEta_InvMass_after_cuts_all_Target",
2644 "WrongEta_InvMass_after_cuts_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2647 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_all_Target",
2648 "TrueEta_AngleBetweenGammas_after_cuts_all_Target (between "
2649 "#gamma#gamma from #eta); angle [deg]",
2653 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_all_Target",
2654 "WrongEta_AngleBetweenGammas_after_cuts_all_Target (between "
2655 "#gamma#gamma); angle [deg]",
2659 new TH1D(
"EMT_Eta_InvMass_all_Target",
"EMT_Eta_InvMass_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2663 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_onetwo_Target",
2664 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Target (between "
2665 "#gamma#gamma from #eta); angle [deg]",
2669 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target",
2670 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target (between "
2671 "#gamma#gamma); angle [deg]",
2675 new TH1D(
"InvMass_GammaCombinations_onetwo_Target",
2676 "InvMass_GammaCombinations_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2679 new TH1D(
"Particles_PDG_onetwo_Target",
"Particles_PDG_onetwo_Target ; PDG", 2300, 0, 2300);
2682 new TH2D(
"TrueEta_pt_vs_rap_onetwo_Target",
"TrueEta_pt_vs_rap_onetwo_Target; rapidity y; p_{t} in GeV/c ", 90, -2.,
2686 new TH2D(
"WrongEta_pt_vs_rap_onetwo_Target",
"WrongEta_pt_vs_rap_onetwo_Target; rapidity y; p_{t} in GeV/c ", 90,
2687 -2., 7., 60, -1., 5.);
2690 new TH2D(
"TrueEta_pt_vs_rap_est_onetwo_Target",
"TrueEta_pt_vs_rap_est_onetwo_Target; rapidity y; p_{t} in GeV/c ",
2691 10, 0., 4., 5, 0., 2.);
2694 new TH2D(
"WrongEta_pt_vs_rap_est_onetwo_Target",
2695 "WrongEta_pt_vs_rap_est_onetwo_Target; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 5, 0., 2.);
2698 new TH1D(
"TrueEta_InvMass_after_cuts_onetwo_Target",
2699 "TrueEta_InvMass_after_cuts_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2702 new TH1D(
"WrongEta_InvMass_after_cuts_onetwo_Target",
2703 "WrongEta_InvMass_after_cuts_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2706 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target",
2707 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target (between "
2708 "#gamma#gamma from #eta); angle [deg]",
2712 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target",
2713 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target (between "
2714 "#gamma#gamma); angle [deg]",
2718 "EMT_Eta_InvMass_onetwo_Target",
"EMT_Eta_InvMass_onetwo_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2722 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_two_Target",
2723 "AngleBetweenGammas_trueEta_before_cuts_two_Target (between "
2724 "#gamma#gamma from #eta); angle [deg]",
2728 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_two_Target",
2729 "AngleBetweenGammas_wrongEta_before_cuts_two_Target (between "
2730 "#gamma#gamma); angle [deg]",
2734 new TH1D(
"InvMass_GammaCombinations_two_Target",
2735 "InvMass_GammaCombinations_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2740 new TH2D(
"TrueEta_pt_vs_rap_two_Target",
"TrueEta_pt_vs_rap_two_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2744 new TH2D(
"WrongEta_pt_vs_rap_two_Target",
"WrongEta_pt_vs_rap_two_Target; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2748 new TH2D(
"TrueEta_pt_vs_rap_est_two_Target",
"TrueEta_pt_vs_rap_est_two_Target; rapidity y; p_{t} in GeV/c ", 10,
2752 new TH2D(
"WrongEta_pt_vs_rap_est_two_Target",
"WrongEta_pt_vs_rap_est_two_Target; rapidity y; p_{t} in GeV/c ", 10,
2756 new TH1D(
"TrueEta_InvMass_after_cuts_two_Target",
2757 "TrueEta_InvMass_after_cuts_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2760 new TH1D(
"WrongEta_InvMass_after_cuts_two_Target",
2761 "WrongEta_InvMass_after_cuts_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2764 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_two_Target",
2765 "TrueEta_AngleBetweenGammas_after_cuts_two_Target (between "
2766 "#gamma#gamma from #eta); angle [deg]",
2770 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_two_Target",
2771 "WrongEta_AngleBetweenGammas_after_cuts_two_Target (between "
2772 "#gamma#gamma); angle [deg]",
2776 new TH1D(
"EMT_Eta_InvMass_two_Target",
"EMT_Eta_InvMass_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2784 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_all_Outside",
2785 "AngleBetweenGammas_trueEta_before_cuts_all_Outside (between "
2786 "#gamma#gamma from #eta); angle [deg]",
2790 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_all_Outside",
2791 "AngleBetweenGammas_wrongEta_before_cuts_all_Outside (between "
2792 "#gamma#gamma); angle [deg]",
2796 new TH1D(
"InvMass_GammaCombinations_all_Outside",
2797 "InvMass_GammaCombinations_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2802 new TH2D(
"TrueEta_pt_vs_rap_all_Outside",
"TrueEta_pt_vs_rap_all_Outside; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2806 new TH2D(
"WrongEta_pt_vs_rap_all_Outside",
"WrongEta_pt_vs_rap_all_Outside; rapidity y; p_{t} in GeV/c ", 90, -2.,
2810 new TH2D(
"TrueEta_pt_vs_rap_est_all_Outside",
"TrueEta_pt_vs_rap_est_all_Outside; rapidity y; p_{t} in GeV/c ", 10,
2814 new TH2D(
"WrongEta_pt_vs_rap_est_all_Outside",
"WrongEta_pt_vs_rap_est_all_Outside; rapidity y; p_{t} in GeV/c ",
2815 10, 0., 4., 5, 0., 2.);
2818 new TH1D(
"TrueEta_InvMass_after_cuts_all_Outside",
2819 "TrueEta_InvMass_after_cuts_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2822 new TH1D(
"WrongEta_InvMass_after_cuts_all_Outside",
2823 "WrongEta_InvMass_after_cuts_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2826 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_all_Outside",
2827 "TrueEta_AngleBetweenGammas_after_cuts_all_Outside (between "
2828 "#gamma#gamma from #eta); angle [deg]",
2832 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_all_Outside",
2833 "WrongEta_AngleBetweenGammas_after_cuts_all_Outside (between "
2834 "#gamma#gamma); angle [deg]",
2838 "EMT_Eta_InvMass_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2842 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside",
2843 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside (between "
2844 "#gamma#gamma from #eta); angle [deg]",
2848 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside",
2849 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside (between "
2850 "#gamma#gamma); angle [deg]",
2854 new TH1D(
"InvMass_GammaCombinations_onetwo_Outside",
2855 "InvMass_GammaCombinations_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2858 new TH1D(
"Particles_PDG_onetwo_Outside",
"Particles_PDG_onetwo_Outside ; PDG", 2300, 0, 2300);
2861 new TH2D(
"TrueEta_pt_vs_rap_onetwo_Outside",
"TrueEta_pt_vs_rap_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 90,
2862 -2., 7., 60, -1., 5.);
2865 new TH2D(
"WrongEta_pt_vs_rap_onetwo_Outside",
"WrongEta_pt_vs_rap_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 90,
2866 -2., 7., 60, -1., 5.);
2869 new TH2D(
"TrueEta_pt_vs_rap_est_onetwo_Outside",
2870 "TrueEta_pt_vs_rap_est_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 5, 0., 2.);
2873 new TH2D(
"WrongEta_pt_vs_rap_est_onetwo_Outside",
2874 "WrongEta_pt_vs_rap_est_onetwo_Outside; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 5, 0., 2.);
2877 new TH1D(
"TrueEta_InvMass_after_cuts_onetwo_Outside",
2878 "TrueEta_InvMass_after_cuts_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2881 new TH1D(
"WrongEta_InvMass_after_cuts_onetwo_Outside",
2882 "WrongEta_InvMass_after_cuts_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2885 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside",
2886 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside (between "
2887 "#gamma#gamma from #eta); angle [deg]",
2891 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside",
2892 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside (between "
2893 "#gamma#gamma); angle [deg]",
2897 "EMT_Eta_InvMass_onetwo_Outside",
"EMT_Eta_InvMass_onetwo_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2901 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_two_Outside",
2902 "AngleBetweenGammas_trueEta_before_cuts_two_Outside (between "
2903 "#gamma#gamma from #eta); angle [deg]",
2907 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_two_Outside",
2908 "AngleBetweenGammas_wrongEta_before_cuts_two_Outside (between "
2909 "#gamma#gamma); angle [deg]",
2913 new TH1D(
"InvMass_GammaCombinations_two_Outside",
2914 "InvMass_GammaCombinations_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2919 new TH2D(
"TrueEta_pt_vs_rap_two_Outside",
"TrueEta_pt_vs_rap_two_Outside; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2923 new TH2D(
"WrongEta_pt_vs_rap_two_Outside",
"WrongEta_pt_vs_rap_two_Outside; rapidity y; p_{t} in GeV/c ", 90, -2.,
2927 new TH2D(
"TrueEta_pt_vs_rap_est_two_Outside",
"TrueEta_pt_vs_rap_est_two_Outside; rapidity y; p_{t} in GeV/c ", 10,
2931 new TH2D(
"WrongEta_pt_vs_rap_est_two_Outside",
"WrongEta_pt_vs_rap_est_two_Outside; rapidity y; p_{t} in GeV/c ",
2932 10, 0., 4., 5, 0., 2.);
2935 new TH1D(
"TrueEta_InvMass_after_cuts_two_Outside",
2936 "TrueEta_InvMass_after_cuts_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2939 new TH1D(
"WrongEta_InvMass_after_cuts_two_Outside",
2940 "WrongEta_InvMass_after_cuts_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2943 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_two_Outside",
2944 "TrueEta_AngleBetweenGammas_after_cuts_two_Outside (between "
2945 "#gamma#gamma from #eta); angle [deg]",
2949 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_two_Outside",
2950 "WrongEta_AngleBetweenGammas_after_cuts_two_Outside (between "
2951 "#gamma#gamma); angle [deg]",
2955 "EMT_Eta_InvMass_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2963 "AngleBetweenGammas_trueEta_before_cuts_all_Both (between "
2964 "#gamma#gamma from #eta); angle [deg]",
2968 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_all_Both",
2969 "AngleBetweenGammas_wrongEta_before_cuts_all_Both (between "
2970 "#gamma#gamma); angle [deg]",
2974 new TH1D(
"InvMass_GammaCombinations_all_Both",
"InvMass_GammaCombinations_all_Both; invariant mass in GeV/c^{2};#",
2977 Particles_PDG_all_Both =
new TH1D(
"Particles_PDG_all_Both",
"Particles_PDG_all_Both ; PDG", 2300, 0, 2300);
2980 "TrueEta_pt_vs_rap_all_Both",
"TrueEta_pt_vs_rap_all_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60, -1., 5.);
2983 new TH2D(
"WrongEta_pt_vs_rap_all_Both",
"WrongEta_pt_vs_rap_all_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60,
2987 new TH2D(
"TrueEta_pt_vs_rap_est_all_Both",
"TrueEta_pt_vs_rap_est_all_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
2991 new TH2D(
"WrongEta_pt_vs_rap_est_all_Both",
"WrongEta_pt_vs_rap_est_all_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
2995 new TH1D(
"TrueEta_InvMass_after_cuts_all_Both",
2996 "TrueEta_InvMass_after_cuts_all_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2999 new TH1D(
"WrongEta_InvMass_after_cuts_all_Both",
3000 "WrongEta_InvMass_after_cuts_all_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3003 "TrueEta_AngleBetweenGammas_after_cuts_all_Both (between "
3004 "#gamma#gamma from #eta); angle [deg]",
3008 "WrongEta_AngleBetweenGammas_after_cuts_all_Both (between "
3009 "#gamma#gamma); angle [deg]",
3013 new TH1D(
"EMT_Eta_InvMass_all_Both",
"EMT_Eta_InvMass_all_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3018 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3019 "invariant mass in GeV/c^{2};#",
3023 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3024 "invariant mass in GeV/c^{2};#",
3028 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3029 "invariant mass in GeV/c^{2};#",
3033 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3034 "invariant mass in GeV/c^{2};#",
3038 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3039 "invariant mass in GeV/c^{2};#",
3043 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3044 "invariant mass in GeV/c^{2};#",
3048 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3049 "invariant mass in GeV/c^{2};#",
3053 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3054 "invariant mass in GeV/c^{2};#",
3058 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3059 "invariant mass in GeV/c^{2};#",
3063 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3064 "invariant mass in GeV/c^{2};#",
3068 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3069 "invariant mass in GeV/c^{2};#",
3073 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3074 "invariant mass in GeV/c^{2};#",
3078 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3079 "invariant mass in GeV/c^{2};#",
3083 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3084 "invariant mass in GeV/c^{2};#",
3088 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3089 "invariant mass in GeV/c^{2};#",
3093 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3094 "invariant mass in GeV/c^{2};#",
3098 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3099 "invariant mass in GeV/c^{2};#",
3103 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3104 "invariant mass in GeV/c^{2};#",
3108 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3109 "invariant mass in GeV/c^{2};#",
3113 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3114 "invariant mass in GeV/c^{2};#",
3118 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3119 "invariant mass in GeV/c^{2};#",
3123 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3124 "invariant mass in GeV/c^{2};#",
3128 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3129 "invariant mass in GeV/c^{2};#",
3133 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3134 "invariant mass in GeV/c^{2};#",
3138 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3139 "invariant mass in GeV/c^{2};#",
3143 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3144 "invariant mass in GeV/c^{2};#",
3148 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3149 "invariant mass in GeV/c^{2};#",
3153 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3154 "invariant mass in GeV/c^{2};#",
3158 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3159 "invariant mass in GeV/c^{2};#",
3163 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3164 "invariant mass in GeV/c^{2};#",
3168 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3169 "invariant mass in GeV/c^{2};#",
3173 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3174 "invariant mass in GeV/c^{2};#",
3178 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3179 "invariant mass in GeV/c^{2};#",
3183 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3184 "invariant mass in GeV/c^{2};#",
3191 new TH1D(
"AngleBetweenGammas_trueEta_before_cuts_onetwo_Both",
3192 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Both (between "
3193 "#gamma#gamma from #eta); angle [deg]",
3197 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both",
3198 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both (between "
3199 "#gamma#gamma); angle [deg]",
3203 new TH1D(
"InvMass_GammaCombinations_onetwo_Both",
3204 "InvMass_GammaCombinations_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3209 new TH2D(
"TrueEta_pt_vs_rap_onetwo_Both",
"TrueEta_pt_vs_rap_onetwo_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
3213 new TH2D(
"WrongEta_pt_vs_rap_onetwo_Both",
"WrongEta_pt_vs_rap_onetwo_Both; rapidity y; p_{t} in GeV/c ", 90, -2.,
3217 new TH2D(
"TrueEta_pt_vs_rap_est_onetwo_Both",
"TrueEta_pt_vs_rap_est_onetwo_Both; rapidity y; p_{t} in GeV/c ", 10,
3221 new TH2D(
"WrongEta_pt_vs_rap_est_onetwo_Both",
"WrongEta_pt_vs_rap_est_onetwo_Both; rapidity y; p_{t} in GeV/c ",
3222 10, 0., 4., 5, 0., 2.);
3225 new TH1D(
"TrueEta_InvMass_after_cuts_onetwo_Both",
3226 "TrueEta_InvMass_after_cuts_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3229 new TH1D(
"WrongEta_InvMass_after_cuts_onetwo_Both",
3230 "WrongEta_InvMass_after_cuts_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3233 new TH1D(
"TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both",
3234 "TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both (between "
3235 "#gamma#gamma from #eta); angle [deg]",
3239 new TH1D(
"WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both",
3240 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both (between "
3241 "#gamma#gamma); angle [deg]",
3245 "EMT_Eta_InvMass_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3250 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3251 "invariant mass in GeV/c^{2};#",
3255 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3256 "invariant mass in GeV/c^{2};#",
3260 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3261 "invariant mass in GeV/c^{2};#",
3265 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3266 "invariant mass in GeV/c^{2};#",
3270 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3271 "invariant mass in GeV/c^{2};#",
3275 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3276 "invariant mass in GeV/c^{2};#",
3280 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3281 "invariant mass in GeV/c^{2};#",
3285 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3286 "invariant mass in GeV/c^{2};#",
3290 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3291 "invariant mass in GeV/c^{2};#",
3295 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3296 "invariant mass in GeV/c^{2};#",
3300 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3301 "invariant mass in GeV/c^{2};#",
3305 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3306 "invariant mass in GeV/c^{2};#",
3310 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3311 "invariant mass in GeV/c^{2};#",
3315 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3316 "invariant mass in GeV/c^{2};#",
3320 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3321 "invariant mass in GeV/c^{2};#",
3325 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3326 "invariant mass in GeV/c^{2};#",
3330 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3331 "invariant mass in GeV/c^{2};#",
3335 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3336 "invariant mass in GeV/c^{2};#",
3340 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3341 "invariant mass in GeV/c^{2};#",
3345 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3346 "invariant mass in GeV/c^{2};#",
3350 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3351 "invariant mass in GeV/c^{2};#",
3355 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3356 "invariant mass in GeV/c^{2};#",
3360 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3361 "invariant mass in GeV/c^{2};#",
3365 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3366 "invariant mass in GeV/c^{2};#",
3370 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3371 "invariant mass in GeV/c^{2};#",
3375 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3376 "invariant mass in GeV/c^{2};#",
3380 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3381 "invariant mass in GeV/c^{2};#",
3385 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3386 "invariant mass in GeV/c^{2};#",
3390 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3391 "invariant mass in GeV/c^{2};#",
3395 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3396 "invariant mass in GeV/c^{2};#",
3400 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3401 "invariant mass in GeV/c^{2};#",
3405 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3406 "invariant mass in GeV/c^{2};#",
3410 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3411 "invariant mass in GeV/c^{2};#",
3415 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3416 "invariant mass in GeV/c^{2};#",
3423 "AngleBetweenGammas_trueEta_before_cuts_two_Both (between "
3424 "#gamma#gamma from #eta); angle [deg]",
3428 new TH1D(
"AngleBetweenGammas_wrongEta_before_cuts_two_Both",
3429 "AngleBetweenGammas_wrongEta_before_cuts_two_Both (between "
3430 "#gamma#gamma); angle [deg]",
3434 new TH1D(
"InvMass_GammaCombinations_two_Both",
"InvMass_GammaCombinations_two_Both; invariant mass in GeV/c^{2};#",
3437 Particles_PDG_two_Both =
new TH1D(
"Particles_PDG_two_Both",
"Particles_PDG_two_Both ; PDG", 2300, 0, 2300);
3440 "TrueEta_pt_vs_rap_two_Both",
"TrueEta_pt_vs_rap_two_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60, -1., 5.);
3443 new TH2D(
"WrongEta_pt_vs_rap_two_Both",
"WrongEta_pt_vs_rap_two_Both; rapidity y; p_{t} in GeV/c ", 90, -2., 7., 60,
3447 new TH2D(
"TrueEta_pt_vs_rap_est_two_Both",
"TrueEta_pt_vs_rap_est_two_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
3451 new TH2D(
"WrongEta_pt_vs_rap_est_two_Both",
"WrongEta_pt_vs_rap_est_two_Both; rapidity y; p_{t} in GeV/c ", 10, 0.,
3455 new TH1D(
"TrueEta_InvMass_after_cuts_two_Both",
3456 "TrueEta_InvMass_after_cuts_two_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3459 new TH1D(
"WrongEta_InvMass_after_cuts_two_Both",
3460 "WrongEta_InvMass_after_cuts_two_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3463 "TrueEta_AngleBetweenGammas_after_cuts_two_Both (between "
3464 "#gamma#gamma from #eta); angle [deg]",
3468 "WrongEta_AngleBetweenGammas_after_cuts_two_Both (between "
3469 "#gamma#gamma); angle [deg]",
3473 new TH1D(
"EMT_Eta_InvMass_two_Both",
"EMT_Eta_InvMass_two_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3478 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3479 "invariant mass in GeV/c^{2};#",
3483 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3484 "invariant mass in GeV/c^{2};#",
3488 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3489 "invariant mass in GeV/c^{2};#",
3493 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3494 "invariant mass in GeV/c^{2};#",
3498 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3499 "invariant mass in GeV/c^{2};#",
3503 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3504 "invariant mass in GeV/c^{2};#",
3508 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3509 "invariant mass in GeV/c^{2};#",
3513 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3514 "invariant mass in GeV/c^{2};#",
3518 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3519 "invariant mass in GeV/c^{2};#",
3523 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3524 "invariant mass in GeV/c^{2};#",
3528 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3529 "invariant mass in GeV/c^{2};#",
3533 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3534 "invariant mass in GeV/c^{2};#",
3538 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3539 "invariant mass in GeV/c^{2};#",
3543 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3544 "invariant mass in GeV/c^{2};#",
3548 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3549 "invariant mass in GeV/c^{2};#",
3553 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3554 "invariant mass in GeV/c^{2};#",
3558 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3559 "invariant mass in GeV/c^{2};#",
3563 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3564 "invariant mass in GeV/c^{2};#",
3568 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3569 "invariant mass in GeV/c^{2};#",
3573 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3574 "invariant mass in GeV/c^{2};#",
3578 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3579 "invariant mass in GeV/c^{2};#",
3583 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3584 "invariant mass in GeV/c^{2};#",
3588 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3589 "invariant mass in GeV/c^{2};#",
3593 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3594 "invariant mass in GeV/c^{2};#",
3598 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3599 "invariant mass in GeV/c^{2};#",
3603 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3604 "invariant mass in GeV/c^{2};#",
3608 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3609 "invariant mass in GeV/c^{2};#",
3613 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3614 "invariant mass in GeV/c^{2};#",
3618 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3619 "invariant mass in GeV/c^{2};#",
3623 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3624 "invariant mass in GeV/c^{2};#",
3628 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3629 "invariant mass in GeV/c^{2};#",
3633 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3634 "invariant mass in GeV/c^{2};#",
3638 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3639 "invariant mass in GeV/c^{2};#",
3643 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3644 "invariant mass in GeV/c^{2};#",
Here the ring is fitted with Taubin algorithm from A. Ayriyan, G. Ososkov, N. Chernov.
Data class for a reconstructed hit in the STS.
Data class for STS tracks.
int32_t GetStsTrackIndex() const
int32_t GetRichRingIndex() const
static void SetKFParticleFromStsTrack(CbmStsTrack *track, KFParticle *particle, Int_t pdg=211, Bool_t firstPoint=kTRUE)
TH1D * multi_EMT_Eta_onetwo_Both_5
std::vector< std::vector< int > > Gammas_MCIndex_all_Outside
TH1D * multi_InvMass_Eta_onetwo_Both_10
TH1D * multi_EMT_Eta_two_Both_5
TH2D * WrongEta_pt_vs_rap_est_all_Outside
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Target
vector< TH1 * > fHistoList_Eta_all_Target
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Target
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Both
TH1D * PlaneAngles_first_fromEta_zeroInRich_Outside
TH1D * EMT_Eta_InvMass_two_Both
std::vector< std::vector< TVector3 > > Gammas_onetwo_Both
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Target
TH1D * multi_EMT_Eta_onetwo_Both_11
TH1D * TrueEta_InvMass_after_cuts_two_Both
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Target
vector< CbmRichRing * > VRichRing_minus_Outside
TH1D * multi_EMT_Eta_onetwo_Both_15
TH1D * multi_InvMass_Eta_two_Both_4
TH1D * Particle_pull_X_Outside
TH1D * multi_InvMass_Eta_all_Both_7
TH1D * multi_EMT_Eta_all_Both_9
std::vector< int > EMT_eta_gg_Event_Target
TH1D * multi_InvMass_Eta_two_Both_9
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Target
TH1D * multi_EMT_Eta_two_Both_2
TH1D * multi_InvMass_Eta_two_Both_6
TH1D * Particles_PDG_two_Target
vector< TH1 * > fHistoList_rap_pt_Eta_all_Both
std::vector< int > frefId
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Both
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Outside
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Target
TH1D * InvMass_GammaCombinations_all_Target
TH1D * Particles_PDG_all_Target
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Outside
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Target
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside
TH1D * Particles_PDG_all_Both
TH1D * multi_EMT_Eta_all_Both_2
std::vector< int > EMT_eta_gg_Event_Both
TH1D * Particle_pull_px_Target
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Target
TH1D * multi_EMT_Eta_two_Both_9
TH1D * PlaneAngles_first_fromEta_oneInRich_Outside
TH2D * TrueEta_pt_vs_rap_est_all_Outside
TH2D * TrueEta_pt_vs_rap_est_two_Both
TH1D * multi_EMT_Eta_all_Both_15
TH1D * Particles_PDG_two_Outside
vector< CbmRichRing * > VRichRing_plus_Target
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both
TH1D * multi_EMT_Eta_onetwo_Both_9
TH1D * WrongEta_InvMass_after_cuts_all_Both
TH1D * Particle_pull_pz_Target
vector< CbmStsTrack * > VStsTrack_minus_Target
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Both
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target
TH1D * InvMass_GammaCombinations_onetwo_Outside
vector< TH1 * > fHistoList_Eta_cuts_Outside
TH1D * multi_InvMass_Eta_onetwo_Both_15
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Outside
TH2D * WrongEta_pt_vs_rap_two_Outside
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Outside
TClonesArray * fStsTracks
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Both
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Target
std::vector< std::vector< int > > Gammas_stsIndex_all_Both
TH1D * multi_EMT_Eta_two_Both_17
TH2D * TrueEta_pt_vs_rap_est_two_Outside
vector< TH1 * > fHistoList_rap_pt_Eta_two_Target
TH1D * PlaneAngles_last_fromEta_afterCuts_oneInRich_Both
std::vector< std::vector< int > > Gammas_stsIndex_two_Target
TH1D * PlaneAngles_last_fromEta_zeroInRich_Target
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Target
TH1D * InvMass_GammaCombinations_all_Both
TH2D * TrueEta_pt_vs_rap_est_onetwo_Target
std::vector< int > VStsIndex_minus_Target
TH1D * multi_InvMass_Eta_all_Both_5
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Outside
TH1D * multi_EMT_Eta_onetwo_Both_6
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside
TH1D * PlaneAngles_last_fromEta_twoInRich_Both
TH1D * Particles_PDG_onetwo_Both
TH1D * WrongEta_InvMass_after_cuts_two_Both
TH1D * multi_EMT_Eta_onetwo_Both_10
std::vector< int > VRings_minus_Outside
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Both
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Both
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target
TH1D * multi_EMT_Eta_all_Both_7
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside
TH1D * PlaneAngles_first_fromEta_oneInRich_Both
vector< TH1 * > fHistoList_Eta_two_Outside
std::vector< int > VStsIndex_plus_Outside
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Outside
TH1D * InvMass_GammaCombinations_two_Outside
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Target
vector< TH1 * > fHistoList_Eta_cuts_Target
TH2D * WrongEta_pt_vs_rap_onetwo_Target
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Outside
TH2D * TrueEta_pt_vs_rap_est_two_Target
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside
TH1D * WrongEta_InvMass_after_cuts_onetwo_Target
TH1D * multi_EMT_Eta_all_Both_17
TH1D * multi_InvMass_Eta_all_Both_12
TH1D * WrongEta_InvMass_after_cuts_all_Outside
std::vector< std::vector< TVector3 > > Gammas_two_Outside
TH1D * multi_EMT_Eta_two_Both_8
std::vector< std::vector< int > > Gammas_stsIndex_two_Outside
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target
TH1D * multi_EMT_Eta_onetwo_Both_2
std::vector< std::vector< TVector3 > > Gammas_onetwo_Outside
TH2D * TrueEta_pt_vs_rap_all_Target
TH1D * PlaneAngles_first_fromEta_zeroInRich_Both
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Both
vector< TH1 * > fHistoList_Eta_onetwo_Target
TH1D * multi_InvMass_Eta_two_Both_12
TH1D * EMT_Eta_InvMass_onetwo_Both
TH1D * multi_InvMass_Eta_two_Both_7
TClonesArray * fArrayStsHit
TH1D * TrueEta_InvMass_after_cuts_onetwo_Outside
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both
void Exec(int fEventNumEta, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
TH1D * multi_InvMass_Eta_two_Both_1
vector< TH1 * > fHistoList_Eta_two_Target
std::vector< int > VRings_plus_Outside
TH1D * multi_EMT_Eta_all_Both_16
TH2D * TrueEta_pt_vs_rap_onetwo_Target
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside
TH1D * PlaneAngles_last_fromEta_zeroInRich_Both
void FindGammasOutside(int EventNumEta, double AngleCut, double InvMassCut, int RealPID, vector< CbmMCTrack * > MCtracks_minus_Outside, vector< CbmMCTrack * > MCtracks_plus_Outside, vector< CbmStsTrack * > StsTrack_minus_Outside, vector< CbmStsTrack * > StsTrack_plus_Outside, std::vector< int > Rings_minus_Outside, std::vector< int > Rings_plus_Outside, std::vector< int > stsIndex_minus_Outside, std::vector< int > stsIndex_plus_Outside, vector< CbmRichRing * > richRing_minus_Outside, vector< CbmRichRing * > richRing_plus_Outside, vector< Int_t > MCIndex_minus_Outside, vector< Int_t > MCIndex_plus_Outside)
void FindEta(TString mod, TString position, vector< vector< TVector3 > > Gammas, vector< vector< int > > StsIndex, vector< vector< int > > MCIndex, vector< vector< CbmMCTrack * > > GammasMC, vector< TH1 * > gg, vector< TH1 * > rap_pt_separation)
TH1D * multi_EMT_Eta_two_Both_12
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Both
TH1D * TrueEta_InvMass_after_cuts_all_Outside
TH1D * InvMass_GammaCombinations_two_Target
TH1D * multi_InvMass_Eta_two_Both_5
TH1D * PlaneAngles_last_fromEta_zeroInRich_Outside
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Outside
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both
TH1D * multi_EMT_Eta_two_Both_4
double CalculatePlaneAngle_last(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
TH1D * multi_EMT_Eta_all_Both_3
TH1D * multi_EMT_Eta_onetwo_Both_7
int FindInRich(int richInd, int stsMcTrackId)
TH1D * multi_EMT_Eta_all_Both_14
TH1D * multi_InvMass_Eta_onetwo_Both_12
TH1D * PlaneAngles_first_fromEta_twoInRich_Target
TH1D * multi_EMT_Eta_onetwo_Both_16
TH2D * TrueEta_pt_vs_rap_two_Both
TH2D * WrongEta_pt_vs_rap_all_Both
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Target
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside
TH1D * multi_InvMass_Eta_onetwo_Both_7
TH1D * multi_InvMass_Eta_all_Both_13
TH1D * multi_InvMass_Eta_all_Both_4
vector< CbmMCTrack * > VMCtracks_plus_Outside
TH1D * multi_EMT_Eta_two_Both_14
TH1D * multi_InvMass_Eta_two_Both_8
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside
TH1D * multi_InvMass_Eta_onetwo_Both_4
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Both
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target
TH1D * multi_InvMass_Eta_two_Both_15
CbmRichRingFitterEllipseTau * fTauFit
TH1D * PlaneAngles_first_fromEta_zeroInRich_Target
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Target
TH2D * WrongEta_pt_vs_rap_all_Outside
TH1D * multi_InvMass_Eta_all_Both_16
TH1D * multi_InvMass_Eta_onetwo_Both_14
vector< CbmMCTrack * > VMCtracks_minus_Outside
std::vector< int > EMT_eta_gg_NofRings_Target
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both
vector< CbmStsTrack * > VStsTrack_minus_Outside
std::vector< std::vector< TVector3 > > Gammas_onetwo_Target
TH2D * WrongEta_pt_vs_rap_est_onetwo_Both
vector< CbmRichRing * > VRichRing_plus_Outside
TH1D * multi_EMT_Eta_all_Both_4
vector< Int_t > VMCIndex_plus_Outside
vector< CbmRichRing * > VRichRing_minus_Target
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Outside
TH2D * TrueEta_pt_vs_rap_all_Outside
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Target
TH1D * PlaneAngles_first_fromEta_oneInRich_Target
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Both
std::vector< std::vector< int > > Gammas_MCIndex_two_Outside
TH1D * multi_InvMass_Eta_two_Both_10
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Target
TH1D * multi_InvMass_Eta_two_Both_13
TH1D * PlaneAngles_last_fromEta_twoInRich_Target
TH1D * multi_EMT_Eta_two_Both_1
TClonesArray * fRichRings
TH1D * multi_InvMass_Eta_all_Both_8
TH1D * Particle_pull_py_Outside
std::vector< int > VRings_plus_Target
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Both
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Both
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Target
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Outside
vector< TVector3 > VMomenta_plus_Target
TH1D * multi_InvMass_Eta_all_Both_10
std::vector< std::vector< TVector3 > > Gammas_all_Both
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Both
TH2D * TrueEta_pt_vs_rap_est_onetwo_Both
TH1D * multi_EMT_Eta_onetwo_Both_17
vector< TH1 * > fHistoList_Eta_two_Both
TH1D * PlaneAngles_last_fromEta_oneInRich_Target
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Both
TClonesArray * fArrayMvdHit
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both
std::vector< std::vector< TVector3 > > Gammas_all_Target
TH1D * multi_InvMass_Eta_all_Both_9
TH1D * Particle_pull_pz_Outside
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Target
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Both
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Both
std::vector< std::vector< int > > Gammas_MCIndex_all_Target
TH2D * WrongEta_pt_vs_rap_est_all_Both
TH1D * multi_InvMass_Eta_all_Both_14
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Outside
TH1D * WrongEta_InvMass_after_cuts_all_Target
vector< TH1 * > fHistoList_rap_pt_Eta_all_Target
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Outside
TH1D * Particle_pull_Y_Outside
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target
TH1D * EMT_Eta_InvMass_onetwo_Target
TH1D * multi_InvMass_Eta_all_Both_2
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Both
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Target
TH1D * multi_EMT_Eta_all_Both_10
TClonesArray * fRichProjections
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Outside
TH1D * TrueEta_InvMass_after_cuts_two_Outside
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target
TH1D * multi_InvMass_Eta_all_Both_11
TH1D * EMT_Eta_InvMass_two_Target
TH1D * multi_InvMass_Eta_two_Both_14
TH2D * TrueEta_pt_vs_rap_onetwo_Outside
void FindGammasTarget(int EventNumEta, double AngleCut, double InvMassCut, int RealPID, vector< CbmMCTrack * > MCtracks_minus, vector< CbmMCTrack * > MCtracks_plus, vector< CbmStsTrack * > StsTrack_minus, vector< CbmStsTrack * > StsTrack_plus, vector< TVector3 > Momenta_minus, vector< TVector3 > Momenta_plus, std::vector< int > Rings_minus, std::vector< int > Rings_plus, std::vector< int > stsIndex_minus, std::vector< int > stsIndex_plus, vector< CbmRichRing * > richRing_minus, vector< CbmRichRing * > richRing_plus, vector< Int_t > MCIndex_minus, vector< Int_t > MCIndex_plus)
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Outside
TH1D * multi_EMT_Eta_all_Both_6
vector< CbmMCTrack * > VMCtracks_plus_Target
TH1D * PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target
TH1D * multi_EMT_Eta_onetwo_Both_12
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Target
TClonesArray * fStsTrackMatches
TH1D * Particles_PDG_onetwo_Target
double CalculatePlaneAngle_first(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
std::vector< std::vector< int > > Gammas_MCIndex_all_Both
TH1D * multi_EMT_Eta_two_Both_16
TH2D * TrueEta_pt_vs_rap_est_all_Both
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside
vector< Int_t > VMCIndex_plus_Target
TH2D * WrongEta_pt_vs_rap_est_onetwo_Outside
TH1D * multi_InvMass_Eta_two_Both_11
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target
TH1D * InvMass_GammaCombinations_onetwo_Both
TH1D * PlaneAngles_last_fromEta_oneInRich_Outside
TH1D * Particle_pull_Z_Outside
vector< TH1 * > fHistoList_Eta_all_Outside
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Target
TH1D * multi_EMT_Eta_two_Both_6
std::vector< int > EMT_eta_gg_NofRings_Both
int CheckIfElectron(CbmRichRing *ring, double momentum)
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside
TH1D * multi_EMT_Eta_onetwo_Both_13
TH1D * InvMass_GammaCombinations_onetwo_Target
TH1D * multi_InvMass_Eta_onetwo_Both_11
TH1D * multi_EMT_Eta_onetwo_Both_4
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Outside
TH1D * multi_InvMass_Eta_onetwo_Both_8
TH2D * WrongEta_pt_vs_rap_two_Target
TH1D * multi_InvMass_Eta_onetwo_Both_1
vector< CbmMCTrack * > VMCtracks_minus_Target
TH1D * EMT_Eta_InvMass_all_Outside
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target
TH2D * WrongEta_pt_vs_rap_est_two_Target
TH1D * PlaneAngles_last_fromEta_oneInRich_Both
TH2D * WrongEta_pt_vs_rap_est_two_Both
TH1D * multi_EMT_Eta_all_Both_12
vector< Int_t > VMCIndex_minus_Target
vector< TH1 * > fHistoList_Eta_onetwo_Outside
std::vector< std::vector< TVector3 > > Gammas_two_Both
std::vector< int > EMT_eta_gg_Event_Outside
TH1D * multi_EMT_Eta_all_Both_8
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside
TH1D * Particle_pull_px_Outside
TH1D * EMT_Eta_InvMass_all_Target
TH1D * WrongEta_InvMass_after_cuts_two_Outside
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside
TH1D * multi_InvMass_Eta_all_Both_17
TH1D * multi_EMT_Eta_onetwo_Both_14
TH1D * TrueEta_InvMass_after_cuts_onetwo_Both
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Target
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Target
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Both
vector< Int_t > VMCIndex_minus_Outside
TH1D * Particles_PDG_onetwo_Outside
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Target
std::vector< std::vector< int > > Gammas_MCIndex_two_Target
TH1D * multi_InvMass_Eta_onetwo_Both_3
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside
TH2D * WrongEta_pt_vs_rap_onetwo_Outside
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Both
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Target
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target
TH2D * TrueEta_pt_vs_rap_est_all_Target
TH1D * multi_InvMass_Eta_all_Both_1
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Outside
vector< TH1 * > fHistoList_Eta_onetwo_Both
TH1D * PlaneAngles_first_fromEta_afterCuts_oneInRich_Both
TH1D * multi_EMT_Eta_all_Both_13
TH1D * Particles_PDG_two_Both
TH1D * PlaneAngles_first_fromEta_twoInRich_Outside
TH1D * PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Outside
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Target
vector< TH1 * > fHistoList_rap_pt_Eta_all_Outside
vector< TVector3 > VMomenta_minus_Target
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Both
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside
vector< CbmStsTrack * > VStsTrack_plus_Target
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Outside
std::vector< CbmMCTrack * > fMCtracks
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Outside
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Both
TH1D * Particle_pull_py_Target
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Both
TH1D * WrongEta_InvMass_after_cuts_two_Target
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Both
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Target
TH1D * multi_InvMass_Eta_all_Both_3
TH2D * TrueEta_pt_vs_rap_onetwo_Both
vector< TH1 * > fHistoList_rap_pt_Eta_two_Both
vector< CbmStsTrack * > VStsTrack_plus_Outside
std::vector< int > EMT_eta_gg_NofRings_Outside
TH1D * InvMass_GammaCombinations_two_Both
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Both
TH1D * multi_InvMass_Eta_all_Both_15
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Both
TH1D * InvMass_GammaCombinations_all_Outside
TH1D * multi_EMT_Eta_two_Both_7
vector< TH1 * > fHistoList_Eta_cuts_Both
TH1D * multi_InvMass_Eta_onetwo_Both_5
std::vector< std::vector< int > > Gammas_stsIndex_all_Outside
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both
TH1D * multi_InvMass_Eta_all_Both_6
void SaveOutsideTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
TH1D * TrueEta_InvMass_after_cuts_onetwo_Target
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Target
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Outside
TH1D * multi_InvMass_Eta_onetwo_Both_17
TH1D * multi_EMT_Eta_onetwo_Both_1
TH1D * WrongEta_InvMass_after_cuts_onetwo_Outside
TH1D * TrueEta_InvMass_after_cuts_all_Target
TH1D * WrongEta_InvMass_after_cuts_onetwo_Both
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Both
TH1D * multi_EMT_Eta_all_Both_5
vector< TH1 * > fHistoList_Eta_all_Both
TH1D * multi_EMT_Eta_onetwo_Both_3
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Outside
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Target
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both
TH1D * multi_InvMass_Eta_onetwo_Both_16
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Outside
std::vector< std::vector< int > > Gammas_stsIndex_two_Both
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Both
std::vector< std::vector< int > > Gammas_MCIndex_two_Both
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside
TH1D * multi_EMT_Eta_two_Both_3
TH1D * Particles_PDG_all_Outside
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both
TH2D * WrongEta_pt_vs_rap_est_two_Outside
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Both
std::vector< std::vector< TVector3 > > Gammas_all_Outside
TH2D * TrueEta_pt_vs_rap_two_Outside
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Target
TH2D * WrongEta_pt_vs_rap_onetwo_Both
vector< TVector3 > frefmomenta
TH1D * multi_InvMass_Eta_onetwo_Both_9
vector< TH1 * > fHistoList_rap_pt_Eta_two_Outside
std::vector< int > VRings_minus_Target
std::vector< std::vector< TVector3 > > Gammas_two_Target
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Outside
TH1D * multi_EMT_Eta_two_Both_10
TH1D * EMT_Eta_InvMass_onetwo_Outside
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Both
TH1D * multi_InvMass_Eta_onetwo_Both_2
TH2D * WrongEta_pt_vs_rap_two_Both
TH1D * multi_EMT_Eta_two_Both_15
void Mixing_Both(vector< TH1 * > rap_pt_separation_all, vector< TH1 * > rap_pt_separation_onetwo, vector< TH1 * > rap_pt_separation_two)
std::vector< std::vector< int > > Gammas_stsIndex_all_Target
TH1D * PlaneAngles_first_fromEta_twoInRich_Both
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Outside
TH1D * multi_EMT_Eta_all_Both_1
void SaveTargetTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, TVector3 refmom, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
std::vector< int > VStsIndex_minus_Outside
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Outside
TH1D * EMT_Eta_InvMass_two_Outside
TH2D * TrueEta_pt_vs_rap_est_onetwo_Outside
TH2D * TrueEta_pt_vs_rap_two_Target
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Target
TH2D * WrongEta_pt_vs_rap_all_Target
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Outside
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Both
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Outside
TH1D * multi_InvMass_Eta_onetwo_Both_13
std::vector< int > VStsIndex_plus_Target
vector< TH1 * > fHistoList_Eta
TH2D * TrueEta_pt_vs_rap_all_Both
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Both
TClonesArray * fGlobalTracks
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Outside
TH1D * multi_InvMass_Eta_two_Both_16
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Target
TH1D * TrueEta_InvMass_after_cuts_all_Both
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Target
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Both
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Outside
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Target
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both
TH2D * WrongEta_pt_vs_rap_est_all_Target
TClonesArray * fRichRingMatches
TH1D * multi_EMT_Eta_two_Both_13
TH1D * multi_InvMass_Eta_two_Both_17
TH1D * multi_InvMass_Eta_two_Both_2
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Both
TH1D * TrueEta_InvMass_after_cuts_two_Target
TH1D * EMT_Eta_InvMass_all_Both
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Outside
TH1D * multi_EMT_Eta_two_Both_11
TH2D * WrongEta_pt_vs_rap_est_onetwo_Target
TH1D * PlaneAngles_last_fromEta_twoInRich_Outside
TH1D * multi_EMT_Eta_onetwo_Both_8
TH1D * multi_InvMass_Eta_two_Both_3
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Target
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Both
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both
TH1D * multi_EMT_Eta_all_Both_11
TH1D * multi_InvMass_Eta_onetwo_Both_6
static Double_t CalculateOpeningAngleBetweenGammas_Reco(TVector3 electron1, TVector3 electron2, TVector3 electron3, TVector3 electron4)
static TVector3 FitToVertexAndGetChi(CbmStsTrack *stsTrack, double x, double y, double z, double &chi)
static Double_t CalculateOpeningAngle_Reco(TVector3 electron1, TVector3 electron2)
static LmvmKinePar CalculateKinematicParams_4particles(const TVector3 part1, const TVector3 part2, const TVector3 part3, const TVector3 part4)
static double Invmass_2particles_RECO(const TVector3 part1, const TVector3 part2)
static LmvmKinePar CalculateKinematicParamsReco(const TVector3 electron1, const TVector3 electron2)
int32_t GetMotherId() const
int32_t GetPdgCode() const
int32_t GetNofLinks() const
const CbmLink & GetMatchedLink() const
Here the ring is fitted with Taubin algorithm from A. Ayriyan, G. Ososkov, N. Chernov.
void AddHit(CbmRichHitLight hit)
Add new hit to the ring.
uint32_t GetHit(int32_t i) const
int32_t GetNofHits() const
data class for a reconstructed 3-d hit in the STS
int32_t GetNofMvdHits() const
int32_t GetMvdHitIndex(int32_t iHit) const
int32_t GetStsHitIndex(int32_t iHit) const
int32_t GetNofStsHits() const
const FairTrackParam * GetParamFirst() const