CbmRoot
Loading...
Searching...
No Matches
CbmKresEta.cxx
Go to the documentation of this file.
1/* Copyright (C) 2018-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
2 SPDX-License-Identifier: GPL-3.0-only
3 Authors: Ievgenii Kres, Florian Uhlig [committer] */
4
20
21#include "CbmKresEta.h"
22
23#include "CbmGlobalTrack.h"
25#include "CbmKresFunctions.h"
26#include "CbmMCTrack.h"
27#include "CbmMvdHit.h"
28#include "CbmRichHit.h"
29#include "CbmRichRing.h"
31#include "CbmRichRingLight.h"
32#include "CbmRichUtil.h"
33#include "CbmStsHit.h"
34#include "CbmStsTrack.h"
35#include "CbmTrackMatchNew.h"
36#include "FairRootManager.h"
37#include "KFParticle/KFParticle.h"
38#include "TDirectory.h"
39
40#include <iostream>
41
42
43using namespace std;
44
46 : fMcTracks(nullptr)
47 , fGlobalTracks(nullptr)
48 , fStsTracks(nullptr)
49 , fStsTrackMatches(nullptr)
50 , fRichProjections(nullptr)
51 , fRichRings(nullptr)
52 , fRichRingMatches(nullptr)
53 , fRichHits(nullptr)
54 , fArrayMvdHit(nullptr)
55 , fArrayStsHit(nullptr)
56 , fPrimVertex(nullptr)
57 , fKFVertex()
58 , fTauFit(nullptr)
85 , frefmomenta()
86 , frefId()
87 , fMCId()
88 , fMCtracks()
135 , Particle_pull_px_Target(nullptr)
136 , Particle_pull_py_Target(nullptr)
137 , Particle_pull_pz_Target(nullptr)
138 , Particle_pull_px_Outside(nullptr)
139 , Particle_pull_py_Outside(nullptr)
140 , Particle_pull_pz_Outside(nullptr)
141 , Particle_pull_X_Outside(nullptr)
142 , Particle_pull_Y_Outside(nullptr)
143 , Particle_pull_Z_Outside(nullptr)
272 , Particles_PDG_all_Target(nullptr)
300 , Particles_PDG_two_Target(nullptr)
362 , Particles_PDG_all_Both(nullptr)
371 , EMT_Eta_InvMass_all_Both(nullptr)
390 , multi_EMT_Eta_all_Both_1(nullptr)
391 , multi_EMT_Eta_all_Both_2(nullptr)
392 , multi_EMT_Eta_all_Both_3(nullptr)
393 , multi_EMT_Eta_all_Both_4(nullptr)
394 , multi_EMT_Eta_all_Both_5(nullptr)
395 , multi_EMT_Eta_all_Both_6(nullptr)
396 , multi_EMT_Eta_all_Both_7(nullptr)
397 , multi_EMT_Eta_all_Both_8(nullptr)
398 , multi_EMT_Eta_all_Both_9(nullptr)
460 , Particles_PDG_two_Both(nullptr)
469 , EMT_Eta_InvMass_two_Both(nullptr)
488 , multi_EMT_Eta_two_Both_1(nullptr)
489 , multi_EMT_Eta_two_Both_2(nullptr)
490 , multi_EMT_Eta_two_Both_3(nullptr)
491 , multi_EMT_Eta_two_Both_4(nullptr)
492 , multi_EMT_Eta_two_Both_5(nullptr)
493 , multi_EMT_Eta_two_Both_6(nullptr)
494 , multi_EMT_Eta_two_Both_7(nullptr)
495 , multi_EMT_Eta_two_Both_8(nullptr)
496 , multi_EMT_Eta_two_Both_9(nullptr)
505{
506}
507
509
511{
512 FairRootManager* ioman = FairRootManager::Instance();
513 if (nullptr == ioman) { Fatal("CbmKresEta::Init", "RootManager not instantised!"); }
514
515 fMcTracks = (TClonesArray*) ioman->GetObject("MCTrack");
516 if (nullptr == fMcTracks) { Fatal("CbmKresEta::Init", "No MCTrack array!"); }
517
518 fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex."));
519 if (nullptr == fPrimVertex) { fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex")); }
520 if (nullptr == fPrimVertex) { LOG(fatal) << "CbmKresEta::Init No PrimaryVertex array!"; }
521
522 fGlobalTracks = (TClonesArray*) ioman->GetObject("GlobalTrack");
523 if (nullptr == fGlobalTracks) { Fatal("CbmKresEta::Init", "No GlobalTrack array!"); }
524
525 fStsTracks = (TClonesArray*) ioman->GetObject("StsTrack");
526 if (nullptr == fStsTracks) { Fatal("CbmKresEta::Init", "No StsTrack array!"); }
527
528 fStsTrackMatches = (TClonesArray*) ioman->GetObject("StsTrackMatch");
529 if (nullptr == fStsTrackMatches) { Fatal("CbmKresEta::Init", "No StsTrackMatch array!"); }
530
531 fRichProjections = (TClonesArray*) ioman->GetObject("RichProjection");
532 if (nullptr == fRichProjections) { Fatal("CbmKresEta::Init", "No RichProjection array!"); }
533
534 fRichRings = (TClonesArray*) ioman->GetObject("RichRing");
535 if (nullptr == fRichRings) { Fatal("CbmKresEta::Init", "No RichRing array!"); }
536
537 fRichRingMatches = (TClonesArray*) ioman->GetObject("RichRingMatch");
538 if (nullptr == fRichRingMatches) { Fatal("CbmKresEta::Init", "No RichRingMatch array!"); }
539
540 fRichHits = (TClonesArray*) ioman->GetObject("RichHit");
541 if (nullptr == fRichHits) { Fatal("CbmKresEta::Init", "No RichHit array!"); }
542
543 fArrayMvdHit = (TClonesArray*) ioman->GetObject("MvdHit");
544 if (nullptr == fArrayMvdHit) { Fatal("CbmKresEta::Init", "No MvdHit array!"); }
545
546 fArrayStsHit = (TClonesArray*) ioman->GetObject("StsHit");
547 if (nullptr == fArrayStsHit) { Fatal("CbmKresEta::Init", "No StsHit array!"); }
548
550
552}
553
554void CbmKresEta::Exec(int fEventNumEta, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
555{
556 // cout << "CbmKresEta, event No. " << fEventNumEta << endl;
557
558 if (fPrimVertex != nullptr) { fKFVertex = CbmKFVertex(*fPrimVertex); }
559 else {
560 Fatal("CbmKresConversionManual::Exec", "No PrimaryVertex array!");
561 }
562
563 Gammas_all_Target.clear();
564 Gammas_two_Target.clear();
565 Gammas_onetwo_Target.clear();
572 Gammas_MC_all_Target.clear();
573 Gammas_MC_two_Target.clear();
575 Gammas_all_Outside.clear();
576 Gammas_two_Outside.clear();
577 Gammas_onetwo_Outside.clear();
584 Gammas_MC_all_Outside.clear();
585 Gammas_MC_two_Outside.clear();
587 Gammas_all_Both.clear();
588 Gammas_two_Both.clear();
589 Gammas_onetwo_Both.clear();
596 Gammas_MC_all_Both.clear();
597 Gammas_MC_two_Both.clear();
598 Gammas_MC_onetwo_Both.clear();
599
602 VRings_minus_Outside.clear();
608 VRings_plus_Outside.clear();
611 VMCIndex_plus_Outside.clear();
612
614 VMomenta_minus_Target.clear();
616 VRings_minus_Target.clear();
619 VMCIndex_minus_Target.clear();
620 VMCtracks_plus_Target.clear();
621 VStsTrack_plus_Target.clear();
622 VMomenta_plus_Target.clear();
623 VRings_plus_Target.clear();
624 VStsIndex_plus_Target.clear();
625 VRichRing_plus_Target.clear();
626 VMCIndex_plus_Target.clear();
627
628 Int_t ngTracks = fGlobalTracks->GetEntriesFast();
629 for (Int_t i = 0; i < ngTracks; i++) {
630 CbmGlobalTrack* gTrack = (CbmGlobalTrack*) fGlobalTracks->At(i);
631 if (nullptr == gTrack) continue;
632 Int_t stsInd = gTrack->GetStsTrackIndex();
633 Int_t richInd = gTrack->GetRichRingIndex();
634
635 if (stsInd < 0) continue;
636 CbmStsTrack* stsTrack = (CbmStsTrack*) fStsTracks->At(stsInd);
637 if (stsTrack == nullptr) continue;
638 CbmTrackMatchNew* stsMatch = (CbmTrackMatchNew*) fStsTrackMatches->At(stsInd);
639 if (stsMatch == nullptr) continue;
640 if (stsMatch->GetNofLinks() <= 0) continue;
641 Int_t stsMcTrackId = stsMatch->GetMatchedLink().GetIndex();
642 if (stsMcTrackId < 0) continue;
643 CbmMCTrack* mcTrack = (CbmMCTrack*) fMcTracks->At(stsMcTrackId);
644 if (mcTrack == nullptr) continue;
645
646 FairTrackParam* proj = (FairTrackParam*) fRichProjections->At(i);
647 if (richInd < 0 && proj->GetX() > -115 && proj->GetX() < 115
648 && ((proj->GetY() < -120 && proj->GetY() > -200) || (proj->GetY() > 120 && proj->GetY() < 200)))
649 continue;
650 CbmRichRing* Ring = nullptr;
651 if (richInd > -1) { Ring = static_cast<CbmRichRing*>(fRichRings->At(richInd)); }
652
653 // Doing fit with Fit To primary Vertex and calculate chi2 to primary vertex
654 double chi2 = 0;
655 TVector3 Momentum = CbmKresFunctions::FitToVertexAndGetChi(stsTrack, fKFVertex.GetRefX(), fKFVertex.GetRefY(),
656 fKFVertex.GetRefZ(), chi2);
657 const FairTrackParam* track_par = stsTrack->GetParamFirst();
658 double charge = track_par->GetQp();
659
660 if (chi2 != chi2) continue;
661 if (chi2 == 0) continue;
662
663 if (chi2 > 3) { SaveOutsideTracks(mcTrack, stsTrack, charge, stsInd, richInd, stsMcTrackId, Ring); }
664 if (chi2 > 3) continue;
665
666 SaveTargetTracks(mcTrack, stsTrack, Momentum, charge, stsInd, richInd, stsMcTrackId, Ring);
667 }
668
669 FindGammasTarget(fEventNumEta, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Target,
674
675 FindGammasOutside(fEventNumEta, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Outside,
679
681
688
695
702
703 int numformix = 500;
704 if (fEventNumEta % numformix == 0) {
718 EMT_eta_gg_Event_Both.clear();
721 }
722
723 if (fEventNumEta % numformix == 0) {
728 }
729
730 if (fEventNumEta % numformix == 0) {
735 }
736}
737
738
739void CbmKresEta::SaveOutsideTracks(CbmMCTrack* mcTrack1, CbmStsTrack* stsTrack, double charge, int stsInd, int richInd,
740 int stsMcTrackId, CbmRichRing* RING)
741{
742 int InRich = FindInRich(richInd, stsMcTrackId);
743 if (charge < 0) {
744 VMCtracks_minus_Outside.push_back(mcTrack1);
745 VStsTrack_minus_Outside.push_back(stsTrack);
746 VRings_minus_Outside.push_back(InRich);
747 VStsIndex_minus_Outside.push_back(stsInd);
748 VRichRing_minus_Outside.push_back(RING);
749 VMCIndex_minus_Outside.push_back(stsMcTrackId);
750 }
751 if (charge > 0) {
752 VMCtracks_plus_Outside.push_back(mcTrack1);
753 VStsTrack_plus_Outside.push_back(stsTrack);
754 VRings_plus_Outside.push_back(InRich);
755 VStsIndex_plus_Outside.push_back(stsInd);
756 VRichRing_plus_Outside.push_back(RING);
757 VMCIndex_plus_Outside.push_back(stsMcTrackId);
758 }
759}
760
761
762void CbmKresEta::SaveTargetTracks(CbmMCTrack* mcTrack1, CbmStsTrack* stsTrack, TVector3 refmom, double charge,
763 int stsInd, int richInd, int stsMcTrackId, CbmRichRing* RING)
764{
765 int InRich = FindInRich(richInd, stsMcTrackId);
766 if (charge < 0) {
767 VMCtracks_minus_Target.push_back(mcTrack1);
768 VStsTrack_minus_Target.push_back(stsTrack);
769 VMomenta_minus_Target.push_back(refmom);
770 VRings_minus_Target.push_back(InRich);
771 VStsIndex_minus_Target.push_back(stsInd);
772 VRichRing_minus_Target.push_back(RING);
773 VMCIndex_minus_Target.push_back(stsMcTrackId);
774 }
775 if (charge > 0) {
776 VMCtracks_plus_Target.push_back(mcTrack1);
777 VStsTrack_plus_Target.push_back(stsTrack);
778 VMomenta_plus_Target.push_back(refmom);
779 VRings_plus_Target.push_back(InRich);
780 VStsIndex_plus_Target.push_back(stsInd);
781 VRichRing_plus_Target.push_back(RING);
782 VMCIndex_plus_Target.push_back(stsMcTrackId);
783 }
784}
785
786int CbmKresEta::FindInRich(int richInd, int stsMcTrackId)
787{
788 int RingsInRich = 0;
789 if (richInd > -1) {
790 CbmTrackMatchNew* richMatch = (CbmTrackMatchNew*) fRichRingMatches->At(richInd);
791 if (richMatch != nullptr && richMatch->GetNofLinks() > 0) {
792 int richMcTrackId = richMatch->GetMatchedLink().GetIndex();
793 if (richMcTrackId > 0) {
794 if (stsMcTrackId == richMcTrackId) { // check that global track was matched correctly for STS and RICH together
795 CbmMCTrack* mcTrack2 = (CbmMCTrack*) fMcTracks->At(richMcTrackId);
796 if (mcTrack2 != nullptr) {
797 int pdgRICH = mcTrack2->GetPdgCode();
798 if (TMath::Abs(pdgRICH) == 11) RingsInRich++;
799 }
800 }
801 }
802 }
803 }
804 return RingsInRich;
805}
806
807int CbmKresEta::CheckIfElectron(CbmRichRing* ring, double momentum)
808{
809 int identified = 0;
810
811 if (nullptr != ring) {
812 CbmRichRingLight ringHit;
813 int nofHits = ring->GetNofHits();
814 for (int i = 0; i < nofHits; i++) {
815 int hitInd = ring->GetHit(i);
816 CbmRichHit* hit = (CbmRichHit*) fRichHits->At(hitInd);
817 if (nullptr == hit) continue;
818 CbmRichHitLight hl(hit->GetX(), hit->GetY());
819 ringHit.AddHit(hl);
820 }
821 fTauFit->DoFit(&ringHit);
822 if (ringHit.GetAaxis() > 4 && ringHit.GetAaxis() < 6 && ringHit.GetBaxis() > 4 && ringHit.GetBaxis() < 6
823 && momentum > 0.2 && momentum < 4.)
824 identified++;
825 }
826
827 return identified;
828}
829
831{
832 double FinalAngle = 400;
833 Int_t hits1sts = Sts_1->GetNofStsHits();
834 Int_t hits2sts = Sts_2->GetNofStsHits();
835 Int_t hits1mvd = Sts_1->GetNofMvdHits();
836 Int_t hits2mvd = Sts_2->GetNofMvdHits();
837
838 double Xpart1 = 0;
839 double Ypart1 = 0;
840 double Zpart1 = 0;
841 double Xpart2 = 0;
842 double Ypart2 = 0;
843 double Zpart2 = 0;
844 if (hits1sts > 0) {
845 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(hits1sts - 1);
846 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
847 Xpart1 = stsHit1->GetX();
848 Ypart1 = stsHit1->GetY();
849 Zpart1 = stsHit1->GetZ();
850 }
851 else {
852 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(hits1mvd - 1);
853 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
854 Xpart1 = mvdHit1->GetX();
855 Ypart1 = mvdHit1->GetY();
856 Zpart1 = mvdHit1->GetZ();
857 }
858
859 if (hits2sts > 0) {
860 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(hits2sts - 1);
861 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
862 Xpart2 = stsHit2->GetX();
863 Ypart2 = stsHit2->GetY();
864 Zpart2 = stsHit2->GetZ();
865 }
866 else {
867 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(hits2mvd - 1);
868 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
869 Xpart2 = mvdHit2->GetX();
870 Ypart2 = mvdHit2->GetY();
871 Zpart2 = mvdHit2->GetZ();
872 }
873
874 // check difference in 2 cm, because of two slices of every STS and Mvd stations
875 if (TMath::Abs(Zpart1 - Zpart2) > 2
876 && Zpart1 > Zpart2) { // if last hits are in different stations --> try to find the latest common station
877 for (int i = hits1sts - 2; i > -1; i--) { // start from second last station
878 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
879 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
880 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
881 Xpart1 = stsHit->GetX();
882 Ypart1 = stsHit->GetY();
883 Zpart1 = stsHit->GetZ();
884 }
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;
888 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
889 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
890 Xpart1 = mvdHit->GetX();
891 Ypart1 = mvdHit->GetY();
892 Zpart1 = mvdHit->GetZ();
893 }
894 }
895 }
896
897 if (TMath::Abs(Zpart1 - Zpart2) > 2
898 && Zpart1 < Zpart2) { // if last hits are in different stations --> try to find the latest common station
899 for (int i = hits2sts - 2; i > -1; i--) { // start from second last station
900 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
901 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
902 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
903 Xpart2 = stsHit->GetX();
904 Ypart2 = stsHit->GetY();
905 Zpart2 = stsHit->GetZ();
906 }
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;
910 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
911 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
912 Xpart2 = mvdHit->GetX();
913 Ypart2 = mvdHit->GetY();
914 Zpart2 = mvdHit->GetZ();
915 }
916 }
917 }
918
919 // calculate angle if we found common station
920 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
921 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
922 }
923
924 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
925}
926
928{
929 double FinalAngle = 400;
930 Int_t hits1sts = Sts_1->GetNofStsHits();
931 Int_t hits2sts = Sts_2->GetNofStsHits();
932 Int_t hits1mvd = Sts_1->GetNofMvdHits();
933 Int_t hits2mvd = Sts_2->GetNofMvdHits();
934
935 double Xpart1 = 0;
936 double Ypart1 = 0;
937 double Zpart1 = 0;
938 double Xpart2 = 0;
939 double Ypart2 = 0;
940 double Zpart2 = 0;
941 if (hits1mvd > 0) {
942 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(0);
943 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
944 Xpart1 = mvdHit1->GetX();
945 Ypart1 = mvdHit1->GetY();
946 Zpart1 = mvdHit1->GetZ();
947 }
948 else {
949 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(0);
950 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
951 Xpart1 = stsHit1->GetX();
952 Ypart1 = stsHit1->GetY();
953 Zpart1 = stsHit1->GetZ();
954 }
955
956 if (hits2mvd > 0) {
957 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(0);
958 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
959 Xpart2 = mvdHit2->GetX();
960 Ypart2 = mvdHit2->GetY();
961 Zpart2 = mvdHit2->GetZ();
962 }
963 else {
964 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(0);
965 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
966 Xpart2 = stsHit2->GetX();
967 Ypart2 = stsHit2->GetY();
968 Zpart2 = stsHit2->GetZ();
969 }
970
971 // check difference in 2 cm, because of two slices of every STS and Mvd stations
972 if (TMath::Abs(Zpart1 - Zpart2) > 2
973 && Zpart1 < Zpart2) { // if first hits are in different stations --> try to find the earliest common station
974 for (int i = 1; i < hits1mvd; i++) { // start from second hit
975 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
976 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
977 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
978 Xpart1 = mvdHit->GetX();
979 Ypart1 = mvdHit->GetY();
980 Zpart1 = mvdHit->GetZ();
981 }
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;
985 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
986 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
987 Xpart1 = stsHit->GetX();
988 Ypart1 = stsHit->GetY();
989 Zpart1 = stsHit->GetZ();
990 }
991 }
992 }
993
994 if (TMath::Abs(Zpart1 - Zpart2) > 2
995 && Zpart1 > Zpart2) { // if first hits are in different stations --> try to find the earliest common station
996 for (int i = 1; i < hits2mvd; i++) { // start from second hit
997 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
998 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
999 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1000 Xpart2 = mvdHit->GetX();
1001 Ypart2 = mvdHit->GetY();
1002 Zpart2 = mvdHit->GetZ();
1003 }
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;
1007 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
1008 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1009 Xpart2 = stsHit->GetX();
1010 Ypart2 = stsHit->GetY();
1011 Zpart2 = stsHit->GetZ();
1012 }
1013 }
1014 }
1015
1016 // calculate angle if we found common station
1017 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1018 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
1019 }
1020
1021 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1022}
1023
1024
1025void CbmKresEta::FindGammasTarget(int EventNumEta, double AngleCut, double InvMassCut, int RealPID,
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)
1033{
1034 for (size_t i = 0; i < Momenta_minus.size(); i++) {
1035 for (size_t j = 0; j < Momenta_plus.size(); j++) {
1036
1037 CbmStsTrack* sts1 = StsTrack_minus.at(i);
1038 CbmStsTrack* sts2 = StsTrack_plus.at(j);
1039
1040 TVector3 part1 = Momenta_minus[i];
1041 TVector3 part2 = Momenta_plus[j];
1042 CbmMCTrack* part1MC = MCtracks_minus[i];
1043 CbmMCTrack* part2MC = MCtracks_plus[j];
1044
1045 int richcheck_0 = 0;
1046 int richcheck_1 = 0;
1047 if (RealPID == 1) {
1048 // Real PID
1049 richcheck_0 = CheckIfElectron(richRing_minus[i], Momenta_minus[i].Mag());
1050 richcheck_1 = CheckIfElectron(richRing_plus[j], Momenta_plus[j].Mag());
1051 }
1052 else {
1053 // MC PID
1054 richcheck_0 = Rings_minus[i];
1055 richcheck_1 = Rings_plus[j];
1056 }
1057 int richcheck = richcheck_0 + richcheck_1;
1058
1059 double InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
1060 double OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
1061 double PlaneAngle_last = CalculatePlaneAngle_last(sts1, sts2);
1062 double PlaneAngle_first = CalculatePlaneAngle_first(sts1, sts2);
1064
1065 int TruePair =
1066 0; //correctly reconstructed photon from eta meson: 0 - means wrong pair combination; 1 - means correct pair
1067 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
1068 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
1069 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
1070 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
1071 if (mcTrackgrmama->GetPdgCode() == 221) { TruePair = 1; }
1072 }
1073 }
1074
1075 if (TruePair == 1) {
1076 if (richcheck == 0) {
1081 PlaneAngles_last_fromEta_zeroInRich_Target->Fill(PlaneAngle_last);
1082 PlaneAngles_first_fromEta_zeroInRich_Target->Fill(PlaneAngle_first);
1083
1088 PlaneAngles_first_fromEta_zeroInRich_Both->Fill(PlaneAngle_first);
1089 PlaneAngles_last_fromEta_zeroInRich_Both->Fill(PlaneAngle_last);
1090 }
1091 if (richcheck == 1) {
1096 PlaneAngles_last_fromEta_oneInRich_Target->Fill(PlaneAngle_last);
1097 PlaneAngles_first_fromEta_oneInRich_Target->Fill(PlaneAngle_first);
1098
1103 PlaneAngles_first_fromEta_oneInRich_Both->Fill(PlaneAngle_first);
1104 PlaneAngles_last_fromEta_oneInRich_Both->Fill(PlaneAngle_last);
1105 }
1106 if (richcheck == 2) {
1111 PlaneAngles_last_fromEta_twoInRich_Target->Fill(PlaneAngle_last);
1112 PlaneAngles_first_fromEta_twoInRich_Target->Fill(PlaneAngle_first);
1113
1118 PlaneAngles_first_fromEta_twoInRich_Both->Fill(PlaneAngle_first);
1119 PlaneAngles_last_fromEta_twoInRich_Both->Fill(PlaneAngle_last);
1120 }
1121 Particle_pull_px_Target->Fill(part1MC->GetPx() - part1.X());
1122 Particle_pull_py_Target->Fill(part1MC->GetPy() - part1.Y());
1123 Particle_pull_pz_Target->Fill(part1MC->GetPz() - part1.Z());
1124 Particle_pull_px_Target->Fill(part2MC->GetPx() - part2.X());
1125 Particle_pull_py_Target->Fill(part2MC->GetPy() - part2.Y());
1126 Particle_pull_pz_Target->Fill(part2MC->GetPz() - part2.Z());
1127 }
1128 else {
1129 if (richcheck == 0) {
1134 PlaneAngles_last_wrongpairs_zeroInRich_Target->Fill(PlaneAngle_last);
1135 PlaneAngles_first_wrongpairs_zeroInRich_Target->Fill(PlaneAngle_first);
1136
1141 PlaneAngles_first_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_first);
1142 PlaneAngles_last_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_last);
1143 }
1144 if (richcheck == 1) {
1149 PlaneAngles_last_wrongpairs_oneInRich_Target->Fill(PlaneAngle_last);
1150 PlaneAngles_first_wrongpairs_oneInRich_Target->Fill(PlaneAngle_first);
1151
1156 PlaneAngles_first_wrongpairs_oneInRich_Both->Fill(PlaneAngle_first);
1157 PlaneAngles_last_wrongpairs_oneInRich_Both->Fill(PlaneAngle_last);
1158 }
1159 if (richcheck == 2) {
1164 PlaneAngles_last_wrongpairs_twoInRich_Target->Fill(PlaneAngle_last);
1165 PlaneAngles_first_wrongpairs_twoInRich_Target->Fill(PlaneAngle_first);
1166
1171 PlaneAngles_first_wrongpairs_twoInRich_Both->Fill(PlaneAngle_first);
1172 PlaneAngles_last_wrongpairs_twoInRich_Both->Fill(PlaneAngle_last);
1173 }
1174 }
1175
1176
1177 // cuts for gamma reconstruction
1178 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1179 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1180
1181 if (richcheck == 1 && TruePair == 1) {
1184 }
1185 if (richcheck == 1 && TruePair == 0) {
1188 }
1189
1190 frefmomenta.clear();
1191 frefmomenta.push_back(part1);
1192 frefmomenta.push_back(part2);
1193 frefId.clear();
1194 frefId.push_back(stsIndex_minus[i]);
1195 frefId.push_back(stsIndex_plus[j]);
1196 fMCId.clear();
1197 fMCId.push_back(MCIndex_minus[i]);
1198 fMCId.push_back(MCIndex_plus[j]);
1199 fMCtracks.clear();
1200 fMCtracks.push_back(part1MC);
1201 fMCtracks.push_back(part2MC);
1202
1203 // for event mixing
1204 EMT_eta_gg_Event_Target.push_back(EventNumEta);
1206 EMT_eta_gg_NofRings_Target.push_back(richcheck);
1207
1208
1209 // everything (RICH == 0, RICH == 1, RICH == 2) together
1210 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1211 Gammas_all_Target.push_back(frefmomenta);
1215 }
1216
1217 // only cases, when RICH == 2
1218 if (richcheck == 2) {
1219 Gammas_two_Target.push_back(frefmomenta);
1223 }
1224
1225 // cases, when RICH == 1 or RICH == 2 together
1226 if (richcheck == 1 || richcheck == 2) {
1231 }
1232 }
1233 }
1234}
1235
1236
1237void CbmKresEta::FindGammasOutside(int EventNumEta, double AngleCut, double InvMassCut, int RealPID,
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)
1246{
1247 for (size_t i = 0; i < StsTrack_minus_Outside.size(); i++) {
1248 for (size_t j = 0; j < StsTrack_plus_Outside.size(); j++) {
1249
1250 CbmStsTrack* sts1 = StsTrack_minus_Outside[i];
1251 CbmStsTrack* sts2 = StsTrack_plus_Outside[j];
1252
1253 CbmMCTrack* part1MC = MCtracks_minus_Outside[i];
1254 CbmMCTrack* part2MC = MCtracks_plus_Outside[j];
1255
1256 KFParticle electron;
1258 KFParticle positron;
1260 const KFParticle* daughters[2] = {&electron, &positron};
1261 KFParticle intersection;
1262 intersection.Construct(daughters, 2);
1263
1264 if (intersection.GetZ() > 75 || intersection.GetZ() < -5) continue; // kick weird intersections
1265
1266 // fit to the vertex fitter
1267 // TVector3 part1 = CbmKresFunctions::FitToVertex(sts1, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1268 // TVector3 part2 = CbmKresFunctions::FitToVertex(sts2, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1269
1270 double chi_electron = 0;
1271 double chi_positron = 0;
1272 TVector3 part1 = CbmKresFunctions::FitToVertexAndGetChi(sts1, intersection.GetX(), intersection.GetY(),
1273 intersection.GetZ(), chi_electron);
1274 TVector3 part2 = CbmKresFunctions::FitToVertexAndGetChi(sts2, intersection.GetX(), intersection.GetY(),
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;
1282
1283 int richcheck = 0;
1284 int richcheck_0 = 0;
1285 int richcheck_1 = 0;
1286 if (RealPID == 1) {
1287 // Real PID
1288 richcheck_0 = CheckIfElectron(richRing_minus_Outside[i], part1.Mag());
1289 richcheck_1 = CheckIfElectron(richRing_plus_Outside[j], part2.Mag());
1290 richcheck = richcheck_0 + richcheck_1;
1291 }
1292 else {
1293 // MC PID
1294 richcheck = Rings_minus_Outside[i] + Rings_plus_Outside[j];
1295 }
1296
1297 double InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
1298 double OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
1299 double PlaneAngle_last = CalculatePlaneAngle_last(sts1, sts2);
1300 double PlaneAngle_first = CalculatePlaneAngle_first(sts1, sts2);
1302
1303 int TruePair =
1304 0; //correctly reconstructed photon from eta meson: 0 - means wrong pair combination; 1 - means correct pair
1305 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
1306 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
1307 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
1308 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
1309 if (mcTrackgrmama->GetPdgCode() == 221) { TruePair = 1; }
1310 }
1311 }
1312
1313 if (TruePair == 1) {
1314 if (richcheck == 0) {
1319 PlaneAngles_last_fromEta_zeroInRich_Outside->Fill(PlaneAngle_last);
1320 PlaneAngles_first_fromEta_zeroInRich_Outside->Fill(PlaneAngle_first);
1321
1326 PlaneAngles_first_fromEta_zeroInRich_Both->Fill(PlaneAngle_first);
1327 PlaneAngles_last_fromEta_zeroInRich_Both->Fill(PlaneAngle_last);
1328 }
1329 if (richcheck == 1) {
1334 PlaneAngles_last_fromEta_oneInRich_Outside->Fill(PlaneAngle_last);
1335 PlaneAngles_first_fromEta_oneInRich_Outside->Fill(PlaneAngle_first);
1336
1341 PlaneAngles_first_fromEta_oneInRich_Both->Fill(PlaneAngle_first);
1342 PlaneAngles_last_fromEta_oneInRich_Both->Fill(PlaneAngle_last);
1343 }
1344 if (richcheck == 2) {
1349 PlaneAngles_last_fromEta_twoInRich_Outside->Fill(PlaneAngle_last);
1350 PlaneAngles_first_fromEta_twoInRich_Outside->Fill(PlaneAngle_first);
1351
1356 PlaneAngles_first_fromEta_twoInRich_Both->Fill(PlaneAngle_first);
1357 PlaneAngles_last_fromEta_twoInRich_Both->Fill(PlaneAngle_last);
1358 }
1359 Particle_pull_px_Outside->Fill(part1MC->GetPx() - part1.X());
1360 Particle_pull_py_Outside->Fill(part1MC->GetPy() - part1.Y());
1361 Particle_pull_pz_Outside->Fill(part1MC->GetPz() - part1.Z());
1362 Particle_pull_px_Outside->Fill(part2MC->GetPx() - part2.X());
1363 Particle_pull_py_Outside->Fill(part2MC->GetPy() - part2.Y());
1364 Particle_pull_pz_Outside->Fill(part2MC->GetPz() - part2.Z());
1365 Particle_pull_X_Outside->Fill(part1MC->GetStartX() - intersection.GetX());
1366 Particle_pull_Y_Outside->Fill(part1MC->GetStartY() - intersection.GetY());
1367 Particle_pull_Z_Outside->Fill(part1MC->GetStartZ() - intersection.GetZ());
1368 Particle_pull_X_Outside->Fill(part2MC->GetStartX() - intersection.GetX());
1369 Particle_pull_Y_Outside->Fill(part2MC->GetStartY() - intersection.GetY());
1370 Particle_pull_Z_Outside->Fill(part2MC->GetStartZ() - intersection.GetZ());
1371 }
1372 else {
1373 if (richcheck == 0) {
1380
1385 PlaneAngles_first_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_first);
1386 PlaneAngles_last_wrongpairs_zeroInRich_Both->Fill(PlaneAngle_last);
1387 }
1388 if (richcheck == 1) {
1393 PlaneAngles_last_wrongpairs_oneInRich_Outside->Fill(PlaneAngle_last);
1394 PlaneAngles_first_wrongpairs_oneInRich_Outside->Fill(PlaneAngle_first);
1395
1400 PlaneAngles_first_wrongpairs_oneInRich_Both->Fill(PlaneAngle_first);
1401 PlaneAngles_last_wrongpairs_oneInRich_Both->Fill(PlaneAngle_last);
1402 }
1403 if (richcheck == 2) {
1408 PlaneAngles_last_wrongpairs_twoInRich_Outside->Fill(PlaneAngle_last);
1409 PlaneAngles_first_wrongpairs_twoInRich_Outside->Fill(PlaneAngle_first);
1410
1415 PlaneAngles_first_wrongpairs_twoInRich_Both->Fill(PlaneAngle_first);
1416 PlaneAngles_last_wrongpairs_twoInRich_Both->Fill(PlaneAngle_last);
1417 }
1418 }
1419
1420
1421 // cuts for gamma reconstruction
1422 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1423 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1424
1425
1426 if (richcheck == 1 && TruePair == 1) {
1429 }
1430 if (richcheck == 1 && TruePair == 0) {
1433 }
1434
1435
1436 frefmomenta.clear();
1437 frefmomenta.push_back(part1);
1438 frefmomenta.push_back(part2);
1439 frefId.clear();
1440 frefId.push_back(stsIndex_minus_Outside[i]);
1441 frefId.push_back(stsIndex_plus_Outside[j]);
1442 fMCId.clear();
1443 fMCId.push_back(MCIndex_minus_Outside[i]);
1444 fMCId.push_back(MCIndex_plus_Outside[j]);
1445 fMCtracks.clear();
1446 fMCtracks.push_back(part1MC);
1447 fMCtracks.push_back(part2MC);
1448
1449 // for event mixing
1450 EMT_eta_gg_Event_Outside.push_back(EventNumEta);
1452 EMT_eta_gg_NofRings_Outside.push_back(richcheck);
1453
1454
1455 // everything (RICH == 0, RICH == 1, RICH == 2) together
1456 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1461 }
1462
1463 // only cases, when RICH == 2
1464 if (richcheck == 2) {
1469 }
1470
1471 // cases, when RICH == 1 or RICH == 2 together
1472 if (richcheck == 1 || richcheck == 2) {
1477 }
1478 }
1479 }
1480}
1481
1483{
1490
1503
1516
1525}
1526
1527
1528void CbmKresEta::FindEta(TString /*mod*/, TString position, vector<vector<TVector3>> Gammas,
1529 vector<vector<int>> StsIndex, vector<vector<int>> MCIndex,
1530 vector<vector<CbmMCTrack*>> GammasMC, vector<TH1*> gg, vector<TH1*> rap_pt_separation)
1531{
1532 if (Gammas.size() < 2) return; // min 2 gammas to form pi0 are required
1533 for (size_t gamma1 = 0; gamma1 < Gammas.size() - 1; gamma1++) {
1534 for (size_t gamma2 = gamma1 + 1; gamma2 < Gammas.size(); gamma2++) {
1535 // 4 reconstructed particles from gammas
1536 TVector3 e11 = Gammas[gamma1][0];
1537 TVector3 e12 = Gammas[gamma1][1];
1538 TVector3 e21 = Gammas[gamma2][0];
1539 TVector3 e22 = Gammas[gamma2][1];
1540
1541 // MC true data for this particles
1542 CbmMCTrack* mcTrack1 = GammasMC[gamma1][0];
1543 CbmMCTrack* mcTrack2 = GammasMC[gamma1][1];
1544 CbmMCTrack* mcTrack3 = GammasMC[gamma2][0];
1545 CbmMCTrack* mcTrack4 = GammasMC[gamma2][1];
1546
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])
1549 continue; // particles are not used twice --> different
1550
1552 double openingAngleBetweenGammasReco =
1554
1555 int TrueEta = 0;
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) {
1562 if (TMath::Abs(mcTrack1->GetPdgCode()) == 11 && TMath::Abs(mcTrack2->GetPdgCode()) == 11
1563 && TMath::Abs(mcTrack3->GetPdgCode()) == 11 && TMath::Abs(mcTrack4->GetPdgCode()) == 11) {
1564 if ((mcTrack1->GetMotherId() == mcTrack2->GetMotherId() && mcTrack3->GetMotherId() == mcTrack4->GetMotherId())
1565 || (mcTrack1->GetMotherId() == mcTrack3->GetMotherId()
1566 && mcTrack2->GetMotherId() == mcTrack4->GetMotherId())
1567 || (mcTrack1->GetMotherId() == mcTrack4->GetMotherId()
1568 && mcTrack2->GetMotherId() == mcTrack3->GetMotherId())) {
1569 CbmMCTrack* MothTrack1 = (CbmMCTrack*) fMcTracks->At(mcTrack1->GetMotherId());
1570 CbmMCTrack* MothTrack2 = (CbmMCTrack*) fMcTracks->At(mcTrack2->GetMotherId());
1571 CbmMCTrack* MothTrack3 = (CbmMCTrack*) fMcTracks->At(mcTrack3->GetMotherId());
1572 CbmMCTrack* MothTrack4 = (CbmMCTrack*) fMcTracks->At(mcTrack4->GetMotherId());
1573 if (MothTrack1->GetPdgCode() == 22 && MothTrack2->GetPdgCode() == 22 && MothTrack3->GetPdgCode() == 22
1574 && MothTrack4->GetPdgCode() == 22 && MothTrack1->GetMotherId() != -1
1575 && MothTrack1->GetMotherId() == MothTrack2->GetMotherId()
1576 && MothTrack1->GetMotherId() == MothTrack3->GetMotherId()
1577 && MothTrack1->GetMotherId() == MothTrack4->GetMotherId()) {
1578 CbmMCTrack* GrandMothTrack1 = (CbmMCTrack*) fMcTracks->At(MothTrack1->GetMotherId());
1579 if (GrandMothTrack1->GetPdgCode() == 221) {
1580 TrueEta = 1;
1581 // cout << " " << mod << "; " << position << " ***** CbmKresEta ***** " << endl;
1582 // cout << "Decay eta -> gamma gamma -> e+e- e+e- detected!\t\t reco mass: " << params.fMinv << endl;
1583 // cout << "motherids: " << mcTrack1->GetMotherId() << "/" << mcTrack2->GetMotherId() << "/" << mcTrack3->GetMotherId() << "/" << mcTrack4->GetMotherId() << endl;
1584 // cout << "grandmotherid: " << MothTrack1->GetMotherId() << "/" << MothTrack2->GetMotherId() << "/" << MothTrack3->GetMotherId() << "/" << MothTrack4->GetMotherId() << endl;
1585 // cout << "pdgs: " << GrandMothTrack1->GetPdgCode() << "-->" << MothTrack1->GetPdgCode() << "/" << MothTrack2->GetPdgCode() << "/" << MothTrack3->GetPdgCode() << "/" << MothTrack4->GetPdgCode() << endl;
1586 // cout << "particle 1: \t" << mcTrack1->GetPdgCode() << ";\t pt = " << mcTrack1->GetPt() << ";\t X = " << mcTrack1->GetStartX() << ";\t Y = " << mcTrack1->GetStartY() << ";\t Z = " << mcTrack1->GetStartZ() << ";\t E = " << mcTrack1->GetEnergy() << endl;
1587 // cout << "particle 2: \t" << mcTrack2->GetPdgCode() << ";\t pt = " << mcTrack2->GetPt() << ";\t X = " << mcTrack2->GetStartX() << ";\t Y = " << mcTrack2->GetStartY() << ";\t Z = " << mcTrack2->GetStartZ() << ";\t E = " << mcTrack2->GetEnergy() << endl;
1588 // cout << "particle 3: \t" << mcTrack3->GetPdgCode() << ";\t pt = " << mcTrack3->GetPt() << ";\t X = " << mcTrack3->GetStartX() << ";\t Y = " << mcTrack3->GetStartY() << ";\t Z = " << mcTrack3->GetStartZ() << ";\t E = " << mcTrack3->GetEnergy() << endl;
1589 // cout << "particle 4: \t" << mcTrack4->GetPdgCode() << ";\t pt = " << mcTrack4->GetPt() << ";\t X = " << mcTrack4->GetStartX() << ";\t Y = " << mcTrack4->GetStartY() << ";\t Z = " << mcTrack4->GetStartZ() << ";\t E = " << mcTrack4->GetEnergy() << endl;
1590 // cout << " ***** CbmKresEta (End) ***** " << endl;
1591 }
1592 }
1593 }
1594 }
1595 }
1596
1597 if (TrueEta == 1) gg[0]->Fill(openingAngleBetweenGammasReco);
1598 if (TrueEta == 0) gg[1]->Fill(openingAngleBetweenGammasReco);
1599
1600 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1601
1602 gg[2]->Fill(params.fMinv);
1603
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()));
1608
1609 // separate by rap and Pt only for case "Both"
1610 if (position == "Both") {
1611 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
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);
1616 }
1617 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
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);
1622 }
1623 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
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);
1628 }
1629 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
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);
1633 }
1634 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
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);
1637 }
1638 }
1639
1640 if (TrueEta == 1) {
1641 gg[4]->Fill(params.fRapidity, params.fPt);
1642 gg[6]->Fill(params.fRapidity, params.fPt);
1643 gg[8]->Fill(params.fMinv);
1644 gg[10]->Fill(openingAngleBetweenGammasReco);
1645 }
1646
1647 if (TrueEta == 0) {
1648 gg[5]->Fill(params.fRapidity, params.fPt);
1649 gg[7]->Fill(params.fRapidity, params.fPt);
1650 gg[9]->Fill(params.fMinv);
1651 gg[11]->Fill(openingAngleBetweenGammasReco);
1652 }
1653 }
1654 }
1655}
1656
1658// TARGET
1659{
1660 int nof_Target = EMT_eta_gg_Event_Target.size();
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++) {
1665 continue; // to make sure that the photons are from two different events
1666 TVector3 e11 = EMT_eta_gg_pair_momenta_Target[a][0];
1667 TVector3 e12 = EMT_eta_gg_pair_momenta_Target[a][1];
1668 TVector3 e21 = EMT_eta_gg_pair_momenta_Target[b][0];
1669 TVector3 e22 = EMT_eta_gg_pair_momenta_Target[b][1];
1671 double openingAngleBetweenGammasReco =
1673 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1674
1675 EMT_Eta_InvMass_all_Target->Fill(params.fMinv);
1677 EMT_Eta_InvMass_two_Target->Fill(params.fMinv);
1681 }
1682 }
1683}
1684
1685
1687// OUTSIDE
1688{
1689 Int_t nof_Outside = EMT_eta_gg_Event_Outside.size();
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++) {
1694 continue; // to make sure that the photons are from two different events
1695 TVector3 e11 = EMT_eta_gg_pair_momenta_Outside[a][0];
1696 TVector3 e12 = EMT_eta_gg_pair_momenta_Outside[a][1];
1697 TVector3 e21 = EMT_eta_gg_pair_momenta_Outside[b][0];
1698 TVector3 e22 = EMT_eta_gg_pair_momenta_Outside[b][1];
1700 double openingAngleBetweenGammasReco =
1702 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1703
1704 EMT_Eta_InvMass_all_Outside->Fill(params.fMinv);
1706 EMT_Eta_InvMass_two_Outside->Fill(params.fMinv);
1710 }
1711 }
1712}
1713
1714void CbmKresEta::Mixing_Both(vector<TH1*> rap_pt_separation_all, vector<TH1*> rap_pt_separation_onetwo,
1715 vector<TH1*> rap_pt_separation_two)
1716// BOTH
1717{
1718 Int_t nof_Both = EMT_eta_gg_Event_Both.size();
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++) {
1723 continue; // to make sure that the photons are from two different events
1724 TVector3 e11 = EMT_eta_gg_pair_momenta_Both[a][0];
1725 TVector3 e12 = EMT_eta_gg_pair_momenta_Both[a][1];
1726 TVector3 e21 = EMT_eta_gg_pair_momenta_Both[b][0];
1727 TVector3 e22 = EMT_eta_gg_pair_momenta_Both[b][1];
1729 double openingAngleBetweenGammasReco =
1731 if (openingAngleBetweenGammasReco < 10 || openingAngleBetweenGammasReco > 40) continue;
1732
1733 EMT_Eta_InvMass_all_Both->Fill(params.fMinv);
1735 EMT_Eta_InvMass_two_Both->Fill(params.fMinv);
1736 if ((EMT_eta_gg_NofRings_Both[a] == 1 || EMT_eta_gg_NofRings_Both[a] == 2)
1737 && (EMT_eta_gg_NofRings_Both[b] == 1 || EMT_eta_gg_NofRings_Both[b] == 2))
1738 EMT_Eta_InvMass_onetwo_Both->Fill(params.fMinv);
1739
1740 // separate by rap and Pt only for all
1741 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
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);
1746 }
1747 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
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);
1752 }
1753 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
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);
1758 }
1759 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
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);
1763 }
1764 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
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);
1767 }
1768
1769 // separate by rap and Pt only for onetwo
1770 if ((EMT_eta_gg_NofRings_Both[a] == 1 || EMT_eta_gg_NofRings_Both[a] == 2)
1771 && (EMT_eta_gg_NofRings_Both[b] == 1 || EMT_eta_gg_NofRings_Both[b] == 2)) {
1772 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
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);
1777 }
1778 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
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);
1783 }
1784 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
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);
1789 }
1790 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
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);
1794 }
1795 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
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);
1798 }
1799 }
1800
1801 // separate by rap and Pt only for two
1802 if (EMT_eta_gg_NofRings_Both[a] == 2 && EMT_eta_gg_NofRings_Both[b] == 2) {
1803 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
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);
1808 }
1809 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
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);
1814 }
1815 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
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);
1820 }
1821 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
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);
1825 }
1826 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
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);
1829 }
1830 }
1831 }
1832 }
1833}
1834
1835
1837{
1838 gDirectory->mkdir("Eta");
1839 gDirectory->cd("Eta");
1840
1841 gDirectory->mkdir("Cuts_for_Eta_gg");
1842 gDirectory->cd("Cuts_for_Eta_gg");
1843
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");
1848 for (UInt_t i = 0; i < fHistoList_Eta_cuts_zeroInRich_Target.size(); i++) {
1850 }
1851 gDirectory->cd("..");
1852 gDirectory->mkdir("Cuts_for_Eta_gg_oneInRich_Target");
1853 gDirectory->cd("Cuts_for_Eta_gg_oneInRich_Target");
1854 for (UInt_t i = 0; i < fHistoList_Eta_cuts_oneInRich_Target.size(); i++) {
1856 }
1857 gDirectory->cd("..");
1858 gDirectory->mkdir("Cuts_for_Eta_gg_twoInRich_Target");
1859 gDirectory->cd("Cuts_for_Eta_gg_twoInRich_Target");
1860 for (UInt_t i = 0; i < fHistoList_Eta_cuts_twoInRich_Target.size(); i++) {
1862 }
1863 gDirectory->cd("..");
1864 for (UInt_t i = 0; i < fHistoList_Eta_cuts_Target.size(); i++) {
1865 fHistoList_Eta_cuts_Target[i]->Write();
1866 }
1867 gDirectory->cd("..");
1868
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");
1873 for (UInt_t i = 0; i < fHistoList_Eta_cuts_zeroInRich_Outside.size(); i++) {
1875 }
1876 gDirectory->cd("..");
1877 gDirectory->mkdir("Cuts_for_Eta_gg_oneInRich_Outside");
1878 gDirectory->cd("Cuts_for_Eta_gg_oneInRich_Outside");
1879 for (UInt_t i = 0; i < fHistoList_Eta_cuts_oneInRich_Outside.size(); i++) {
1881 }
1882 gDirectory->cd("..");
1883 gDirectory->mkdir("Cuts_for_Eta_gg_twoInRich_Outside");
1884 gDirectory->cd("Cuts_for_Eta_gg_twoInRich_Outside");
1885 for (UInt_t i = 0; i < fHistoList_Eta_cuts_twoInRich_Outside.size(); i++) {
1887 }
1888 gDirectory->cd("..");
1889 for (UInt_t i = 0; i < fHistoList_Eta_cuts_Outside.size(); i++) {
1890 fHistoList_Eta_cuts_Outside[i]->Write();
1891 }
1892 gDirectory->cd("..");
1893
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");
1898 for (UInt_t i = 0; i < fHistoList_Eta_cuts_zeroInRich_Both.size(); i++) {
1900 }
1901 gDirectory->cd("..");
1902 gDirectory->mkdir("Cuts_for_Eta_gg_oneInRich_Both");
1903 gDirectory->cd("Cuts_for_Eta_gg_oneInRich_Both");
1904 for (UInt_t i = 0; i < fHistoList_Eta_cuts_oneInRich_Both.size(); i++) {
1906 }
1907 gDirectory->cd("..");
1908 gDirectory->mkdir("Cuts_for_Eta_gg_twoInRich_Both");
1909 gDirectory->cd("Cuts_for_Eta_gg_twoInRich_Both");
1910 for (UInt_t i = 0; i < fHistoList_Eta_cuts_twoInRich_Both.size(); i++) {
1912 }
1913 gDirectory->cd("..");
1914 for (UInt_t i = 0; i < fHistoList_Eta_cuts_Both.size(); i++) {
1915 fHistoList_Eta_cuts_Both[i]->Write();
1916 }
1917 gDirectory->cd("..");
1918
1919 gDirectory->cd("..");
1920
1921
1922 gDirectory->mkdir("Eta_gg_Target");
1923 gDirectory->cd("Eta_gg_Target");
1924
1925 gDirectory->mkdir("all");
1926 gDirectory->cd("all");
1927 for (UInt_t i = 0; i < fHistoList_Eta_all_Target.size(); i++) {
1928 fHistoList_Eta_all_Target[i]->Write();
1929 }
1930 gDirectory->cd("..");
1931
1932 gDirectory->mkdir("onetwo");
1933 gDirectory->cd("onetwo");
1934 for (UInt_t i = 0; i < fHistoList_Eta_onetwo_Target.size(); i++) {
1935 fHistoList_Eta_onetwo_Target[i]->Write();
1936 }
1937 gDirectory->cd("..");
1938
1939 gDirectory->mkdir("two");
1940 gDirectory->cd("two");
1941 for (UInt_t i = 0; i < fHistoList_Eta_two_Target.size(); i++) {
1942 fHistoList_Eta_two_Target[i]->Write();
1943 }
1944 gDirectory->cd("..");
1945
1946 gDirectory->cd("..");
1947
1948
1949 gDirectory->mkdir("Eta_gg_Outside");
1950 gDirectory->cd("Eta_gg_Outside");
1951
1952 gDirectory->mkdir("all");
1953 gDirectory->cd("all");
1954 for (UInt_t i = 0; i < fHistoList_Eta_all_Outside.size(); i++) {
1955 fHistoList_Eta_all_Outside[i]->Write();
1956 }
1957 gDirectory->cd("..");
1958
1959 gDirectory->mkdir("onetwo");
1960 gDirectory->cd("onetwo");
1961 for (UInt_t i = 0; i < fHistoList_Eta_onetwo_Outside.size(); i++) {
1963 }
1964 gDirectory->cd("..");
1965
1966 gDirectory->mkdir("two");
1967 gDirectory->cd("two");
1968 for (UInt_t i = 0; i < fHistoList_Eta_two_Outside.size(); i++) {
1969 fHistoList_Eta_two_Outside[i]->Write();
1970 }
1971 gDirectory->cd("..");
1972
1973 gDirectory->cd("..");
1974
1975
1976 gDirectory->mkdir("Eta_gg_Both");
1977 gDirectory->cd("Eta_gg_Both");
1978
1979 gDirectory->mkdir("all");
1980 gDirectory->cd("all");
1981 gDirectory->mkdir("multi_all_Both");
1982 gDirectory->cd("multi_all_Both");
1983 for (UInt_t i = 0; i < fHistoList_rap_pt_Eta_all_Both.size(); i++) {
1985 }
1986 gDirectory->cd("..");
1987 for (UInt_t i = 0; i < fHistoList_Eta_all_Both.size(); i++) {
1988 fHistoList_Eta_all_Both[i]->Write();
1989 }
1990 gDirectory->cd("..");
1991
1992 gDirectory->mkdir("onetwo");
1993 gDirectory->cd("onetwo");
1994 gDirectory->mkdir("multi_onetwo_Both");
1995 gDirectory->cd("multi_onetwo_Both");
1996 for (UInt_t i = 0; i < fHistoList_rap_pt_Eta_onetwo_Both.size(); i++) {
1998 }
1999 gDirectory->cd("..");
2000 for (UInt_t i = 0; i < fHistoList_Eta_onetwo_Both.size(); i++) {
2001 fHistoList_Eta_onetwo_Both[i]->Write();
2002 }
2003 gDirectory->cd("..");
2004
2005 gDirectory->mkdir("two");
2006 gDirectory->cd("two");
2007 gDirectory->mkdir("multi_two_Both");
2008 gDirectory->cd("multi_two_Both");
2009 for (UInt_t i = 0; i < fHistoList_rap_pt_Eta_two_Both.size(); i++) {
2011 }
2012 gDirectory->cd("..");
2013 for (UInt_t i = 0; i < fHistoList_Eta_two_Both.size(); i++) {
2014 fHistoList_Eta_two_Both[i]->Write();
2015 }
2016 gDirectory->cd("..");
2017
2018 gDirectory->cd("..");
2019
2020
2021 for (UInt_t i = 0; i < fHistoList_Eta.size(); i++) {
2022 fHistoList_Eta[i]->Write();
2023 }
2024
2025
2026 gDirectory->cd("..");
2027}
2028
2030{
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);
2058
2061 new TH1D("GammasInvMass_fromEta_beforeCuts_zeroInRich_Target",
2062 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Target; invariant "
2063 "mass in GeV/c^{2}",
2064 510, -0.01, 0.5);
2067 new TH1D("GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target",
2068 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target; opening "
2069 "angle in degree",
2070 300, -0.1, 29.9);
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}",
2092 510, -0.01, 0.5);
2095 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target",
2096 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target; opening "
2097 "angle in degree",
2098 300, -0.1, 29.9);
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}",
2121 510, -0.01, 0.5);
2124 new TH1D("GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target",
2125 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target; opening "
2126 "angle in degree",
2127 300, -0.1, 29.9);
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}",
2149 510, -0.01, 0.5);
2152 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target",
2153 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target; opening "
2154 "angle in degree",
2155 300, -0.1, 29.9);
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}",
2178 510, -0.01, 0.5);
2181 new TH1D("GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target",
2182 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target; opening "
2183 "angle in degree",
2184 300, -0.1, 29.9);
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}",
2206 510, -0.01, 0.5);
2209 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target",
2210 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target; opening "
2211 "angle in degree",
2212 300, -0.1, 29.9);
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.);
2230
2231
2234 new TH1D("GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside",
2235 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside; invariant "
2236 "mass in GeV/c^{2}",
2237 510, -0.01, 0.5);
2240 new TH1D("GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside",
2241 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside; opening "
2242 "angle in degree",
2243 300, -0.1, 29.9);
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}",
2265 510, -0.01, 0.5);
2268 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside",
2269 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside; "
2270 "opening angle in degree",
2271 300, -0.1, 29.9);
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}",
2294 510, -0.01, 0.5);
2297 new TH1D("GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside",
2298 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside; opening "
2299 "angle in degree",
2300 300, -0.1, 29.9);
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}",
2322 510, -0.01, 0.5);
2325 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside",
2326 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside; opening "
2327 "angle in degree",
2328 300, -0.1, 29.9);
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}",
2351 510, -0.01, 0.5);
2354 new TH1D("GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside",
2355 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside; opening "
2356 "angle in degree",
2357 300, -0.1, 29.9);
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}",
2379 510, -0.01, 0.5);
2382 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside",
2383 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside; opening "
2384 "angle in degree",
2385 300, -0.1, 29.9);
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.);
2403
2404
2407 new TH1D("GammasInvMass_fromEta_beforeCuts_zeroInRich_Both",
2408 "GammasInvMass_fromEta_beforeCuts_zeroInRich_Both; invariant mass "
2409 "in GeV/c^{2}",
2410 510, -0.01, 0.5);
2413 new TH1D("GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both",
2414 "GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both; opening "
2415 "angle in degree",
2416 300, -0.1, 29.9);
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}",
2438 510, -0.01, 0.5);
2441 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both",
2442 "GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both; opening "
2443 "angle in degree",
2444 300, -0.1, 29.9);
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 "
2466 "in GeV/c^{2}",
2467 510, -0.01, 0.5);
2470 new TH1D("GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both",
2471 "GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both; opening angle "
2472 "in degree",
2473 300, -0.1, 29.9);
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}",
2495 510, -0.01, 0.5);
2498 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both",
2499 "GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both; opening "
2500 "angle in degree",
2501 300, -0.1, 29.9);
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 "
2522 "in degree;#",
2523 720, -1., 179.);
2526 new TH1D("PlaneAngles_first_fromEta_afterCuts_oneInRich_Both",
2527 "PlaneAngles_first_fromEta_afterCuts_oneInRich_Both; #theta angle "
2528 "in degree;#",
2529 720, -1., 179.);
2532 new TH1D("PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both",
2533 "PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both; #theta "
2534 "angle in degree;#",
2535 720, -1., 179.);
2538 new TH1D("PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both",
2539 "PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both; #theta "
2540 "angle in degree;#",
2541 720, -1., 179.);
2545 new TH1D("GammasInvMass_fromEta_beforeCuts_twoInRich_Both",
2546 "GammasInvMass_fromEta_beforeCuts_twoInRich_Both; invariant mass "
2547 "in GeV/c^{2}",
2548 510, -0.01, 0.5);
2551 new TH1D("GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both",
2552 "GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both; opening angle "
2553 "in degree",
2554 300, -0.1, 29.9);
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}",
2576 510, -0.01, 0.5);
2579 new TH1D("GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both",
2580 "GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both; opening "
2581 "angle in degree",
2582 300, -0.1, 29.9);
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.);
2600
2601
2605 new TH1D("AngleBetweenGammas_trueEta_before_cuts_all_Target",
2606 "AngleBetweenGammas_trueEta_before_cuts_all_Target (between "
2607 "#gamma#gamma from #eta); angle [deg]",
2608 900, -0.1, 89.9);
2611 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_all_Target",
2612 "AngleBetweenGammas_wrongEta_before_cuts_all_Target (between "
2613 "#gamma#gamma); angle [deg]",
2614 900, -0.1, 89.9);
2617 new TH1D("InvMass_GammaCombinations_all_Target",
2618 "InvMass_GammaCombinations_all_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2620 Particles_PDG_all_Target = new TH1D("Particles_PDG_all_Target", "Particles_PDG_all_Target ; PDG", 2300, 0, 2300);
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.,
2624 60, -1., 5.);
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.,
2628 60, -1., 5.);
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,
2632 0., 4., 5, 0., 2.);
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,
2636 0., 4., 5, 0., 2.);
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]",
2650 900, -0.1, 89.9);
2653 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_all_Target",
2654 "WrongEta_AngleBetweenGammas_after_cuts_all_Target (between "
2655 "#gamma#gamma); angle [deg]",
2656 900, -0.1, 89.9);
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]",
2666 900, -0.1, 89.9);
2669 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target",
2670 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target (between "
2671 "#gamma#gamma); angle [deg]",
2672 900, -0.1, 89.9);
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.,
2683 7., 60, -1., 5.);
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]",
2709 900, -0.1, 89.9);
2712 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target",
2713 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target (between "
2714 "#gamma#gamma); angle [deg]",
2715 900, -0.1, 89.9);
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]",
2725 900, -0.1, 89.9);
2728 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_two_Target",
2729 "AngleBetweenGammas_wrongEta_before_cuts_two_Target (between "
2730 "#gamma#gamma); angle [deg]",
2731 900, -0.1, 89.9);
2734 new TH1D("InvMass_GammaCombinations_two_Target",
2735 "InvMass_GammaCombinations_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2737 Particles_PDG_two_Target = new TH1D("Particles_PDG_two_Target", "Particles_PDG_two_Target ; PDG", 2300, 0, 2300);
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.,
2741 60, -1., 5.);
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.,
2745 60, -1., 5.);
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,
2749 0., 4., 5, 0., 2.);
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,
2753 0., 4., 5, 0., 2.);
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]",
2767 900, -0.1, 89.9);
2770 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_two_Target",
2771 "WrongEta_AngleBetweenGammas_after_cuts_two_Target (between "
2772 "#gamma#gamma); angle [deg]",
2773 900, -0.1, 89.9);
2776 new TH1D("EMT_Eta_InvMass_two_Target", "EMT_Eta_InvMass_two_Target; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2779
2780
2784 new TH1D("AngleBetweenGammas_trueEta_before_cuts_all_Outside",
2785 "AngleBetweenGammas_trueEta_before_cuts_all_Outside (between "
2786 "#gamma#gamma from #eta); angle [deg]",
2787 900, -0.1, 89.9);
2790 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_all_Outside",
2791 "AngleBetweenGammas_wrongEta_before_cuts_all_Outside (between "
2792 "#gamma#gamma); angle [deg]",
2793 900, -0.1, 89.9);
2796 new TH1D("InvMass_GammaCombinations_all_Outside",
2797 "InvMass_GammaCombinations_all_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2799 Particles_PDG_all_Outside = new TH1D("Particles_PDG_all_Outside", "Particles_PDG_all_Outside ; PDG", 2300, 0, 2300);
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.,
2803 60, -1., 5.);
2806 new TH2D("WrongEta_pt_vs_rap_all_Outside", "WrongEta_pt_vs_rap_all_Outside; rapidity y; p_{t} in GeV/c ", 90, -2.,
2807 7., 60, -1., 5.);
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,
2811 0., 4., 5, 0., 2.);
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]",
2829 900, -0.1, 89.9);
2832 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_all_Outside",
2833 "WrongEta_AngleBetweenGammas_after_cuts_all_Outside (between "
2834 "#gamma#gamma); angle [deg]",
2835 900, -0.1, 89.9);
2837 EMT_Eta_InvMass_all_Outside = new TH1D("EMT_Eta_InvMass_all_Outside",
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]",
2845 900, -0.1, 89.9);
2848 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside",
2849 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside (between "
2850 "#gamma#gamma); angle [deg]",
2851 900, -0.1, 89.9);
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]",
2888 900, -0.1, 89.9);
2891 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside",
2892 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside (between "
2893 "#gamma#gamma); angle [deg]",
2894 900, -0.1, 89.9);
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]",
2904 900, -0.1, 89.9);
2907 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_two_Outside",
2908 "AngleBetweenGammas_wrongEta_before_cuts_two_Outside (between "
2909 "#gamma#gamma); angle [deg]",
2910 900, -0.1, 89.9);
2913 new TH1D("InvMass_GammaCombinations_two_Outside",
2914 "InvMass_GammaCombinations_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2916 Particles_PDG_two_Outside = new TH1D("Particles_PDG_two_Outside", "Particles_PDG_two_Outside ; PDG", 2300, 0, 2300);
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.,
2920 60, -1., 5.);
2923 new TH2D("WrongEta_pt_vs_rap_two_Outside", "WrongEta_pt_vs_rap_two_Outside; rapidity y; p_{t} in GeV/c ", 90, -2.,
2924 7., 60, -1., 5.);
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,
2928 0., 4., 5, 0., 2.);
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]",
2946 900, -0.1, 89.9);
2949 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_two_Outside",
2950 "WrongEta_AngleBetweenGammas_after_cuts_two_Outside (between "
2951 "#gamma#gamma); angle [deg]",
2952 900, -0.1, 89.9);
2954 EMT_Eta_InvMass_two_Outside = new TH1D("EMT_Eta_InvMass_two_Outside",
2955 "EMT_Eta_InvMass_two_Outside; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
2958
2959
2962 AngleBetweenGammas_trueEta_before_cuts_all_Both = new TH1D("AngleBetweenGammas_trueEta_before_cuts_all_Both",
2963 "AngleBetweenGammas_trueEta_before_cuts_all_Both (between "
2964 "#gamma#gamma from #eta); angle [deg]",
2965 900, -0.1, 89.9);
2968 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_all_Both",
2969 "AngleBetweenGammas_wrongEta_before_cuts_all_Both (between "
2970 "#gamma#gamma); angle [deg]",
2971 900, -0.1, 89.9);
2974 new TH1D("InvMass_GammaCombinations_all_Both", "InvMass_GammaCombinations_all_Both; invariant mass in GeV/c^{2};#",
2975 200, 0.0, 2.0);
2977 Particles_PDG_all_Both = new TH1D("Particles_PDG_all_Both", "Particles_PDG_all_Both ; PDG", 2300, 0, 2300);
2979 TrueEta_pt_vs_rap_all_Both = new TH2D(
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,
2984 -1., 5.);
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.,
2988 4., 5, 0., 2.);
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.,
2992 4., 5, 0., 2.);
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);
3002 TrueEta_AngleBetweenGammas_after_cuts_all_Both = new TH1D("TrueEta_AngleBetweenGammas_after_cuts_all_Both",
3003 "TrueEta_AngleBetweenGammas_after_cuts_all_Both (between "
3004 "#gamma#gamma from #eta); angle [deg]",
3005 900, -0.1, 89.9);
3007 WrongEta_AngleBetweenGammas_after_cuts_all_Both = new TH1D("WrongEta_AngleBetweenGammas_after_cuts_all_Both",
3008 "WrongEta_AngleBetweenGammas_after_cuts_all_Both (between "
3009 "#gamma#gamma); angle [deg]",
3010 900, -0.1, 89.9);
3013 new TH1D("EMT_Eta_InvMass_all_Both", "EMT_Eta_InvMass_all_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3017 multi_InvMass_Eta_all_Both_1 = new TH1D("multi_InvMass_Eta_all_Both_1",
3018 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3019 "invariant mass in GeV/c^{2};#",
3020 200, 0.0, 2.0);
3022 multi_InvMass_Eta_all_Both_2 = new TH1D("multi_InvMass_Eta_all_Both_2",
3023 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3024 "invariant mass in GeV/c^{2};#",
3025 200, 0.0, 2.0);
3027 multi_InvMass_Eta_all_Both_3 = new TH1D("multi_InvMass_Eta_all_Both_3",
3028 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3029 "invariant mass in GeV/c^{2};#",
3030 200, 0.0, 2.0);
3032 multi_InvMass_Eta_all_Both_4 = new TH1D("multi_InvMass_Eta_all_Both_4",
3033 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3034 "invariant mass in GeV/c^{2};#",
3035 200, 0.0, 2.0);
3037 multi_InvMass_Eta_all_Both_5 = new TH1D("multi_InvMass_Eta_all_Both_5",
3038 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3039 "invariant mass in GeV/c^{2};#",
3040 200, 0.0, 2.0);
3042 multi_InvMass_Eta_all_Both_6 = new TH1D("multi_InvMass_Eta_all_Both_6",
3043 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3044 "invariant mass in GeV/c^{2};#",
3045 200, 0.0, 2.0);
3047 multi_InvMass_Eta_all_Both_7 = new TH1D("multi_InvMass_Eta_all_Both_7",
3048 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3049 "invariant mass in GeV/c^{2};#",
3050 200, 0.0, 2.0);
3052 multi_InvMass_Eta_all_Both_8 = new TH1D("multi_InvMass_Eta_all_Both_8",
3053 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3054 "invariant mass in GeV/c^{2};#",
3055 200, 0.0, 2.0);
3057 multi_InvMass_Eta_all_Both_9 = new TH1D("multi_InvMass_Eta_all_Both_9",
3058 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3059 "invariant mass in GeV/c^{2};#",
3060 200, 0.0, 2.0);
3062 multi_InvMass_Eta_all_Both_10 = new TH1D("multi_InvMass_Eta_all_Both_10",
3063 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3064 "invariant mass in GeV/c^{2};#",
3065 200, 0.0, 2.0);
3067 multi_InvMass_Eta_all_Both_11 = new TH1D("multi_InvMass_Eta_all_Both_11",
3068 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3069 "invariant mass in GeV/c^{2};#",
3070 200, 0.0, 2.0);
3072 multi_InvMass_Eta_all_Both_12 = new TH1D("multi_InvMass_Eta_all_Both_12",
3073 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3074 "invariant mass in GeV/c^{2};#",
3075 200, 0.0, 2.0);
3077 multi_InvMass_Eta_all_Both_13 = new TH1D("multi_InvMass_Eta_all_Both_13",
3078 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3079 "invariant mass in GeV/c^{2};#",
3080 200, 0.0, 2.0);
3082 multi_InvMass_Eta_all_Both_14 = new TH1D("multi_InvMass_Eta_all_Both_14",
3083 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3084 "invariant mass in GeV/c^{2};#",
3085 200, 0.0, 2.0);
3087 multi_InvMass_Eta_all_Both_15 = new TH1D("multi_InvMass_Eta_all_Both_15",
3088 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3089 "invariant mass in GeV/c^{2};#",
3090 200, 0.0, 2.0);
3092 multi_InvMass_Eta_all_Both_16 = new TH1D("multi_InvMass_Eta_all_Both_16",
3093 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3094 "invariant mass in GeV/c^{2};#",
3095 200, 0.0, 2.0);
3097 multi_InvMass_Eta_all_Both_17 = new TH1D("multi_InvMass_Eta_all_Both_17",
3098 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3099 "invariant mass in GeV/c^{2};#",
3100 200, 0.0, 2.0);
3102 multi_EMT_Eta_all_Both_1 = new TH1D("multi_EMT_Eta_all_Both_1",
3103 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3104 "invariant mass in GeV/c^{2};#",
3105 200, 0.0, 2.0);
3107 multi_EMT_Eta_all_Both_2 = new TH1D("multi_EMT_Eta_all_Both_2",
3108 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3109 "invariant mass in GeV/c^{2};#",
3110 200, 0.0, 2.0);
3112 multi_EMT_Eta_all_Both_3 = new TH1D("multi_EMT_Eta_all_Both_3",
3113 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3114 "invariant mass in GeV/c^{2};#",
3115 200, 0.0, 2.0);
3117 multi_EMT_Eta_all_Both_4 = new TH1D("multi_EMT_Eta_all_Both_4",
3118 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3119 "invariant mass in GeV/c^{2};#",
3120 200, 0.0, 2.0);
3122 multi_EMT_Eta_all_Both_5 = new TH1D("multi_EMT_Eta_all_Both_5",
3123 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3124 "invariant mass in GeV/c^{2};#",
3125 200, 0.0, 2.0);
3127 multi_EMT_Eta_all_Both_6 = new TH1D("multi_EMT_Eta_all_Both_6",
3128 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3129 "invariant mass in GeV/c^{2};#",
3130 200, 0.0, 2.0);
3132 multi_EMT_Eta_all_Both_7 = new TH1D("multi_EMT_Eta_all_Both_7",
3133 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3134 "invariant mass in GeV/c^{2};#",
3135 200, 0.0, 2.0);
3137 multi_EMT_Eta_all_Both_8 = new TH1D("multi_EMT_Eta_all_Both_8",
3138 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3139 "invariant mass in GeV/c^{2};#",
3140 200, 0.0, 2.0);
3142 multi_EMT_Eta_all_Both_9 = new TH1D("multi_EMT_Eta_all_Both_9",
3143 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3144 "invariant mass in GeV/c^{2};#",
3145 200, 0.0, 2.0);
3147 multi_EMT_Eta_all_Both_10 = new TH1D("multi_EMT_Eta_all_Both_10",
3148 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3149 "invariant mass in GeV/c^{2};#",
3150 200, 0.0, 2.0);
3152 multi_EMT_Eta_all_Both_11 = new TH1D("multi_EMT_Eta_all_Both_11",
3153 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3154 "invariant mass in GeV/c^{2};#",
3155 200, 0.0, 2.0);
3157 multi_EMT_Eta_all_Both_12 = new TH1D("multi_EMT_Eta_all_Both_12",
3158 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3159 "invariant mass in GeV/c^{2};#",
3160 200, 0.0, 2.0);
3162 multi_EMT_Eta_all_Both_13 = new TH1D("multi_EMT_Eta_all_Both_13",
3163 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3164 "invariant mass in GeV/c^{2};#",
3165 200, 0.0, 2.0);
3167 multi_EMT_Eta_all_Both_14 = new TH1D("multi_EMT_Eta_all_Both_14",
3168 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3169 "invariant mass in GeV/c^{2};#",
3170 200, 0.0, 2.0);
3172 multi_EMT_Eta_all_Both_15 = new TH1D("multi_EMT_Eta_all_Both_15",
3173 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3174 "invariant mass in GeV/c^{2};#",
3175 200, 0.0, 2.0);
3177 multi_EMT_Eta_all_Both_16 = new TH1D("multi_EMT_Eta_all_Both_16",
3178 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3179 "invariant mass in GeV/c^{2};#",
3180 200, 0.0, 2.0);
3182 multi_EMT_Eta_all_Both_17 = new TH1D("multi_EMT_Eta_all_Both_17",
3183 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3184 "invariant mass in GeV/c^{2};#",
3185 200, 0.0, 2.0);
3187
3188
3191 new TH1D("AngleBetweenGammas_trueEta_before_cuts_onetwo_Both",
3192 "AngleBetweenGammas_trueEta_before_cuts_onetwo_Both (between "
3193 "#gamma#gamma from #eta); angle [deg]",
3194 900, -0.1, 89.9);
3197 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both",
3198 "AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both (between "
3199 "#gamma#gamma); angle [deg]",
3200 900, -0.1, 89.9);
3203 new TH1D("InvMass_GammaCombinations_onetwo_Both",
3204 "InvMass_GammaCombinations_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3206 Particles_PDG_onetwo_Both = new TH1D("Particles_PDG_onetwo_Both", "Particles_PDG_onetwo_Both ; PDG", 2300, 0, 2300);
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.,
3210 60, -1., 5.);
3213 new TH2D("WrongEta_pt_vs_rap_onetwo_Both", "WrongEta_pt_vs_rap_onetwo_Both; rapidity y; p_{t} in GeV/c ", 90, -2.,
3214 7., 60, -1., 5.);
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,
3218 0., 4., 5, 0., 2.);
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]",
3236 900, -0.1, 89.9);
3239 new TH1D("WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both",
3240 "WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both (between "
3241 "#gamma#gamma); angle [deg]",
3242 900, -0.1, 89.9);
3244 EMT_Eta_InvMass_onetwo_Both = new TH1D("EMT_Eta_InvMass_onetwo_Both",
3245 "EMT_Eta_InvMass_onetwo_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3249 multi_InvMass_Eta_onetwo_Both_1 = new TH1D("multi_InvMass_Eta_onetwo_Both_1",
3250 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3251 "invariant mass in GeV/c^{2};#",
3252 200, 0.0, 2.0);
3254 multi_InvMass_Eta_onetwo_Both_2 = new TH1D("multi_InvMass_Eta_onetwo_Both_2",
3255 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3256 "invariant mass in GeV/c^{2};#",
3257 200, 0.0, 2.0);
3259 multi_InvMass_Eta_onetwo_Both_3 = new TH1D("multi_InvMass_Eta_onetwo_Both_3",
3260 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3261 "invariant mass in GeV/c^{2};#",
3262 200, 0.0, 2.0);
3264 multi_InvMass_Eta_onetwo_Both_4 = new TH1D("multi_InvMass_Eta_onetwo_Both_4",
3265 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3266 "invariant mass in GeV/c^{2};#",
3267 200, 0.0, 2.0);
3269 multi_InvMass_Eta_onetwo_Both_5 = new TH1D("multi_InvMass_Eta_onetwo_Both_5",
3270 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3271 "invariant mass in GeV/c^{2};#",
3272 200, 0.0, 2.0);
3274 multi_InvMass_Eta_onetwo_Both_6 = new TH1D("multi_InvMass_Eta_onetwo_Both_6",
3275 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3276 "invariant mass in GeV/c^{2};#",
3277 200, 0.0, 2.0);
3279 multi_InvMass_Eta_onetwo_Both_7 = new TH1D("multi_InvMass_Eta_onetwo_Both_7",
3280 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3281 "invariant mass in GeV/c^{2};#",
3282 200, 0.0, 2.0);
3284 multi_InvMass_Eta_onetwo_Both_8 = new TH1D("multi_InvMass_Eta_onetwo_Both_8",
3285 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3286 "invariant mass in GeV/c^{2};#",
3287 200, 0.0, 2.0);
3289 multi_InvMass_Eta_onetwo_Both_9 = new TH1D("multi_InvMass_Eta_onetwo_Both_9",
3290 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3291 "invariant mass in GeV/c^{2};#",
3292 200, 0.0, 2.0);
3294 multi_InvMass_Eta_onetwo_Both_10 = new TH1D("multi_InvMass_Eta_onetwo_Both_10",
3295 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3296 "invariant mass in GeV/c^{2};#",
3297 200, 0.0, 2.0);
3299 multi_InvMass_Eta_onetwo_Both_11 = new TH1D("multi_InvMass_Eta_onetwo_Both_11",
3300 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3301 "invariant mass in GeV/c^{2};#",
3302 200, 0.0, 2.0);
3304 multi_InvMass_Eta_onetwo_Both_12 = new TH1D("multi_InvMass_Eta_onetwo_Both_12",
3305 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3306 "invariant mass in GeV/c^{2};#",
3307 200, 0.0, 2.0);
3309 multi_InvMass_Eta_onetwo_Both_13 = new TH1D("multi_InvMass_Eta_onetwo_Both_13",
3310 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3311 "invariant mass in GeV/c^{2};#",
3312 200, 0.0, 2.0);
3314 multi_InvMass_Eta_onetwo_Both_14 = new TH1D("multi_InvMass_Eta_onetwo_Both_14",
3315 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3316 "invariant mass in GeV/c^{2};#",
3317 200, 0.0, 2.0);
3319 multi_InvMass_Eta_onetwo_Both_15 = new TH1D("multi_InvMass_Eta_onetwo_Both_15",
3320 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3321 "invariant mass in GeV/c^{2};#",
3322 200, 0.0, 2.0);
3324 multi_InvMass_Eta_onetwo_Both_16 = new TH1D("multi_InvMass_Eta_onetwo_Both_16",
3325 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3326 "invariant mass in GeV/c^{2};#",
3327 200, 0.0, 2.0);
3329 multi_InvMass_Eta_onetwo_Both_17 = new TH1D("multi_InvMass_Eta_onetwo_Both_17",
3330 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3331 "invariant mass in GeV/c^{2};#",
3332 200, 0.0, 2.0);
3334 multi_EMT_Eta_onetwo_Both_1 = new TH1D("multi_EMT_Eta_onetwo_Both_1",
3335 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3336 "invariant mass in GeV/c^{2};#",
3337 200, 0.0, 2.0);
3339 multi_EMT_Eta_onetwo_Both_2 = new TH1D("multi_EMT_Eta_onetwo_Both_2",
3340 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3341 "invariant mass in GeV/c^{2};#",
3342 200, 0.0, 2.0);
3344 multi_EMT_Eta_onetwo_Both_3 = new TH1D("multi_EMT_Eta_onetwo_Both_3",
3345 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3346 "invariant mass in GeV/c^{2};#",
3347 200, 0.0, 2.0);
3349 multi_EMT_Eta_onetwo_Both_4 = new TH1D("multi_EMT_Eta_onetwo_Both_4",
3350 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3351 "invariant mass in GeV/c^{2};#",
3352 200, 0.0, 2.0);
3354 multi_EMT_Eta_onetwo_Both_5 = new TH1D("multi_EMT_Eta_onetwo_Both_5",
3355 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3356 "invariant mass in GeV/c^{2};#",
3357 200, 0.0, 2.0);
3359 multi_EMT_Eta_onetwo_Both_6 = new TH1D("multi_EMT_Eta_onetwo_Both_6",
3360 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3361 "invariant mass in GeV/c^{2};#",
3362 200, 0.0, 2.0);
3364 multi_EMT_Eta_onetwo_Both_7 = new TH1D("multi_EMT_Eta_onetwo_Both_7",
3365 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3366 "invariant mass in GeV/c^{2};#",
3367 200, 0.0, 2.0);
3369 multi_EMT_Eta_onetwo_Both_8 = new TH1D("multi_EMT_Eta_onetwo_Both_8",
3370 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3371 "invariant mass in GeV/c^{2};#",
3372 200, 0.0, 2.0);
3374 multi_EMT_Eta_onetwo_Both_9 = new TH1D("multi_EMT_Eta_onetwo_Both_9",
3375 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3376 "invariant mass in GeV/c^{2};#",
3377 200, 0.0, 2.0);
3379 multi_EMT_Eta_onetwo_Both_10 = new TH1D("multi_EMT_Eta_onetwo_Both_10",
3380 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3381 "invariant mass in GeV/c^{2};#",
3382 200, 0.0, 2.0);
3384 multi_EMT_Eta_onetwo_Both_11 = new TH1D("multi_EMT_Eta_onetwo_Both_11",
3385 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3386 "invariant mass in GeV/c^{2};#",
3387 200, 0.0, 2.0);
3389 multi_EMT_Eta_onetwo_Both_12 = new TH1D("multi_EMT_Eta_onetwo_Both_12",
3390 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3391 "invariant mass in GeV/c^{2};#",
3392 200, 0.0, 2.0);
3394 multi_EMT_Eta_onetwo_Both_13 = new TH1D("multi_EMT_Eta_onetwo_Both_13",
3395 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3396 "invariant mass in GeV/c^{2};#",
3397 200, 0.0, 2.0);
3399 multi_EMT_Eta_onetwo_Both_14 = new TH1D("multi_EMT_Eta_onetwo_Both_14",
3400 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3401 "invariant mass in GeV/c^{2};#",
3402 200, 0.0, 2.0);
3404 multi_EMT_Eta_onetwo_Both_15 = new TH1D("multi_EMT_Eta_onetwo_Both_15",
3405 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3406 "invariant mass in GeV/c^{2};#",
3407 200, 0.0, 2.0);
3409 multi_EMT_Eta_onetwo_Both_16 = new TH1D("multi_EMT_Eta_onetwo_Both_16",
3410 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3411 "invariant mass in GeV/c^{2};#",
3412 200, 0.0, 2.0);
3414 multi_EMT_Eta_onetwo_Both_17 = new TH1D("multi_EMT_Eta_onetwo_Both_17",
3415 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3416 "invariant mass in GeV/c^{2};#",
3417 200, 0.0, 2.0);
3419
3420
3422 AngleBetweenGammas_trueEta_before_cuts_two_Both = new TH1D("AngleBetweenGammas_trueEta_before_cuts_two_Both",
3423 "AngleBetweenGammas_trueEta_before_cuts_two_Both (between "
3424 "#gamma#gamma from #eta); angle [deg]",
3425 900, -0.1, 89.9);
3428 new TH1D("AngleBetweenGammas_wrongEta_before_cuts_two_Both",
3429 "AngleBetweenGammas_wrongEta_before_cuts_two_Both (between "
3430 "#gamma#gamma); angle [deg]",
3431 900, -0.1, 89.9);
3434 new TH1D("InvMass_GammaCombinations_two_Both", "InvMass_GammaCombinations_two_Both; invariant mass in GeV/c^{2};#",
3435 200, 0.0, 2.0);
3437 Particles_PDG_two_Both = new TH1D("Particles_PDG_two_Both", "Particles_PDG_two_Both ; PDG", 2300, 0, 2300);
3439 TrueEta_pt_vs_rap_two_Both = new TH2D(
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,
3444 -1., 5.);
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.,
3448 4., 5, 0., 2.);
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.,
3452 4., 5, 0., 2.);
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);
3462 TrueEta_AngleBetweenGammas_after_cuts_two_Both = new TH1D("TrueEta_AngleBetweenGammas_after_cuts_two_Both",
3463 "TrueEta_AngleBetweenGammas_after_cuts_two_Both (between "
3464 "#gamma#gamma from #eta); angle [deg]",
3465 900, -0.1, 89.9);
3467 WrongEta_AngleBetweenGammas_after_cuts_two_Both = new TH1D("WrongEta_AngleBetweenGammas_after_cuts_two_Both",
3468 "WrongEta_AngleBetweenGammas_after_cuts_two_Both (between "
3469 "#gamma#gamma); angle [deg]",
3470 900, -0.1, 89.9);
3473 new TH1D("EMT_Eta_InvMass_two_Both", "EMT_Eta_InvMass_two_Both; invariant mass in GeV/c^{2};#", 200, 0.0, 2.0);
3477 multi_InvMass_Eta_two_Both_1 = new TH1D("multi_InvMass_Eta_two_Both_1",
3478 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3479 "invariant mass in GeV/c^{2};#",
3480 200, 0.0, 2.0);
3482 multi_InvMass_Eta_two_Both_2 = new TH1D("multi_InvMass_Eta_two_Both_2",
3483 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3484 "invariant mass in GeV/c^{2};#",
3485 200, 0.0, 2.0);
3487 multi_InvMass_Eta_two_Both_3 = new TH1D("multi_InvMass_Eta_two_Both_3",
3488 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3489 "invariant mass in GeV/c^{2};#",
3490 200, 0.0, 2.0);
3492 multi_InvMass_Eta_two_Both_4 = new TH1D("multi_InvMass_Eta_two_Both_4",
3493 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3494 "invariant mass in GeV/c^{2};#",
3495 200, 0.0, 2.0);
3497 multi_InvMass_Eta_two_Both_5 = new TH1D("multi_InvMass_Eta_two_Both_5",
3498 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3499 "invariant mass in GeV/c^{2};#",
3500 200, 0.0, 2.0);
3502 multi_InvMass_Eta_two_Both_6 = new TH1D("multi_InvMass_Eta_two_Both_6",
3503 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3504 "invariant mass in GeV/c^{2};#",
3505 200, 0.0, 2.0);
3507 multi_InvMass_Eta_two_Both_7 = new TH1D("multi_InvMass_Eta_two_Both_7",
3508 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3509 "invariant mass in GeV/c^{2};#",
3510 200, 0.0, 2.0);
3512 multi_InvMass_Eta_two_Both_8 = new TH1D("multi_InvMass_Eta_two_Both_8",
3513 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3514 "invariant mass in GeV/c^{2};#",
3515 200, 0.0, 2.0);
3517 multi_InvMass_Eta_two_Both_9 = new TH1D("multi_InvMass_Eta_two_Both_9",
3518 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3519 "invariant mass in GeV/c^{2};#",
3520 200, 0.0, 2.0);
3522 multi_InvMass_Eta_two_Both_10 = new TH1D("multi_InvMass_Eta_two_Both_10",
3523 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3524 "invariant mass in GeV/c^{2};#",
3525 200, 0.0, 2.0);
3527 multi_InvMass_Eta_two_Both_11 = new TH1D("multi_InvMass_Eta_two_Both_11",
3528 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3529 "invariant mass in GeV/c^{2};#",
3530 200, 0.0, 2.0);
3532 multi_InvMass_Eta_two_Both_12 = new TH1D("multi_InvMass_Eta_two_Both_12",
3533 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3534 "invariant mass in GeV/c^{2};#",
3535 200, 0.0, 2.0);
3537 multi_InvMass_Eta_two_Both_13 = new TH1D("multi_InvMass_Eta_two_Both_13",
3538 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3539 "invariant mass in GeV/c^{2};#",
3540 200, 0.0, 2.0);
3542 multi_InvMass_Eta_two_Both_14 = new TH1D("multi_InvMass_Eta_two_Both_14",
3543 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3544 "invariant mass in GeV/c^{2};#",
3545 200, 0.0, 2.0);
3547 multi_InvMass_Eta_two_Both_15 = new TH1D("multi_InvMass_Eta_two_Both_15",
3548 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3549 "invariant mass in GeV/c^{2};#",
3550 200, 0.0, 2.0);
3552 multi_InvMass_Eta_two_Both_16 = new TH1D("multi_InvMass_Eta_two_Both_16",
3553 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3554 "invariant mass in GeV/c^{2};#",
3555 200, 0.0, 2.0);
3557 multi_InvMass_Eta_two_Both_17 = new TH1D("multi_InvMass_Eta_two_Both_17",
3558 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3559 "invariant mass in GeV/c^{2};#",
3560 200, 0.0, 2.0);
3562 multi_EMT_Eta_two_Both_1 = new TH1D("multi_EMT_Eta_two_Both_1",
3563 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3564 "invariant mass in GeV/c^{2};#",
3565 200, 0.0, 2.0);
3567 multi_EMT_Eta_two_Both_2 = new TH1D("multi_EMT_Eta_two_Both_2",
3568 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3569 "invariant mass in GeV/c^{2};#",
3570 200, 0.0, 2.0);
3572 multi_EMT_Eta_two_Both_3 = new TH1D("multi_EMT_Eta_two_Both_3",
3573 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3574 "invariant mass in GeV/c^{2};#",
3575 200, 0.0, 2.0);
3577 multi_EMT_Eta_two_Both_4 = new TH1D("multi_EMT_Eta_two_Both_4",
3578 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3579 "invariant mass in GeV/c^{2};#",
3580 200, 0.0, 2.0);
3582 multi_EMT_Eta_two_Both_5 = new TH1D("multi_EMT_Eta_two_Both_5",
3583 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3584 "invariant mass in GeV/c^{2};#",
3585 200, 0.0, 2.0);
3587 multi_EMT_Eta_two_Both_6 = new TH1D("multi_EMT_Eta_two_Both_6",
3588 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3589 "invariant mass in GeV/c^{2};#",
3590 200, 0.0, 2.0);
3592 multi_EMT_Eta_two_Both_7 = new TH1D("multi_EMT_Eta_two_Both_7",
3593 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3594 "invariant mass in GeV/c^{2};#",
3595 200, 0.0, 2.0);
3597 multi_EMT_Eta_two_Both_8 = new TH1D("multi_EMT_Eta_two_Both_8",
3598 "rapidity = (1.6-2.0) P_{t} = (1.2-2.0 GeV/c^{2}) ; "
3599 "invariant mass in GeV/c^{2};#",
3600 200, 0.0, 2.0);
3602 multi_EMT_Eta_two_Both_9 = new TH1D("multi_EMT_Eta_two_Both_9",
3603 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3604 "invariant mass in GeV/c^{2};#",
3605 200, 0.0, 2.0);
3607 multi_EMT_Eta_two_Both_10 = new TH1D("multi_EMT_Eta_two_Both_10",
3608 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3609 "invariant mass in GeV/c^{2};#",
3610 200, 0.0, 2.0);
3612 multi_EMT_Eta_two_Both_11 = new TH1D("multi_EMT_Eta_two_Both_11",
3613 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3614 "invariant mass in GeV/c^{2};#",
3615 200, 0.0, 2.0);
3617 multi_EMT_Eta_two_Both_12 = new TH1D("multi_EMT_Eta_two_Both_12",
3618 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ; "
3619 "invariant mass in GeV/c^{2};#",
3620 200, 0.0, 2.0);
3622 multi_EMT_Eta_two_Both_13 = new TH1D("multi_EMT_Eta_two_Both_13",
3623 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3624 "invariant mass in GeV/c^{2};#",
3625 200, 0.0, 2.0);
3627 multi_EMT_Eta_two_Both_14 = new TH1D("multi_EMT_Eta_two_Both_14",
3628 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3629 "invariant mass in GeV/c^{2};#",
3630 200, 0.0, 2.0);
3632 multi_EMT_Eta_two_Both_15 = new TH1D("multi_EMT_Eta_two_Both_15",
3633 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ; "
3634 "invariant mass in GeV/c^{2};#",
3635 200, 0.0, 2.0);
3637 multi_EMT_Eta_two_Both_16 = new TH1D("multi_EMT_Eta_two_Both_16",
3638 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ; "
3639 "invariant mass in GeV/c^{2};#",
3640 200, 0.0, 2.0);
3642 multi_EMT_Eta_two_Both_17 = new TH1D("multi_EMT_Eta_two_Both_17",
3643 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ; "
3644 "invariant mass in GeV/c^{2};#",
3645 200, 0.0, 2.0);
3648}
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.
int Int_t
int32_t GetStsTrackIndex() const
int32_t GetRichRingIndex() const
double GetZ() const
Definition CbmHit.h:74
static void SetKFParticleFromStsTrack(CbmStsTrack *track, KFParticle *particle, Int_t pdg=211, Bool_t firstPoint=kTRUE)
TClonesArray * fRichHits
Definition CbmKresEta.h:85
TH1D * multi_EMT_Eta_onetwo_Both_5
Definition CbmKresEta.h:553
std::vector< std::vector< int > > Gammas_MCIndex_all_Outside
Definition CbmKresEta.h:176
TH1D * multi_InvMass_Eta_onetwo_Both_10
Definition CbmKresEta.h:541
TH1D * multi_EMT_Eta_two_Both_5
Definition CbmKresEta.h:606
TH2D * WrongEta_pt_vs_rap_est_all_Outside
Definition CbmKresEta.h:419
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Target
Definition CbmKresEta.h:135
vector< TH1 * > fHistoList_Eta_all_Target
Definition CbmKresEta.h:359
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:251
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Both
Definition CbmKresEta.h:336
TH1D * PlaneAngles_first_fromEta_zeroInRich_Outside
Definition CbmKresEta.h:271
TH1D * EMT_Eta_InvMass_two_Both
Definition CbmKresEta.h:582
std::vector< std::vector< TVector3 > > Gammas_onetwo_Both
Definition CbmKresEta.h:188
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Target
Definition CbmKresEta.h:361
TH1D * multi_EMT_Eta_onetwo_Both_11
Definition CbmKresEta.h:559
TH1D * TrueEta_InvMass_after_cuts_two_Both
Definition CbmKresEta.h:578
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:237
vector< CbmRichRing * > VRichRing_minus_Outside
Definition CbmKresEta.h:100
TH1D * multi_EMT_Eta_onetwo_Both_15
Definition CbmKresEta.h:563
TH1D * multi_InvMass_Eta_two_Both_4
Definition CbmKresEta.h:588
TH1D * Particle_pull_X_Outside
Definition CbmKresEta.h:211
TH1D * multi_InvMass_Eta_all_Both_7
Definition CbmKresEta.h:485
void Mixing_Outside()
TH1D * multi_EMT_Eta_all_Both_9
Definition CbmKresEta.h:504
std::vector< int > EMT_eta_gg_Event_Target
Definition CbmKresEta.h:134
TH1D * multi_InvMass_Eta_two_Both_9
Definition CbmKresEta.h:593
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:333
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Target
Definition CbmKresEta.h:160
TH1D * multi_EMT_Eta_two_Both_2
Definition CbmKresEta.h:603
TH1D * multi_InvMass_Eta_two_Both_6
Definition CbmKresEta.h:590
TH1D * Particles_PDG_two_Target
Definition CbmKresEta.h:393
vector< TH1 * > fHistoList_rap_pt_Eta_all_Both
Definition CbmKresEta.h:478
std::vector< int > frefId
Definition CbmKresEta.h:129
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:313
void Mixing_Target()
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:282
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:223
TH1D * InvMass_GammaCombinations_all_Target
Definition CbmKresEta.h:362
TH1D * Particles_PDG_all_Target
Definition CbmKresEta.h:363
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Outside
Definition CbmKresEta.h:182
TClonesArray * fMcTracks
Definition CbmKresEta.h:78
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:326
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:221
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:256
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:286
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Target
Definition CbmKresEta.h:375
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:273
TH1D * Particles_PDG_all_Both
Definition CbmKresEta.h:467
TH1D * multi_EMT_Eta_all_Both_2
Definition CbmKresEta.h:497
std::vector< int > EMT_eta_gg_Event_Both
Definition CbmKresEta.h:144
TH1D * Particle_pull_px_Target
Definition CbmKresEta.h:205
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Target
Definition CbmKresEta.h:390
TH1D * multi_EMT_Eta_two_Both_9
Definition CbmKresEta.h:610
TH1D * PlaneAngles_first_fromEta_oneInRich_Outside
Definition CbmKresEta.h:285
TH2D * TrueEta_pt_vs_rap_est_all_Outside
Definition CbmKresEta.h:418
TH2D * TrueEta_pt_vs_rap_est_two_Both
Definition CbmKresEta.h:576
TH1D * multi_EMT_Eta_all_Both_15
Definition CbmKresEta.h:510
TH1D * Particles_PDG_two_Outside
Definition CbmKresEta.h:445
vector< CbmRichRing * > VRichRing_plus_Target
Definition CbmKresEta.h:125
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:318
TH1D * multi_EMT_Eta_onetwo_Both_9
Definition CbmKresEta.h:557
TH1D * WrongEta_InvMass_after_cuts_all_Both
Definition CbmKresEta.h:473
TH1D * Particle_pull_pz_Target
Definition CbmKresEta.h:207
vector< CbmStsTrack * > VStsTrack_minus_Target
Definition CbmKresEta.h:113
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Both
Definition CbmKresEta.h:196
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:229
TH1D * InvMass_GammaCombinations_onetwo_Outside
Definition CbmKresEta.h:429
vector< TH1 * > fHistoList_Eta_cuts_Outside
Definition CbmKresEta.h:263
TH1D * multi_InvMass_Eta_onetwo_Both_15
Definition CbmKresEta.h:546
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:296
TH2D * WrongEta_pt_vs_rap_two_Outside
Definition CbmKresEta.h:447
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Outside
Definition CbmKresEta.h:180
TClonesArray * fStsTracks
Definition CbmKresEta.h:80
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Both
Definition CbmKresEta.h:322
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Target
Definition CbmKresEta.h:360
std::vector< std::vector< int > > Gammas_stsIndex_all_Both
Definition CbmKresEta.h:190
TH1D * multi_EMT_Eta_two_Both_17
Definition CbmKresEta.h:618
TH2D * TrueEta_pt_vs_rap_est_two_Outside
Definition CbmKresEta.h:448
vector< TH1 * > fHistoList_rap_pt_Eta_two_Target
Definition CbmKresEta.h:406
TH1D * PlaneAngles_last_fromEta_afterCuts_oneInRich_Both
Definition CbmKresEta.h:337
std::vector< std::vector< int > > Gammas_stsIndex_two_Target
Definition CbmKresEta.h:155
TH1D * PlaneAngles_last_fromEta_zeroInRich_Target
Definition CbmKresEta.h:224
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:288
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Target
Definition CbmKresEta.h:163
TH1D * InvMass_GammaCombinations_all_Both
Definition CbmKresEta.h:466
TH2D * TrueEta_pt_vs_rap_est_onetwo_Target
Definition CbmKresEta.h:381
std::vector< int > VStsIndex_minus_Target
Definition CbmKresEta.h:116
TH1D * multi_InvMass_Eta_all_Both_5
Definition CbmKresEta.h:483
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:300
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Outside
Definition CbmKresEta.h:291
TH1D * multi_EMT_Eta_onetwo_Both_6
Definition CbmKresEta.h:554
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:274
TH1D * PlaneAngles_last_fromEta_twoInRich_Both
Definition CbmKresEta.h:347
TH1D * Particles_PDG_onetwo_Both
Definition CbmKresEta.h:520
TH1D * WrongEta_InvMass_after_cuts_two_Both
Definition CbmKresEta.h:579
TH1D * multi_EMT_Eta_onetwo_Both_10
Definition CbmKresEta.h:558
std::vector< int > VRings_minus_Outside
Definition CbmKresEta.h:98
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Both
Definition CbmKresEta.h:531
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Both
Definition CbmKresEta.h:192
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:249
TH1D * multi_EMT_Eta_all_Both_7
Definition CbmKresEta.h:502
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:266
TH1D * PlaneAngles_first_fromEta_oneInRich_Both
Definition CbmKresEta.h:330
vector< TH1 * > fHistoList_Eta_two_Outside
Definition CbmKresEta.h:441
std::vector< int > VStsIndex_plus_Outside
Definition CbmKresEta.h:106
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Outside
Definition CbmKresEta.h:423
TH1D * InvMass_GammaCombinations_two_Outside
Definition CbmKresEta.h:444
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Target
Definition CbmKresEta.h:391
vector< TH1 * > fHistoList_Eta_cuts_Target
Definition CbmKresEta.h:217
TH2D * WrongEta_pt_vs_rap_onetwo_Target
Definition CbmKresEta.h:380
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Outside
Definition CbmKresEta.h:442
TH2D * TrueEta_pt_vs_rap_est_two_Target
Definition CbmKresEta.h:396
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:289
TH1D * WrongEta_InvMass_after_cuts_onetwo_Target
Definition CbmKresEta.h:384
TH1D * multi_EMT_Eta_all_Both_17
Definition CbmKresEta.h:512
TH1D * multi_InvMass_Eta_all_Both_12
Definition CbmKresEta.h:490
TH1D * WrongEta_InvMass_after_cuts_all_Outside
Definition CbmKresEta.h:421
std::vector< std::vector< TVector3 > > Gammas_two_Outside
Definition CbmKresEta.h:169
TH1D * multi_EMT_Eta_two_Both_8
Definition CbmKresEta.h:609
std::vector< std::vector< int > > Gammas_stsIndex_two_Outside
Definition CbmKresEta.h:173
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:228
TH1D * multi_EMT_Eta_onetwo_Both_2
Definition CbmKresEta.h:550
std::vector< std::vector< TVector3 > > Gammas_onetwo_Outside
Definition CbmKresEta.h:170
TH2D * TrueEta_pt_vs_rap_all_Target
Definition CbmKresEta.h:364
TH1D * PlaneAngles_first_fromEta_zeroInRich_Both
Definition CbmKresEta.h:316
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:255
std::vector< int > fMCId
Definition CbmKresEta.h:130
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Both
Definition CbmKresEta.h:353
vector< TH1 * > fHistoList_Eta_onetwo_Target
Definition CbmKresEta.h:374
TH1D * multi_InvMass_Eta_two_Both_12
Definition CbmKresEta.h:596
TH1D * EMT_Eta_InvMass_onetwo_Both
Definition CbmKresEta.h:529
TH1D * multi_InvMass_Eta_two_Both_7
Definition CbmKresEta.h:591
TClonesArray * fArrayStsHit
Definition CbmKresEta.h:87
TH1D * TrueEta_InvMass_after_cuts_onetwo_Outside
Definition CbmKresEta.h:435
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:331
void Exec(int fEventNumEta, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
TH1D * multi_InvMass_Eta_two_Both_1
Definition CbmKresEta.h:585
vector< TH1 * > fHistoList_Eta_two_Target
Definition CbmKresEta.h:389
std::vector< int > VRings_plus_Outside
Definition CbmKresEta.h:105
TH1D * multi_EMT_Eta_all_Both_16
Definition CbmKresEta.h:511
TH2D * TrueEta_pt_vs_rap_onetwo_Target
Definition CbmKresEta.h:379
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Outside
Definition CbmKresEta.h:437
TH1D * PlaneAngles_last_fromEta_zeroInRich_Both
Definition CbmKresEta.h:315
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
Definition CbmKresEta.h:613
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Both
Definition CbmKresEta.h:199
TH1D * TrueEta_InvMass_after_cuts_all_Outside
Definition CbmKresEta.h:420
TH1D * InvMass_GammaCombinations_two_Target
Definition CbmKresEta.h:392
TH1D * multi_InvMass_Eta_two_Both_5
Definition CbmKresEta.h:589
TH1D * PlaneAngles_last_fromEta_zeroInRich_Outside
Definition CbmKresEta.h:270
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Outside
Definition CbmKresEta.h:452
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:351
TH1D * multi_EMT_Eta_two_Both_4
Definition CbmKresEta.h:605
double CalculatePlaneAngle_last(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
TH1D * multi_EMT_Eta_all_Both_3
Definition CbmKresEta.h:498
TH1D * multi_EMT_Eta_onetwo_Both_7
Definition CbmKresEta.h:555
int FindInRich(int richInd, int stsMcTrackId)
TH1D * multi_EMT_Eta_all_Both_14
Definition CbmKresEta.h:509
TH1D * multi_InvMass_Eta_onetwo_Both_12
Definition CbmKresEta.h:543
TH1D * PlaneAngles_first_fromEta_twoInRich_Target
Definition CbmKresEta.h:253
TH1D * multi_EMT_Eta_onetwo_Both_16
Definition CbmKresEta.h:564
TH2D * TrueEta_pt_vs_rap_two_Both
Definition CbmKresEta.h:574
TH2D * WrongEta_pt_vs_rap_all_Both
Definition CbmKresEta.h:469
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Target
Definition CbmKresEta.h:401
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:275
TH1D * multi_InvMass_Eta_onetwo_Both_7
Definition CbmKresEta.h:538
TH1D * multi_InvMass_Eta_all_Both_13
Definition CbmKresEta.h:491
void InitHistograms()
TH1D * multi_InvMass_Eta_all_Both_4
Definition CbmKresEta.h:482
vector< CbmMCTrack * > VMCtracks_plus_Outside
Definition CbmKresEta.h:104
TH1D * multi_EMT_Eta_two_Both_14
Definition CbmKresEta.h:615
TH1D * multi_InvMass_Eta_two_Both_8
Definition CbmKresEta.h:592
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:295
TH1D * multi_InvMass_Eta_onetwo_Both_4
Definition CbmKresEta.h:535
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Both
Definition CbmKresEta.h:198
TH1D * GammasMomentum_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:242
TH1D * multi_InvMass_Eta_two_Both_15
Definition CbmKresEta.h:599
CbmRichRingFitterEllipseTau * fTauFit
Definition CbmKresEta.h:92
TH1D * PlaneAngles_first_fromEta_zeroInRich_Target
Definition CbmKresEta.h:225
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Target
Definition CbmKresEta.h:230
TH2D * WrongEta_pt_vs_rap_all_Outside
Definition CbmKresEta.h:417
TH1D * multi_InvMass_Eta_all_Both_16
Definition CbmKresEta.h:494
TH1D * multi_InvMass_Eta_onetwo_Both_14
Definition CbmKresEta.h:545
vector< CbmMCTrack * > VMCtracks_minus_Outside
Definition CbmKresEta.h:97
std::vector< int > EMT_eta_gg_NofRings_Target
Definition CbmKresEta.h:136
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:352
vector< CbmStsTrack * > VStsTrack_minus_Outside
Definition CbmKresEta.h:96
std::vector< std::vector< TVector3 > > Gammas_onetwo_Target
Definition CbmKresEta.h:152
TH2D * WrongEta_pt_vs_rap_est_onetwo_Both
Definition CbmKresEta.h:524
vector< CbmRichRing * > VRichRing_plus_Outside
Definition CbmKresEta.h:107
TH1D * multi_EMT_Eta_all_Both_4
Definition CbmKresEta.h:499
vector< Int_t > VMCIndex_plus_Outside
Definition CbmKresEta.h:108
vector< CbmRichRing * > VRichRing_minus_Target
Definition CbmKresEta.h:117
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Outside
Definition CbmKresEta.h:279
TH2D * TrueEta_pt_vs_rap_all_Outside
Definition CbmKresEta.h:416
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Target
Definition CbmKresEta.h:233
TH1D * PlaneAngles_first_fromEta_oneInRich_Target
Definition CbmKresEta.h:239
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:311
std::vector< std::vector< int > > Gammas_MCIndex_two_Outside
Definition CbmKresEta.h:177
TH1D * multi_InvMass_Eta_two_Both_10
Definition CbmKresEta.h:594
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Target
Definition CbmKresEta.h:371
TH1D * multi_InvMass_Eta_two_Both_13
Definition CbmKresEta.h:597
TH1D * PlaneAngles_last_fromEta_twoInRich_Target
Definition CbmKresEta.h:252
TH1D * multi_EMT_Eta_two_Both_1
Definition CbmKresEta.h:602
TClonesArray * fRichRings
Definition CbmKresEta.h:83
TH1D * multi_InvMass_Eta_all_Both_8
Definition CbmKresEta.h:486
TH1D * Particle_pull_py_Outside
Definition CbmKresEta.h:209
std::vector< int > VRings_plus_Target
Definition CbmKresEta.h:123
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Both
Definition CbmKresEta.h:464
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Outside
Definition CbmKresEta.h:428
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Both
Definition CbmKresEta.h:517
TH1D * PlaneAngles_first_wrongpairs_oneInRich_Target
Definition CbmKresEta.h:245
TH1D * AngleBetweenGammas_trueEta_before_cuts_all_Outside
Definition CbmKresEta.h:412
vector< TVector3 > VMomenta_plus_Target
Definition CbmKresEta.h:122
TH1D * multi_InvMass_Eta_all_Both_10
Definition CbmKresEta.h:488
std::vector< std::vector< TVector3 > > Gammas_all_Both
Definition CbmKresEta.h:186
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:303
vector< TH1 * > fHistoList_Eta_cuts_oneInRich_Both
Definition CbmKresEta.h:324
TH2D * TrueEta_pt_vs_rap_est_onetwo_Both
Definition CbmKresEta.h:523
TH1D * multi_EMT_Eta_onetwo_Both_17
Definition CbmKresEta.h:565
vector< TH1 * > fHistoList_Eta_two_Both
Definition CbmKresEta.h:569
TH1D * PlaneAngles_last_fromEta_oneInRich_Target
Definition CbmKresEta.h:238
virtual ~CbmKresEta()
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Both
Definition CbmKresEta.h:354
TClonesArray * fArrayMvdHit
Definition CbmKresEta.h:86
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:317
std::vector< std::vector< TVector3 > > Gammas_all_Target
Definition CbmKresEta.h:150
TH1D * multi_InvMass_Eta_all_Both_9
Definition CbmKresEta.h:487
TH1D * Particle_pull_pz_Outside
Definition CbmKresEta.h:210
TH1D * GammasInvMass_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:220
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Both
Definition CbmKresEta.h:527
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:235
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Both
Definition CbmKresEta.h:474
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:345
std::vector< std::vector< int > > Gammas_MCIndex_all_Target
Definition CbmKresEta.h:158
TH2D * WrongEta_pt_vs_rap_est_all_Both
Definition CbmKresEta.h:471
TH1D * multi_InvMass_Eta_all_Both_14
Definition CbmKresEta.h:492
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Outside
Definition CbmKresEta.h:293
TH1D * WrongEta_InvMass_after_cuts_all_Target
Definition CbmKresEta.h:369
vector< TH1 * > fHistoList_rap_pt_Eta_all_Target
Definition CbmKresEta.h:404
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Outside
Definition CbmKresEta.h:276
TH1D * Particle_pull_Y_Outside
Definition CbmKresEta.h:212
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Both
Definition CbmKresEta.h:528
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:241
TH1D * EMT_Eta_InvMass_onetwo_Target
Definition CbmKresEta.h:387
TH1D * multi_InvMass_Eta_all_Both_2
Definition CbmKresEta.h:480
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Both
Definition CbmKresEta.h:580
TH1D * TrueEta_AngleBetweenGammas_after_cuts_two_Target
Definition CbmKresEta.h:400
TH1D * multi_EMT_Eta_all_Both_10
Definition CbmKresEta.h:505
TClonesArray * fRichProjections
Definition CbmKresEta.h:82
TH1D * GammasMomentum_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:319
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Outside
Definition CbmKresEta.h:443
TH1D * TrueEta_InvMass_after_cuts_two_Outside
Definition CbmKresEta.h:450
TH1D * TrueEta_AngleBetweenGammas_after_cuts_onetwo_Target
Definition CbmKresEta.h:385
TH1D * multi_InvMass_Eta_all_Both_11
Definition CbmKresEta.h:489
TH1D * EMT_Eta_InvMass_two_Target
Definition CbmKresEta.h:402
TH1D * multi_InvMass_Eta_two_Both_14
Definition CbmKresEta.h:598
TH2D * TrueEta_pt_vs_rap_onetwo_Outside
Definition CbmKresEta.h:431
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
Definition CbmKresEta.h:280
TH1D * multi_EMT_Eta_all_Both_6
Definition CbmKresEta.h:501
vector< CbmMCTrack * > VMCtracks_plus_Target
Definition CbmKresEta.h:120
TH1D * PlaneAngles_last_wrongpairs_afterCuts_oneInRich_Both
Definition CbmKresEta.h:339
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:226
TH1D * multi_EMT_Eta_onetwo_Both_12
Definition CbmKresEta.h:560
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:234
TClonesArray * fStsTrackMatches
Definition CbmKresEta.h:81
TH1D * Particles_PDG_onetwo_Target
Definition CbmKresEta.h:378
double CalculatePlaneAngle_first(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
std::vector< std::vector< int > > Gammas_MCIndex_all_Both
Definition CbmKresEta.h:194
TH1D * multi_EMT_Eta_two_Both_16
Definition CbmKresEta.h:617
TH2D * TrueEta_pt_vs_rap_est_all_Both
Definition CbmKresEta.h:470
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:301
vector< Int_t > VMCIndex_plus_Target
Definition CbmKresEta.h:126
TH2D * WrongEta_pt_vs_rap_est_onetwo_Outside
Definition CbmKresEta.h:434
TH1D * multi_InvMass_Eta_two_Both_11
Definition CbmKresEta.h:595
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:227
TH1D * InvMass_GammaCombinations_onetwo_Both
Definition CbmKresEta.h:519
TH1D * PlaneAngles_last_fromEta_oneInRich_Outside
Definition CbmKresEta.h:284
TH1D * Particle_pull_Z_Outside
Definition CbmKresEta.h:213
vector< TH1 * > fHistoList_Eta_all_Outside
Definition CbmKresEta.h:411
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:248
TH1D * multi_EMT_Eta_two_Both_6
Definition CbmKresEta.h:607
std::vector< int > EMT_eta_gg_NofRings_Both
Definition CbmKresEta.h:146
int CheckIfElectron(CbmRichRing *ring, double momentum)
CbmKFVertex fKFVertex
Definition CbmKresEta.h:90
TH1D * AngleBetweenGammas_trueEta_before_cuts_onetwo_Outside
Definition CbmKresEta.h:427
TH1D * multi_EMT_Eta_onetwo_Both_13
Definition CbmKresEta.h:561
TH1D * InvMass_GammaCombinations_onetwo_Target
Definition CbmKresEta.h:377
TH1D * multi_InvMass_Eta_onetwo_Both_11
Definition CbmKresEta.h:542
TH1D * multi_EMT_Eta_onetwo_Both_4
Definition CbmKresEta.h:552
TH1D * GammasRapidity_wrongpairs_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:320
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Outside
Definition CbmKresEta.h:277
TH1D * multi_InvMass_Eta_onetwo_Both_8
Definition CbmKresEta.h:539
TH2D * WrongEta_pt_vs_rap_two_Target
Definition CbmKresEta.h:395
TH1D * multi_InvMass_Eta_onetwo_Both_1
Definition CbmKresEta.h:532
vector< CbmMCTrack * > VMCtracks_minus_Target
Definition CbmKresEta.h:112
TH1D * EMT_Eta_InvMass_all_Outside
Definition CbmKresEta.h:424
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:254
TH2D * WrongEta_pt_vs_rap_est_two_Target
Definition CbmKresEta.h:397
TH1D * PlaneAngles_last_fromEta_oneInRich_Both
Definition CbmKresEta.h:329
TH2D * WrongEta_pt_vs_rap_est_two_Both
Definition CbmKresEta.h:577
TH1D * multi_EMT_Eta_all_Both_12
Definition CbmKresEta.h:507
vector< Int_t > VMCIndex_minus_Target
Definition CbmKresEta.h:118
vector< TH1 * > fHistoList_Eta_onetwo_Outside
Definition CbmKresEta.h:426
std::vector< std::vector< TVector3 > > Gammas_two_Both
Definition CbmKresEta.h:187
std::vector< int > EMT_eta_gg_Event_Outside
Definition CbmKresEta.h:139
TH1D * multi_EMT_Eta_all_Both_8
Definition CbmKresEta.h:503
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:268
TH1D * Particle_pull_px_Outside
Definition CbmKresEta.h:208
TH1D * EMT_Eta_InvMass_all_Target
Definition CbmKresEta.h:372
TH1D * WrongEta_InvMass_after_cuts_two_Outside
Definition CbmKresEta.h:451
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:287
TH1D * multi_InvMass_Eta_all_Both_17
Definition CbmKresEta.h:495
TH1D * multi_EMT_Eta_onetwo_Both_14
Definition CbmKresEta.h:562
CbmVertex * fPrimVertex
Definition CbmKresEta.h:89
TH1D * TrueEta_InvMass_after_cuts_onetwo_Both
Definition CbmKresEta.h:525
TH1D * GammasMomentum_fromEta_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:250
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Target
Definition CbmKresEta.h:258
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:328
vector< Int_t > VMCIndex_minus_Outside
Definition CbmKresEta.h:101
TH1D * Particles_PDG_onetwo_Outside
Definition CbmKresEta.h:430
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Target
Definition CbmKresEta.h:219
std::vector< std::vector< int > > Gammas_MCIndex_two_Target
Definition CbmKresEta.h:159
TH1D * multi_InvMass_Eta_onetwo_Both_3
Definition CbmKresEta.h:534
TH1D * GammasOpenAngle_fromEta_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:281
TH2D * WrongEta_pt_vs_rap_onetwo_Outside
Definition CbmKresEta.h:432
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:314
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Target
Definition CbmKresEta.h:244
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Target
Definition CbmKresEta.h:376
TH2D * TrueEta_pt_vs_rap_est_all_Target
Definition CbmKresEta.h:366
TH1D * multi_InvMass_Eta_all_Both_1
Definition CbmKresEta.h:479
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:297
vector< TH1 * > fHistoList_Eta_onetwo_Both
Definition CbmKresEta.h:516
TH1D * PlaneAngles_first_fromEta_afterCuts_oneInRich_Both
Definition CbmKresEta.h:338
TH1D * multi_EMT_Eta_all_Both_13
Definition CbmKresEta.h:508
TH1D * Particles_PDG_two_Both
Definition CbmKresEta.h:573
TH1D * PlaneAngles_first_fromEta_twoInRich_Outside
Definition CbmKresEta.h:299
TH1D * PlaneAngles_first_wrongpairs_afterCuts_oneInRich_Both
Definition CbmKresEta.h:340
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:243
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Outside
Definition CbmKresEta.h:140
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Target
Definition CbmKresEta.h:370
vector< TH1 * > fHistoList_rap_pt_Eta_all_Outside
Definition CbmKresEta.h:456
vector< TVector3 > VMomenta_minus_Target
Definition CbmKresEta.h:114
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:343
TH1D * GammasRapidity_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:269
vector< CbmStsTrack * > VStsTrack_plus_Target
Definition CbmKresEta.h:121
TH1D * GammasRapidity_wrongpairs_beforeCuts_twoInRich_Target
Definition CbmKresEta.h:257
TH1D * PlaneAngles_last_wrongpairs_twoInRich_Outside
Definition CbmKresEta.h:304
std::vector< CbmMCTrack * > fMCtracks
Definition CbmKresEta.h:131
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Outside
Definition CbmKresEta.h:181
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Both
Definition CbmKresEta.h:310
TH1D * Particle_pull_py_Target
Definition CbmKresEta.h:206
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Both
Definition CbmKresEta.h:200
TH1D * WrongEta_InvMass_after_cuts_two_Target
Definition CbmKresEta.h:399
TH1D * GammasInvMass_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:325
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Target
Definition CbmKresEta.h:259
TH1D * multi_InvMass_Eta_all_Both_3
Definition CbmKresEta.h:481
TH2D * TrueEta_pt_vs_rap_onetwo_Both
Definition CbmKresEta.h:521
vector< TH1 * > fHistoList_rap_pt_Eta_two_Both
Definition CbmKresEta.h:584
vector< CbmStsTrack * > VStsTrack_plus_Outside
Definition CbmKresEta.h:103
std::vector< int > EMT_eta_gg_NofRings_Outside
Definition CbmKresEta.h:141
TH1D * InvMass_GammaCombinations_two_Both
Definition CbmKresEta.h:572
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Both
Definition CbmKresEta.h:581
TH1D * multi_InvMass_Eta_all_Both_15
Definition CbmKresEta.h:493
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Outside
Definition CbmKresEta.h:438
std::vector< std::vector< TVector3 > > EMT_eta_gg_pair_momenta_Both
Definition CbmKresEta.h:145
TH1D * InvMass_GammaCombinations_all_Outside
Definition CbmKresEta.h:414
TH1D * multi_EMT_Eta_two_Both_7
Definition CbmKresEta.h:608
vector< TH1 * > fHistoList_Eta_cuts_Both
Definition CbmKresEta.h:308
TH1D * multi_InvMass_Eta_onetwo_Both_5
Definition CbmKresEta.h:536
std::vector< std::vector< int > > Gammas_stsIndex_all_Outside
Definition CbmKresEta.h:172
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:332
TH1D * multi_InvMass_Eta_all_Both_6
Definition CbmKresEta.h:484
void SaveOutsideTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
TH1D * TrueEta_InvMass_after_cuts_onetwo_Target
Definition CbmKresEta.h:383
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Target
Definition CbmKresEta.h:247
TH1D * TrueEta_AngleBetweenGammas_after_cuts_all_Outside
Definition CbmKresEta.h:422
TH1D * multi_InvMass_Eta_onetwo_Both_17
Definition CbmKresEta.h:548
TH1D * multi_EMT_Eta_onetwo_Both_1
Definition CbmKresEta.h:549
TH1D * WrongEta_InvMass_after_cuts_onetwo_Outside
Definition CbmKresEta.h:436
TH1D * TrueEta_InvMass_after_cuts_all_Target
Definition CbmKresEta.h:368
TH1D * WrongEta_InvMass_after_cuts_onetwo_Both
Definition CbmKresEta.h:526
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Both
Definition CbmKresEta.h:465
TH1D * multi_EMT_Eta_all_Both_5
Definition CbmKresEta.h:500
vector< TH1 * > fHistoList_Eta_all_Both
Definition CbmKresEta.h:463
TH1D * multi_EMT_Eta_onetwo_Both_3
Definition CbmKresEta.h:551
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Outside
Definition CbmKresEta.h:174
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Target
Definition CbmKresEta.h:405
TH1D * GammasRapidity_wrongpairs_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:334
TH1D * multi_InvMass_Eta_onetwo_Both_16
Definition CbmKresEta.h:547
TH1D * WrongEta_AngleBetweenGammas_after_cuts_two_Outside
Definition CbmKresEta.h:453
std::vector< std::vector< int > > Gammas_stsIndex_two_Both
Definition CbmKresEta.h:191
TH1D * AngleBetweenGammas_wrongEta_before_cuts_two_Both
Definition CbmKresEta.h:571
std::vector< std::vector< int > > Gammas_MCIndex_two_Both
Definition CbmKresEta.h:195
TH1D * GammasInvMass_wrongpairs_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:272
TH1D * multi_EMT_Eta_two_Both_3
Definition CbmKresEta.h:604
TH1D * Particles_PDG_all_Outside
Definition CbmKresEta.h:415
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Both
Definition CbmKresEta.h:312
TH2D * WrongEta_pt_vs_rap_est_two_Outside
Definition CbmKresEta.h:449
TH1D * GammasOpenAngle_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:350
TH1D * GammasRapidity_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:346
std::vector< std::vector< TVector3 > > Gammas_all_Outside
Definition CbmKresEta.h:168
TH2D * TrueEta_pt_vs_rap_two_Outside
Definition CbmKresEta.h:446
TH1D * PlaneAngles_first_wrongpairs_zeroInRich_Target
Definition CbmKresEta.h:231
TH2D * WrongEta_pt_vs_rap_onetwo_Both
Definition CbmKresEta.h:522
vector< TVector3 > frefmomenta
Definition CbmKresEta.h:128
TH1D * multi_InvMass_Eta_onetwo_Both_9
Definition CbmKresEta.h:540
vector< TH1 * > fHistoList_rap_pt_Eta_two_Outside
Definition CbmKresEta.h:458
std::vector< int > VRings_minus_Target
Definition CbmKresEta.h:115
std::vector< std::vector< TVector3 > > Gammas_two_Target
Definition CbmKresEta.h:151
vector< TH1 * > fHistoList_Eta_cuts_zeroInRich_Outside
Definition CbmKresEta.h:265
TH1D * multi_EMT_Eta_two_Both_10
Definition CbmKresEta.h:611
TH1D * EMT_Eta_InvMass_onetwo_Outside
Definition CbmKresEta.h:439
TH1D * WrongEta_AngleBetweenGammas_after_cuts_all_Both
Definition CbmKresEta.h:475
TH1D * multi_InvMass_Eta_onetwo_Both_2
Definition CbmKresEta.h:533
TH2D * WrongEta_pt_vs_rap_two_Both
Definition CbmKresEta.h:575
TH1D * multi_EMT_Eta_two_Both_15
Definition CbmKresEta.h:616
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
Definition CbmKresEta.h:154
TH1D * PlaneAngles_first_fromEta_twoInRich_Both
Definition CbmKresEta.h:348
TH1D * GammasMomentum_wrongpairs_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:302
TH1D * GammasRapidity_fromEta_beforeCuts_oneInRich_Outside
Definition CbmKresEta.h:283
TH1D * multi_EMT_Eta_all_Both_1
Definition CbmKresEta.h:496
void SaveTargetTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, TVector3 refmom, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
std::vector< int > VStsIndex_minus_Outside
Definition CbmKresEta.h:99
TH1D * PlaneAngles_first_wrongpairs_twoInRich_Outside
Definition CbmKresEta.h:305
TH1D * EMT_Eta_InvMass_two_Outside
Definition CbmKresEta.h:454
TH2D * TrueEta_pt_vs_rap_est_onetwo_Outside
Definition CbmKresEta.h:433
TH2D * TrueEta_pt_vs_rap_two_Target
Definition CbmKresEta.h:394
void FindGammasBoth()
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Target
Definition CbmKresEta.h:156
TH2D * WrongEta_pt_vs_rap_all_Target
Definition CbmKresEta.h:365
TH1D * AngleBetweenGammas_wrongEta_before_cuts_all_Outside
Definition CbmKresEta.h:413
vector< TH1 * > fHistoList_Eta_cuts_twoInRich_Both
Definition CbmKresEta.h:342
std::vector< std::vector< int > > Gammas_MCIndex_onetwo_Outside
Definition CbmKresEta.h:178
TH1D * multi_InvMass_Eta_onetwo_Both_13
Definition CbmKresEta.h:544
std::vector< int > VStsIndex_plus_Target
Definition CbmKresEta.h:124
vector< TH1 * > fHistoList_Eta
Definition CbmKresEta.h:204
TH2D * TrueEta_pt_vs_rap_all_Both
Definition CbmKresEta.h:468
TH1D * PlaneAngles_last_wrongpairs_zeroInRich_Both
Definition CbmKresEta.h:321
TClonesArray * fGlobalTracks
Definition CbmKresEta.h:79
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Outside
Definition CbmKresEta.h:290
TH1D * multi_InvMass_Eta_two_Both_16
Definition CbmKresEta.h:600
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Target
Definition CbmKresEta.h:164
TH1D * TrueEta_InvMass_after_cuts_all_Both
Definition CbmKresEta.h:472
TH1D * GammasInvMass_wrongpairs_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:240
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Target
Definition CbmKresEta.h:162
TH1D * PlaneAngles_last_wrongpairs_oneInRich_Both
Definition CbmKresEta.h:335
vector< TH1 * > fHistoList_rap_pt_Eta_onetwo_Outside
Definition CbmKresEta.h:457
TH1D * GammasMomentum_fromEta_beforeCuts_zeroInRich_Target
Definition CbmKresEta.h:222
TH1D * GammasOpenAngle_fromEta_beforeCuts_zeroInRich_Outside
Definition CbmKresEta.h:267
TH1D * GammasOpenAngle_fromEta_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:344
TH2D * WrongEta_pt_vs_rap_est_all_Target
Definition CbmKresEta.h:367
TClonesArray * fRichRingMatches
Definition CbmKresEta.h:84
TH1D * multi_EMT_Eta_two_Both_13
Definition CbmKresEta.h:614
TH1D * multi_InvMass_Eta_two_Both_17
Definition CbmKresEta.h:601
TH1D * multi_InvMass_Eta_two_Both_2
Definition CbmKresEta.h:586
TH1D * WrongEta_AngleBetweenGammas_after_cuts_onetwo_Target
Definition CbmKresEta.h:386
TH1D * AngleBetweenGammas_trueEta_before_cuts_two_Both
Definition CbmKresEta.h:570
TH1D * TrueEta_InvMass_after_cuts_two_Target
Definition CbmKresEta.h:398
TH1D * EMT_Eta_InvMass_all_Both
Definition CbmKresEta.h:476
TH1D * GammasInvMass_fromEta_beforeCuts_twoInRich_Outside
Definition CbmKresEta.h:294
TH1D * multi_EMT_Eta_two_Both_11
Definition CbmKresEta.h:612
TH2D * WrongEta_pt_vs_rap_est_onetwo_Target
Definition CbmKresEta.h:382
TH1D * PlaneAngles_last_fromEta_twoInRich_Outside
Definition CbmKresEta.h:298
TH1D * multi_EMT_Eta_onetwo_Both_8
Definition CbmKresEta.h:556
TH1D * multi_InvMass_Eta_two_Both_3
Definition CbmKresEta.h:587
TH1D * GammasInvMass_wrongpairs_beforeCuts_twoInRich_Both
Definition CbmKresEta.h:349
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Target
Definition CbmKresEta.h:236
TH1D * GammasMomentum_fromEta_beforeCuts_oneInRich_Both
Definition CbmKresEta.h:327
TH1D * AngleBetweenGammas_wrongEta_before_cuts_onetwo_Both
Definition CbmKresEta.h:518
TH1D * multi_EMT_Eta_all_Both_11
Definition CbmKresEta.h:506
TH1D * multi_InvMass_Eta_onetwo_Both_6
Definition CbmKresEta.h:537
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)
double GetPz() const
Definition CbmMCTrack.h:71
double GetPx() const
Definition CbmMCTrack.h:69
double GetStartZ() const
Definition CbmMCTrack.h:74
int32_t GetMotherId() const
Definition CbmMCTrack.h:68
double GetStartX() const
Definition CbmMCTrack.h:72
int32_t GetPdgCode() const
Definition CbmMCTrack.h:67
double GetStartY() const
Definition CbmMCTrack.h:73
double GetPy() const
Definition CbmMCTrack.h:70
int32_t GetNofLinks() const
Definition CbmMatch.h:42
const CbmLink & GetMatchedLink() const
Definition CbmMatch.h:41
double GetY() const
Definition CbmPixelHit.h:74
double GetX() const
Definition CbmPixelHit.h:73
Here the ring is fitted with Taubin algorithm from A. Ayriyan, G. Ososkov, N. Chernov.
float GetAaxis() const
float GetBaxis() const
void AddHit(CbmRichHitLight hit)
Add new hit to the ring.
uint32_t GetHit(int32_t i) const
Definition CbmRichRing.h:39
int32_t GetNofHits() const
Definition CbmRichRing.h:37
data class for a reconstructed 3-d hit in the STS
Definition CbmStsHit.h:35
int32_t GetNofMvdHits() const
Definition CbmStsTrack.h:87
int32_t GetMvdHitIndex(int32_t iHit) const
Definition CbmStsTrack.h:75
int32_t GetStsHitIndex(int32_t iHit) const
int32_t GetNofStsHits() const
Definition CbmStsTrack.h:93
const FairTrackParam * GetParamFirst() const
Definition CbmTrack.h:68
Double_t fMinv
Definition LmvmKinePar.h:20
Double_t fPt
Definition LmvmKinePar.h:18
Double_t fMomentumMag
Definition LmvmKinePar.h:17
Double_t fRapidity
Definition LmvmKinePar.h:19
Hash for CbmL1LinkKey.