CbmRoot
Loading...
Searching...
No Matches
CbmKresConversionManualmbias2.cxx
Go to the documentation of this file.
1/* Copyright (C) 2017-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
2 SPDX-License-Identifier: GPL-3.0-only
3 Authors: Ievgenii Kres, Florian Uhlig [committer] */
4
22
24
25#include "CbmGlobalTrack.h"
27#include "CbmKresConversionBG.h"
28#include "CbmKresFunctions.h"
29#include "CbmMCTrack.h"
30#include "CbmMvdHit.h"
31#include "CbmRichHit.h"
32#include "CbmRichRing.h"
34#include "CbmRichRingLight.h"
35#include "CbmStsHit.h"
36#include "CbmStsTrack.h"
37#include "CbmTrackMatchNew.h"
38#include "FairRootManager.h"
39#include "KFParticle/KFParticle.h"
40#include "TDirectory.h"
41
42#include <iostream>
43
44
45using namespace std;
46
48 : fMcTracks(nullptr)
49 , fGlobalTracks(nullptr)
50 , fStsTracks(nullptr)
51 , fStsTrackMatches(nullptr)
52 , fRichProjections(nullptr)
53 , fRichRings(nullptr)
54 , fRichRingMatches(nullptr)
55 , fRichHits(nullptr)
56 , fArrayMvdHit(nullptr)
57 , fArrayStsHit(nullptr)
58 , fArrayCentrality(nullptr)
59 , fPrimVertex(nullptr)
60 , fKFVertex()
61 , fTauFit(nullptr)
62 , fAnaBG(nullptr)
85 , frefmomenta()
86 , frefId()
87 , fMCtracks()
178 , Pdg_all_Target_mbias2(nullptr)
179 , P_reco_all_Target_mbias2(nullptr)
188 , Pdg_zero_Target_mbias2(nullptr)
198 , Pdg_one_Target_mbias2(nullptr)
199 , P_reco_one_Target_mbias2(nullptr)
208 , Pdg_two_Target_mbias2(nullptr)
209 , P_reco_two_Target_mbias2(nullptr)
218 , Pdg_onetwo_Target_mbias2(nullptr)
228 , Pdg_all_Outside_mbias2(nullptr)
238 , Pdg_zero_Outside_mbias2(nullptr)
248 , Pdg_one_Outside_mbias2(nullptr)
258 , Pdg_two_Outside_mbias2(nullptr)
278 , Pdg_all_Both_mbias2(nullptr)
279 , P_reco_all_Both_mbias2(nullptr)
280 , Pt_reco_all_Both_mbias2(nullptr)
288 , Pdg_zero_Both_mbias2(nullptr)
289 , P_reco_zero_Both_mbias2(nullptr)
290 , Pt_reco_zero_Both_mbias2(nullptr)
298 , Pdg_one_Both_mbias2(nullptr)
299 , P_reco_one_Both_mbias2(nullptr)
300 , Pt_reco_one_Both_mbias2(nullptr)
308 , Pdg_two_Both_mbias2(nullptr)
309 , P_reco_two_Both_mbias2(nullptr)
310 , Pt_reco_two_Both_mbias2(nullptr)
318 , Pdg_onetwo_Both_mbias2(nullptr)
326 , Pdg_vs_Distance_mbias2(nullptr)
327 , P_vs_Distance_mbias2(nullptr)
362 , rap_vs_Pt_InM_1_mbias2(nullptr)
363 , rap_vs_Pt_InM_2_mbias2(nullptr)
364 , rap_vs_Pt_InM_3_mbias2(nullptr)
365 , rap_vs_Pt_InM_4_mbias2(nullptr)
366 , rap_vs_Pt_InM_5_mbias2(nullptr)
367 , rap_vs_Pt_InM_6_mbias2(nullptr)
368 , rap_vs_Pt_InM_7_mbias2(nullptr)
369 , rap_vs_Pt_InM_8_mbias2(nullptr)
370 , rap_vs_Pt_InM_9_mbias2(nullptr)
371 , rap_vs_Pt_InM_10_mbias2(nullptr)
372 , rap_vs_Pt_InM_11_mbias2(nullptr)
373 , rap_vs_Pt_InM_12_mbias2(nullptr)
374 , rap_vs_Pt_InM_13_mbias2(nullptr)
375 , rap_vs_Pt_InM_14_mbias2(nullptr)
376 , rap_vs_Pt_InM_15_mbias2(nullptr)
377 , rap_vs_Pt_InM_16_mbias2(nullptr)
378 , rap_vs_Pt_InM_17_mbias2(nullptr)
379 , rap_vs_Pt_InM_18_mbias2(nullptr)
380 , rap_vs_Pt_InM_19_mbias2(nullptr)
381 , rap_vs_Pt_InM_20_mbias2(nullptr)
382 , rap_vs_Pt_InM_21_mbias2(nullptr)
383 , rap_vs_Pt_InM_22_mbias2(nullptr)
384 , rap_vs_Pt_InM_23_mbias2(nullptr)
385 , rap_vs_Pt_InM_24_mbias2(nullptr)
386 , rap_vs_Pt_InM_25_mbias2(nullptr)
387 , rap_vs_Pt_InM_26_mbias2(nullptr)
388 , rap_vs_Pt_InM_27_mbias2(nullptr)
389 , rap_vs_Pt_InM_28_mbias2(nullptr)
390 , rap_vs_Pt_InM_29_mbias2(nullptr)
391 , rap_vs_Pt_InM_30_mbias2(nullptr)
433 , BG1_InM_all_Both_mbias2(nullptr)
434 , BG2_InM_all_Both_mbias2(nullptr)
435 , BG3_InM_all_Both_mbias2(nullptr)
436 , BG4_InM_all_Both_mbias2(nullptr)
437 , BG5_InM_all_Both_mbias2(nullptr)
438 , BG6_InM_all_Both_mbias2(nullptr)
439 , BG7_InM_all_Both_mbias2(nullptr)
440 , BG8_InM_all_Both_mbias2(nullptr)
441 , BG9_InM_all_Both_mbias2(nullptr)
442 , BG10_InM_all_Both_mbias2(nullptr)
464 , BG1_InM_zero_Both_mbias2(nullptr)
465 , BG2_InM_zero_Both_mbias2(nullptr)
466 , BG3_InM_zero_Both_mbias2(nullptr)
467 , BG4_InM_zero_Both_mbias2(nullptr)
468 , BG5_InM_zero_Both_mbias2(nullptr)
469 , BG6_InM_zero_Both_mbias2(nullptr)
470 , BG7_InM_zero_Both_mbias2(nullptr)
471 , BG8_InM_zero_Both_mbias2(nullptr)
472 , BG9_InM_zero_Both_mbias2(nullptr)
494 , BG1_InM_one_Both_mbias2(nullptr)
495 , BG2_InM_one_Both_mbias2(nullptr)
496 , BG3_InM_one_Both_mbias2(nullptr)
497 , BG4_InM_one_Both_mbias2(nullptr)
498 , BG5_InM_one_Both_mbias2(nullptr)
499 , BG6_InM_one_Both_mbias2(nullptr)
500 , BG7_InM_one_Both_mbias2(nullptr)
501 , BG8_InM_one_Both_mbias2(nullptr)
502 , BG9_InM_one_Both_mbias2(nullptr)
503 , BG10_InM_one_Both_mbias2(nullptr)
525 , BG1_InM_two_Both_mbias2(nullptr)
526 , BG2_InM_two_Both_mbias2(nullptr)
527 , BG3_InM_two_Both_mbias2(nullptr)
528 , BG4_InM_two_Both_mbias2(nullptr)
529 , BG5_InM_two_Both_mbias2(nullptr)
530 , BG6_InM_two_Both_mbias2(nullptr)
531 , BG7_InM_two_Both_mbias2(nullptr)
532 , BG8_InM_two_Both_mbias2(nullptr)
533 , BG9_InM_two_Both_mbias2(nullptr)
534 , BG10_InM_two_Both_mbias2(nullptr)
586{
587}
588
590
592{
593 FairRootManager* ioman = FairRootManager::Instance();
594 if (nullptr == ioman) { Fatal("CbmKresConversionManualmbias2::Init", "RootManager not instantised!"); }
595
596 fMcTracks = (TClonesArray*) ioman->GetObject("MCTrack");
597 if (nullptr == fMcTracks) { Fatal("CbmKresConversionManualmbias2::Init", "No MCTrack array!"); }
598
599 fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex."));
600 if (nullptr == fPrimVertex) { fPrimVertex = dynamic_cast<CbmVertex*>(ioman->GetObject("PrimaryVertex")); }
601 if (nullptr == fPrimVertex) { LOG(fatal) << "CbmKresConversionManualmbias2::Init No PrimaryVertex array!"; }
602
603 fGlobalTracks = (TClonesArray*) ioman->GetObject("GlobalTrack");
604 if (nullptr == fGlobalTracks) { Fatal("CbmKresConversionManualmbias2::Init", "No GlobalTrack array!"); }
605
606 fStsTracks = (TClonesArray*) ioman->GetObject("StsTrack");
607 if (nullptr == fStsTracks) { Fatal("CbmKresConversionManualmbias2::Init", "No StsTrack array!"); }
608
609 fStsTrackMatches = (TClonesArray*) ioman->GetObject("StsTrackMatch");
610 if (nullptr == fStsTrackMatches) { Fatal("CbmKresConversionManualmbias2::Init", "No StsTrackMatch array!"); }
611
612 fRichProjections = (TClonesArray*) ioman->GetObject("RichProjection");
613 if (nullptr == fRichProjections) { Fatal("CbmKresConversionManualmbias2::Init", "No RichProjection array!"); }
614
615 fRichRings = (TClonesArray*) ioman->GetObject("RichRing");
616 if (nullptr == fRichRings) { Fatal("CbmKresConversionManualmbias2::Init", "No RichRing array!"); }
617
618 fRichRingMatches = (TClonesArray*) ioman->GetObject("RichRingMatch");
619 if (nullptr == fRichRingMatches) { Fatal("CbmKresConversionManualmbias2::Init", "No RichRingMatch array!"); }
620
621 fRichHits = (TClonesArray*) ioman->GetObject("RichHit");
622 if (nullptr == fRichHits) { Fatal("CbmKresConversionManualmbias2::Init", "No RichHit array!"); }
623
624 fArrayMvdHit = (TClonesArray*) ioman->GetObject("MvdHit");
625 if (nullptr == fArrayMvdHit) { Fatal("CbmKresConversionManualmbias2::Init", "No MvdHit array!"); }
626
627 fArrayStsHit = (TClonesArray*) ioman->GetObject("StsHit");
628 if (nullptr == fArrayStsHit) { Fatal("CbmKresConversionManualmbias2::Init", "No StsHit array!"); }
629
630 fArrayCentrality = (FairMCEventHeader*) ioman->GetObject("MCEventHeader.");
631 if (nullptr == fArrayCentrality) { Fatal("CbmAnaConversion2::Init", "No fArrayCentrality array!"); }
632
633
635
637
639 fAnaBG->Init();
640}
641
642void CbmKresConversionManualmbias2::Exec(int fEventNumMan, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
643{
644 cout << "CbmKresConversionManualmbias2, event No. " << fEventNumMan << endl;
645
646 double centrality_mbias2 = fArrayCentrality->GetB();
647
648 if (fEventNumMan == 1000) {
663 Mixing_Both();
667
672
677 }
678 if (centrality_mbias2 <= 2 || centrality_mbias2 > 6) return;
679
680
681 if (fPrimVertex != nullptr) { fKFVertex = CbmKFVertex(*fPrimVertex); }
682 else {
683 return;
684 }
685
686
732
733
744
757
758
759 Int_t ngTracks = fGlobalTracks->GetEntriesFast();
760 for (Int_t i = 0; i < ngTracks; i++) {
761 CbmGlobalTrack* gTrack = (CbmGlobalTrack*) fGlobalTracks->At(i);
762 if (nullptr == gTrack) continue;
763 Int_t stsInd = gTrack->GetStsTrackIndex();
764 Int_t richInd = gTrack->GetRichRingIndex();
765
766 if (stsInd < 0) continue;
767 CbmStsTrack* stsTrack = (CbmStsTrack*) fStsTracks->At(stsInd);
768 if (stsTrack == nullptr) continue;
769 CbmTrackMatchNew* stsMatch = (CbmTrackMatchNew*) fStsTrackMatches->At(stsInd);
770 if (stsMatch == nullptr) continue;
771 if (stsMatch->GetNofLinks() <= 0) continue;
772 Int_t stsMcTrackId = stsMatch->GetMatchedLink().GetIndex();
773 if (stsMcTrackId < 0) continue;
774 CbmMCTrack* mcTrack = (CbmMCTrack*) fMcTracks->At(stsMcTrackId);
775 if (mcTrack == nullptr) continue;
776
777
778 FairTrackParam* proj = (FairTrackParam*) fRichProjections->At(i);
779 if (richInd < 0 && proj->GetX() > -115 && proj->GetX() < 115
780 && ((proj->GetY() < -120 && proj->GetY() > -200) || (proj->GetY() > 120 && proj->GetY() < 200)))
781 continue;
782 CbmRichRing* Ring = nullptr;
783 if (richInd > -1) {
784 Ring = static_cast<CbmRichRing*>(fRichRings->At(richInd));
785 //Pdg_vs_Distance_mbias2->Fill(TMath::Abs(mcTrack->GetPdgCode()), Ring->GetDistance());
786 //if (TMath::Abs(mcTrack->GetPdgCode()) == 11) P_vs_Distance_mbias2->Fill(mcTrack->GetP(), Ring->GetDistance());
787 }
788
789
790 // Doing fit with Fit To primary Vertex and calculate chi2 to primary vertex
791 double chi2 = 0;
792 TVector3 Momentum = CbmKresFunctions::FitToVertexAndGetChi(stsTrack, fKFVertex.GetRefX(), fKFVertex.GetRefY(),
793 fKFVertex.GetRefZ(), chi2);
794 const FairTrackParam* track_par = stsTrack->GetParamFirst();
795 double charge = track_par->GetQp();
796
797
798 if (chi2 != chi2) continue;
799 if (chi2 == 0) continue;
800
801 if (chi2 > 3) { SaveOutsideTracks(mcTrack, stsTrack, charge, stsInd, richInd, stsMcTrackId, Ring); }
802 if (chi2 > 3) continue;
803
804 SaveTargetTracks(mcTrack, stsTrack, Momentum, charge, stsInd, richInd, stsMcTrackId, Ring);
805 }
806
807
808 FindGammasTarget(fEventNumMan, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Target_mbias2,
813
814 FindGammasOutside(fEventNumMan, OpeningAngleCut, GammaInvMassCut, RealPID, VMCtracks_minus_Outside_mbias2,
818
820
821
842
863
884}
885
886
888 int stsInd, int richInd, int stsMcTrackId, CbmRichRing* RING)
889{
890 int InRich = FindInRich(richInd, stsMcTrackId);
891 if (charge < 0) {
892 VMCtracks_minus_Outside_mbias2.push_back(mcTrack1);
893 VStsTrack_minus_Outside_mbias2.push_back(stsTrack);
894 VRings_minus_Outside_mbias2.push_back(InRich);
895 VStsIndex_minus_Outside_mbias2.push_back(stsInd);
896 VRichRing_minus_Outside_mbias2.push_back(RING);
897 }
898 if (charge > 0) {
899 VMCtracks_plus_Outside_mbias2.push_back(mcTrack1);
900 VStsTrack_plus_Outside_mbias2.push_back(stsTrack);
901 VRings_plus_Outside_mbias2.push_back(InRich);
902 VStsIndex_plus_Outside_mbias2.push_back(stsInd);
903 VRichRing_plus_Outside_mbias2.push_back(RING);
904 }
905}
906
908 double charge, int stsInd, int richInd, int stsMcTrackId,
909 CbmRichRing* RING)
910{
911 int InRich = FindInRich(richInd, stsMcTrackId);
912 if (charge < 0) {
913 VMCtracks_minus_Target_mbias2.push_back(mcTrack1);
914 VStsTrack_minus_Target_mbias2.push_back(stsTrack);
915 VMomenta_minus_Target_mbias2.push_back(refmom);
916 VRings_minus_Target_mbias2.push_back(InRich);
917 VStsIndex_minus_Target_mbias2.push_back(stsInd);
918 VRichRing_minus_Target_mbias2.push_back(RING);
919 }
920 if (charge > 0) {
921 VMCtracks_plus_Target_mbias2.push_back(mcTrack1);
922 VStsTrack_plus_Target_mbias2.push_back(stsTrack);
923 VMomenta_plus_Target_mbias2.push_back(refmom);
924 VRings_plus_Target_mbias2.push_back(InRich);
925 VStsIndex_plus_Target_mbias2.push_back(stsInd);
926 VRichRing_plus_Target_mbias2.push_back(RING);
927 }
928}
929
930
932 int EventNumMan, double AngleCut, double InvMassCut, int RealPID, vector<CbmMCTrack*> MCtracks_minus,
933 vector<CbmMCTrack*> MCtracks_plus, vector<CbmStsTrack*> StsTrack_minus, vector<CbmStsTrack*> StsTrack_plus,
934 vector<TVector3> Momenta_minus, vector<TVector3> Momenta_plus, std::vector<int> Rings_minus,
935 std::vector<int> Rings_plus, std::vector<int> stsIndex_minus, std::vector<int> stsIndex_plus,
936 vector<CbmRichRing*> richRing_minus, vector<CbmRichRing*> richRing_plus)
937{
938 for (size_t i = 0; i < Momenta_minus.size(); i++) {
939 for (size_t j = 0; j < Momenta_plus.size(); j++) {
940
941 CbmStsTrack* sts1 = StsTrack_minus.at(i);
942 CbmStsTrack* sts2 = StsTrack_plus.at(j);
943
944 TVector3 part1 = Momenta_minus[i];
945 TVector3 part2 = Momenta_plus[j];
946 CbmMCTrack* part1MC = MCtracks_minus[i];
947 CbmMCTrack* part2MC = MCtracks_plus[j];
948
949 int richcheck_0 = 0;
950 int richcheck_1 = 0;
951 if (RealPID == 1) {
952 // Real PID
953 richcheck_0 = CheckIfElectron(richRing_minus[i], Momenta_minus[i].Mag());
954 richcheck_1 = CheckIfElectron(richRing_plus[j], Momenta_plus[j].Mag());
955 }
956 else {
957 // MC PID
958 richcheck_0 = Rings_minus[i];
959 richcheck_1 = Rings_plus[j];
960 }
961 int richcheck = richcheck_0 + richcheck_1;
962
963 Double_t InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
964 Double_t OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
965
966 InvMass_vs_OA_candidates_Target_mbias2->Fill(InvmassReco, OpeningAngle);
967 InvMass_vs_OA_candidates_Both_mbias2->Fill(InvmassReco, OpeningAngle);
969 GammasInvMass_candidates_Both_mbias2->Fill(InvmassReco);
970 GammasOA_candidates_Target_mbias2->Fill(OpeningAngle);
971 GammasOA_candidates_Both_mbias2->Fill(OpeningAngle);
972
973 double PlaneAngle_last = CalculatePlaneAngle_last(sts1, sts2);
974 double PlaneAngle_first = CalculatePlaneAngle_first(sts1, sts2);
975
976 //if (PlaneAngle_last > 20) continue;
977
978 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
979 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
980 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
981 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
982 if (mcTrackgrmama->GetPdgCode() == 111) {
983 GammasInvMass_fromPi0_Target_mbias2->Fill(InvmassReco);
984 GammasInvMass_fromPi0_Both_mbias2->Fill(InvmassReco);
985 GammasOA_fromPi0_Target_mbias2->Fill(OpeningAngle);
986 GammasOA_fromPi0_Both_mbias2->Fill(OpeningAngle);
987 InvMass_vs_OA_fromPi0_Target_mbias2->Fill(InvmassReco, OpeningAngle);
988 InvMass_vs_OA_fromPi0_Both_mbias2->Fill(InvmassReco, OpeningAngle);
989 PlaneAngles_last_fromPi0_Target_mbias2->Fill(PlaneAngle_last);
990 PlaneAngles_last_fromPi0_Both_mbias2->Fill(PlaneAngle_last);
991 PlaneAngles_first_fromPi0_Target_mbias2->Fill(PlaneAngle_first);
992 PlaneAngles_first_fromPi0_Both_mbias2->Fill(PlaneAngle_first);
993 }
994 }
995 }
996
997 PlaneAngles_last_candidates_Target_mbias2->Fill(PlaneAngle_last);
998 PlaneAngles_last_candidates_Both_mbias2->Fill(PlaneAngle_last);
999 PlaneAngles_first_candidates_Target_mbias2->Fill(PlaneAngle_first);
1000 PlaneAngles_first_candidates_Both_mbias2->Fill(PlaneAngle_first);
1001
1002 // cuts
1003 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1004 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1005
1006 frefmomenta.clear();
1007 frefmomenta.push_back(part1);
1008 frefmomenta.push_back(part2);
1009 frefId.clear();
1010 frefId.push_back(stsIndex_minus[i]);
1011 frefId.push_back(stsIndex_plus[j]);
1012 fMCtracks.clear();
1013 fMCtracks.push_back(part1MC);
1014 fMCtracks.push_back(part2MC);
1015
1016 // for event mixing
1017 EMT_man_Event_Target_mbias2.push_back(EventNumMan);
1019 EMT_man_NofRings_Target_mbias2.push_back(richcheck);
1020
1021
1022 // everything (RICH == 0, RICH == 1, RICH == 2) together
1023 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1027 GammaInvMassReco_all_Target_mbias2->Fill(InvmassReco);
1028 GammaOpeningAngleReco_all_Target_mbias2->Fill(OpeningAngle);
1029 Pdg_all_Target_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1030 Pdg_all_Target_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1031 P_reco_all_Target_mbias2->Fill(part1.Mag());
1032 P_reco_all_Target_mbias2->Fill(part2.Mag());
1033 Pt_reco_all_Target_mbias2->Fill(part1.Perp());
1034 Pt_reco_all_Target_mbias2->Fill(part2.Perp());
1035 GammaInvMassReco_all_Both_mbias2->Fill(InvmassReco);
1036 GammaOpeningAngleReco_all_Both_mbias2->Fill(OpeningAngle);
1037 Pdg_all_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1038 Pdg_all_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1039 P_reco_all_Both_mbias2->Fill(part1.Mag());
1040 P_reco_all_Both_mbias2->Fill(part2.Mag());
1041 Pt_reco_all_Both_mbias2->Fill(part1.Perp());
1042 Pt_reco_all_Both_mbias2->Fill(part2.Perp());
1043 }
1044
1045 // only cases, when RICH == 0
1046 if (richcheck == 0) {
1050 GammaInvMassReco_zero_Target_mbias2->Fill(InvmassReco);
1051 GammaOpeningAngleReco_zero_Target_mbias2->Fill(OpeningAngle);
1052 Pdg_zero_Target_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1053 Pdg_zero_Target_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1054 P_reco_zero_Target_mbias2->Fill(part1.Mag());
1055 P_reco_zero_Target_mbias2->Fill(part2.Mag());
1056 Pt_reco_zero_Target_mbias2->Fill(part1.Perp());
1057 Pt_reco_zero_Target_mbias2->Fill(part2.Perp());
1058 GammaInvMassReco_zero_Both_mbias2->Fill(InvmassReco);
1059 GammaOpeningAngleReco_zero_Both_mbias2->Fill(OpeningAngle);
1060 Pdg_zero_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1061 Pdg_zero_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1062 P_reco_zero_Both_mbias2->Fill(part1.Mag());
1063 P_reco_zero_Both_mbias2->Fill(part2.Mag());
1064 Pt_reco_zero_Both_mbias2->Fill(part1.Perp());
1065 Pt_reco_zero_Both_mbias2->Fill(part2.Perp());
1066 }
1067
1068 // only cases, when RICH == 1
1069 if (richcheck == 1) {
1073 GammaInvMassReco_one_Target_mbias2->Fill(InvmassReco);
1074 GammaOpeningAngleReco_one_Target_mbias2->Fill(OpeningAngle);
1075 Pdg_one_Target_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1076 Pdg_one_Target_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1077 P_reco_one_Target_mbias2->Fill(part1.Mag());
1078 P_reco_one_Target_mbias2->Fill(part2.Mag());
1079 Pt_reco_one_Target_mbias2->Fill(part1.Perp());
1080 Pt_reco_one_Target_mbias2->Fill(part2.Perp());
1081 GammaInvMassReco_one_Both_mbias2->Fill(InvmassReco);
1082 GammaOpeningAngleReco_one_Both_mbias2->Fill(OpeningAngle);
1083 Pdg_one_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1084 Pdg_one_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1085 P_reco_one_Both_mbias2->Fill(part1.Mag());
1086 P_reco_one_Both_mbias2->Fill(part2.Mag());
1087 Pt_reco_one_Both_mbias2->Fill(part1.Perp());
1088 Pt_reco_one_Both_mbias2->Fill(part2.Perp());
1089 }
1090
1091 // only cases, when RICH == 2
1092 if (richcheck == 2) {
1096 GammaInvMassReco_two_Target_mbias2->Fill(InvmassReco);
1097 GammaOpeningAngleReco_two_Target_mbias2->Fill(OpeningAngle);
1098 Pdg_two_Target_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1099 Pdg_two_Target_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1100 P_reco_two_Target_mbias2->Fill(part1.Mag());
1101 P_reco_two_Target_mbias2->Fill(part2.Mag());
1102 Pt_reco_two_Target_mbias2->Fill(part1.Perp());
1103 Pt_reco_two_Target_mbias2->Fill(part2.Perp());
1104 GammaInvMassReco_two_Both_mbias2->Fill(InvmassReco);
1105 GammaOpeningAngleReco_two_Both_mbias2->Fill(OpeningAngle);
1106 Pdg_two_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1107 Pdg_two_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1108 P_reco_two_Both_mbias2->Fill(part1.Mag());
1109 P_reco_two_Both_mbias2->Fill(part2.Mag());
1110 Pt_reco_two_Both_mbias2->Fill(part1.Perp());
1111 Pt_reco_two_Both_mbias2->Fill(part2.Perp());
1112 }
1113
1114 // cases, when RICH == 1 or RICH == 2 together
1115 if (richcheck == 1 || richcheck == 2) {
1119 GammaInvMassReco_onetwo_Target_mbias2->Fill(InvmassReco);
1121 Pdg_onetwo_Target_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1122 Pdg_onetwo_Target_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1123 P_reco_onetwo_Target_mbias2->Fill(part1.Mag());
1124 P_reco_onetwo_Target_mbias2->Fill(part2.Mag());
1125 Pt_reco_onetwo_Target_mbias2->Fill(part1.Perp());
1126 Pt_reco_onetwo_Target_mbias2->Fill(part2.Perp());
1127 GammaInvMassReco_onetwo_Both_mbias2->Fill(InvmassReco);
1128 GammaOpeningAngleReco_onetwo_Both_mbias2->Fill(OpeningAngle);
1129 Pdg_onetwo_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1130 Pdg_onetwo_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1131 P_reco_onetwo_Both_mbias2->Fill(part1.Mag());
1132 P_reco_onetwo_Both_mbias2->Fill(part2.Mag());
1133 Pt_reco_onetwo_Both_mbias2->Fill(part1.Perp());
1134 Pt_reco_onetwo_Both_mbias2->Fill(part2.Perp());
1135 }
1136 }
1137 }
1138
1139 //cout << "number of gammas Target with 0-2 electron identified in RICH = " << Gammas_all_Target_mbias2.size() << endl;
1140 //cout << "number of gammas Target with 0 electron identified in RICH = " << Gammas_zero_Target_mbias2.size() << endl;
1141 //cout << "number of gammas Target with 1 electron identified in RICH = " << Gammas_one_Target_mbias2.size() << endl;
1142 //cout << "number of gammas Target with 2 electron identified in RICH = " << Gammas_two_Target_mbias2.size() << endl;
1143 //cout << "number of gammas Target with 1-2 electron identified in RICH = " << Gammas_onetwo_Target_mbias2.size() << endl;
1144}
1145
1146
1148 int EventNumMan, double AngleCut, double InvMassCut, int RealPID, vector<CbmMCTrack*> MCtracks_minus_Outside_mbias2,
1149 vector<CbmMCTrack*> MCtracks_plus_Outside_mbias2, vector<CbmStsTrack*> StsTrack_minus_Outside_mbias2,
1150 vector<CbmStsTrack*> StsTrack_plus_Outside_mbias2, std::vector<int> Rings_minus_Outside_mbias2,
1151 std::vector<int> Rings_plus_Outside_mbias2, std::vector<int> stsIndex_minus_Outside_mbias2,
1152 std::vector<int> stsIndex_plus_Outside_mbias2, vector<CbmRichRing*> richRing_minus_Outside_mbias2,
1153 vector<CbmRichRing*> richRing_plus_Outside_mbias2)
1154{
1155 for (size_t i = 0; i < StsTrack_minus_Outside_mbias2.size(); i++) {
1156 for (size_t j = 0; j < StsTrack_plus_Outside_mbias2.size(); j++) {
1157
1158 CbmStsTrack* part1STS = StsTrack_minus_Outside_mbias2[i];
1159 CbmStsTrack* part2STS = StsTrack_plus_Outside_mbias2[j];
1160 CbmMCTrack* part1MC = MCtracks_minus_Outside_mbias2[i];
1161 CbmMCTrack* part2MC = MCtracks_plus_Outside_mbias2[j];
1162
1163 KFParticle electron;
1165 KFParticle positron;
1166 CbmKFParticleInterface::SetKFParticleFromStsTrack(part2STS, &positron, -11);
1167 const KFParticle* daughters[2] = {&electron, &positron};
1168 KFParticle intersection;
1169 intersection.Construct(daughters, 2);
1170
1171 if (intersection.GetZ() > 75 || intersection.GetZ() < -5) continue; // kick weird intersections
1172
1173
1174 TVector3 part1 =
1175 CbmKresFunctions::FitToVertex(part1STS, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1176 TVector3 part2 =
1177 CbmKresFunctions::FitToVertex(part2STS, intersection.GetX(), intersection.GetY(), intersection.GetZ());
1178
1179
1180 //cout << "=================" << endl;
1181 //cout << "FitToVertex part1: px = " << part1.X() << "; \t py = " << part1.Y() << "; \t pz = " << part1.Z() << endl;
1182 //cout << "FitToVertex part2: px = " << part2.X() << "; \t py = " << part2.Y() << "; \t pz = " << part2.Z() << endl;
1183 //cout << "KFParticle part1: px = " << electron.KFParticleBase::Px() << "; \t py = " << electron.KFParticleBase::Py() << "; \t pz = " << electron.KFParticleBase::Pz() << endl;
1184 //cout << "KFParticle part2: px = " << positron.KFParticleBase::Px() << "; \t py = " << positron.KFParticleBase::Py() << "; \t pz = " << positron.KFParticleBase::Pz() << endl;
1185 //TVector3 part1(electron.KFParticleBase::Px(), electron.KFParticleBase::Py(), electron.KFParticleBase::Pz());
1186 //TVector3 part2(positron.KFParticleBase::Px(), positron.KFParticleBase::Py(), positron.KFParticleBase::Pz());
1187
1188
1189 int richcheck = 0;
1190 int richcheck_0 = 0;
1191 int richcheck_1 = 0;
1192 if (RealPID == 1) {
1193 // Real PID
1194 richcheck_0 = CheckIfElectron(richRing_minus_Outside_mbias2[i], part1.Mag());
1195 richcheck_1 = CheckIfElectron(richRing_plus_Outside_mbias2[j], part2.Mag());
1196 richcheck = richcheck_0 + richcheck_1;
1197 }
1198 else {
1199 // MC PID
1200 richcheck = Rings_minus_Outside_mbias2[i] + Rings_plus_Outside_mbias2[j];
1201 }
1202
1203
1204 Double_t InvmassReco = CbmKresFunctions::Invmass_2particles_RECO(part1, part2);
1205 Double_t OpeningAngle = CbmKresFunctions::CalculateOpeningAngle_Reco(part1, part2);
1206
1207 InvMass_vs_OA_candidates_Outside_mbias2->Fill(InvmassReco, OpeningAngle);
1208 InvMass_vs_OA_candidates_Both_mbias2->Fill(InvmassReco, OpeningAngle);
1209 GammasInvMass_candidates_Outside_mbias2->Fill(InvmassReco);
1210 GammasInvMass_candidates_Both_mbias2->Fill(InvmassReco);
1211 GammasOA_candidates_Outside_mbias2->Fill(OpeningAngle);
1212 GammasOA_candidates_Both_mbias2->Fill(OpeningAngle);
1213
1214 double PlaneAngle_last = CalculatePlaneAngle_last(part1STS, part2STS);
1215 double PlaneAngle_first = CalculatePlaneAngle_first(part1STS, part2STS);
1216
1217 // if (PlaneAngle_last > 20) continue;
1218
1219 if (part1MC->GetMotherId() == part2MC->GetMotherId() && part1MC->GetMotherId() != -1) {
1220 CbmMCTrack* mcTrackmama = (CbmMCTrack*) fMcTracks->At(part1MC->GetMotherId());
1221 if (mcTrackmama->GetMotherId() != -1 && mcTrackmama->GetPdgCode() == 22) {
1222 CbmMCTrack* mcTrackgrmama = (CbmMCTrack*) fMcTracks->At(mcTrackmama->GetMotherId());
1223 if (mcTrackgrmama->GetPdgCode() == 111) {
1224 GammasInvMass_fromPi0_Outside_mbias2->Fill(InvmassReco);
1225 GammasInvMass_fromPi0_Both_mbias2->Fill(InvmassReco);
1226 GammasOA_fromPi0_Outside_mbias2->Fill(OpeningAngle);
1227 GammasOA_fromPi0_Both_mbias2->Fill(OpeningAngle);
1228 InvMass_vs_OA_fromPi0_Outside_mbias2->Fill(InvmassReco, OpeningAngle);
1229 InvMass_vs_OA_fromPi0_Both_mbias2->Fill(InvmassReco, OpeningAngle);
1230 PlaneAngles_last_fromPi0_Outside_mbias2->Fill(PlaneAngle_last);
1231 PlaneAngles_last_fromPi0_Both_mbias2->Fill(PlaneAngle_last);
1232 PlaneAngles_first_fromPi0_Outside_mbias2->Fill(PlaneAngle_first);
1233 PlaneAngles_first_fromPi0_Both_mbias2->Fill(PlaneAngle_first);
1234 }
1235 }
1236 }
1237
1238 PlaneAngles_last_candidates_Outside_mbias2->Fill(PlaneAngle_last);
1239 PlaneAngles_last_candidates_Both_mbias2->Fill(PlaneAngle_last);
1240 PlaneAngles_first_candidates_Outside_mbias2->Fill(PlaneAngle_first);
1241 PlaneAngles_first_candidates_Both_mbias2->Fill(PlaneAngle_first);
1242
1243 // cuts
1244 if (TMath::Abs(OpeningAngle) > AngleCut) continue;
1245 if (TMath::Abs(InvmassReco) > InvMassCut) continue;
1246
1247 frefmomenta.clear();
1248 frefmomenta.push_back(part1);
1249 frefmomenta.push_back(part2);
1250
1251 frefId.clear();
1252 frefId.push_back(stsIndex_minus_Outside_mbias2[i]);
1253 frefId.push_back(stsIndex_plus_Outside_mbias2[j]);
1254 fMCtracks.clear();
1255 fMCtracks.push_back(part1MC);
1256 fMCtracks.push_back(part2MC);
1257
1258
1259 // for event mixing
1260 EMT_man_Event_Outside_mbias2.push_back(EventNumMan);
1262 EMT_man_NofRings_Outside_mbias2.push_back(richcheck);
1263
1264
1265 // everything (RICH == 0, RICH == 1, RICH == 2) together
1266 if (richcheck == 0 || richcheck == 1 || richcheck == 2) {
1270 GammaInvMassReco_all_Outside_mbias2->Fill(InvmassReco);
1271 GammaOpeningAngleReco_all_Outside_mbias2->Fill(OpeningAngle);
1272 Pdg_all_Outside_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1273 Pdg_all_Outside_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1274 P_reco_all_Outside_mbias2->Fill(part1.Mag());
1275 P_reco_all_Outside_mbias2->Fill(part2.Mag());
1276 Pt_reco_all_Outside_mbias2->Fill(part1.Perp());
1277 Pt_reco_all_Outside_mbias2->Fill(part2.Perp());
1278 GammaInvMassReco_all_Both_mbias2->Fill(InvmassReco);
1279 GammaOpeningAngleReco_all_Both_mbias2->Fill(OpeningAngle);
1280 Pdg_all_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1281 Pdg_all_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1282 P_reco_all_Both_mbias2->Fill(part1.Mag());
1283 P_reco_all_Both_mbias2->Fill(part2.Mag());
1284 Pt_reco_all_Both_mbias2->Fill(part1.Perp());
1285 Pt_reco_all_Both_mbias2->Fill(part2.Perp());
1286 }
1287
1288 // only cases, when RICH == 0
1289 if (richcheck == 0) {
1293 GammaInvMassReco_zero_Outside_mbias2->Fill(InvmassReco);
1295 Pdg_zero_Outside_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1296 Pdg_zero_Outside_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1297 P_reco_zero_Outside_mbias2->Fill(part1.Mag());
1298 P_reco_zero_Outside_mbias2->Fill(part2.Mag());
1299 Pt_reco_zero_Outside_mbias2->Fill(part1.Perp());
1300 Pt_reco_zero_Outside_mbias2->Fill(part2.Perp());
1301 GammaInvMassReco_zero_Both_mbias2->Fill(InvmassReco);
1302 GammaOpeningAngleReco_zero_Both_mbias2->Fill(OpeningAngle);
1303 Pdg_zero_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1304 Pdg_zero_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1305 P_reco_zero_Both_mbias2->Fill(part1.Mag());
1306 P_reco_zero_Both_mbias2->Fill(part2.Mag());
1307 Pt_reco_zero_Both_mbias2->Fill(part1.Perp());
1308 Pt_reco_zero_Both_mbias2->Fill(part2.Perp());
1309 }
1310
1311 // only cases, when RICH == 1
1312 if (richcheck == 1) {
1316 GammaInvMassReco_one_Outside_mbias2->Fill(InvmassReco);
1317 GammaOpeningAngleReco_one_Outside_mbias2->Fill(OpeningAngle);
1318 Pdg_one_Outside_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1319 Pdg_one_Outside_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1320 P_reco_one_Outside_mbias2->Fill(part1.Mag());
1321 P_reco_one_Outside_mbias2->Fill(part2.Mag());
1322 Pt_reco_one_Outside_mbias2->Fill(part1.Perp());
1323 Pt_reco_one_Outside_mbias2->Fill(part2.Perp());
1324 GammaInvMassReco_one_Both_mbias2->Fill(InvmassReco);
1325 GammaOpeningAngleReco_one_Both_mbias2->Fill(OpeningAngle);
1326 Pdg_one_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1327 Pdg_one_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1328 P_reco_one_Both_mbias2->Fill(part1.Mag());
1329 P_reco_one_Both_mbias2->Fill(part2.Mag());
1330 Pt_reco_one_Both_mbias2->Fill(part1.Perp());
1331 Pt_reco_one_Both_mbias2->Fill(part2.Perp());
1332 }
1333
1334 // only cases, when RICH == 2
1335 if (richcheck == 2) {
1339 GammaInvMassReco_two_Outside_mbias2->Fill(InvmassReco);
1340 GammaOpeningAngleReco_two_Outside_mbias2->Fill(OpeningAngle);
1341 Pdg_two_Outside_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1342 Pdg_two_Outside_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1343 P_reco_two_Outside_mbias2->Fill(part1.Mag());
1344 P_reco_two_Outside_mbias2->Fill(part2.Mag());
1345 Pt_reco_two_Outside_mbias2->Fill(part1.Perp());
1346 Pt_reco_two_Outside_mbias2->Fill(part2.Perp());
1347 GammaInvMassReco_two_Both_mbias2->Fill(InvmassReco);
1348 GammaOpeningAngleReco_two_Both_mbias2->Fill(OpeningAngle);
1349 Pdg_two_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1350 Pdg_two_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1351 P_reco_two_Both_mbias2->Fill(part1.Mag());
1352 P_reco_two_Both_mbias2->Fill(part2.Mag());
1353 Pt_reco_two_Both_mbias2->Fill(part1.Perp());
1354 Pt_reco_two_Both_mbias2->Fill(part2.Perp());
1355 }
1356
1357 // cases, when RICH == 1 or RICH == 2 together
1358 if (richcheck == 1 || richcheck == 2) {
1362 GammaInvMassReco_onetwo_Outside_mbias2->Fill(InvmassReco);
1364 Pdg_onetwo_Outside_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1365 Pdg_onetwo_Outside_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1366 P_reco_onetwo_Outside_mbias2->Fill(part1.Mag());
1367 P_reco_onetwo_Outside_mbias2->Fill(part2.Mag());
1368 Pt_reco_onetwo_Outside_mbias2->Fill(part1.Perp());
1369 Pt_reco_onetwo_Outside_mbias2->Fill(part2.Perp());
1370 GammaInvMassReco_onetwo_Both_mbias2->Fill(InvmassReco);
1371 GammaOpeningAngleReco_onetwo_Both_mbias2->Fill(OpeningAngle);
1372 Pdg_onetwo_Both_mbias2->Fill(TMath::Abs(part1MC->GetPdgCode()));
1373 Pdg_onetwo_Both_mbias2->Fill(TMath::Abs(part2MC->GetPdgCode()));
1374 P_reco_onetwo_Both_mbias2->Fill(part1.Mag());
1375 P_reco_onetwo_Both_mbias2->Fill(part2.Mag());
1376 Pt_reco_onetwo_Both_mbias2->Fill(part1.Perp());
1377 Pt_reco_onetwo_Both_mbias2->Fill(part2.Perp());
1378 }
1379 }
1380 }
1381
1382 //cout << "number of gammas Outside with 0-2 electron identified in RICH = " << Gammas_all_Outside_mbias2.size() << endl;
1383 //cout << "number of gammas Outside with 0 electron identified in RICH = " << Gammas_zero_Outside_mbias2.size() << endl;
1384 //cout << "number of gammas Outside with 1 electron identified in RICH = " << Gammas_one_Outside_mbias2.size() << endl;
1385 //cout << "number of gammas Outside with 2 electron identified in RICH = " << Gammas_two_Outside_mbias2.size() << endl;
1386 //cout << "number of gammas Outside with 1-2 electron identified in RICH = " << Gammas_onetwo_Outside_mbias2.size() << endl;
1387}
1388
1389
1391{
1412
1443
1464
1465 //cout << "number of gammas Both with 0-2 electron identified in RICH = " << Gammas_all_Both_mbias2.size() << endl;
1466 //cout << "number of gammas Both with 0 electron identified in RICH = " << Gammas_zero_Both_mbias2.size() << endl;
1467 //cout << "number of gammas Both with 1 electron identified in RICH = " << Gammas_one_Both_mbias2.size() << endl;
1468 //cout << "number of gammas Both with 2 electron identified in RICH = " << Gammas_two_Both_mbias2.size() << endl;
1469 //cout << "number of gammas Both with 1-2 electron identified in RICH = " << Gammas_onetwo_Both_mbias2.size() << endl;
1470}
1471
1472
1473void CbmKresConversionManualmbias2::FindPi0(TString mod, TString position, vector<vector<TVector3>> Gammas,
1474 vector<vector<int>> StsIndex, vector<vector<CbmMCTrack*>> GammasMC,
1475 TH1D* Pi0InvMassReco, TH2D* Pi0_pt_vs_rap, TH2D* Pi0_pt_vs_rap_est,
1476 TH2D* MultiplicityGamma, TH2D* MultiplicityChargedParticles,
1477 vector<TH1*> BGCases)
1478{
1479 // combine all gamma in pi0 --> calculate inv mass for gammas and pi0 // not the case, when one particle used twice for different gammas
1480 if (Gammas.size() < 2) return; // min 2 gammas to form pi0 are required
1481 for (size_t gamma1 = 0; gamma1 < Gammas.size() - 1; gamma1++) {
1482 for (size_t gamma2 = gamma1 + 1; gamma2 < Gammas.size(); gamma2++) {
1483 // 4 reconstructed particles from gammas
1484 TVector3 e11 = Gammas[gamma1][0];
1485 TVector3 e12 = Gammas[gamma1][1];
1486 TVector3 e21 = Gammas[gamma2][0];
1487 TVector3 e22 = Gammas[gamma2][1];
1488
1489 // MC true data for this particles
1490 CbmMCTrack* mcTrack1 = GammasMC[gamma1][0];
1491 CbmMCTrack* mcTrack2 = GammasMC[gamma1][1];
1492 CbmMCTrack* mcTrack3 = GammasMC[gamma2][0];
1493 CbmMCTrack* mcTrack4 = GammasMC[gamma2][1];
1494
1495 if (StsIndex[gamma1][0] == StsIndex[gamma2][0] || StsIndex[gamma1][0] == StsIndex[gamma2][1]
1496 || StsIndex[gamma1][1] == StsIndex[gamma2][0] || StsIndex[gamma1][1] == StsIndex[gamma2][1])
1497 continue; // particles not used twice --> different
1498
1500
1501 //Double_t OA1 = CbmKresFunctions::CalculateOpeningAngle_Reco(e11, e12);
1502 //Double_t OA2 = CbmKresFunctions::CalculateOpeningAngle_Reco(e21, e22);
1503 //if (params.fRapidity > 2.8 && params.fRapidity <= 3.6 && params.fPt > 0.0 && params.fPt <= 0.4 && OA1 > 1 && OA2 > 1) continue;
1504
1505 Pi0InvMassReco->Fill(params.fMinv);
1506 MultiplicityGamma->Fill(Gammas.size(), params.fMinv);
1507 MultiplicityChargedParticles->Fill(fGlobalTracks->GetEntriesFast(), params.fMinv);
1508
1509 // separate by rap and Pt only for onetwo
1510 if (mod == "OneTwo" && position == "Both") {
1511 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1512 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_1_mbias2->Fill(params.fMinv);
1513 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_2_mbias2->Fill(params.fMinv);
1514 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_3_mbias2->Fill(params.fMinv);
1515 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_4_mbias2->Fill(params.fMinv);
1516 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_5_mbias2->Fill(params.fMinv);
1517 }
1518 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1519 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_6_mbias2->Fill(params.fMinv);
1520 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_7_mbias2->Fill(params.fMinv);
1521 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_8_mbias2->Fill(params.fMinv);
1522 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_9_mbias2->Fill(params.fMinv);
1523 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_10_mbias2->Fill(params.fMinv);
1524 }
1525 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1526 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_11_mbias2->Fill(params.fMinv);
1527 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_12_mbias2->Fill(params.fMinv);
1528 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_13_mbias2->Fill(params.fMinv);
1529 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_14_mbias2->Fill(params.fMinv);
1530 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_15_mbias2->Fill(params.fMinv);
1531 }
1532 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1533 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_16_mbias2->Fill(params.fMinv);
1534 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_17_mbias2->Fill(params.fMinv);
1535 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_18_mbias2->Fill(params.fMinv);
1536 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_19_mbias2->Fill(params.fMinv);
1537 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_20_mbias2->Fill(params.fMinv);
1538 }
1539 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1540 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_21_mbias2->Fill(params.fMinv);
1541 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_22_mbias2->Fill(params.fMinv);
1542 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_23_mbias2->Fill(params.fMinv);
1543 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_24_mbias2->Fill(params.fMinv);
1544 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_25_mbias2->Fill(params.fMinv);
1545 }
1546 if (params.fRapidity > 3.2 && params.fRapidity <= 3.6) {
1547 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_26_mbias2->Fill(params.fMinv);
1548 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_27_mbias2->Fill(params.fMinv);
1549 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_28_mbias2->Fill(params.fMinv);
1550 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_29_mbias2->Fill(params.fMinv);
1551 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_30_mbias2->Fill(params.fMinv);
1552 }
1553 }
1554
1555
1556 if (position == "Both") fAnaBG->Exec(mcTrack1, mcTrack2, mcTrack3, mcTrack4, params.fMinv, BGCases);
1557
1558
1559 // fill histos with rapidity and Pt for correctly reconstructed pi0(dalitz) and pi0(g+g)
1560 if (nullptr == mcTrack1 || nullptr == mcTrack2 || nullptr == mcTrack3 || nullptr == mcTrack4) continue;
1561 if (TMath::Abs(mcTrack1->GetPdgCode()) != 11 || TMath::Abs(mcTrack2->GetPdgCode()) != 11
1562 || TMath::Abs(mcTrack3->GetPdgCode()) != 11 || TMath::Abs(mcTrack4->GetPdgCode()) != 11)
1563 continue;
1564 if (mcTrack1->GetPdgCode() + mcTrack2->GetPdgCode() != 0) continue;
1565 if (mcTrack3->GetPdgCode() + mcTrack4->GetPdgCode() != 0) continue;
1566 int motherId1 = mcTrack1->GetMotherId();
1567 int motherId2 = mcTrack2->GetMotherId();
1568 int motherId3 = mcTrack3->GetMotherId();
1569 int motherId4 = mcTrack4->GetMotherId();
1570 if (motherId1 == -1 || motherId2 == -1 || motherId3 == -1 || motherId4 == -1) continue;
1571 if (motherId1 != motherId2 || motherId3 != motherId4) continue;
1572 CbmMCTrack* mother1 = (CbmMCTrack*) fMcTracks->At(motherId1);
1573 CbmMCTrack* mother2 = (CbmMCTrack*) fMcTracks->At(motherId2);
1574 CbmMCTrack* mother3 = (CbmMCTrack*) fMcTracks->At(motherId3);
1575 CbmMCTrack* mother4 = (CbmMCTrack*) fMcTracks->At(motherId4);
1576 if (nullptr == mother1 || nullptr == mother2 || nullptr == mother3 || nullptr == mother4) continue;
1577 int mcMotherPdg1 = mother1->GetPdgCode();
1578 int mcMotherPdg2 = mother2->GetPdgCode();
1579 int mcMotherPdg3 = mother3->GetPdgCode();
1580 int mcMotherPdg4 = mother4->GetPdgCode();
1581 int grandmotherId1 = mother1->GetMotherId();
1582 int grandmotherId2 = mother2->GetMotherId();
1583 int grandmotherId3 = mother3->GetMotherId();
1584 int grandmotherId4 = mother4->GetMotherId();
1585
1586 if (mcMotherPdg1 == 22 && mcMotherPdg2 == 22 && mcMotherPdg3 == 111 && mcMotherPdg4 == 111) {
1587 if (grandmotherId1 != motherId3) continue;
1588 Pi0_pt_vs_rap->Fill(params.fRapidity, params.fPt);
1589 Pi0_pt_vs_rap_est->Fill(params.fRapidity, params.fPt);
1590 }
1591
1592 if (mcMotherPdg1 == 111 && mcMotherPdg2 == 111 && mcMotherPdg3 == 22 && mcMotherPdg4 == 22) {
1593 if (grandmotherId3 != motherId1) continue;
1594 Pi0_pt_vs_rap->Fill(params.fRapidity, params.fPt);
1595 Pi0_pt_vs_rap_est->Fill(params.fRapidity, params.fPt);
1596 }
1597
1598 if (mcMotherPdg1 == 22 && mcMotherPdg2 == 22 && mcMotherPdg3 == 22 && mcMotherPdg4 == 22) {
1599 if (grandmotherId1 != grandmotherId2 || grandmotherId3 != grandmotherId4 || grandmotherId1 != grandmotherId3)
1600 continue;
1601 if (grandmotherId1 == -1) continue;
1602 CbmMCTrack* grandmother1 = (CbmMCTrack*) fMcTracks->At(grandmotherId1);
1603 if (nullptr == grandmother1) continue;
1604 int mcGrandMotherPdg1 = grandmother1->GetPdgCode();
1605 if (mcGrandMotherPdg1 != 111) continue;
1606 Pi0_pt_vs_rap->Fill(params.fRapidity, params.fPt);
1607 Pi0_pt_vs_rap_est->Fill(params.fRapidity, params.fPt);
1608 }
1609 }
1610 }
1611}
1612
1613
1615// TARGET
1616{
1617 Int_t nof_Target = EMT_man_Event_Target_mbias2.size();
1618 cout << "Mixing in Manual(Target_mbias2) - nof entries " << nof_Target << endl;
1619 for (Int_t a = 0; a < nof_Target - 1; a++) {
1620 for (Int_t b = a + 1; b < nof_Target; b++) {
1622 continue; // to make sure that the photons are from two different events
1623 TVector3 e11 = EMT_man_pair_momenta_Target_mbias2[a][0];
1624 TVector3 e12 = EMT_man_pair_momenta_Target_mbias2[a][1];
1625 TVector3 e21 = EMT_man_pair_momenta_Target_mbias2[b][0];
1626 TVector3 e22 = EMT_man_pair_momenta_Target_mbias2[b][1];
1628
1639 }
1640 }
1641}
1642
1643
1645// OUTSIDE
1646{
1647 Int_t nof_Outside = EMT_man_Event_Outside_mbias2.size();
1648 cout << "Mixing in Manual(Outside_mbias2) - nof entries " << nof_Outside << endl;
1649 for (Int_t a = 0; a < nof_Outside - 1; a++) {
1650 for (Int_t b = a + 1; b < nof_Outside; b++) {
1652 continue; // to make sure that the photons are from two different events
1653 TVector3 e11 = EMT_man_pair_momenta_Outside_mbias2[a][0];
1654 TVector3 e12 = EMT_man_pair_momenta_Outside_mbias2[a][1];
1655 TVector3 e21 = EMT_man_pair_momenta_Outside_mbias2[b][0];
1656 TVector3 e22 = EMT_man_pair_momenta_Outside_mbias2[b][1];
1658
1669 }
1670 }
1671}
1672
1673
1675// BOTH
1676{
1677 Int_t nof_Both = EMT_man_Event_Both_mbias2.size();
1678 cout << "Mixing in Manual(Both_mbias2) - nof entries " << nof_Both << endl;
1679 for (Int_t a = 0; a < nof_Both - 1; a++) {
1680 for (Int_t b = a + 1; b < nof_Both; b++) {
1682 continue; // to make sure that the photons are from two different events
1683 TVector3 e11 = EMT_man_pair_momenta_Both_mbias2[a][0];
1684 TVector3 e12 = EMT_man_pair_momenta_Both_mbias2[a][1];
1685 TVector3 e21 = EMT_man_pair_momenta_Both_mbias2[b][0];
1686 TVector3 e22 = EMT_man_pair_momenta_Both_mbias2[b][1];
1688
1689 EMT_InvMass_all_Both_mbias2->Fill(params.fMinv);
1691 EMT_InvMass_zero_Both_mbias2->Fill(params.fMinv);
1693 EMT_InvMass_one_Both_mbias2->Fill(params.fMinv);
1695 EMT_InvMass_two_Both_mbias2->Fill(params.fMinv);
1699
1700 // separate by rap and Pt only for onetwo
1703 if (params.fRapidity > 1.2 && params.fRapidity <= 1.6) {
1704 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_1_mbias2->Fill(params.fMinv);
1705 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_2_mbias2->Fill(params.fMinv);
1706 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_3_mbias2->Fill(params.fMinv);
1707 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_4_mbias2->Fill(params.fMinv);
1708 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_5_mbias2->Fill(params.fMinv);
1709 }
1710 if (params.fRapidity > 1.6 && params.fRapidity <= 2.0) {
1711 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_6_mbias2->Fill(params.fMinv);
1712 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_7_mbias2->Fill(params.fMinv);
1713 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_8_mbias2->Fill(params.fMinv);
1714 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_9_mbias2->Fill(params.fMinv);
1715 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_10_mbias2->Fill(params.fMinv);
1716 }
1717 if (params.fRapidity > 2.0 && params.fRapidity <= 2.4) {
1718 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_11_mbias2->Fill(params.fMinv);
1719 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_12_mbias2->Fill(params.fMinv);
1720 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_13_mbias2->Fill(params.fMinv);
1721 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_14_mbias2->Fill(params.fMinv);
1722 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_15_mbias2->Fill(params.fMinv);
1723 }
1724 if (params.fRapidity > 2.4 && params.fRapidity <= 2.8) {
1725 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_16_mbias2->Fill(params.fMinv);
1726 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_17_mbias2->Fill(params.fMinv);
1727 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_18_mbias2->Fill(params.fMinv);
1728 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_19_mbias2->Fill(params.fMinv);
1729 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_20_mbias2->Fill(params.fMinv);
1730 }
1731 if (params.fRapidity > 2.8 && params.fRapidity <= 3.2) {
1732 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_21_mbias2->Fill(params.fMinv);
1733 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_22_mbias2->Fill(params.fMinv);
1734 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_23_mbias2->Fill(params.fMinv);
1735 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_24_mbias2->Fill(params.fMinv);
1736 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_25_mbias2->Fill(params.fMinv);
1737 }
1738 if (params.fRapidity > 3.2 && params.fRapidity <= 3.6) {
1739 if (params.fPt > 0.0 && params.fPt <= 0.4) rap_vs_Pt_InM_mixing_26_mbias2->Fill(params.fMinv);
1740 if (params.fPt > 0.4 && params.fPt <= 0.8) rap_vs_Pt_InM_mixing_27_mbias2->Fill(params.fMinv);
1741 if (params.fPt > 0.8 && params.fPt <= 1.2) rap_vs_Pt_InM_mixing_28_mbias2->Fill(params.fMinv);
1742 if (params.fPt > 1.2 && params.fPt <= 1.6) rap_vs_Pt_InM_mixing_29_mbias2->Fill(params.fMinv);
1743 if (params.fPt > 1.6 && params.fPt <= 2.0) rap_vs_Pt_InM_mixing_30_mbias2->Fill(params.fMinv);
1744 }
1745 }
1746 }
1747 }
1748}
1749
1750
1751int CbmKresConversionManualmbias2::FindInRich(int richInd, int stsMcTrackId)
1752{
1753 int RingsInRich = 0;
1754 if (richInd > 0) {
1755 CbmTrackMatchNew* richMatch = (CbmTrackMatchNew*) fRichRingMatches->At(richInd);
1756 if (richMatch != nullptr && richMatch->GetNofLinks() > 0) {
1757 int richMcTrackId = richMatch->GetMatchedLink().GetIndex();
1758 if (richMcTrackId > 0) {
1759 if (stsMcTrackId == richMcTrackId) { // check that global track was matched correctly for STS and RICH together
1760 CbmMCTrack* mcTrack2 = (CbmMCTrack*) fMcTracks->At(richMcTrackId);
1761 if (mcTrack2 != nullptr) {
1762 int pdgRICH = mcTrack2->GetPdgCode();
1763 if (TMath::Abs(pdgRICH) == 11) RingsInRich++;
1764 }
1765 }
1766 }
1767 }
1768 }
1769 return RingsInRich;
1770}
1771
1772
1774{
1775 int identified = 0;
1776
1777 if (nullptr != ring) {
1778 CbmRichRingLight ringHit;
1779 int nofHits = ring->GetNofHits();
1780 for (int i = 0; i < nofHits; i++) {
1781 int hitInd = ring->GetHit(i);
1782 CbmRichHit* hit = (CbmRichHit*) fRichHits->At(hitInd);
1783 if (nullptr == hit) continue;
1784 CbmRichHitLight hl(hit->GetX(), hit->GetY());
1785 ringHit.AddHit(hl);
1786 }
1787 fTauFit->DoFit(&ringHit);
1788 if (ringHit.GetAaxis() > 4 && ringHit.GetAaxis() < 6 && ringHit.GetBaxis() > 4 && ringHit.GetBaxis() < 6
1789 && momentum > 0.2 && momentum < 4.)
1790 identified++;
1791 //if (ring->GetDistance() < 2 && ringHit.GetAaxis() > 4 && ringHit.GetAaxis() < 6 && ringHit.GetBaxis() > 4 && ringHit.GetBaxis() < 6 && momentum > 0.2 && momentum < 4.) identified ++;
1792 }
1793
1794 return identified;
1795}
1796
1797
1799{
1800 double FinalAngle = 400;
1801 Int_t hits1sts = Sts_1->GetNofStsHits();
1802 Int_t hits2sts = Sts_2->GetNofStsHits();
1803 Int_t hits1mvd = Sts_1->GetNofMvdHits();
1804 Int_t hits2mvd = Sts_2->GetNofMvdHits();
1805
1806 double Xpart1 = 0;
1807 double Ypart1 = 0;
1808 double Zpart1 = 0;
1809 double Xpart2 = 0;
1810 double Ypart2 = 0;
1811 double Zpart2 = 0;
1812 if (hits1sts > 0) {
1813 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(hits1sts - 1);
1814 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
1815 Xpart1 = stsHit1->GetX();
1816 Ypart1 = stsHit1->GetY();
1817 Zpart1 = stsHit1->GetZ();
1818 }
1819 else {
1820 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(hits1mvd - 1);
1821 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
1822 Xpart1 = mvdHit1->GetX();
1823 Ypart1 = mvdHit1->GetY();
1824 Zpart1 = mvdHit1->GetZ();
1825 }
1826
1827 if (hits2sts > 0) {
1828 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(hits2sts - 1);
1829 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
1830 Xpart2 = stsHit2->GetX();
1831 Ypart2 = stsHit2->GetY();
1832 Zpart2 = stsHit2->GetZ();
1833 }
1834 else {
1835 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(hits2mvd - 1);
1836 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
1837 Xpart2 = mvdHit2->GetX();
1838 Ypart2 = mvdHit2->GetY();
1839 Zpart2 = mvdHit2->GetZ();
1840 }
1841
1842 // check difference in 2 cm, because of two slices of every STS and Mvd stations
1843 if (TMath::Abs(Zpart1 - Zpart2) > 2
1844 && Zpart1 > Zpart2) { // if last hits are in different stations --> try to find the latest common station
1845 for (int i = hits1sts - 2; i > -1; i--) { // start from second last station
1846 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
1847 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1848 Xpart1 = stsHit->GetX();
1849 Ypart1 = stsHit->GetY();
1850 Zpart1 = stsHit->GetZ();
1851 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1852 }
1853 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1854 for (int i = hits1mvd - 2; i > -1; i--) {
1855 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
1856 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1857 Xpart1 = mvdHit->GetX();
1858 Ypart1 = mvdHit->GetY();
1859 Zpart1 = mvdHit->GetZ();
1860 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1861 }
1862 }
1863 }
1864
1865 if (TMath::Abs(Zpart1 - Zpart2) > 2
1866 && Zpart1 < Zpart2) { // if last hits are in different stations --> try to find the latest common station
1867 for (int i = hits2sts - 2; i > -1; i--) { // start from second last station
1868 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
1869 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1870 Xpart2 = stsHit->GetX();
1871 Ypart2 = stsHit->GetY();
1872 Zpart2 = stsHit->GetZ();
1873 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1874 }
1875 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
1876 for (int i = hits2mvd - 2; i > -1; i--) {
1877 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
1878 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1879 Xpart2 = mvdHit->GetX();
1880 Ypart2 = mvdHit->GetY();
1881 Zpart2 = mvdHit->GetZ();
1882 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1883 }
1884 }
1885 }
1886
1887 // calculate angle if we found common station
1888 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1889 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
1890 }
1891
1892 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1893}
1894
1896{
1897 double FinalAngle = 400;
1898 Int_t hits1sts = Sts_1->GetNofStsHits();
1899 Int_t hits2sts = Sts_2->GetNofStsHits();
1900 Int_t hits1mvd = Sts_1->GetNofMvdHits();
1901 Int_t hits2mvd = Sts_2->GetNofMvdHits();
1902
1903 double Xpart1 = 0;
1904 double Ypart1 = 0;
1905 double Zpart1 = 0;
1906 double Xpart2 = 0;
1907 double Ypart2 = 0;
1908 double Zpart2 = 0;
1909 if (hits1mvd > 0) {
1910 Int_t mvdHitIndex1 = Sts_1->GetMvdHitIndex(0);
1911 CbmMvdHit* mvdHit1 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex1);
1912 Xpart1 = mvdHit1->GetX();
1913 Ypart1 = mvdHit1->GetY();
1914 Zpart1 = mvdHit1->GetZ();
1915 }
1916 else {
1917 Int_t stsHitIndex1 = Sts_1->GetStsHitIndex(0);
1918 CbmStsHit* stsHit1 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex1);
1919 Xpart1 = stsHit1->GetX();
1920 Ypart1 = stsHit1->GetY();
1921 Zpart1 = stsHit1->GetZ();
1922 }
1923
1924 if (hits2mvd > 0) {
1925 Int_t mvdHitIndex2 = Sts_2->GetMvdHitIndex(0);
1926 CbmMvdHit* mvdHit2 = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex2);
1927 Xpart2 = mvdHit2->GetX();
1928 Ypart2 = mvdHit2->GetY();
1929 Zpart2 = mvdHit2->GetZ();
1930 }
1931 else {
1932 Int_t stsHitIndex2 = Sts_2->GetStsHitIndex(0);
1933 CbmStsHit* stsHit2 = (CbmStsHit*) fArrayStsHit->At(stsHitIndex2);
1934 Xpart2 = stsHit2->GetX();
1935 Ypart2 = stsHit2->GetY();
1936 Zpart2 = stsHit2->GetZ();
1937 }
1938
1939 // check difference in 2 cm, because of two slices of every STS and Mvd stations
1940 if (TMath::Abs(Zpart1 - Zpart2) > 2
1941 && Zpart1 < Zpart2) { // if first hits are in different stations --> try to find the earliest common station
1942 for (int i = 1; i < hits1mvd; i++) { // start from second hit
1943 Int_t mvdHitIndex = Sts_1->GetMvdHitIndex(i);
1944 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1945 Xpart1 = mvdHit->GetX();
1946 Ypart1 = mvdHit->GetY();
1947 Zpart1 = mvdHit->GetZ();
1948 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1949 }
1950 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 < Zpart2) {
1951 for (int i = 0; i < hits1sts; i++) {
1952 Int_t stsHitIndex = Sts_1->GetStsHitIndex(i);
1953 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1954 Xpart1 = stsHit->GetX();
1955 Ypart1 = stsHit->GetY();
1956 Zpart1 = stsHit->GetZ();
1957 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1958 }
1959 }
1960 }
1961
1962 if (TMath::Abs(Zpart1 - Zpart2) > 2
1963 && Zpart1 > Zpart2) { // if first hits are in different stations --> try to find the earliest common station
1964 for (int i = 1; i < hits2mvd; i++) { // start from second hit
1965 Int_t mvdHitIndex = Sts_2->GetMvdHitIndex(i);
1966 CbmMvdHit* mvdHit = (CbmMvdHit*) fArrayMvdHit->At(mvdHitIndex);
1967 Xpart2 = mvdHit->GetX();
1968 Ypart2 = mvdHit->GetY();
1969 Zpart2 = mvdHit->GetZ();
1970 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1971 }
1972 if (TMath::Abs(Zpart1 - Zpart2) > 2 && Zpart1 > Zpart2) {
1973 for (int i = 0; i < hits2sts; i++) {
1974 Int_t stsHitIndex = Sts_2->GetStsHitIndex(i);
1975 CbmStsHit* stsHit = (CbmStsHit*) fArrayStsHit->At(stsHitIndex);
1976 Xpart2 = stsHit->GetX();
1977 Ypart2 = stsHit->GetY();
1978 Zpart2 = stsHit->GetZ();
1979 if (TMath::Abs(Zpart1 - Zpart2) < 2) continue;
1980 }
1981 }
1982 }
1983
1984 // calculate angle if we found common station
1985 if (TMath::Abs(Zpart1 - Zpart2) < 2 && Zpart1 != 0 && Zpart2 != 0) {
1986 FinalAngle = TMath::ATan2(Ypart1 - Ypart2, Xpart1 - Xpart2) * (180 / TMath::Pi());
1987 }
1988
1989 return TMath::Abs(TMath::Abs(FinalAngle) - 180);
1990}
1991
1992
1994{
1995 gDirectory->mkdir("Target_mbias2");
1996 gDirectory->cd("Target_mbias2");
1997 gDirectory->mkdir("CheckCuts_Target_mbias2");
1998 gDirectory->cd("CheckCuts_Target_mbias2");
1999 for (UInt_t i = 0; i < fHistoList_man_cuts_Target_mbias2.size(); i++) {
2001 }
2002 gDirectory->cd("..");
2003 gDirectory->mkdir("multiplicity");
2004 gDirectory->cd("multiplicity");
2005 for (UInt_t i = 0; i < fHistoList_multiplicity_man_Target_mbias2.size(); i++) {
2007 }
2008 gDirectory->cd("..");
2009 gDirectory->mkdir("all");
2010 gDirectory->cd("all");
2011 for (UInt_t i = 0; i < fHistoList_man_all_Target_mbias2.size(); i++) {
2013 }
2014 gDirectory->cd("..");
2015 gDirectory->mkdir("zero");
2016 gDirectory->cd("zero");
2017 for (UInt_t i = 0; i < fHistoList_man_zero_Target_mbias2.size(); i++) {
2019 }
2020 gDirectory->cd("..");
2021 gDirectory->mkdir("one");
2022 gDirectory->cd("one");
2023 for (UInt_t i = 0; i < fHistoList_man_one_Target_mbias2.size(); i++) {
2025 }
2026 gDirectory->cd("..");
2027 gDirectory->mkdir("two");
2028 gDirectory->cd("two");
2029 for (UInt_t i = 0; i < fHistoList_man_two_Target_mbias2.size(); i++) {
2031 }
2032 gDirectory->cd("..");
2033 gDirectory->mkdir("onetwo");
2034 gDirectory->cd("onetwo");
2035 for (UInt_t i = 0; i < fHistoList_man_onetwo_Target_mbias2.size(); i++) {
2037 }
2038 gDirectory->cd("..");
2039
2040 gDirectory->cd("..");
2041
2042
2043 gDirectory->mkdir("Outside_mbias2");
2044 gDirectory->cd("Outside_mbias2");
2045 gDirectory->mkdir("CheckCuts_Outside_mbias2");
2046 gDirectory->cd("CheckCuts_Outside_mbias2");
2047 for (UInt_t i = 0; i < fHistoList_man_cuts_Outside_mbias2.size(); i++) {
2049 }
2050 gDirectory->cd("..");
2051 gDirectory->mkdir("multiplicity");
2052 gDirectory->cd("multiplicity");
2053 for (UInt_t i = 0; i < fHistoList_multiplicity_man_Outside_mbias2.size(); i++) {
2055 }
2056 gDirectory->cd("..");
2057 gDirectory->mkdir("all");
2058 gDirectory->cd("all");
2059 for (UInt_t i = 0; i < fHistoList_man_all_Outside_mbias2.size(); i++) {
2061 }
2062 gDirectory->cd("..");
2063 gDirectory->mkdir("zero");
2064 gDirectory->cd("zero");
2065 for (UInt_t i = 0; i < fHistoList_man_zero_Outside_mbias2.size(); i++) {
2067 }
2068 gDirectory->cd("..");
2069 gDirectory->mkdir("one");
2070 gDirectory->cd("one");
2071 for (UInt_t i = 0; i < fHistoList_man_one_Outside_mbias2.size(); i++) {
2073 }
2074 gDirectory->cd("..");
2075 gDirectory->mkdir("two");
2076 gDirectory->cd("two");
2077 for (UInt_t i = 0; i < fHistoList_man_two_Outside_mbias2.size(); i++) {
2079 }
2080 gDirectory->cd("..");
2081 gDirectory->mkdir("onetwo");
2082 gDirectory->cd("onetwo");
2083 for (UInt_t i = 0; i < fHistoList_man_onetwo_Outside_mbias2.size(); i++) {
2085 }
2086 gDirectory->cd("..");
2087 gDirectory->cd("..");
2088
2089
2090 gDirectory->mkdir("Both_mbias2");
2091 gDirectory->cd("Both_mbias2");
2092 gDirectory->mkdir("CheckCuts_Both_mbias2");
2093 gDirectory->cd("CheckCuts_Both_mbias2");
2094 for (UInt_t i = 0; i < fHistoList_man_cuts_Both_mbias2.size(); i++) {
2096 }
2097 gDirectory->cd("..");
2098 gDirectory->mkdir("multiplicity");
2099 gDirectory->cd("multiplicity");
2100 for (UInt_t i = 0; i < fHistoList_multiplicity_man_Both_mbias2.size(); i++) {
2102 }
2103 gDirectory->cd("..");
2104 gDirectory->mkdir("rap_vs_Pt_onetwo_Both_mbias2");
2105 gDirectory->cd("rap_vs_Pt_onetwo_Both_mbias2");
2106 for (UInt_t i = 0; i < fHistoList_rap_vs_pt_InM_mbias2.size(); i++) {
2108 }
2109 gDirectory->cd("..");
2110 gDirectory->mkdir("BG_Both_mbias2");
2111 gDirectory->cd("BG_Both_mbias2");
2112 gDirectory->mkdir("BG_Both_mbias2_all");
2113 gDirectory->cd("BG_Both_mbias2_all");
2114 for (UInt_t i = 0; i < fHistoList_bg_InM_all_Both_mbias2.size(); i++) {
2116 }
2117 gDirectory->cd("..");
2118 gDirectory->mkdir("BG_Both_mbias2_zero");
2119 gDirectory->cd("BG_Both_mbias2_zero");
2120 for (UInt_t i = 0; i < fHistoList_bg_InM_zero_Both_mbias2.size(); i++) {
2122 }
2123 gDirectory->cd("..");
2124 gDirectory->mkdir("BG_Both_mbias2_one");
2125 gDirectory->cd("BG_Both_mbias2_one");
2126 for (UInt_t i = 0; i < fHistoList_bg_InM_one_Both_mbias2.size(); i++) {
2128 }
2129 gDirectory->cd("..");
2130 gDirectory->mkdir("BG_Both_mbias2_two");
2131 gDirectory->cd("BG_Both_mbias2_two");
2132 for (UInt_t i = 0; i < fHistoList_bg_InM_two_Both_mbias2.size(); i++) {
2134 }
2135 gDirectory->cd("..");
2136 gDirectory->mkdir("BG_Both_mbias2_onetwo");
2137 gDirectory->cd("BG_Both_mbias2_onetwo");
2138 for (UInt_t i = 0; i < fHistoList_bg_InM_onetwo_Both_mbias2.size(); i++) {
2140 }
2141 gDirectory->cd("..");
2142 gDirectory->cd("..");
2143 gDirectory->mkdir("all");
2144 gDirectory->cd("all");
2145 for (UInt_t i = 0; i < fHistoList_man_all_Both_mbias2.size(); i++) {
2147 }
2148 gDirectory->cd("..");
2149 gDirectory->mkdir("zero");
2150 gDirectory->cd("zero");
2151 for (UInt_t i = 0; i < fHistoList_man_zero_Both_mbias2.size(); i++) {
2153 }
2154 gDirectory->cd("..");
2155 gDirectory->mkdir("one");
2156 gDirectory->cd("one");
2157 for (UInt_t i = 0; i < fHistoList_man_one_Both_mbias2.size(); i++) {
2159 }
2160 gDirectory->cd("..");
2161 gDirectory->mkdir("two");
2162 gDirectory->cd("two");
2163 for (UInt_t i = 0; i < fHistoList_man_two_Both_mbias2.size(); i++) {
2165 }
2166 gDirectory->cd("..");
2167 gDirectory->mkdir("onetwo");
2168 gDirectory->cd("onetwo");
2169 for (UInt_t i = 0; i < fHistoList_man_onetwo_Both_mbias2.size(); i++) {
2171 }
2172 gDirectory->cd("..");
2173 for (UInt_t i = 0; i < fHistoList_man_Both_mbias2.size(); i++) {
2174 fHistoList_man_Both_mbias2[i]->Write();
2175 }
2176 gDirectory->cd("..");
2177}
2178
2180{
2182 InvMass_vs_OA_candidates_Both_mbias2 = new TH2D("InvMass_vs_OA_candidates_Both_mbias2",
2183 "InvMass_vs_OA_candidates_Both_mbias2; invariant mass in "
2184 "GeV/c^{2}; opening angle in degree",
2185 500, 0, 0.5, 500, 0, 50);
2187 InvMass_vs_OA_fromPi0_Both_mbias2 = new TH2D("InvMass_vs_OA_fromPi0_Both_mbias2",
2188 "InvMass_vs_OA_fromPi0_Both_mbias2; invariant mass in GeV/c^{2}; "
2189 "opening angle in degree",
2190 500, 0, 0.5, 500, 0, 50);
2193 new TH1D("GammasInvMass_candidates_Both_mbias2",
2194 "GammasInvMass_candidates_Both_mbias2; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2197 "GammasOA_candidates_Both_mbias2", "GammasOA_candidates_Both_mbias2; opening angle in degree;#", 300, -0.1, 29.9);
2200 new TH1D("GammasInvMass_fromPi0_Both_mbias2", "GammasInvMass_fromPi0_Both_mbias2; invariant mass in GeV/c^{2};#",
2201 510, -0.01, 0.5);
2203 GammasOA_fromPi0_Both_mbias2 = new TH1D("GammasOA_fromPi0_Both_mbias2",
2204 "GammasOA_fromPi0_Both_mbias2; opening angle in degree;#", 300, -0.1, 29.9);
2207 new TH1D("PlaneAngles_last_candidates_Both_mbias2",
2208 "PlaneAngles_last_candidates_Both_mbias2; #theta angle in degree;#", 720, -1., 179.);
2211 new TH1D("PlaneAngles_last_fromPi0_Both_mbias2", "PlaneAngles_last_fromPi0_Both_mbias2; #theta angle in degree;#",
2212 720, -1., 179.);
2215 new TH1D("PlaneAngles_first_candidates_Both_mbias2",
2216 "PlaneAngles_first_candidates_Both_mbias2; #theta angle in degree;#", 720, -1., 179.);
2219 new TH1D("PlaneAngles_first_fromPi0_Both_mbias2", "PlaneAngles_first_fromPi0_Both_mbias2; #theta angle in degree;#",
2220 720, -1., 179.);
2222
2223
2224 InvMass_vs_OA_candidates_Target_mbias2 = new TH2D("InvMass_vs_OA_candidates_Target_mbias2",
2225 "InvMass_vs_OA_candidates_Target_mbias2; invariant mass in "
2226 "GeV/c^{2}; opening angle in degree",
2227 500, 0, 0.5, 500, 0, 50);
2229 InvMass_vs_OA_fromPi0_Target_mbias2 = new TH2D("InvMass_vs_OA_fromPi0_Target_mbias2",
2230 "InvMass_vs_OA_fromPi0_Target_mbias2; invariant mass in "
2231 "GeV/c^{2}; opening angle in degree",
2232 500, 0, 0.5, 500, 0, 50);
2235 new TH1D("GammasInvMass_candidates_Target_mbias2",
2236 "GammasInvMass_candidates_Target_mbias2; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2239 new TH1D("GammasOA_candidates_Target_mbias2", "GammasOA_candidates_Target_mbias2; opening angle in degree;#", 300,
2240 -0.1, 29.9);
2243 new TH1D("GammasInvMass_fromPi0_Target_mbias2",
2244 "GammasInvMass_fromPi0_Target_mbias2; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2247 "GammasOA_fromPi0_Target_mbias2", "GammasOA_fromPi0_Target_mbias2; opening angle in degree;#", 300, -0.1, 29.9);
2250 new TH1D("PlaneAngles_last_candidates_Target_mbias2",
2251 "PlaneAngles_last_candidates_Target_mbias2; #theta angle in degree;#", 720, -1., 179.);
2254 new TH1D("PlaneAngles_last_fromPi0_Target_mbias2",
2255 "PlaneAngles_last_fromPi0_Target_mbias2; #theta angle in degree;#", 720, -1., 179.);
2258 new TH1D("PlaneAngles_first_candidates_Target_mbias2",
2259 "PlaneAngles_first_candidates_Target_mbias2; #theta angle in degree;#", 720, -1., 179.);
2262 new TH1D("PlaneAngles_first_fromPi0_Target_mbias2",
2263 "PlaneAngles_first_fromPi0_Target_mbias2; #theta angle in degree;#", 720, -1., 179.);
2265
2266
2267 InvMass_vs_OA_candidates_Outside_mbias2 = new TH2D("InvMass_vs_OA_candidates_Outside_mbias2",
2268 "InvMass_vs_OA_candidates_Outside_mbias2; invariant mass in "
2269 "GeV/c^{2}; opening angle in degree",
2270 500, 0, 0.5, 500, 0, 50);
2272 InvMass_vs_OA_fromPi0_Outside_mbias2 = new TH2D("InvMass_vs_OA_fromPi0_Outside_mbias2",
2273 "InvMass_vs_OA_fromPi0_Outside_mbias2; invariant mass in "
2274 "GeV/c^{2}; opening angle in degree",
2275 500, 0, 0.5, 500, 0, 50);
2278 new TH1D("GammasInvMass_candidates_Outside_mbias2",
2279 "GammasInvMass_candidates_Outside_mbias2; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2282 new TH1D("GammasOA_candidates_Outside_mbias2", "GammasOA_candidates_Outside_mbias2; opening angle in degree;#", 300,
2283 -0.1, 29.9);
2286 new TH1D("GammasInvMass_fromPi0_Outside_mbias2",
2287 "GammasInvMass_fromPi0_Outside_mbias2; invariant mass in GeV/c^{2};#", 510, -0.01, 0.5);
2290 "GammasOA_fromPi0_Outside_mbias2", "GammasOA_fromPi0_Outside_mbias2; opening angle in degree;#", 300, -0.1, 29.9);
2293 new TH1D("PlaneAngles_last_candidates_Outside_mbias2",
2294 "PlaneAngles_last_candidates_Outside_mbias2; #theta angle in degree;#", 720, -1., 179.);
2297 new TH1D("PlaneAngles_last_fromPi0_Outside_mbias2",
2298 "PlaneAngles_last_fromPi0_Outside_mbias2; #theta angle in degree;#", 720, -1., 179.);
2301 new TH1D("PlaneAngles_first_candidates_Outside_mbias2",
2302 "PlaneAngles_first_candidates_Outside_mbias2; #theta angle in degree;#", 720, -1., 179.);
2305 new TH1D("PlaneAngles_first_fromPi0_Outside_mbias2",
2306 "PlaneAngles_first_fromPi0_Outside_mbias2; #theta angle in degree;#", 720, -1., 179.);
2308
2309
2310 // Target_mbias2 => all
2312 new TH1D("GammaInvMassReco_all_Target_mbias2", "GammaInvMassReco_all_Target_mbias2; invariant mass in GeV/c^{2};#",
2313 110, -0.01, 0.1);
2316 new TH1D("GammaOpeningAngleReco_all_Target_mbias2", "GammaOpeningAngleReco_all_Target_mbias2; angle [deg];#", 200,
2317 -0.1, 19.9);
2319 Pdg_all_Target_mbias2 = new TH1D("Pdg_all_Target_mbias2", "Pdg_all_Target_mbias2; Id;#", 800, 0, 400);
2322 new TH1D("P_reco_all_Target_mbias2", "P_reco_all_Target_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2325 new TH1D("Pt_reco_all_Target_mbias2", "Pt_reco_all_Target_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2328 new TH1D("Pi0InvMassReco_all_Target_mbias2", "Pi0InvMassReco_all_Target_mbias2; invariant mass in GeV/c^{2};#",
2329 1000, 0, 2.0);
2332 "EMT_InvMass_all_Target_mbias2", "EMT_InvMass_all_Target_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2335 new TH2D("Pi0_pt_vs_rap_all_Target_mbias2", "Pi0_pt_vs_rap_all_Target_mbias2; rapidity y; p_{t} in GeV/c ", 90, -2.,
2336 7., 60, -1., 5.);
2339 new TH2D("Pi0_pt_vs_rap_est_all_Target_mbias2", "Pi0_pt_vs_rap_est_all_Target_mbias2; rapidity y; p_{t} in GeV/c ",
2340 10, 0., 4., 10, 0., 4.);
2342
2343
2344 // Target_mbias2 => zero
2346 new TH1D("GammaInvMassReco_zero_Target_mbias2",
2347 "GammaInvMassReco_zero_Target_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2350 new TH1D("GammaOpeningAngleReco_zero_Target_mbias2", "GammaOpeningAngleReco_zero_Target_mbias2; angle [deg];#", 200,
2351 -0.1, 19.9);
2353 Pdg_zero_Target_mbias2 = new TH1D("Pdg_zero_Target_mbias2", "Pdg_zero_Target_mbias2; Id;#", 800, 0, 400);
2356 new TH1D("P_reco_zero_Target_mbias2", "P_reco_zero_Target_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2359 new TH1D("Pt_reco_zero_Target_mbias2", "Pt_reco_zero_Target_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2362 new TH1D("Pi0InvMassReco_zero_Target_mbias2", "Pi0InvMassReco_zero_Target_mbias2; invariant mass in GeV/c^{2};#",
2363 1000, 0, 2.0);
2366 "EMT_InvMass_zero_Target_mbias2", "EMT_InvMass_zero_Target_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2369 new TH2D("Pi0_pt_vs_rap_zero_Target_mbias2", "Pi0_pt_vs_rap_zero_Target_mbias2; rapidity y; p_{t} in GeV/c ", 90,
2370 -2., 7., 60, -1., 5.);
2373 new TH2D("Pi0_pt_vs_rap_est_zero_Target_mbias2",
2374 "Pi0_pt_vs_rap_est_zero_Target_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2376
2377
2378 // Target_mbias2 => one
2380 new TH1D("GammaInvMassReco_one_Target_mbias2", "GammaInvMassReco_one_Target_mbias2; invariant mass in GeV/c^{2};#",
2381 110, -0.01, 0.1);
2384 new TH1D("GammaOpeningAngleReco_one_Target_mbias2", "GammaOpeningAngleReco_one_Target_mbias2; angle [deg];#", 200,
2385 -0.1, 19.9);
2387 Pdg_one_Target_mbias2 = new TH1D("Pdg_one_Target_mbias2", "Pdg_one_Target_mbias2; Id;#", 800, 0, 400);
2390 new TH1D("P_reco_one_Target_mbias2", "P_reco_one_Target_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2393 new TH1D("Pt_reco_one_Target_mbias2", "Pt_reco_one_Target_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2396 new TH1D("Pi0InvMassReco_one_Target_mbias2", "Pi0InvMassReco_one_Target_mbias2; invariant mass in GeV/c^{2};#",
2397 1000, 0, 2.0);
2400 "EMT_InvMass_one_Target_mbias2", "EMT_InvMass_one_Target_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2403 new TH2D("Pi0_pt_vs_rap_one_Target_mbias2", "Pi0_pt_vs_rap_one_Target_mbias2; rapidity y; p_{t} in GeV/c ", 90, -2.,
2404 7., 60, -1., 5.);
2407 new TH2D("Pi0_pt_vs_rap_est_one_Target_mbias2", "Pi0_pt_vs_rap_est_one_Target_mbias2; rapidity y; p_{t} in GeV/c ",
2408 10, 0., 4., 10, 0., 4.);
2410
2411
2412 // Target_mbias2 => two
2414 new TH1D("GammaInvMassReco_two_Target_mbias2", "GammaInvMassReco_two_Target_mbias2; invariant mass in GeV/c^{2};#",
2415 110, -0.01, 0.1);
2418 new TH1D("GammaOpeningAngleReco_two_Target_mbias2", "GammaOpeningAngleReco_two_Target_mbias2; angle [deg];#", 200,
2419 -0.1, 19.9);
2421 Pdg_two_Target_mbias2 = new TH1D("Pdg_two_Target_mbias2", "Pdg_two_Target_mbias2; Id;#", 800, 0, 400);
2424 new TH1D("P_reco_two_Target_mbias2", "P_reco_two_Target_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2427 new TH1D("Pt_reco_two_Target_mbias2", "Pt_reco_two_Target_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2430 new TH1D("Pi0InvMassReco_two_Target_mbias2", "Pi0InvMassReco_two_Target_mbias2; invariant mass in GeV/c^{2};#",
2431 1000, 0, 2.0);
2434 "EMT_InvMass_two_Target_mbias2", "EMT_InvMass_two_Target_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2437 new TH2D("Pi0_pt_vs_rap_two_Target_mbias2", "Pi0_pt_vs_rap_two_Target_mbias2; rapidity y; p_{t} in GeV/c ", 90, -2.,
2438 7., 60, -1., 5.);
2441 new TH2D("Pi0_pt_vs_rap_est_two_Target_mbias2", "Pi0_pt_vs_rap_est_two_Target_mbias2; rapidity y; p_{t} in GeV/c ",
2442 10, 0., 4., 10, 0., 4.);
2444
2445
2446 // Target_mbias2 => onetwo
2448 new TH1D("GammaInvMassReco_onetwo_Target_mbias2",
2449 "GammaInvMassReco_onetwo_Target_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2452 new TH1D("GammaOpeningAngleReco_onetwo_Target_mbias2", "GammaOpeningAngleReco_onetwo_Target_mbias2; angle [deg];#",
2453 200, -0.1, 19.9);
2455 Pdg_onetwo_Target_mbias2 = new TH1D("Pdg_onetwo_Target_mbias2", "Pdg_onetwo_Target_mbias2; Id;#", 800, 0, 400);
2458 new TH1D("P_reco_onetwo_Target_mbias2", "P_reco_onetwo_Target_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2461 new TH1D("Pt_reco_onetwo_Target_mbias2", "Pt_reco_onetwo_Target_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2464 new TH1D("Pi0InvMassReco_onetwo_Target_mbias2",
2465 "Pi0InvMassReco_onetwo_Target_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2468 new TH1D("EMT_InvMass_onetwo_Target_mbias2", "EMT_InvMass_onetwo_Target_mbias2; invariant mass in GeV/c^{2};#",
2469 1000, 0, 2.0);
2472 new TH2D("Pi0_pt_vs_rap_onetwo_Target_mbias2", "Pi0_pt_vs_rap_onetwo_Target_mbias2; rapidity y; p_{t} in GeV/c ",
2473 90, -2., 7., 60, -1., 5.);
2476 new TH2D("Pi0_pt_vs_rap_est_onetwo_Target_mbias2",
2477 "Pi0_pt_vs_rap_est_onetwo_Target_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2479
2480
2481 // Outside_mbias2 => all
2483 new TH1D("GammaInvMassReco_all_Outside_mbias2",
2484 "GammaInvMassReco_all_Outside_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2487 new TH1D("GammaOpeningAngleReco_all_Outside_mbias2", "GammaOpeningAngleReco_all_Outside_mbias2; angle [deg];#", 200,
2488 -0.1, 19.9);
2490 Pdg_all_Outside_mbias2 = new TH1D("Pdg_all_Outside_mbias2", "Pdg_all_Outside_mbias2; Id;#", 800, 0, 400);
2493 new TH1D("P_reco_all_Outside_mbias2", "P_reco_all_Outside_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2496 new TH1D("Pt_reco_all_Outside_mbias2", "Pt_reco_all_Outside_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2499 new TH1D("Pi0InvMassReco_all_Outside_mbias2", "Pi0InvMassReco_all_Outside_mbias2; invariant mass in GeV/c^{2};#",
2500 1000, 0, 2.0);
2503 "EMT_InvMass_all_Outside_mbias2", "EMT_InvMass_all_Outside_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2506 new TH2D("Pi0_pt_vs_rap_all_Outside_mbias2", "Pi0_pt_vs_rap_all_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 90,
2507 -2., 7., 60, -1., 5.);
2510 new TH2D("Pi0_pt_vs_rap_est_all_Outside_mbias2",
2511 "Pi0_pt_vs_rap_est_all_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2513
2514
2515 // Outside_mbias2 => zero
2517 new TH1D("GammaInvMassReco_zero_Outside_mbias2",
2518 "GammaInvMassReco_zero_Outside_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2521 new TH1D("GammaOpeningAngleReco_zero_Outside_mbias2", "GammaOpeningAngleReco_zero_Outside_mbias2; angle [deg];#",
2522 200, -0.1, 19.9);
2524 Pdg_zero_Outside_mbias2 = new TH1D("Pdg_zero_Outside_mbias2", "Pdg_zero_Outside_mbias2; Id;#", 800, 0, 400);
2527 new TH1D("P_reco_zero_Outside_mbias2", "P_reco_zero_Outside_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2530 new TH1D("Pt_reco_zero_Outside_mbias2", "Pt_reco_zero_Outside_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2533 new TH1D("Pi0InvMassReco_zero_Outside_mbias2", "Pi0InvMassReco_zero_Outside_mbias2; invariant mass in GeV/c^{2};#",
2534 1000, 0, 2.0);
2537 "EMT_InvMass_zero_Outside_mbias2", "EMT_InvMass_zero_Outside_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2540 new TH2D("Pi0_pt_vs_rap_zero_Outside_mbias2", "Pi0_pt_vs_rap_zero_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 90,
2541 -2., 7., 60, -1., 5.);
2544 new TH2D("Pi0_pt_vs_rap_est_zero_Outside_mbias2",
2545 "Pi0_pt_vs_rap_est_zero_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2547
2548
2549 // Outside_mbias2 => one
2551 new TH1D("GammaInvMassReco_one_Outside_mbias2",
2552 "GammaInvMassReco_one_Outside_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2555 new TH1D("GammaOpeningAngleReco_one_Outside_mbias2", "GammaOpeningAngleReco_one_Outside_mbias2; angle [deg];#", 200,
2556 -0.1, 19.9);
2558 Pdg_one_Outside_mbias2 = new TH1D("Pdg_one_Outside_mbias2", "Pdg_one_Outside_mbias2; Id;#", 800, 0, 400);
2561 new TH1D("P_reco_one_Outside_mbias2", "P_reco_one_Outside_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2564 new TH1D("Pt_reco_one_Outside_mbias2", "Pt_reco_one_Outside_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2567 new TH1D("Pi0InvMassReco_one_Outside_mbias2", "Pi0InvMassReco_one_Outside_mbias2; invariant mass in GeV/c^{2};#",
2568 1000, 0, 2.0);
2571 "EMT_InvMass_one_Outside_mbias2", "EMT_InvMass_one_Outside_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2574 new TH2D("Pi0_pt_vs_rap_one_Outside_mbias2", "Pi0_pt_vs_rap_one_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 90,
2575 -2., 7., 60, -1., 5.);
2578 new TH2D("Pi0_pt_vs_rap_est_one_Outside_mbias2",
2579 "Pi0_pt_vs_rap_est_one_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2581
2582
2583 // Outside_mbias2 => two
2585 new TH1D("GammaInvMassReco_two_Outside_mbias2",
2586 "GammaInvMassReco_two_Outside_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2589 new TH1D("GammaOpeningAngleReco_two_Outside_mbias2", "GammaOpeningAngleReco_two_Outside_mbias2; angle [deg];#", 200,
2590 -0.1, 19.9);
2592 Pdg_two_Outside_mbias2 = new TH1D("Pdg_two_Outside_mbias2", "Pdg_two_Outside_mbias2; Id;#", 800, 0, 400);
2595 new TH1D("P_reco_two_Outside_mbias2", "P_reco_two_Outside_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2598 new TH1D("Pt_reco_two_Outside_mbias2", "Pt_reco_two_Outside_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2601 new TH1D("Pi0InvMassReco_two_Outside_mbias2", "Pi0InvMassReco_two_Outside_mbias2; invariant mass in GeV/c^{2};#",
2602 1000, 0, 2.0);
2605 "EMT_InvMass_two_Outside_mbias2", "EMT_InvMass_two_Outside_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2608 new TH2D("Pi0_pt_vs_rap_two_Outside_mbias2", "Pi0_pt_vs_rap_two_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 90,
2609 -2., 7., 60, -1., 5.);
2612 new TH2D("Pi0_pt_vs_rap_est_two_Outside_mbias2",
2613 "Pi0_pt_vs_rap_est_two_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2615
2616
2617 // Outside_mbias2 => onetwo
2619 new TH1D("GammaInvMassReco_onetwo_Outside_mbias2",
2620 "GammaInvMassReco_onetwo_Outside_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2623 new TH1D("GammaOpeningAngleReco_onetwo_Outside_mbias2",
2624 "GammaOpeningAngleReco_onetwo_Outside_mbias2; angle [deg];#", 200, -0.1, 19.9);
2626 Pdg_onetwo_Outside_mbias2 = new TH1D("Pdg_onetwo_Outside_mbias2", "Pdg_onetwo_Outside_mbias2; Id;#", 800, 0, 400);
2629 new TH1D("P_reco_onetwo_Outside_mbias2", "P_reco_onetwo_Outside_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2632 new TH1D("Pt_reco_onetwo_Outside_mbias2", "Pt_reco_onetwo_Outside_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2635 new TH1D("Pi0InvMassReco_onetwo_Outside_mbias2",
2636 "Pi0InvMassReco_onetwo_Outside_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2639 new TH1D("EMT_InvMass_onetwo_Outside_mbias2", "EMT_InvMass_onetwo_Outside_mbias2; invariant mass in GeV/c^{2};#",
2640 1000, 0, 2.0);
2643 new TH2D("Pi0_pt_vs_rap_onetwo_Outside_mbias2", "Pi0_pt_vs_rap_onetwo_Outside_mbias2; rapidity y; p_{t} in GeV/c ",
2644 90, -2., 7., 60, -1., 5.);
2647 new TH2D("Pi0_pt_vs_rap_est_onetwo_Outside_mbias2",
2648 "Pi0_pt_vs_rap_est_onetwo_Outside_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2650
2651
2652 // Both_mbias2 => all
2654 new TH1D("GammaInvMassReco_all_Both_mbias2", "GammaInvMassReco_all_Both_mbias2; invariant mass in GeV/c^{2};#", 110,
2655 -0.01, 0.1);
2658 "GammaOpeningAngleReco_all_Both_mbias2", "GammaOpeningAngleReco_all_Both_mbias2; angle [deg];#", 200, -0.1, 19.9);
2660 Pdg_all_Both_mbias2 = new TH1D("Pdg_all_Both_mbias2", "Pdg_all_Both_mbias2; Id;#", 800, 0, 400);
2662 P_reco_all_Both_mbias2 = new TH1D("P_reco_all_Both_mbias2", "P_reco_all_Both_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2665 new TH1D("Pt_reco_all_Both_mbias2", "Pt_reco_all_Both_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2668 "Pi0InvMassReco_all_Both_mbias2", "Pi0InvMassReco_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2671 new TH1D("EMT_InvMass_all_Both_mbias2", "EMT_InvMass_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2674 new TH2D("Pi0_pt_vs_rap_all_Both_mbias2", "Pi0_pt_vs_rap_all_Both_mbias2; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2675 60, -1., 5.);
2678 new TH2D("Pi0_pt_vs_rap_est_all_Both_mbias2", "Pi0_pt_vs_rap_est_all_Both_mbias2; rapidity y; p_{t} in GeV/c ", 10,
2679 0., 4., 10, 0., 4.);
2681
2682
2683 // Both_mbias2 => zero
2685 new TH1D("GammaInvMassReco_zero_Both_mbias2", "GammaInvMassReco_zero_Both_mbias2; invariant mass in GeV/c^{2};#",
2686 110, -0.01, 0.1);
2689 "GammaOpeningAngleReco_zero_Both_mbias2", "GammaOpeningAngleReco_zero_Both_mbias2; angle [deg];#", 200, -0.1, 19.9);
2691 Pdg_zero_Both_mbias2 = new TH1D("Pdg_zero_Both_mbias2", "Pdg_zero_Both_mbias2; Id;#", 800, 0, 400);
2693 P_reco_zero_Both_mbias2 = new TH1D("P_reco_zero_Both_mbias2", "P_reco_zero_Both_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2696 new TH1D("Pt_reco_zero_Both_mbias2", "Pt_reco_zero_Both_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2699 "Pi0InvMassReco_zero_Both_mbias2", "Pi0InvMassReco_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2701 EMT_InvMass_zero_Both_mbias2 = new TH1D("EMT_InvMass_zero_Both_mbias2",
2702 "EMT_InvMass_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2705 new TH2D("Pi0_pt_vs_rap_zero_Both_mbias2", "Pi0_pt_vs_rap_zero_Both_mbias2; rapidity y; p_{t} in GeV/c ", 90, -2.,
2706 7., 60, -1., 5.);
2709 new TH2D("Pi0_pt_vs_rap_est_zero_Both_mbias2", "Pi0_pt_vs_rap_est_zero_Both_mbias2; rapidity y; p_{t} in GeV/c ",
2710 10, 0., 4., 10, 0., 4.);
2712
2713
2714 // Both_mbias2 => one
2716 new TH1D("GammaInvMassReco_one_Both_mbias2", "GammaInvMassReco_one_Both_mbias2; invariant mass in GeV/c^{2};#", 110,
2717 -0.01, 0.1);
2720 "GammaOpeningAngleReco_one_Both_mbias2", "GammaOpeningAngleReco_one_Both_mbias2; angle [deg];#", 200, -0.1, 19.9);
2722 Pdg_one_Both_mbias2 = new TH1D("Pdg_one_Both_mbias2", "Pdg_one_Both_mbias2; Id;#", 800, 0, 400);
2724 P_reco_one_Both_mbias2 = new TH1D("P_reco_one_Both_mbias2", "P_reco_one_Both_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2727 new TH1D("Pt_reco_one_Both_mbias2", "Pt_reco_one_Both_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2730 "Pi0InvMassReco_one_Both_mbias2", "Pi0InvMassReco_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2733 new TH1D("EMT_InvMass_one_Both_mbias2", "EMT_InvMass_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2736 new TH2D("Pi0_pt_vs_rap_one_Both_mbias2", "Pi0_pt_vs_rap_one_Both_mbias2; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2737 60, -1., 5.);
2740 new TH2D("Pi0_pt_vs_rap_est_one_Both_mbias2", "Pi0_pt_vs_rap_est_one_Both_mbias2; rapidity y; p_{t} in GeV/c ", 10,
2741 0., 4., 10, 0., 4.);
2743
2744
2745 // Both_mbias2 => two
2747 new TH1D("GammaInvMassReco_two_Both_mbias2", "GammaInvMassReco_two_Both_mbias2; invariant mass in GeV/c^{2};#", 110,
2748 -0.01, 0.1);
2751 "GammaOpeningAngleReco_two_Both_mbias2", "GammaOpeningAngleReco_two_Both_mbias2; angle [deg];#", 200, -0.1, 19.9);
2753 Pdg_two_Both_mbias2 = new TH1D("Pdg_two_Both_mbias2", "Pdg_two_Both_mbias2; Id;#", 800, 0, 400);
2755 P_reco_two_Both_mbias2 = new TH1D("P_reco_two_Both_mbias2", "P_reco_two_Both_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2758 new TH1D("Pt_reco_two_Both_mbias2", "Pt_reco_two_Both_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2761 "Pi0InvMassReco_two_Both_mbias2", "Pi0InvMassReco_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2764 new TH1D("EMT_InvMass_two_Both_mbias2", "EMT_InvMass_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2767 new TH2D("Pi0_pt_vs_rap_two_Both_mbias2", "Pi0_pt_vs_rap_two_Both_mbias2; rapidity y; p_{t} in GeV/c ", 90, -2., 7.,
2768 60, -1., 5.);
2771 new TH2D("Pi0_pt_vs_rap_est_two_Both_mbias2", "Pi0_pt_vs_rap_est_two_Both_mbias2; rapidity y; p_{t} in GeV/c ", 10,
2772 0., 4., 10, 0., 4.);
2774
2775
2776 // Both_mbias2 => onetwo
2778 new TH1D("GammaInvMassReco_onetwo_Both_mbias2",
2779 "GammaInvMassReco_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 110, -0.01, 0.1);
2782 new TH1D("GammaOpeningAngleReco_onetwo_Both_mbias2", "GammaOpeningAngleReco_onetwo_Both_mbias2; angle [deg];#", 200,
2783 -0.1, 19.9);
2785 Pdg_onetwo_Both_mbias2 = new TH1D("Pdg_onetwo_Both_mbias2", "Pdg_onetwo_Both_mbias2; Id;#", 800, 0, 400);
2788 new TH1D("P_reco_onetwo_Both_mbias2", "P_reco_onetwo_Both_mbias2; P in GeV/c^{2};#", 600, 0, 6);
2791 new TH1D("Pt_reco_onetwo_Both_mbias2", "Pt_reco_onetwo_Both_mbias2; P_{t} in GeV/c^{2};#", 300, 0, 3);
2794 new TH1D("Pi0InvMassReco_onetwo_Both_mbias2", "Pi0InvMassReco_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#",
2795 1000, 0, 2.0);
2798 "EMT_InvMass_onetwo_Both_mbias2", "EMT_InvMass_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
2801 new TH2D("Pi0_pt_vs_rap_onetwo_Both_mbias2", "Pi0_pt_vs_rap_onetwo_Both_mbias2; rapidity y; p_{t} in GeV/c ", 90,
2802 -2., 7., 60, -1., 5.);
2805 new TH2D("Pi0_pt_vs_rap_est_onetwo_Both_mbias2",
2806 "Pi0_pt_vs_rap_est_onetwo_Both_mbias2; rapidity y; p_{t} in GeV/c ", 10, 0., 4., 10, 0., 4.);
2808
2809
2810 // Both_mbias2
2812 new TH2D("Pdg_vs_Distance_mbias2", "Pdg_vs_Distance_mbias2; pdg; distance in cm", 2500, 0, 2499, 500, 0, 50);
2814 P_vs_Distance_mbias2 = new TH2D("P_vs_Distance_mbias2",
2815 "Distance between projected track and center of the ring (for e+ "
2816 "and e-); P in GeV/c^{2}; distance in cm",
2817 300, 0, 3, 300, 0, 15);
2819
2820
2821 // Multiplicity Target_mbias2
2822 MultiplicityGamma_all_Target_mbias2 = new TH2D("MultiplicityGamma_all_Target_mbias2",
2823 "MultiplicityGamma_all_Target_mbias2; Nof gammas in event; "
2824 "invariant mass in GeV/c^{2};#",
2825 400, 0, 30, 1000, 0, 2.0);
2827 MultiplicityGamma_zero_Target_mbias2 = new TH2D("MultiplicityGamma_zero_Target_mbias2",
2828 "MultiplicityGamma_zero_Target_mbias2; Nof gammas in event; "
2829 "invariant mass in GeV/c^{2};#",
2830 400, 0, 30, 1000, 0, 2.0);
2832 MultiplicityGamma_one_Target_mbias2 = new TH2D("MultiplicityGamma_one_Target_mbias2",
2833 "MultiplicityGamma_one_Target_mbias2; Nof gammas in event; "
2834 "invariant mass in GeV/c^{2};#",
2835 400, 0, 30, 1000, 0, 2.0);
2837 MultiplicityGamma_two_Target_mbias2 = new TH2D("MultiplicityGamma_two_Target_mbias2",
2838 "MultiplicityGamma_two_Target_mbias2; Nof gammas in event; "
2839 "invariant mass in GeV/c^{2};#",
2840 400, 0, 30, 1000, 0, 2.0);
2842 MultiplicityGamma_onetwo_Target_mbias2 = new TH2D("MultiplicityGamma_onetwo_Target_mbias2",
2843 "MultiplicityGamma_onetwo_Target_mbias2; Nof gammas in event; "
2844 "invariant mass in GeV/c^{2};#",
2845 400, 0, 30, 1000, 0, 2.0);
2847
2849 new TH2D("MultiplicityChargedParticles_all_Target_mbias2",
2850 "MultiplicityChargedParticles_all_Target_mbias2; Nof charged "
2851 "particles in event; invariant mass in GeV/c^{2};#",
2852 1000, 0, 1000, 1000, 0, 2.0);
2855 new TH2D("MultiplicityChargedParticles_zero_Target_mbias2",
2856 "MultiplicityChargedParticles_zero_Target_mbias2; Nof charged "
2857 "particles in event; invariant mass in GeV/c^{2};#",
2858 1000, 0, 1000, 1000, 0, 2.0);
2861 new TH2D("MultiplicityChargedParticles_one_Target_mbias2",
2862 "MultiplicityChargedParticles_one_Target_mbias2; Nof charged "
2863 "particles in event; invariant mass in GeV/c^{2};#",
2864 1000, 0, 1000, 1000, 0, 2.0);
2867 new TH2D("MultiplicityChargedParticles_two_Target_mbias2",
2868 "MultiplicityChargedParticles_two_Target_mbias2; Nof charged "
2869 "particles in event; invariant mass in GeV/c^{2};#",
2870 1000, 0, 1000, 1000, 0, 2.0);
2873 new TH2D("MultiplicityChargedParticles_onetwo_Target_mbias2",
2874 "MultiplicityChargedParticles_onetwo_Target_mbias2; Nof charged "
2875 "particles in event; invariant mass in GeV/c^{2};#",
2876 1000, 0, 1000, 1000, 0, 2.0);
2878
2879 // Multiplicity Outside_mbias2
2880 MultiplicityGamma_all_Outside_mbias2 = new TH2D("MultiplicityGamma_all_Outside_mbias2",
2881 "MultiplicityGamma_all_Outside_mbias2; Nof gammas in event; "
2882 "invariant mass in GeV/c^{2};#",
2883 400, 0, 30, 1000, 0, 2.0);
2885 MultiplicityGamma_zero_Outside_mbias2 = new TH2D("MultiplicityGamma_zero_Outside_mbias2",
2886 "MultiplicityGamma_zero_Outside_mbias2; Nof gammas in event; "
2887 "invariant mass in GeV/c^{2};#",
2888 400, 0, 30, 1000, 0, 2.0);
2890 MultiplicityGamma_one_Outside_mbias2 = new TH2D("MultiplicityGamma_one_Outside_mbias2",
2891 "MultiplicityGamma_one_Outside_mbias2; Nof gammas in event; "
2892 "invariant mass in GeV/c^{2};#",
2893 400, 0, 30, 1000, 0, 2.0);
2895 MultiplicityGamma_two_Outside_mbias2 = new TH2D("MultiplicityGamma_two_Outside_mbias2",
2896 "MultiplicityGamma_two_Outside_mbias2; Nof gammas in event; "
2897 "invariant mass in GeV/c^{2};#",
2898 400, 0, 30, 1000, 0, 2.0);
2900 MultiplicityGamma_onetwo_Outside_mbias2 = new TH2D("MultiplicityGamma_onetwo_Outside_mbias2",
2901 "MultiplicityGamma_onetwo_Outside_mbias2; Nof gammas in event; "
2902 "invariant mass in GeV/c^{2};#",
2903 400, 0, 30, 1000, 0, 2.0);
2905
2907 new TH2D("MultiplicityChargedParticles_all_Outside_mbias2",
2908 "MultiplicityChargedParticles_all_Outside_mbias2; Nof charged "
2909 "particles in event; invariant mass in GeV/c^{2};#",
2910 1000, 0, 1000, 1000, 0, 2.0);
2913 new TH2D("MultiplicityChargedParticles_zero_Outside_mbias2",
2914 "MultiplicityChargedParticles_zero_Outside_mbias2; Nof charged "
2915 "particles in event; invariant mass in GeV/c^{2};#",
2916 1000, 0, 1000, 1000, 0, 2.0);
2919 new TH2D("MultiplicityChargedParticles_one_Outside_mbias2",
2920 "MultiplicityChargedParticles_one_Outside_mbias2; Nof charged "
2921 "particles in event; invariant mass in GeV/c^{2};#",
2922 1000, 0, 1000, 1000, 0, 2.0);
2925 new TH2D("MultiplicityChargedParticles_two_Outside_mbias2",
2926 "MultiplicityChargedParticles_two_Outside_mbias2; Nof charged "
2927 "particles in event; invariant mass in GeV/c^{2};#",
2928 1000, 0, 1000, 1000, 0, 2.0);
2931 new TH2D("MultiplicityChargedParticles_onetwo_Outside_mbias2",
2932 "MultiplicityChargedParticles_onetwo_Outside_mbias2; Nof charged "
2933 "particles in event; invariant mass in GeV/c^{2};#",
2934 1000, 0, 1000, 1000, 0, 2.0);
2936
2937
2938 // Multiplicity Both_mbias2
2939 MultiplicityGamma_all_Both_mbias2 = new TH2D("MultiplicityGamma_all_Both_mbias2",
2940 "MultiplicityGamma_all_Both_mbias2; Nof gammas in event; "
2941 "invariant mass in GeV/c^{2};#",
2942 400, 0, 30, 1000, 0, 2.0);
2944 MultiplicityGamma_zero_Both_mbias2 = new TH2D("MultiplicityGamma_zero_Both_mbias2",
2945 "MultiplicityGamma_zero_Both_mbias2; Nof gammas in event; "
2946 "invariant mass in GeV/c^{2};#",
2947 400, 0, 30, 1000, 0, 2.0);
2949 MultiplicityGamma_one_Both_mbias2 = new TH2D("MultiplicityGamma_one_Both_mbias2",
2950 "MultiplicityGamma_one_Both_mbias2; Nof gammas in event; "
2951 "invariant mass in GeV/c^{2};#",
2952 400, 0, 30, 1000, 0, 2.0);
2954 MultiplicityGamma_two_Both_mbias2 = new TH2D("MultiplicityGamma_two_Both_mbias2",
2955 "MultiplicityGamma_two_Both_mbias2; Nof gammas in event; "
2956 "invariant mass in GeV/c^{2};#",
2957 400, 0, 30, 1000, 0, 2.0);
2959 MultiplicityGamma_onetwo_Both_mbias2 = new TH2D("MultiplicityGamma_onetwo_Both_mbias2",
2960 "MultiplicityGamma_onetwo_Both_mbias2; Nof gammas in event; "
2961 "invariant mass in GeV/c^{2};#",
2962 400, 0, 30, 1000, 0, 2.0);
2964
2965 MultiplicityChargedParticles_all_Both_mbias2 = new TH2D("MultiplicityChargedParticles_all_Both_mbias2",
2966 "MultiplicityChargedParticles_all_Both_mbias2; Nof charged "
2967 "particles in event; invariant mass in GeV/c^{2};#",
2968 1000, 0, 1000, 1000, 0, 2.0);
2970 MultiplicityChargedParticles_zero_Both_mbias2 = new TH2D("MultiplicityChargedParticles_zero_Both_mbias2",
2971 "MultiplicityChargedParticles_zero_Both_mbias2; Nof charged "
2972 "particles in event; invariant mass in GeV/c^{2};#",
2973 1000, 0, 1000, 1000, 0, 2.0);
2975 MultiplicityChargedParticles_one_Both_mbias2 = new TH2D("MultiplicityChargedParticles_one_Both_mbias2",
2976 "MultiplicityChargedParticles_one_Both_mbias2; Nof charged "
2977 "particles in event; invariant mass in GeV/c^{2};#",
2978 1000, 0, 1000, 1000, 0, 2.0);
2980 MultiplicityChargedParticles_two_Both_mbias2 = new TH2D("MultiplicityChargedParticles_two_Both_mbias2",
2981 "MultiplicityChargedParticles_two_Both_mbias2; Nof charged "
2982 "particles in event; invariant mass in GeV/c^{2};#",
2983 1000, 0, 1000, 1000, 0, 2.0);
2986 new TH2D("MultiplicityChargedParticles_onetwo_Both_mbias2",
2987 "MultiplicityChargedParticles_onetwo_Both_mbias2; Nof charged "
2988 "particles in event; invariant mass in GeV/c^{2};#",
2989 1000, 0, 1000, 1000, 0, 2.0);
2991
2992
2993 // rap_vs_Pt for "OneTwo" and "Both_mbias2"
2995 rap_vs_Pt_InM_1_mbias2 = new TH1D("rap_vs_Pt_InM_1_mbias2",
2996 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
2997 "mass in GeV/c^{2};#",
2998 1000, 0, 2.0);
3000 rap_vs_Pt_InM_2_mbias2 = new TH1D("rap_vs_Pt_InM_2_mbias2",
3001 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3002 "mass in GeV/c^{2};#",
3003 1000, 0, 2.0);
3005 rap_vs_Pt_InM_3_mbias2 = new TH1D("rap_vs_Pt_InM_3_mbias2",
3006 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3007 "mass in GeV/c^{2};#",
3008 1000, 0, 2.0);
3010 rap_vs_Pt_InM_4_mbias2 = new TH1D("rap_vs_Pt_InM_4_mbias2",
3011 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3012 "mass in GeV/c^{2};#",
3013 1000, 0, 2.0);
3015 rap_vs_Pt_InM_5_mbias2 = new TH1D("rap_vs_Pt_InM_5_mbias2",
3016 "rapidity = (1.2-1.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3017 "mass in GeV/c^{2};#",
3018 1000, 0, 2.0);
3020 rap_vs_Pt_InM_6_mbias2 = new TH1D("rap_vs_Pt_InM_6_mbias2",
3021 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3022 "mass in GeV/c^{2};#",
3023 1000, 0, 2.0);
3025 rap_vs_Pt_InM_7_mbias2 = new TH1D("rap_vs_Pt_InM_7_mbias2",
3026 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3027 "mass in GeV/c^{2};#",
3028 1000, 0, 2.0);
3030 rap_vs_Pt_InM_8_mbias2 = new TH1D("rap_vs_Pt_InM_8_mbias2",
3031 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3032 "mass in GeV/c^{2};#",
3033 1000, 0, 2.0);
3035 rap_vs_Pt_InM_9_mbias2 = new TH1D("rap_vs_Pt_InM_9_mbias2",
3036 "rapidity = (1.6-2.0) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3037 "mass in GeV/c^{2};#",
3038 1000, 0, 2.0);
3040 rap_vs_Pt_InM_10_mbias2 = new TH1D("rap_vs_Pt_InM_10_mbias2",
3041 "rapidity = (1.6-2.0) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3042 "mass in GeV/c^{2};#",
3043 1000, 0, 2.0);
3045 rap_vs_Pt_InM_11_mbias2 = new TH1D("rap_vs_Pt_InM_11_mbias2",
3046 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3047 "mass in GeV/c^{2};#",
3048 1000, 0, 2.0);
3050 rap_vs_Pt_InM_12_mbias2 = new TH1D("rap_vs_Pt_InM_12_mbias2",
3051 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3052 "mass in GeV/c^{2};#",
3053 1000, 0, 2.0);
3055 rap_vs_Pt_InM_13_mbias2 = new TH1D("rap_vs_Pt_InM_13_mbias2",
3056 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3057 "mass in GeV/c^{2};#",
3058 1000, 0, 2.0);
3060 rap_vs_Pt_InM_14_mbias2 = new TH1D("rap_vs_Pt_InM_14_mbias2",
3061 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3062 "mass in GeV/c^{2};#",
3063 1000, 0, 2.0);
3065 rap_vs_Pt_InM_15_mbias2 = new TH1D("rap_vs_Pt_InM_15_mbias2",
3066 "rapidity = (2.0-2.4) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3067 "mass in GeV/c^{2};#",
3068 1000, 0, 2.0);
3070 rap_vs_Pt_InM_16_mbias2 = new TH1D("rap_vs_Pt_InM_16_mbias2",
3071 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3072 "mass in GeV/c^{2};#",
3073 1000, 0, 2.0);
3075 rap_vs_Pt_InM_17_mbias2 = new TH1D("rap_vs_Pt_InM_17_mbias2",
3076 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3077 "mass in GeV/c^{2};#",
3078 1000, 0, 2.0);
3080 rap_vs_Pt_InM_18_mbias2 = new TH1D("rap_vs_Pt_InM_18_mbias2",
3081 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3082 "mass in GeV/c^{2};#",
3083 1000, 0, 2.0);
3085 rap_vs_Pt_InM_19_mbias2 = new TH1D("rap_vs_Pt_InM_19_mbias2",
3086 "rapidity = (2.4-2.8) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3087 "mass in GeV/c^{2};#",
3088 1000, 0, 2.0);
3090 rap_vs_Pt_InM_20_mbias2 = new TH1D("rap_vs_Pt_InM_20_mbias2",
3091 "rapidity = (2.4-2.8) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3092 "mass in GeV/c^{2};#",
3093 1000, 0, 2.0);
3095 rap_vs_Pt_InM_21_mbias2 = new TH1D("rap_vs_Pt_InM_21_mbias2",
3096 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3097 "mass in GeV/c^{2};#",
3098 1000, 0, 2.0);
3100 rap_vs_Pt_InM_22_mbias2 = new TH1D("rap_vs_Pt_InM_22_mbias2",
3101 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3102 "mass in GeV/c^{2};#",
3103 1000, 0, 2.0);
3105 rap_vs_Pt_InM_23_mbias2 = new TH1D("rap_vs_Pt_InM_23_mbias2",
3106 "rapidity = (2.8-3.2) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3107 "mass in GeV/c^{2};#",
3108 1000, 0, 2.0);
3110 rap_vs_Pt_InM_24_mbias2 = new TH1D("rap_vs_Pt_InM_24_mbias2",
3111 "rapidity = (2.8-3.2) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3112 "mass in GeV/c^{2};#",
3113 1000, 0, 2.0);
3115 rap_vs_Pt_InM_25_mbias2 = new TH1D("rap_vs_Pt_InM_25_mbias2",
3116 "rapidity = (2.8-3.2) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3117 "mass in GeV/c^{2};#",
3118 1000, 0, 2.0);
3120 rap_vs_Pt_InM_26_mbias2 = new TH1D("rap_vs_Pt_InM_26_mbias2",
3121 "rapidity = (3.2-3.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3122 "mass in GeV/c^{2};#",
3123 1000, 0, 2.0);
3125 rap_vs_Pt_InM_27_mbias2 = new TH1D("rap_vs_Pt_InM_27_mbias2",
3126 "rapidity = (3.2-3.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3127 "mass in GeV/c^{2};#",
3128 1000, 0, 2.0);
3130 rap_vs_Pt_InM_28_mbias2 = new TH1D("rap_vs_Pt_InM_28_mbias2",
3131 "rapidity = (3.2-3.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3132 "mass in GeV/c^{2};#",
3133 1000, 0, 2.0);
3135 rap_vs_Pt_InM_29_mbias2 = new TH1D("rap_vs_Pt_InM_29_mbias2",
3136 "rapidity = (3.2-3.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3137 "mass in GeV/c^{2};#",
3138 1000, 0, 2.0);
3140 rap_vs_Pt_InM_30_mbias2 = new TH1D("rap_vs_Pt_InM_30_mbias2",
3141 "rapidity = (3.2-3.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3142 "mass in GeV/c^{2};#",
3143 1000, 0, 2.0);
3145
3146 rap_vs_Pt_InM_mixing_1_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_1_mbias2",
3147 "rapidity = (1.2-1.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3148 "mass in GeV/c^{2};#",
3149 1000, 0, 2.0);
3151 rap_vs_Pt_InM_mixing_2_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_2_mbias2",
3152 "rapidity = (1.2-1.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3153 "mass in GeV/c^{2};#",
3154 1000, 0, 2.0);
3156 rap_vs_Pt_InM_mixing_3_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_3_mbias2",
3157 "rapidity = (1.2-1.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3158 "mass in GeV/c^{2};#",
3159 1000, 0, 2.0);
3161 rap_vs_Pt_InM_mixing_4_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_4_mbias2",
3162 "rapidity = (1.2-1.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3163 "mass in GeV/c^{2};#",
3164 1000, 0, 2.0);
3166 rap_vs_Pt_InM_mixing_5_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_5_mbias2",
3167 "rapidity = (1.2-1.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3168 "mass in GeV/c^{2};#",
3169 1000, 0, 2.0);
3171 rap_vs_Pt_InM_mixing_6_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_6_mbias2",
3172 "rapidity = (1.6-2.0) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3173 "mass in GeV/c^{2};#",
3174 1000, 0, 2.0);
3176 rap_vs_Pt_InM_mixing_7_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_7_mbias2",
3177 "rapidity = (1.6-2.0) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3178 "mass in GeV/c^{2};#",
3179 1000, 0, 2.0);
3181 rap_vs_Pt_InM_mixing_8_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_8_mbias2",
3182 "rapidity = (1.6-2.0) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3183 "mass in GeV/c^{2};#",
3184 1000, 0, 2.0);
3186 rap_vs_Pt_InM_mixing_9_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_9_mbias2",
3187 "rapidity = (1.6-2.0) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3188 "mass in GeV/c^{2};#",
3189 1000, 0, 2.0);
3191 rap_vs_Pt_InM_mixing_10_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_10_mbias2",
3192 "rapidity = (1.6-2.0) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3193 "mass in GeV/c^{2};#",
3194 1000, 0, 2.0);
3196 rap_vs_Pt_InM_mixing_11_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_11_mbias2",
3197 "rapidity = (2.0-2.4) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3198 "mass in GeV/c^{2};#",
3199 1000, 0, 2.0);
3201 rap_vs_Pt_InM_mixing_12_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_12_mbias2",
3202 "rapidity = (2.0-2.4) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3203 "mass in GeV/c^{2};#",
3204 1000, 0, 2.0);
3206 rap_vs_Pt_InM_mixing_13_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_13_mbias2",
3207 "rapidity = (2.0-2.4) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3208 "mass in GeV/c^{2};#",
3209 1000, 0, 2.0);
3211 rap_vs_Pt_InM_mixing_14_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_14_mbias2",
3212 "rapidity = (2.0-2.4) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3213 "mass in GeV/c^{2};#",
3214 1000, 0, 2.0);
3216 rap_vs_Pt_InM_mixing_15_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_15_mbias2",
3217 "rapidity = (2.0-2.4) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3218 "mass in GeV/c^{2};#",
3219 1000, 0, 2.0);
3221 rap_vs_Pt_InM_mixing_16_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_16_mbias2",
3222 "rapidity = (2.4-2.8) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3223 "mass in GeV/c^{2};#",
3224 1000, 0, 2.0);
3226 rap_vs_Pt_InM_mixing_17_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_17_mbias2",
3227 "rapidity = (2.4-2.8) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3228 "mass in GeV/c^{2};#",
3229 1000, 0, 2.0);
3231 rap_vs_Pt_InM_mixing_18_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_18_mbias2",
3232 "rapidity = (2.4-2.8) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3233 "mass in GeV/c^{2};#",
3234 1000, 0, 2.0);
3236 rap_vs_Pt_InM_mixing_19_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_19_mbias2",
3237 "rapidity = (2.4-2.8) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3238 "mass in GeV/c^{2};#",
3239 1000, 0, 2.0);
3241 rap_vs_Pt_InM_mixing_20_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_20_mbias2",
3242 "rapidity = (2.4-2.8) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3243 "mass in GeV/c^{2};#",
3244 1000, 0, 2.0);
3246 rap_vs_Pt_InM_mixing_21_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_21_mbias2",
3247 "rapidity = (2.8-3.2) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3248 "mass in GeV/c^{2};#",
3249 1000, 0, 2.0);
3251 rap_vs_Pt_InM_mixing_22_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_22_mbias2",
3252 "rapidity = (2.8-3.2) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3253 "mass in GeV/c^{2};#",
3254 1000, 0, 2.0);
3256 rap_vs_Pt_InM_mixing_23_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_23_mbias2",
3257 "rapidity = (2.8-3.2) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3258 "mass in GeV/c^{2};#",
3259 1000, 0, 2.0);
3261 rap_vs_Pt_InM_mixing_24_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_24_mbias2",
3262 "rapidity = (2.8-3.2) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3263 "mass in GeV/c^{2};#",
3264 1000, 0, 2.0);
3266 rap_vs_Pt_InM_mixing_25_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_25_mbias2",
3267 "rapidity = (2.8-3.2) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3268 "mass in GeV/c^{2};#",
3269 1000, 0, 2.0);
3271 rap_vs_Pt_InM_mixing_26_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_26_mbias2",
3272 "rapidity = (3.2-3.6) P_{t} = (0.0-0.4 GeV/c^{2}) ;invariant "
3273 "mass in GeV/c^{2};#",
3274 1000, 0, 2.0);
3276 rap_vs_Pt_InM_mixing_27_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_27_mbias2",
3277 "rapidity = (3.2-3.6) P_{t} = (0.4-0.8 GeV/c^{2}) ;invariant "
3278 "mass in GeV/c^{2};#",
3279 1000, 0, 2.0);
3281 rap_vs_Pt_InM_mixing_28_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_28_mbias2",
3282 "rapidity = (3.2-3.6) P_{t} = (0.8-1.2 GeV/c^{2}) ;invariant "
3283 "mass in GeV/c^{2};#",
3284 1000, 0, 2.0);
3286 rap_vs_Pt_InM_mixing_29_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_29_mbias2",
3287 "rapidity = (3.2-3.6) P_{t} = (1.2-1.6 GeV/c^{2}) ;invariant "
3288 "mass in GeV/c^{2};#",
3289 1000, 0, 2.0);
3291 rap_vs_Pt_InM_mixing_30_mbias2 = new TH1D("rap_vs_Pt_InM_mixing_30_mbias2",
3292 "rapidity = (3.2-3.6) P_{t} = (1.6-2.0 GeV/c^{2}) ;invariant "
3293 "mass in GeV/c^{2};#",
3294 1000, 0, 2.0);
3296
3297
3298 // BG cases
3299 //Both_mbias2 all
3301 new TH1D("BG1_InM_all_Both_mbias2", "BG1_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3304 new TH1D("BG2_InM_all_Both_mbias2", "BG2_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3307 new TH1D("BG3_InM_all_Both_mbias2", "BG3_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3310 new TH1D("BG4_InM_all_Both_mbias2", "BG4_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3313 new TH1D("BG5_InM_all_Both_mbias2", "BG5_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3316 new TH1D("BG6_InM_all_Both_mbias2", "BG6_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3319 new TH1D("BG7_InM_all_Both_mbias2", "BG7_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3322 new TH1D("BG8_InM_all_Both_mbias2", "BG8_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3325 new TH1D("BG9_InM_all_Both_mbias2", "BG9_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3328 new TH1D("BG10_InM_all_Both_mbias2", "BG10_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3331 new TH1D("PdgCase8_InM_all_Both_mbias2", "PdgCase8_InM_all_Both_mbias2; Id ;#", 5000, -2500, 2500);
3334 new TH1D("PdgCase8mothers_InM_all_Both_mbias2", "PdgCase8mothers_InM_all_Both_mbias2; Id ;#", 5000, -2500, 2500);
3337 new TH1D("sameMIDcase8_InM_all_Both_mbias2", "sameMIDcase8_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#",
3338 1000, 0, 2.0);
3341 new TH1D("sameGRIDcase8_InM_all_Both_mbias2", "sameGRIDcase8_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#",
3342 1000, 0, 2.0);
3344 Case1ZYPos_InM_all_Both_mbias2 = new TH2D("Case1ZYPos_InM_all_Both_mbias2",
3345 "Case1ZYPos_InM_all_Both_mbias2; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3348 new TH1D("sameMIDcase8_mothedPDG_InM_all_Both_mbias2", "sameMIDcase8_mothedPDG_InM_all_Both_mbias2; Id ;#", 5000,
3349 -2500, 2500);
3352 new TH1D("PdgCase8NonEComeFromTarget_mbias2_InM_all_Both_mbias2",
3353 "PdgCase8NonEComeFromTarget_mbias2_InM_all_Both_mbias2; Id ;#", 5000, -2500, 2500);
3356 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias2_InM_all_Both_mbias2",
3357 "PdgCase8NonE_NOT_FromTarget_mbias2_InM_all_Both_mbias2; Id ;#", 5000, -2500, 2500);
3359 PdgCase8motherNonE_InM_all_Both_mbias2 = new TH1D("PdgCase8motherNonE_InM_all_Both_mbias2",
3360 "PdgCase8motherNonE_InM_all_Both_mbias2; Id ;#", 5000, -2500, 2500);
3363 new TH1D("Case8ElFromDalitz_InM_all_Both_mbias2",
3364 "Case8ElFromDalitz_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3367 new TH1D("Case8NonElFrom_pn_InM_all_Both_mbias2",
3368 "Case8NonElFrom_pn_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3371 new TH1D("Case8NonElFrom_eta_InM_all_Both_mbias2",
3372 "Case8NonElFrom_eta_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3375 new TH1D("Case8NonElFrom_kaon_InM_all_Both_mbias2",
3376 "Case8NonElFrom_kaon_InM_all_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3379 "sameMIDcase8NonEPdg_InM_all_Both_mbias2", "sameMIDcase8NonEPdg_InM_all_Both_mbias2; Id ;#", 5000, -2500, 2500);
3382 new TH1D("sameMIDcase8NonEMotherPdg_InM_all_Both_mbias2", "sameMIDcase8NonEMotherPdg_InM_all_Both_mbias2; Id ;#",
3383 5000, -2500, 2500);
3386 new TH1D("sameMIDcase8NonEMotherIM_InM_all_Both_mbias2",
3387 "sameMIDcase8NonEMotherIM_InM_all_Both_mbias2; invariant mass in "
3388 "GeV/c^{2};#",
3389 1000, 0, 2.0);
3392 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias2_InM_all_Both_mbias2",
3393 "sameMIDcase8NonEPdgFromTarget_mbias2_InM_all_Both_mbias2; Id ;#", 5000, -2500, 2500);
3396 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2IM_InM_all_Both_mbias2",
3397 "sameMIDcase8NonEComeFromTarget_mbias2IM_InM_all_Both_mbias2; "
3398 "invariant mass in GeV/c^{2};#",
3399 1000, 0, 2.0);
3402 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2P_InM_all_Both_mbias2",
3403 "sameMIDcase8NonEComeFromTarget_mbias2P_InM_all_Both_mbias2; P in "
3404 "GeV/c^{2} ;#",
3405 200, 0, 10);
3408 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_all_Both_mbias2",
3409 "sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_all_Both_mbias2; Pt "
3410 "in GeV/c^{2} ;#",
3411 200, 0, 10);
3413 //Both_mbias2 zero
3415 new TH1D("BG1_InM_zero_Both_mbias2", "BG1_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3418 new TH1D("BG2_InM_zero_Both_mbias2", "BG2_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3421 new TH1D("BG3_InM_zero_Both_mbias2", "BG3_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3424 new TH1D("BG4_InM_zero_Both_mbias2", "BG4_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3427 new TH1D("BG5_InM_zero_Both_mbias2", "BG5_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3430 new TH1D("BG6_InM_zero_Both_mbias2", "BG6_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3433 new TH1D("BG7_InM_zero_Both_mbias2", "BG7_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3436 new TH1D("BG8_InM_zero_Both_mbias2", "BG8_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3439 new TH1D("BG9_InM_zero_Both_mbias2", "BG9_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3442 new TH1D("BG10_InM_zero_Both_mbias2", "BG10_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3445 new TH1D("PdgCase8_InM_zero_Both_mbias2", "PdgCase8_InM_zero_Both_mbias2; Id ;#", 5000, -2500, 2500);
3448 new TH1D("PdgCase8mothers_InM_zero_Both_mbias2", "PdgCase8mothers_InM_zero_Both_mbias2; Id ;#", 5000, -2500, 2500);
3451 new TH1D("sameMIDcase8_InM_zero_Both_mbias2", "sameMIDcase8_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#",
3452 1000, 0, 2.0);
3455 new TH1D("sameGRIDcase8_InM_zero_Both_mbias2", "sameGRIDcase8_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#",
3456 1000, 0, 2.0);
3459 "Case1ZYPos_InM_zero_Both_mbias2", "Case1ZYPos_InM_zero_Both_mbias2; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3462 new TH1D("sameMIDcase8_mothedPDG_InM_zero_Both_mbias2", "sameMIDcase8_mothedPDG_InM_zero_Both_mbias2; Id ;#", 5000,
3463 -2500, 2500);
3466 new TH1D("PdgCase8NonEComeFromTarget_mbias2_InM_zero_Both_mbias2",
3467 "PdgCase8NonEComeFromTarget_mbias2_InM_zero_Both_mbias2; Id ;#", 5000, -2500, 2500);
3470 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias2_InM_zero_Both_mbias2",
3471 "PdgCase8NonE_NOT_FromTarget_mbias2_InM_zero_Both_mbias2; Id ;#", 5000, -2500, 2500);
3474 "PdgCase8motherNonE_InM_zero_Both_mbias2", "PdgCase8motherNonE_InM_zero_Both_mbias2; Id ;#", 5000, -2500, 2500);
3477 new TH1D("Case8ElFromDalitz_InM_zero_Both_mbias2",
3478 "Case8ElFromDalitz_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3481 new TH1D("Case8NonElFrom_pn_InM_zero_Both_mbias2",
3482 "Case8NonElFrom_pn_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3485 new TH1D("Case8NonElFrom_eta_InM_zero_Both_mbias2",
3486 "Case8NonElFrom_eta_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3489 new TH1D("Case8NonElFrom_kaon_InM_zero_Both_mbias2",
3490 "Case8NonElFrom_kaon_InM_zero_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3493 "sameMIDcase8NonEPdg_InM_zero_Both_mbias2", "sameMIDcase8NonEPdg_InM_zero_Both_mbias2; Id ;#", 5000, -2500, 2500);
3496 new TH1D("sameMIDcase8NonEMotherPdg_InM_zero_Both_mbias2", "sameMIDcase8NonEMotherPdg_InM_zero_Both_mbias2; Id ;#",
3497 5000, -2500, 2500);
3500 new TH1D("sameMIDcase8NonEMotherIM_InM_zero_Both_mbias2",
3501 "sameMIDcase8NonEMotherIM_InM_zero_Both_mbias2; invariant mass in "
3502 "GeV/c^{2};#",
3503 1000, 0, 2.0);
3506 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias2_InM_zero_Both_mbias2",
3507 "sameMIDcase8NonEPdgFromTarget_mbias2_InM_zero_Both_mbias2; Id ;#", 5000, -2500, 2500);
3510 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2IM_InM_zero_Both_mbias2",
3511 "sameMIDcase8NonEComeFromTarget_mbias2IM_InM_zero_Both_mbias2; "
3512 "invariant mass in GeV/c^{2};#",
3513 1000, 0, 2.0);
3516 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2P_InM_zero_Both_mbias2",
3517 "sameMIDcase8NonEComeFromTarget_mbias2P_InM_zero_Both_mbias2; P "
3518 "in GeV/c^{2} ;#",
3519 200, 0, 10);
3522 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_zero_Both_mbias2",
3523 "sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_zero_Both_mbias2; Pt "
3524 "in GeV/c^{2} ;#",
3525 200, 0, 10);
3527 //Both_mbias2 one
3529 new TH1D("BG1_InM_one_Both_mbias2", "BG1_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3532 new TH1D("BG2_InM_one_Both_mbias2", "BG2_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3535 new TH1D("BG3_InM_one_Both_mbias2", "BG3_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3538 new TH1D("BG4_InM_one_Both_mbias2", "BG4_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3541 new TH1D("BG5_InM_one_Both_mbias2", "BG5_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3544 new TH1D("BG6_InM_one_Both_mbias2", "BG6_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3547 new TH1D("BG7_InM_one_Both_mbias2", "BG7_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3550 new TH1D("BG8_InM_one_Both_mbias2", "BG8_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3553 new TH1D("BG9_InM_one_Both_mbias2", "BG9_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3556 new TH1D("BG10_InM_one_Both_mbias2", "BG10_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3559 new TH1D("PdgCase8_InM_one_Both_mbias2", "PdgCase8_InM_one_Both_mbias2; Id ;#", 5000, -2500, 2500);
3562 new TH1D("PdgCase8mothers_InM_one_Both_mbias2", "PdgCase8mothers_InM_one_Both_mbias2; Id ;#", 5000, -2500, 2500);
3565 new TH1D("sameMIDcase8_InM_one_Both_mbias2", "sameMIDcase8_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#",
3566 1000, 0, 2.0);
3569 new TH1D("sameGRIDcase8_InM_one_Both_mbias2", "sameGRIDcase8_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#",
3570 1000, 0, 2.0);
3572 Case1ZYPos_InM_one_Both_mbias2 = new TH2D("Case1ZYPos_InM_one_Both_mbias2",
3573 "Case1ZYPos_InM_one_Both_mbias2; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3576 new TH1D("sameMIDcase8_mothedPDG_InM_one_Both_mbias2", "sameMIDcase8_mothedPDG_InM_one_Both_mbias2; Id ;#", 5000,
3577 -2500, 2500);
3580 new TH1D("PdgCase8NonEComeFromTarget_mbias2_InM_one_Both_mbias2",
3581 "PdgCase8NonEComeFromTarget_mbias2_InM_one_Both_mbias2; Id ;#", 5000, -2500, 2500);
3584 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias2_InM_one_Both_mbias2",
3585 "PdgCase8NonE_NOT_FromTarget_mbias2_InM_one_Both_mbias2; Id ;#", 5000, -2500, 2500);
3587 PdgCase8motherNonE_InM_one_Both_mbias2 = new TH1D("PdgCase8motherNonE_InM_one_Both_mbias2",
3588 "PdgCase8motherNonE_InM_one_Both_mbias2; Id ;#", 5000, -2500, 2500);
3591 new TH1D("Case8ElFromDalitz_InM_one_Both_mbias2",
3592 "Case8ElFromDalitz_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3595 new TH1D("Case8NonElFrom_pn_InM_one_Both_mbias2",
3596 "Case8NonElFrom_pn_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3599 new TH1D("Case8NonElFrom_eta_InM_one_Both_mbias2",
3600 "Case8NonElFrom_eta_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3603 new TH1D("Case8NonElFrom_kaon_InM_one_Both_mbias2",
3604 "Case8NonElFrom_kaon_InM_one_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3607 "sameMIDcase8NonEPdg_InM_one_Both_mbias2", "sameMIDcase8NonEPdg_InM_one_Both_mbias2; Id ;#", 5000, -2500, 2500);
3610 new TH1D("sameMIDcase8NonEMotherPdg_InM_one_Both_mbias2", "sameMIDcase8NonEMotherPdg_InM_one_Both_mbias2; Id ;#",
3611 5000, -2500, 2500);
3614 new TH1D("sameMIDcase8NonEMotherIM_InM_one_Both_mbias2",
3615 "sameMIDcase8NonEMotherIM_InM_one_Both_mbias2; invariant mass in "
3616 "GeV/c^{2};#",
3617 1000, 0, 2.0);
3620 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias2_InM_one_Both_mbias2",
3621 "sameMIDcase8NonEPdgFromTarget_mbias2_InM_one_Both_mbias2; Id ;#", 5000, -2500, 2500);
3624 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2IM_InM_one_Both_mbias2",
3625 "sameMIDcase8NonEComeFromTarget_mbias2IM_InM_one_Both_mbias2; "
3626 "invariant mass in GeV/c^{2};#",
3627 1000, 0, 2.0);
3630 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2P_InM_one_Both_mbias2",
3631 "sameMIDcase8NonEComeFromTarget_mbias2P_InM_one_Both_mbias2; P in "
3632 "GeV/c^{2} ;#",
3633 200, 0, 10);
3636 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_one_Both_mbias2",
3637 "sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_one_Both_mbias2; Pt "
3638 "in GeV/c^{2} ;#",
3639 200, 0, 10);
3641 //Both_mbias2 two
3643 new TH1D("BG1_InM_two_Both_mbias2", "BG1_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3646 new TH1D("BG2_InM_two_Both_mbias2", "BG2_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3649 new TH1D("BG3_InM_two_Both_mbias2", "BG3_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3652 new TH1D("BG4_InM_two_Both_mbias2", "BG4_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3655 new TH1D("BG5_InM_two_Both_mbias2", "BG5_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3658 new TH1D("BG6_InM_two_Both_mbias2", "BG6_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3661 new TH1D("BG7_InM_two_Both_mbias2", "BG7_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3664 new TH1D("BG8_InM_two_Both_mbias2", "BG8_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3667 new TH1D("BG9_InM_two_Both_mbias2", "BG9_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3670 new TH1D("BG10_InM_two_Both_mbias2", "BG10_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3673 new TH1D("PdgCase8_InM_two_Both_mbias2", "PdgCase8_InM_two_Both_mbias2; Id ;#", 5000, -2500, 2500);
3676 new TH1D("PdgCase8mothers_InM_two_Both_mbias2", "PdgCase8mothers_InM_two_Both_mbias2; Id ;#", 5000, -2500, 2500);
3679 new TH1D("sameMIDcase8_InM_two_Both_mbias2", "sameMIDcase8_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#",
3680 1000, 0, 2.0);
3683 new TH1D("sameGRIDcase8_InM_two_Both_mbias2", "sameGRIDcase8_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#",
3684 1000, 0, 2.0);
3686 Case1ZYPos_InM_two_Both_mbias2 = new TH2D("Case1ZYPos_InM_two_Both_mbias2",
3687 "Case1ZYPos_InM_two_Both_mbias2; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3690 new TH1D("sameMIDcase8_mothedPDG_InM_two_Both_mbias2", "sameMIDcase8_mothedPDG_InM_two_Both_mbias2; Id ;#", 5000,
3691 -2500, 2500);
3694 new TH1D("PdgCase8NonEComeFromTarget_mbias2_InM_two_Both_mbias2",
3695 "PdgCase8NonEComeFromTarget_mbias2_InM_two_Both_mbias2; Id ;#", 5000, -2500, 2500);
3698 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias2_InM_two_Both_mbias2",
3699 "PdgCase8NonE_NOT_FromTarget_mbias2_InM_two_Both_mbias2; Id ;#", 5000, -2500, 2500);
3701 PdgCase8motherNonE_InM_two_Both_mbias2 = new TH1D("PdgCase8motherNonE_InM_two_Both_mbias2",
3702 "PdgCase8motherNonE_InM_two_Both_mbias2; Id ;#", 5000, -2500, 2500);
3705 new TH1D("Case8ElFromDalitz_InM_two_Both_mbias2",
3706 "Case8ElFromDalitz_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3709 new TH1D("Case8NonElFrom_pn_InM_two_Both_mbias2",
3710 "Case8NonElFrom_pn_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3713 new TH1D("Case8NonElFrom_eta_InM_two_Both_mbias2",
3714 "Case8NonElFrom_eta_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3717 new TH1D("Case8NonElFrom_kaon_InM_two_Both_mbias2",
3718 "Case8NonElFrom_kaon_InM_two_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3721 "sameMIDcase8NonEPdg_InM_two_Both_mbias2", "sameMIDcase8NonEPdg_InM_two_Both_mbias2; Id ;#", 5000, -2500, 2500);
3724 new TH1D("sameMIDcase8NonEMotherPdg_InM_two_Both_mbias2", "sameMIDcase8NonEMotherPdg_InM_two_Both_mbias2; Id ;#",
3725 5000, -2500, 2500);
3728 new TH1D("sameMIDcase8NonEMotherIM_InM_two_Both_mbias2",
3729 "sameMIDcase8NonEMotherIM_InM_two_Both_mbias2; invariant mass in "
3730 "GeV/c^{2};#",
3731 1000, 0, 2.0);
3734 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias2_InM_two_Both_mbias2",
3735 "sameMIDcase8NonEPdgFromTarget_mbias2_InM_two_Both_mbias2; Id ;#", 5000, -2500, 2500);
3738 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2IM_InM_two_Both_mbias2",
3739 "sameMIDcase8NonEComeFromTarget_mbias2IM_InM_two_Both_mbias2; "
3740 "invariant mass in GeV/c^{2};#",
3741 1000, 0, 2.0);
3744 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2P_InM_two_Both_mbias2",
3745 "sameMIDcase8NonEComeFromTarget_mbias2P_InM_two_Both_mbias2; P in "
3746 "GeV/c^{2} ;#",
3747 200, 0, 10);
3750 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_two_Both_mbias2",
3751 "sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_two_Both_mbias2; Pt "
3752 "in GeV/c^{2} ;#",
3753 200, 0, 10);
3755 //Both_mbias2 onetwo
3757 new TH1D("BG1_InM_onetwo_Both_mbias2", "BG1_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3760 new TH1D("BG2_InM_onetwo_Both_mbias2", "BG2_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3763 new TH1D("BG3_InM_onetwo_Both_mbias2", "BG3_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3766 new TH1D("BG4_InM_onetwo_Both_mbias2", "BG4_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3769 new TH1D("BG5_InM_onetwo_Both_mbias2", "BG5_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3772 new TH1D("BG6_InM_onetwo_Both_mbias2", "BG6_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3775 new TH1D("BG7_InM_onetwo_Both_mbias2", "BG7_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3778 new TH1D("BG8_InM_onetwo_Both_mbias2", "BG8_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3781 new TH1D("BG9_InM_onetwo_Both_mbias2", "BG9_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3784 new TH1D("BG10_InM_onetwo_Both_mbias2", "BG10_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3787 new TH1D("PdgCase8_InM_onetwo_Both_mbias2", "PdgCase8_InM_onetwo_Both_mbias2; Id ;#", 5000, -2500, 2500);
3789 PdgCase8mothers_InM_onetwo_Both_mbias2 = new TH1D("PdgCase8mothers_InM_onetwo_Both_mbias2",
3790 "PdgCase8mothers_InM_onetwo_Both_mbias2; Id ;#", 5000, -2500, 2500);
3793 new TH1D("sameMIDcase8_InM_onetwo_Both_mbias2",
3794 "sameMIDcase8_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3797 new TH1D("sameGRIDcase8_InM_onetwo_Both_mbias2",
3798 "sameGRIDcase8_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3801 "Case1ZYPos_InM_onetwo_Both_mbias2", "Case1ZYPos_InM_onetwo_Both_mbias2; z[cm]; y[cm]", 400, -1, 200, 200, -50, 50);
3804 new TH1D("sameMIDcase8_mothedPDG_InM_onetwo_Both_mbias2", "sameMIDcase8_mothedPDG_InM_onetwo_Both_mbias2; Id ;#",
3805 5000, -2500, 2500);
3808 new TH1D("PdgCase8NonEComeFromTarget_mbias2_InM_onetwo_Both_mbias2",
3809 "PdgCase8NonEComeFromTarget_mbias2_InM_onetwo_Both_mbias2; Id ;#", 5000, -2500, 2500);
3812 new TH1D("PdgCase8NonE_NOT_FromTarget_mbias2_InM_onetwo_Both_mbias2",
3813 "PdgCase8NonE_NOT_FromTarget_mbias2_InM_onetwo_Both_mbias2; Id ;#", 5000, -2500, 2500);
3816 "PdgCase8motherNonE_InM_onetwo_Both_mbias2", "PdgCase8motherNonE_InM_onetwo_Both_mbias2; Id ;#", 5000, -2500, 2500);
3819 new TH1D("Case8ElFromDalitz_InM_onetwo_Both_mbias2",
3820 "Case8ElFromDalitz_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3823 new TH1D("Case8NonElFrom_pn_InM_onetwo_Both_mbias2",
3824 "Case8NonElFrom_pn_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3827 new TH1D("Case8NonElFrom_eta_InM_onetwo_Both_mbias2",
3828 "Case8NonElFrom_eta_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3831 new TH1D("Case8NonElFrom_kaon_InM_onetwo_Both_mbias2",
3832 "Case8NonElFrom_kaon_InM_onetwo_Both_mbias2; invariant mass in GeV/c^{2};#", 1000, 0, 2.0);
3835 new TH1D("sameMIDcase8NonEPdg_InM_onetwo_Both_mbias2", "sameMIDcase8NonEPdg_InM_onetwo_Both_mbias2; Id ;#", 5000,
3836 -2500, 2500);
3839 new TH1D("sameMIDcase8NonEMotherPdg_InM_onetwo_Both_mbias2",
3840 "sameMIDcase8NonEMotherPdg_InM_onetwo_Both_mbias2; Id ;#", 5000, -2500, 2500);
3843 new TH1D("sameMIDcase8NonEMotherIM_InM_onetwo_Both_mbias2",
3844 "sameMIDcase8NonEMotherIM_InM_onetwo_Both_mbias2; invariant mass "
3845 "in GeV/c^{2};#",
3846 1000, 0, 2.0);
3849 new TH1D("sameMIDcase8NonEPdgFromTarget_mbias2_InM_onetwo_Both_mbias2",
3850 "sameMIDcase8NonEPdgFromTarget_mbias2_InM_onetwo_Both_mbias2; Id ;#", 5000, -2500, 2500);
3853 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2IM_InM_onetwo_Both_mbias2",
3854 "sameMIDcase8NonEComeFromTarget_mbias2IM_InM_onetwo_Both_mbias2; "
3855 "invariant mass in GeV/c^{2};#",
3856 1000, 0, 2.0);
3859 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2P_InM_onetwo_Both_mbias2",
3860 "sameMIDcase8NonEComeFromTarget_mbias2P_InM_onetwo_Both_mbias2; P "
3861 "in GeV/c^{2} ;#",
3862 200, 0, 10);
3865 new TH1D("sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_onetwo_Both_mbias2",
3866 "sameMIDcase8NonEComeFromTarget_mbias2Pt_InM_onetwo_Both_mbias2; "
3867 "Pt in GeV/c^{2} ;#",
3868 200, 0, 10);
3870}
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)
std::vector< std::vector< TVector3 > > Gammas_onetwo_Both_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_one_Both_mbias2
double CalculatePlaneAngle_first(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Target_mbias2
std::vector< std::vector< TVector3 > > EMT_man_pair_momenta_Both_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Target_mbias2
std::vector< std::vector< TVector3 > > Gammas_all_Outside_mbias2
double CalculatePlaneAngle_last(CbmStsTrack *Sts_1, CbmStsTrack *Sts_2)
vector< CbmRichRing * > VRichRing_plus_Target_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_two_Outside_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_one_Outside_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Target_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_zero_Outside_mbias2
CbmRichRingFitterEllipseTau * fTauFit
std::vector< std::vector< int > > Gammas_stsIndex_all_Both_mbias2
std::vector< std::vector< TVector3 > > Gammas_zero_Outside_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Outside_mbias2
void FindGammasTarget(int EventNumMan, 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)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_one_Outside_mbias2
std::vector< std::vector< TVector3 > > Gammas_two_Both_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_zero_Both_mbias2
int FindInRich(int richInd, int stsMcTrackId)
std::vector< std::vector< TVector3 > > Gammas_two_Outside_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Both_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_zero_Target_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_all_Target_mbias2
vector< CbmRichRing * > VRichRing_minus_Target_mbias2
std::vector< std::vector< TVector3 > > EMT_man_pair_momenta_Target_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Both_mbias2
void FindGammasOutside(int EventNumMan, double AngleCut, double InvMassCut, int RealPID, vector< CbmMCTrack * > MCtracks_minus_Outside_mbias2, vector< CbmMCTrack * > MCtracks_plus_Outside_mbias2, vector< CbmStsTrack * > StsTrack_minus_Outside_mbias2, vector< CbmStsTrack * > StsTrack_plus_Outside_mbias2, std::vector< int > Rings_minus_Outside_mbias2, std::vector< int > Rings_plus_Outside_mbias2, std::vector< int > stsIndex_minus_Outside_mbias2, std::vector< int > stsIndex_plus_Outside_mbias2, vector< CbmRichRing * > richRing_minus_Outside_mbias2, vector< CbmRichRing * > richRing_plus_Outside_mbias2)
std::vector< std::vector< TVector3 > > Gammas_zero_Target_mbias2
void SaveOutsideTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
vector< CbmStsTrack * > VStsTrack_minus_Outside_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_all_Outside_mbias2
std::vector< std::vector< TVector3 > > EMT_man_pair_momenta_Outside_mbias2
std::vector< std::vector< TVector3 > > Gammas_onetwo_Outside_mbias2
vector< CbmStsTrack * > VStsTrack_minus_Target_mbias2
int CheckIfElectron(CbmRichRing *ring, double momentum)
vector< CbmRichRing * > VRichRing_minus_Outside_mbias2
std::vector< std::vector< TVector3 > > Gammas_all_Both_mbias2
std::vector< std::vector< TVector3 > > Gammas_one_Both_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Both_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_zero_Target_mbias2
std::vector< std::vector< TVector3 > > Gammas_zero_Both_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_two_Target_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_one_Both_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_zero_Outside_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Target_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_two_Both_mbias2
vector< CbmStsTrack * > VStsTrack_plus_Target_mbias2
std::vector< std::vector< TVector3 > > Gammas_two_Target_mbias2
vector< CbmStsTrack * > VStsTrack_plus_Outside_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_zero_Both_mbias2
vector< CbmMCTrack * > VMCtracks_minus_Outside_mbias2
void Exec(int fEventNumMan, double OpeningAngleCut, double GammaInvMassCut, int RealPID)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_all_Both_mbias2
void SaveTargetTracks(CbmMCTrack *mcTrack1, CbmStsTrack *stsTrack, TVector3 refmom, double charge, int stsInd, int richInd, int stsMcTrackId, CbmRichRing *RING)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_two_Outside_mbias2
std::vector< std::vector< TVector3 > > Gammas_onetwo_Target_mbias2
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_onetwo_Outside_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_one_Target_mbias2
std::vector< std::vector< int > > Gammas_stsIndex_onetwo_Outside_mbias2
vector< CbmRichRing * > VRichRing_plus_Outside_mbias2
std::vector< std::vector< TVector3 > > Gammas_all_Target_mbias2
std::vector< std::vector< TVector3 > > Gammas_one_Outside_mbias2
std::vector< std::vector< TVector3 > > Gammas_one_Target_mbias2
void FindPi0(TString mod, TString position, vector< vector< TVector3 > > Gammas, vector< vector< int > > StsIndex, vector< vector< CbmMCTrack * > > GammasMC, TH1D *Pi0InvMassReco, TH2D *Pi0_pt_vs_rap, TH2D *Pi0_pt_vs_rap_est, TH2D *MultiplicityGamma, TH2D *MultiplicityChargedParticles, vector< TH1 * > BGCases)
std::vector< std::vector< CbmMCTrack * > > Gammas_MC_one_Target_mbias2
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 TVector3 FitToVertex(CbmStsTrack *stsTrack, double x, double y, double z)
int32_t GetMotherId() const
Definition CbmMCTrack.h:68
int32_t GetPdgCode() const
Definition CbmMCTrack.h:67
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 fRapidity
Definition LmvmKinePar.h:19
Hash for CbmL1LinkKey.