/* * PndEventAnalysisLambdaLambdabar.cpp * * Created on: 25.06.2013 * Author: simone */ #include "PndEventAnalysisLambdaLambdabar.h" #include "TClonesArray.h" #include "FairRootManager.h" #include "FairRunAna.h" #include "FairRuntimeDb.h" #include "TVector3.h" #include "TH1F.h" #include "TH1I.h" #include "TH2F.h" #include "TStopwatch.h" #include "FairRun.h" #include "FairRuntimeDb.h" #include #include // RHO stuff #include "RhoCandidate.h" #include "RhoCandList.h" // Analysis Tools #include "PndAnalysis.h" // Fitters #include "Pnd4CFitter.h" #include "PndKinVtxFitter.h" #include "PndKinFitter.h" #include "PndVtxPRG.h" #include "PndChiVtxFitter.h" #include "PndVtxPoca.h" // Others #include "TDatabasePDG.h" PndEventAnalysisLambdaLambdabar::PndEventAnalysisLambdaLambdabar() : fPndVtxPRGFitter_switch(kTRUE),fInputBranchName(" "),fPndChiVtxFitter_switch(kTRUE),fPndKinVtxFitter_switch(kTRUE) { } PndEventAnalysisLambdaLambdabar::~PndEventAnalysisLambdaLambdabar() { } void PndEventAnalysisLambdaLambdabar::SetParContainers() { FairRun* run = FairRun::Instance(); if ( ! run ) Fatal("SetParContainers", "No analysis run"); } InitStatus PndEventAnalysisLambdaLambdabar::Init() { FairRootManager* ioman = FairRootManager::Instance(); if (!ioman) { std::cout << "-E- PndAnaLambdaLambdaBar::Init: " << "RootManager not instantiated!" << std::endl; return kFATAL; } if(fInputBranchName== " ") { std::cout << "-E- PndEventAnalysisLambdaLambdabar::Init: No InputBranchName Set" << std::endl; return kERROR; } else { theAnalysis=new PndAnalysis(fInputBranchName); } // theAnalysis=new PndAnalysis("SttMvdGemFtsGenIdealTrackSelection"); //theAnalysis=new PndAnalysis(fInputBranchName); TDatabasePDG::Instance()->AddParticle("pbarpSystem","pbarpSystem",1.9,kFALSE,0.1,0,"",88888); std::cout << "-I- PndEventAnalysisLambdaLambdabar: Initialization successful" << std::endl; timer=new TStopwatch(); timer->Start(); eventcounter=0; fErrorCounter=0; multipleLambdacounter=0; fInvMass_2D_Distribution = new TH2F("fInvMass_2D_Distribution","fInvMass_2D_Distribution",1000,1,1.5 ,1000,1,1.5 ); fInvMass_2D_Distribution_MCMatch = new TH2F("fInvMass_2D_Distribution_MCMatch","fInvMass_2D_Distribution_MCMatch",1000,1,1.5 ,1000,1,1.5 ); fAntiProton_Theta_Phi = new TH2F("fAntiProton_Theta_Phi","fAntiProton_Theta_Phi; #Phi ; #Theta", 720,-180,180,320,0,180); fAntiProton_Theta_Phi_Cut = new TH2F("fAntiProton_Theta_Phi_Cut","fAntiProton_Theta_Phi_Cut; #Phi ; #Theta", 720,-180,180,320,0,180); fAntiProton_CosTheta_Phi = new TH2F("fAntiProton_CosTheta_Phi","fAntiProton_CosTheta_Phi; #Phi ;cos #Theta", 720,-180,180,320,0,1); fAntiProton_CosTheta_Phi_Cut = new TH2F("fAntiProton_CosTheta_Phi_Cut","fAntiProton_CosTheta_Phi_Cut; #Phi ; cos #Theta", 720,-180,320,180,0,1); fAntiProton_Theta_Phi_CMS = new TH2F("fAntiProton_Theta_Phi_CMS","fAntiProton_Theta_Phi_CMS; #Phi ; #Theta", 720,-180,180,320,0,180); fAntiProton_Theta_Phi_Cut_CMS = new TH2F("fAntiProton_Theta_Phi_Cut_CMS","fAntiProton_Theta_Phi_Cut_CMS; #Phi ; #Theta", 720,-180,180,320,0,180); fAntiProton_CosTheta_Phi_CMS = new TH2F("fAntiProton_CosTheta_Phi_CMS","fAntiProton_CosTheta_Phi_CMS; #Phi ;cos #Theta", 720,-180,180,320,0,1); fAntiProton_CosTheta_Phi_Cut_CMS = new TH2F("fAntiProton_CosTheta_Phi_Cut_CMS","fAntiProton_CosTheta_Phi_Cut_CMS; #Phi ; cos #Theta", 720,-180,320,180,0,1); fDeltaAngleVsDecayLength = new TH2F("fDeltaAngleVsDecayLength","fDeltaAngleVsDecayLength",1000,0,100, 100,0,5); fReconstructedParticle_MomentaDistributionPionMinus_Lab = new TH1F("fReconstructedParticle_MomentaDistributionPionMinus_Lab" ,"fReconstructedParticle_MomentaDistributionPionMinus_Lab; Reconstructed PionMinus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionPionPlus_Lab = new TH1F("fReconstructedParticle_MomentaDistributionPionPlus_Lab" ,"fReconstructedParticle_MomentaDistributionPionPlus_Lab; Reconstructed PionPlus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionProton_Lab = new TH1F("fReconstructedParticle_MomentaDistributionProton_Lab" ,"fReconstructedParticle_MomentaDistributionProton_Lab; Reconstructed Proton Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionAntiProton_Lab= new TH1F("fReconstructedParticle_MomentaDistributionAntiProton_Lab" ,"fReconstructedParticle_MomentaDistributionAntiProton_Lab; Reconstructed AntiProton Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionPionMinus_CMS = new TH1F("fReconstructedParticle_MomentaDistributionPionMinus_CMS" ,"fReconstructedParticle_MomentaDistributionPionMinus_CMS; Reconstructed PionMinus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionPionPlus_CMS = new TH1F("fReconstructedParticle_MomentaDistributionPionPlus_CMS" ,"fReconstructedParticle_MomentaDistributionPionPlus_CMS; Reconstructed PionPlus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionProton_CMS = new TH1F("fReconstructedParticle_MomentaDistributionProton_CMS" ,"fReconstructedParticle_MomentaDistributionProton_CMS; Reconstructed Proton Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionAntiProton_CMS= new TH1F("fReconstructedParticle_MomentaDistributionAntiProton_CMS" ,"fReconstructedParticle_MomentaDistributionAntiProton_CMS; Reconstructed AntiProton Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionPionMinus_Lab_MCMatch = new TH1F("fReconstructedParticle_MomentaDistributionPionMinus_Lab_MCMatch" ,"fReconstructedParticle_MomentaDistributionPionMinus_Lab_MCMatch; Reconstructed PionMinus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionPionPlus_Lab_MCMatch = new TH1F("fReconstructedParticle_MomentaDistributionPionPlus_Lab_MCMatch" ,"fReconstructedParticle_MomentaDistributionPionPlus_Lab_MCMatch; Reconstructed PionPlus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionProton_Lab_MCMatch = new TH1F("fReconstructedParticle_MomentaDistributionProton_Lab_MCMatch" ,"fReconstructedParticle_MomentaDistributionProton_Lab_MCMatch; Reconstructed Proton Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionAntiProton_Lab_MCMatch= new TH1F("fReconstructedParticle_MomentaDistributionAntiProton_Lab_MCMatch" ,"fReconstructedParticle_MomentaDistributionAntiProton_Lab_MCMatch; Reconstructed AntiProton Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionPionMinus_CMS_MCMatch = new TH1F("fReconstructedParticle_MomentaDistributionPionMinus_CMS_MCMatch" ,"fReconstructedParticle_MomentaDistributionPionMinus_CMS_MCMatch; Reconstructed PionMinus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionPionPlus_CMS_MCMatch = new TH1F("fReconstructedParticle_MomentaDistributionPionPlus_CMS_MCMatch" ,"fReconstructedParticle_MomentaDistributionPionPlus_CMS_MCMatch; Reconstructed PionPlus Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionProton_CMS_MCMatch = new TH1F("fReconstructedParticle_MomentaDistributionProton_CMS_MCMatch" ,"fReconstructedParticle_MomentaDistributionProton_CMS_MCMatch; Reconstructed Proton Momentum [GeV/c];#Events",800,0,6); fReconstructedParticle_MomentaDistributionAntiProton_CMS_MCMatch= new TH1F("fReconstructedParticle_MomentaDistributionAntiProton_CMS_MCMatch" ,"fReconstructedParticle_MomentaDistributionAntiProton_CMS_MCMatch; Reconstructed AntiProton Momentum [GeV/c];#Events",800,0,6); fReconstructedArmenterosPodolanski_PndVtxPRGFitter = new TH2F("fReconstructedArmenterosPodolanski_PndVtxPRGFitter","Reconstructed Armenteros-Podolanski-Plot; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch= new TH2F("fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch","fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut= new TH2F("fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut","Reconstructed Armenteros-Podolanski-Plot _ProbCut; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut_MCMatch= new TH2F("fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut_MCMatch","Reconstructed Armenteros-Podolanski-Plot _ProbCut_MCMatch; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter = new TH2F("fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter","Reconstructed Armenteros-Podolanski-Plot Lambda_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_MCMatch= new TH2F("fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_MCMatch","fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch Lambda_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_ProbCut= new TH2F("fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_ProbCut","Reconstructed Armenteros-Podolanski-Plot _ProbCut L ambda_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_ProbCut_MCMatch= new TH2F("fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_ProbCut_MCMatch","Reconstructed Armenteros-Podolanski-Plot _ProbCut_MCMatch Lambda_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter = new TH2F("fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter","Reconstructed Armenteros-Podolanski-Plot LambdaBar_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_MCMatch= new TH2F("fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_MCMatch","fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch LambdaBar_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut= new TH2F("fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut","Reconstructed Armenteros-Podolanski-Plot _ProbCut LambdaBar_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut_MCMatch= new TH2F("fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut_MCMatch","Reconstructed Armenteros-Podolanski-Plot _ProbCut_MCMatch LambdaBar_; #alpha; P_t [GeV/c]",400,-1.5,1.5,400,0,0.12); fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison = new TH1F("fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison","fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison",800, -2,2); fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison= new TH1F("fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison","fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison",800, -2,2); fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison = new TH1F("fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison","fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison",800, -2,2); fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison= new TH1F("fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison","fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison",800, -2,2); fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison = new TH1F("fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison","fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison",800, -2,2); fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison= new TH1F("fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison","fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison",800, -2,2); fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison = new TH1F("fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison","fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison",800, -2,2); fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison= new TH1F("fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison","fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison",800, -2,2); fReconstructedParticle_AngularDistribution_Lambda=new TH1I("fReconstructedParticle_AngularDistribution_Lambda","fReconstructedParticle_AngularDistribution_Lambda;cos(#Theta); #Events",320,-1,1); fReconstructedParticle_AngularDistribution_LambdaBar=new TH1I("fReconstructedParticle_AngularDistribution_LambdaBar","fReconstructedParticle_AngularDistribution_LambdaBar;cos(#Theta); #Events",320,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X;cos(#Theta_{x});#Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X;cos(#Theta_{x}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_Both_Lambda=new TH1I("fReconstructedParticle_AngularDistribution_Both_Lambda","fReconstructedParticle_AngularDistribution_Both_Lambda;cos(#Theta); #Events",320,-1,1); fReconstructedParticle_AngularDistribution_Both_LambdaBar=new TH1I("fReconstructedParticle_AngularDistribution_Both_LambdaBar","fReconstructedParticle_AngularDistribution_Both_LambdaBar;cos(#Theta); #Events",320,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X;cos(#Theta_{x});#Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X;cos(#Theta_{x}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_Lambda_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistribution_Lambda_Pnd4CFitter","fReconstructedParticle_AngularDistribution_Lambda_Pnd4CFitter;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_LambdaBar_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistribution_LambdaBar_Pnd4CFitter","fReconstructedParticle_AngularDistribution_LambdaBar_Pnd4CFitter;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_Pnd4CFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_Pnd4CFitter;cos(#Theta_{x});#Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_Pnd4CFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_Pnd4CFitter;cos(#Theta_{x}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_Pnd4CFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_Pnd4CFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_Pnd4CFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_Pnd4CFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_Pnd4CFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_Pnd4CFitter;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_Pnd4CFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_Pnd4CFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_Pnd4CFitter;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_Lambda_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistribution_Lambda_PndKinVtxFitter","fReconstructedParticle_AngularDistribution_Lambda;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_LambdaBar_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistribution_LambdaBar_PndKinVtxFitter","fReconstructedParticle_AngularDistribution_LambdaBar;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndKinVtxFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndKinVtxFitter;cos(#Theta_{x});#Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndKinVtxFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndKinVtxFitter;cos(#Theta_{x}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndKinVtxFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndKinVtxFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndKinVtxFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndKinVtxFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndKinVtxFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndKinVtxFitter;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndKinVtxFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndKinVtxFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndKinVtxFitter;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_Lambda_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistribution_Lambda_PndVtxPRGFitter","fReconstructedParticle_AngularDistribution_Lambda_PndVtxPRGFitter;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_LambdaBar_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistribution_LambdaBar_PndVtxPRGFitter","fReconstructedParticle_AngularDistribution_LambdaBar_PndVtxPRGFitter;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndVtxPRGFitter;cos(#Theta_{x});#Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndVtxPRGFitter;cos(#Theta_{x}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndVtxPRGFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndVtxPRGFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndVtxPRGFitter;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndVtxPRGFitter;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_Lambda_Both_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistribution_Lambda_Both_PndVtxPRGFitter","fReconstructedParticle_AngularDistribution_Lambda_Both_PndVtxPRGFitter;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistribution_LambdaBar_Both_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistribution_LambdaBar_Both_PndVtxPRGFitter","fReconstructedParticle_AngularDistribution_LambdaBar_Both_PndVtxPRGFitter;cos(#Theta); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X_PndVtxPRGFitter;cos(#Theta_{x});#Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X_PndVtxPRGFitter;cos(#Theta_{x}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y_PndVtxPRGFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y_PndVtxPRGFitter;cos(#Theta_{y}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z_PndVtxPRGFitter;cos(#Theta_{z}); #Events",80,-1,1); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z_PndVtxPRGFitter=new TH1I("fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z_PndVtxPRGFitter","fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z_PndVtxPRGFitter;cos(#Theta_{z}); #Events",80,-1,1); fAntiProton_DecayPoint =new TH2F("fAntiProton_DecayPoint","fAntiProton_DecayPoint",300,0,300,150,0,150 ); fAntiProton_NumberOfDaughters = new TH1I("fAntiProton_NumberOfDaughters","fAntiProton_NumberOfDaughters",20,-0.5,19.5); fPnd4CFitter_Lambda_AngularDistribution = new TH1F("fPnd4CFitter_Lambda_AngularDistribution","fPnd4CFitter_Lambda_AngularDistribution;cos(#Theta)",40,-1,1 ); fPnd4CFitter_LambdaBar_AngularDistribution = new TH1F("fPnd4CFitter_LambdaBar_AngularDistribution","fPnd4CFitter_LambdaBar_AngularDistribution;cos(#Theta)",40,-1,1 ); fKalmanFitterFlag = new TH1I("fKalmanFitterFlag","fKalmanFitterFlag",10,-5,5); fComparison_Lambda_InvMass_Before_After_PndVtxPRGFitter=new TH2F("fComparison_Lambda_InvMass_Before_After_PndVtxPRGFitter","fComparison_Lambda_InvMass_Before_After_PndVtxPRGFitter",300,1,4,300,1,4); fTracksPerEvent= new TH1I("fTracksPerEvent","fTracksPerEvent", 6,-0.5,5.5); fLambda_InvMass= new TH1F("fLambda_InvMass","All p #pi- combinations ; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fLambdaBar_InvMass= new TH1F("fLambdaBar_InvMass","All #bar{p} #pi+ combinations; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fLambdasPerMCTruthMatchStage= new TH1I("fLambdasPerMCTruthMatchStage", "Lambdas Per MCTruthMatch stage", 6,-0.5,5.5); fLambdaBarsPerMCTruthMatchStage= new TH1I("fLambdaBarsPerMCTruthMatchStage", "LambdaBars Per MCTruthMatch stage", 6,-0.5,5.5); fLambda_InvMass_MCMatch_Rejected= new TH1F("fLambda_InvMass_MCMatch_Rejected","#Lambda mass (MCTruth, rejected) ; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fLambdaBar_InvMass_MCMatch_Rejected= new TH1F("fLambdaBar_InvMass_MCMatch_Rejected","#bar{#Lambda} mass (MCTruth, rejected) ; Inv. Mass [GeV/c]; #Events",1000,1,1.5); // PndVtxPoca fPndVtxPoca_Lambda_QualityNumber = new TH1F("fPndVtxPoca_Lambda_QualityNumber","fPndVtxPoca_Lambda_QualityNumber",400,0,40); fPndVtxPoca_Lambda_VertexPosition = new TH2F("fPndVtxPoca_Lambda_VertexPosition","fPndVtxPoca_Lambda_VertexPosition",100,0,100,30,0,30); fPndVtxPoca_Lambda_QualityNumber_MCMatch = new TH1F("fPndVtxPoca_Lambda_QualityNumber_MCMatch","fPndVtxPoca_Lambda_QualityNumber_MCMatch",400,0,40); fPndVtxPoca_Lambda_QualityNumber_Rejected = new TH1F("fPndVtxPoca_Lambda_QualityNumber_Rejected","fPndVtxPoca_Lambda_QualityNumber_Rejected",400,0,40); fPndVtxPoca_Lambda_VertexPosition_MCMatch = new TH2F("fPndVtxPoca_Lambda_VertexPosition_MCMatch","fPndVtxPoca_Lambda_VertexPosition_MCMatch",100,0,100,30,0,30); fPndVtxPoca_Lambda_VertexPosition_Rejected = new TH2F("fPndVtxPoca_Lambda_VertexPosition_Rejected","fPndVtxPoca_Lambda_VertexPosition_Rejected",100,0,100,30,0,30); fPndVtxPoca_LambdaBar_QualityNumber = new TH1F("fPndVtxPoca_LambdaBar_QualityNumber","fPndVtxPoca_LambdaBar_QualityNumber",400,0,40); fPndVtxPoca_LambdaBar_VertexPosition = new TH2F("fPndVtxPoca_LambdaBar_VertexPosition","fPndVtxPoca_LambdaBar_VertexPosition",100,0,100,30,0,30); fPndVtxPoca_LambdaBar_QualityNumber_MCMatch = new TH1F("fPndVtxPoca_LambdaBar_QualityNumber_MCMatch","fPndVtxPoca_LambdaBar_QualityNumber_MCMatch",400,0,40); fPndVtxPoca_LambdaBar_QualityNumber_Rejected = new TH1F("fPndVtxPoca_LambdaBar_QualityNumber_Rejected","fPndVtxPoca_LambdaBar_QualityNumber_Rejected",400,0,40); fPndVtxPoca_LambdaBar_VertexPosition_MCMatch = new TH2F("fPndVtxPoca_LambdaBar_VertexPosition_MCMatch","fPndVtxPoca_LambdaBar_VertexPosition_MCMatch",100,0,100,30,0,30); fPndVtxPoca_LambdaBar_VertexPosition_Rejected = new TH2F("fPndVtxPoca_LambdaBar_VertexPosition_Rejected","fPndVtxPoca_LambdaBar_VertexPosition_Rejected",100,0,100,30,0,30); // PndVtxPRGFitter fPndVtxPRGFitter_Lambda_VertexRes_X = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X","fPndVtxPRGFitter_Lambda_VertexRes_X; x_{Reco} - x_{MC} [cm]; #Events ",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y","fPndVtxPRGFitter_Lambda_VertexRes_Y; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z","fPndVtxPRGFitter_Lambda_VertexRes_Z; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch_Rejected; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch_Rejected; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch_Rejected; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut","fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut","fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut","fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut","fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut","fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut","fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_Rejected; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_Rejected; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_Rejected; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_Rejected; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch","fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_Both_VertexRes_X_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_Both_VertexRes_X_ProbCut_MCMatch","fPndVtxPRGFitter_Lambda_Both_VertexRes_X_ProbCut_MCMatch; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_Both_VertexRes_Y_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_Both_VertexRes_Y_ProbCut_MCMatch","fPndVtxPRGFitter_Lambda_Both_VertexRes_Y_ProbCut_MCMatch; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_Both_VertexRes_Z_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_Both_VertexRes_Z_ProbCut_MCMatch","fPndVtxPRGFitter_Lambda_Both_VertexRes_Z_ProbCut_MCMatch; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X","fPndVtxPRGFitter_LambdaBar_VertexRes_X; x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y","fPndVtxPRGFitter_LambdaBar_VertexRes_Y ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z","fPndVtxPRGFitter_LambdaBar_VertexRes_Z ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut","fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut","fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_X_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_Both_VertexRes_X_ProbCut_MCMatch","fPndVtxPRGFitter_LambdaBar_Both_VertexRes_X_ProbCut_MCMatch; x_{Reco} - x_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Y_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Y_ProbCut_MCMatch","fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Y_ProbCut_MCMatch; y_{Reco} - y_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Z_ProbCut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Z_ProbCut_MCMatch","fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Z_ProbCut_MCMatch; z_{Reco} - z_{MC} [cm]; #Events",200, -0.5,0.5); fPndVtxPRGFitter_Lambda_InvMass= new TH1F("fPndVtxPRGFitter_Lambda_InvMass","fPndVtxPRGFitter_Lambda_InvMass; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_InvMass_MCMatch","fPndVtxPRGFitter_Lambda_InvMass_MCMatch; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_InvMass_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_InvMass_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_Lambda_Chi2 = new TH1F("fPndVtxPRGFitter_Lambda_Chi2","fPndVtxPRGFitter_Lambda_Chi2",1000,-25,25); fPndVtxPRGFitter_Lambda_Chi2_MCMatch = new TH1F("fPndVtxPRGFitter_Lambda_Chi2_MCMatch","fPndVtxPRGFitter_Lambda_Chi2_MCMatch",1000,-25,25); fPndVtxPRGFitter_Lambda_Chi2_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_Lambda_Chi2_MCMatchRejected","fPndVtxPRGFitter_Lambda_Chi2_MCMatch_Rejected",1000,-25,25); fPndVtxPRGFitter_Lambda_Prob = new TH1F("fPndVtxPRGFitter_Lambda_Prob","fPndVtxPRGFitter_Lambda_Prob",260,-0.15,1.15); fPndVtxPRGFitter_Lambda_DecayLength = new TH1F("fPndVtxPRGFitter_Lambda_DecayLength","fPndVtxPRGFitter_Lambda_DecayLength",1000,0,100); fPndVtxPRGFitter_Lambda_DecayLength_MCMatch = new TH1F("fPndVtxPRGFitter_Lambda_DecayLength_MCMatch","fPndVtxPRGFitter_Lambda_DecayLength_MCMatch",1000,0,100); fPndVtxPRGFitter_Lambda_DecayLength_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_Lambda_DecayLength_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_DecayLength_MCMatch_Rejected",1000,0,100); fPndVtxPRGFitter_Lambda_PocaComparison = new TH2F("fPndVtxPRGFitter_Lambda_PocaComparison","fPndVtxPRGFitter_Lambda_PocaComparison",400,0,40,200,0,20); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut= new TH1F("fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut","fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_Rejected","fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_Rejected; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch","fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_Lambda_DecayLength_Chi2Cut = new TH1F("fPndVtxPRGFitter_Lambda_DecayLength_Chi2Cut","fPndVtxPRGFitter_Lambda_DecayLength_Chi2Cut",1000,0,100); fPndVtxPRGFitter_Lambda_InvMass_ProbCut = new TH1F("fPndVtxPRGFitter_Lambda_InvMass_ProbCut", "fPndVtxPRGFitter_Lambda_InvMass_ProbCut; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_ProbCut_Rejected = new TH1F("fPndVtxPRGFitter_Lambda_InvMass_ProbCut_Rejected", "fPndVtxPRGFitter_Lambda_InvMass_ProbCut_Rejecteds; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch = new TH1F("fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch", "fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected", "fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_Lambda_DecayLength_ProbCut = new TH1F("fPndVtxPRGFitter_Lambda_DecayLength_ProbCut","fPndVtxPRGFitter_Lambda_DecayLength_ProbCut; Inv. Mass [GeV/c]; #Events",1000,0,100); fPndVtxPRGFitter_Lambda_Both_InvMass_ProbCut_MCMatch = new TH1F("fPndVtxPRGFitter_Lambda_Both_InvMass_ProbCut_MCMatch", "fPndVtxPRGFitter_Lambda_Both_InvMass_ProbCut_MCMatch; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass= new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass","fPndVtxPRGFitter_LambdaBar_InvMass; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch","fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_LambdaBar_Chi2 = new TH1F("fPndVtxPRGFitter_LambdaBar_Chi2","fPndVtxPRGFitter_LambdaBar_Chi2",1000,-25,25); fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch = new TH1F("fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch","fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch",1000,-25,25); fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch_Rejected",1000,-25,25); fPndVtxPRGFitter_LambdaBar_Prob = new TH1F("fPndVtxPRGFitter_LambdaBar_Prob","fPndVtxPRGFitter_LambdaBar_Prob",260,-0.15,1.15); fPndVtxPRGFitter_LambdaBar_DecayLength = new TH1F("fPndVtxPRGFitter_LambdaBar_DecayLength","fPndVtxPRGFitter_LambdaBar_DecayLength",1000,0,100); fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch = new TH1F("fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch","fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch",1000,0,100); fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch_Rejected",1000,0,100); fPndVtxPRGFitter_LambdaBar_PocaComparison = new TH2F("fPndVtxPRGFitter_LambdaBar_PocaComparison","fPndVtxPRGFitter_LambdaBar_PocaComparison",400,0,40,200,0,20); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut= new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut","fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_Rejected","fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_Rejected; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch= new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch","fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected= new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected","fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events",1000,1,1.5); fPndVtxPRGFitter_LambdaBar_DecayLength_Chi2Cut = new TH1F("fPndVtxPRGFitter_LambdaBar_DecayLength_Chi2Cut","fPndVtxPRGFitter_LambdaBar_DecayLength_Chi2Cut",1000,0,100); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut = new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut", "fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_Rejected = new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_Rejected", "fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_Rejected; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch = new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch", "fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected = new TH1F("fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected", "fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPndVtxPRGFitter_LambdaBar_DecayLength_ProbCut = new TH1F("fPndVtxPRGFitter_LambdaBar_DecayLength_ProbCut","fPndVtxPRGFitter_LambdaBar_DecayLength_ProbCut",1000,0,100); fPndVtxPRGFitter_LambdaBar_Both_InvMass_ProbCut_MCMatch = new TH1F("fPndVtxPRGFitter_LambdaBar_Both_InvMass_ProbCut_MCMatch", "fPndVtxPRGFitter_LambdaBar_Both_InvMass_ProbCut_MCMatch; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fLambda_InvMass_MCMatch = new TH1F("fLambda_InvMass_MCMatch","#Lambda mass (MCTruth)",1000,1,1.5); fLambdaBar_InvMass_MCMatch = new TH1F("fLambdaBar_InvMass_MCMatch","#bar{#Lambda} mass (MCTruth)",1000,1,1.5); fNumberOfFoundLambda = new TH1I("fNumberOfFoundLambda","Number of found #Lambda",11, -0.5,10.5); fNumberOfFoundLambdaBar = new TH1I("fNumberOfFoundLambdaBar","Number of found #bar{#Lambda}",11, -0.5,10.5); fLambdaPerStage = new TH1I("fLambdaPerStage","Number of #Lambda per Stage", 10, 0.5,10.5); fLambdaBarPerStage = new TH1I("fLambdaBarPerStage","Number of #bar{#Lambda} per Stage", 10, 0.5,10.5); fCompleteSystemEnergie = new TH1F("fCompleteSystemEnergie", "Combined Lambda LambdaBar", 1000,1,1.5); fHitAnalysisTracksPerEvent2D = new TH2I("fHitAnalysisTracksPerEvent2D", "Hit Analysis Tracks/Event total 2D;Number of Tracks/Event", 6,-0.5,5.5,16,-0.5,15.5); fPnd4CFitter_Chi2=new TH1F("fPnd4CFitter_Chi2","fPnd4CFitter_Chi2",1000,-25,25); fPnd4CFitter_Prob= new TH1F("fPnd4CFitter_Prob","fPnd4CFitter_Prob",260,-0.15,1.15); fPnd4CFitter_PbarPSystem_InvMass= new TH1F("fPnd4CFitter_PbarPSystem_InvMass","PbarPSystem mass after 4CFitter ",5000,1.5,6.5); fPnd4CFitter_PbarPSystem_InvMass_MCMatch=new TH1F("fPnd4CFitter_PbarPSystem_InvMass_MCMatch","fPnd4CFitter_PbarPSystem_InvMass_MCMatch", 5000,1.5,6.5); fPnd4CFitter_PbarPSystem_InvMass_MCMatch_Rejected=new TH1F("fPnd4CFitter_PbarPSystem_InvMass_MCMatch_Rejected","fPnd4CFitter_PbarPSystem_InvMass_MCMatch_Rejected", 5000,1.5,6.5); fPnd4CFitter_Lambda_InvMass= new TH1F("fPnd4CFitter_Lambda_InvMass","#Lambda mass after 4CFitter ",1000,1,1.5); fPnd4CFitter_Lambda_InvMass_MCMatch=new TH1F("fPnd4CFitter_Lambda_InvMass_MCMatch","fPnd4CFitter_Lambda_InvMass_MCMatch", 1000,1,1.5); fPnd4CFitter_Lambda_InvMass_MCMatch_Rejected=new TH1F("fPnd4CFitter_Lambda_InvMass_MCMatch_Rejected","fPnd4CFitter_Lambda_InvMass_MCMatch_Rejected", 1000,1,1.5); fPnd4CFitter_Lambda_Chi2=new TH1F("fPnd4CFitter_Lambda_Chi2","fPnd4CFitter_Lambda_Chi2", 1000,-25,25); fPnd4CFitter_Lambda_Chi2_MCMatch=new TH1F("fPnd4CFitter_Lambda_Chi2_MCMatch","fPnd4CFitter_Lambda_Chi2_MCMatch", 1000,-25,25); fPnd4CFitter_Lambda_Chi2_MCMatch_Rejected=new TH1F("fPnd4CFitter_Lambda_Chi2_MCMatch_Rejected","fPnd4CFitter_Lambda_Chi2_MCMatch_Rejected", 1000,-25,25); fPnd4CFitter_Lambda_InvMass_ProbCut = new TH1F("fPnd4CFitter_Lambda_InvMass_ProbCut", "fPnd4CFitter_Lambda_InvMass_ProbCut; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_Lambda_InvMass_ProbCut_Rejected = new TH1F("fPnd4CFitter_Lambda_InvMass_ProbCut_Rejected", "fPnd4CFitter_Lambda_InvMass_ProbCut_Rejecteds; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch = new TH1F("fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch", "fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected = new TH1F("fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected", "fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_LambdaBar_InvMass=new TH1F("fPnd4CFitter_LambdaBar_InvMass","#Lambda mass 4CFitter",1000,1,1.5); fPnd4CFitter_LambdaBar_InvMass_MCMatch=new TH1F("fPnd4CFitter_LambdaBar_InvMass_MCMatch","#Lambda bar mass 4CFitter MCTruth", 1000,1,1.5); fPnd4CFitter_LambdaBar_InvMass_MCMatch_Rejected=new TH1F("fPnd4CFitter_LambdaBar_InvMass_MCMatch_Rejected","fPnd4CFitter_LambdaBar_InvMass_MCMatch_Rejected", 1000,1,1.5); fPnd4CFitter_LambdaBar_Chi2=new TH1F("fPnd4CFitter_LambdaBar_Chi","fPnd4CFitter_LambdaBar_Chi", 1000,-25,25); fPnd4CFitter_LambdaBar_Chi2_MCMatch=new TH1F("fPnd4CFitter_LambdaBar_Chi2_MCMatch","fPnd4CFitter_LambdaBar_Chi2_MCMatch", 1000,-25,25); fPnd4CFitter_LambdaBar_Chi2_MCMatch_Rejected=new TH1F("fPnd4CFitter_LambdaBar_Chi2_MCMatch_Rejected","fPnd4CFitter_LambdaBar_Chi2_MCMatch_Rejected", 1000,-25,25); fPnd4CFitter_LambdaBar_InvMass_ProbCut = new TH1F("fPnd4CFitter_LambdaBar_InvMass_ProbCut", "fPnd4CFitter_LambdaBar_InvMass_ProbCut; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_LambdaBar_InvMass_ProbCut_Rejected = new TH1F("fPnd4CFitter_LambdaBar_InvMass_ProbCut_Rejected", "fPnd4CFitter_LambdaBar_InvMass_ProbCut_Rejecteds; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch = new TH1F("fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch", "fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected = new TH1F("fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected", "fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected; Inv. Mass [GeV/c]; #Events", 1000,1,1.5); fPnd4CFitter_Lambda_VertexRes_X = new TH1F("fPnd4CFitter_Lambda_VertexRes_X","fPnd4CFitter_Lambda_VertexRes_X",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_Y = new TH1F("fPnd4CFitter_Lambda_VertexRes_Y","fPnd4CFitter_Lambda_VertexRes_Y",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_Z = new TH1F("fPnd4CFitter_Lambda_VertexRes_Z","fPnd4CFitter_Lambda_VertexRes_Z",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_X = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_X","fPnd4CFitter_LambdaBar_VertexRes_X",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_Y = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_Y","fPnd4CFitter_LambdaBar_VertexRes_Y",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_Z = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_Z","fPnd4CFitter_LambdaBar_VertexRes_Z",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_X_MCMatch = new TH1F("fPnd4CFitter_Lambda_VertexRes_X_MCMatch","fPnd4CFitter_Lambda_VertexRes_X_MCMatch",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_Y_MCMatch = new TH1F("fPnd4CFitter_Lambda_VertexRes_Y_MCMatch","fPnd4CFitter_Lambda_VertexRes_Y_MCMatch",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_Z_MCMatch = new TH1F("fPnd4CFitter_Lambda_VertexRes_Z_MCMatch","fPnd4CFitter_Lambda_VertexRes_Z_MCMatch",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch","fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch","fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch","fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_X_MCMatch_Rejected = new TH1F("fPnd4CFitter_Lambda_VertexRes_X_MCMatch_Rejected","fPnd4CFitter_Lambda_VertexRes_X_MCMatch_Rejected",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_Y_MCMatch_Rejected = new TH1F("fPnd4CFitter_Lambda_VertexRes_Y_MCMatch_Rejected","fPnd4CFitter_Lambda_VertexRes_Y_MCMatch_Rejected",100, -0.5,0.5); fPnd4CFitter_Lambda_VertexRes_Z_MCMatch_Rejected = new TH1F("fPnd4CFitter_Lambda_VertexRes_Z_MCMatch_Rejected","fPnd4CFitter_Lambda_VertexRes_Z_MCMatch_Rejected",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch_Rejected = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch_Rejected","fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch_Rejected",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected","fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected",100, -0.5,0.5); fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected = new TH1F("fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected","fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected",100, -0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X","fPndKinVtxFitter_Lambda_VertexRes_X;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y","fPndKinVtxFitter_Lambda_VertexRes_Y;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z","fPndKinVtxFitter_Lambda_VertexRes_Z;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut","fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut","fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut","fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut","fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut","fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut","fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_Rejected","fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_Rejected","fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch","fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected","fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X","fPndKinVtxFitter_LambdaBar_VertexRes_X ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y","fPndKinVtxFitter_LambdaBar_VertexRes_Y ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z","fPndKinVtxFitter_LambdaBar_VertexRes_Z ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut","fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut","fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut","fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut","fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut","fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut","fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch","fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected ;x_{Reco} - x_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected ;y_{Reco} - y_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected ;z_{Reco} - z_{MC} [cm]; #Events ",200,-0.5,0.5); fPndKinVtxFitter_Lambda_InvMass= new TH1F("fPndKinVtxFitter_Lambda_InvMass","fPndKinVtxFitter_Lambda_InvMass",1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_InvMass_MCMatch","fPndKinVtxFitter_Lambda_InvMass_MCMatch",1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_InvMass_MCMatch_Rejected","fPndKinVtxFitter_Lambda_InvMass_MCMatch_Rejected",1000,1,1.5); fPndKinVtxFitter_Lambda_Chi2 = new TH1F("fPndKinVtxFitter_Lambda_Chi2","fPndKinVtxFitter_Lambda_Chi2",1000,-25,25); fPndKinVtxFitter_Lambda_Chi2_MCMatch = new TH1F("fPndKinVtxFitter_Lambda_Chi2_MCMatch","fPndKinVtxFitter_Lambda_Chi2_MCMatch",1000,-25,25); fPndKinVtxFitter_Lambda_Chi2_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_Lambda_Chi2_MCMatchRejected","fPndKinVtxFitter_Lambda_Chi2_MCMatch_Rejected",1000,-25,25); fPndKinVtxFitter_Lambda_Prob = new TH1F("fPndKinVtxFitter_Lambda_Prob","fPndKinVtxFitter_Lambda_Prob",260,-0.15,1.15); fPndKinVtxFitter_Lambda_DecayLength = new TH1F("fPndKinVtxFitter_Lambda_DecayLength","fPndKinVtxFitter_Lambda_DecayLength",1000,0,100); fPndKinVtxFitter_Lambda_DecayLength_MCMatch = new TH1F("fPndKinVtxFitter_Lambda_DecayLength_MCMatch","fPndKinVtxFitter_Lambda_DecayLength_MCMatch",1000,0,100); fPndKinVtxFitter_Lambda_DecayLength_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_Lambda_DecayLength_MCMatch_Rejected","fPndKinVtxFitter_Lambda_DecayLength_MCMatch_Rejected",1000,0,100); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut= new TH1F("fPndKinVtxFitter_Lambda_InvMass_Chi2Cut","fPndKinVtxFitter_Lambda_InvMass_Chi2Cut",1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_Rejected","fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_Rejected",1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch","fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch",1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected",1000,1,1.5); fPndKinVtxFitter_Lambda_DecayLength_Chi2Cut = new TH1F("fPndKinVtxFitter_Lambda_DecayLength_Chi2Cut","fPndKinVtxFitter_Lambda_DecayLength_Chi2Cut",1000,0,100); fPndKinVtxFitter_Lambda_InvMass_ProbCut = new TH1F("fPndKinVtxFitter_Lambda_InvMass_ProbCut", "fPndKinVtxFitter_Lambda_InvMass_ProbCut", 1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_ProbCut_Rejected = new TH1F("fPndKinVtxFitter_Lambda_InvMass_ProbCut_Rejected", "fPndKinVtxFitter_Lambda_InvMass_ProbCut_Rejecteds", 1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch = new TH1F("fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch", "fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch", 1000,1,1.5); fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected", "fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected", 1000,1,1.5); fPndKinVtxFitter_Lambda_DecayLength_ProbCut = new TH1F("fPndKinVtxFitter_Lambda_DecayLength_ProbCut","fPndKinVtxFitter_Lambda_DecayLength_ProbCut",1000,0,100); fPndKinVtxFitter_LambdaBar_InvMass= new TH1F("fPndKinVtxFitter_LambdaBar_InvMass","fPndKinVtxFitter_LambdaBar_InvMass",1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_MCMatch","fPndKinVtxFitter_LambdaBar_InvMass_MCMatch",1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_InvMass_MCMatch_Rejected",1000,1,1.5); fPndKinVtxFitter_LambdaBar_Chi2 = new TH1F("fPndKinVtxFitter_LambdaBar_Chi2","fPndKinVtxFitter_LambdaBar_Chi2",1000,-25,25); fPndKinVtxFitter_LambdaBar_Chi2_MCMatch = new TH1F("fPndKinVtxFitter_LambdaBar_Chi2_MCMatch","fPndKinVtxFitter_LambdaBar_Chi2_MCMatch",1000,-25,25); fPndKinVtxFitter_LambdaBar_Chi2_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_LambdaBar_Chi2_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_Chi2_MCMatch_Rejected",1000,-25,25); fPndKinVtxFitter_LambdaBar_Prob = new TH1F("fPndKinVtxFitter_LambdaBar_Prob","fPndKinVtxFitter_LambdaBar_Prob",260,-0.15,1.15); fPndKinVtxFitter_LambdaBar_DecayLength = new TH1F("fPndKinVtxFitter_LambdaBar_DecayLength","fPndKinVtxFitter_LambdaBar_DecayLength",1000,0,100); fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch = new TH1F("fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch","fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch",1000,0,100); fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch_Rejected",1000,0,100); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut= new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut","fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut",1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_Rejected","fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_Rejected",1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch= new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch","fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch",1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected= new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected","fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected",1000,1,1.5); fPndKinVtxFitter_LambdaBar_DecayLength_Chi2Cut = new TH1F("fPndKinVtxFitter_LambdaBar_DecayLength_Chi2Cut","fPndKinVtxFitter_LambdaBar_DecayLength_Chi2Cut",1000,0,100); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut = new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_ProbCut", "fPndKinVtxFitter_LambdaBar_InvMass_ProbCut", 1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_Rejected = new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_Rejected", "fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_Rejected", 1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch = new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch", "fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch", 1000,1,1.5); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected = new TH1F("fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected", "fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected", 1000,1,1.5); fPndKinVtxFitter_LambdaBar_DecayLength_ProbCut = new TH1F("fPndKinVtxFitter_LambdaBar_DecayLength_ProbCut","fPndKinVtxFitter_LambdaBar_DecayLength_ProbCut",1000,0,100); fWrongAssignedList = new TH2I("fWrongAssignedList","fWrongAssignedList", 4, 0.5,4.5,20,0.5,20.5); fWrongAssignedList2 = new TH2I("fWrongAssignedList2","fWrongAssignedList2", 4, 0.5,4.5,20,0.5,20.5); fWrongAssignedList->GetXaxis()->SetBinLabel(1,"Proton List"); fWrongAssignedList->GetXaxis()->SetBinLabel(2,"Antiproton List"); fWrongAssignedList->GetXaxis()->SetBinLabel(3,"Piplus List"); fWrongAssignedList->GetXaxis()->SetBinLabel(4,"Piminus List"); fWrongAssignedList->GetYaxis()->SetBinLabel(1,"p right"); fWrongAssignedList->GetYaxis()->SetBinLabel(2,"pbar right"); fWrongAssignedList->GetYaxis()->SetBinLabel(3,"pi plus right"); fWrongAssignedList->GetYaxis()->SetBinLabel(4,"pi minus right"); fWrongAssignedList->GetYaxis()->SetBinLabel(5,"sum wrong"); fWrongAssignedList->GetYaxis()->SetBinLabel(6,"neutron "); fWrongAssignedList->GetYaxis()->SetBinLabel(7,"gamma"); fWrongAssignedList->GetYaxis()->SetBinLabel(8,"e+"); fWrongAssignedList->GetYaxis()->SetBinLabel(9,"e- "); fWrongAssignedList->GetYaxis()->SetBinLabel(10,"lambda"); fWrongAssignedList->GetYaxis()->SetBinLabel(11,"lambdabar"); fWrongAssignedList->GetYaxis()->SetBinLabel(12,"p"); fWrongAssignedList->GetYaxis()->SetBinLabel(13,"pbar "); fWrongAssignedList->GetYaxis()->SetBinLabel(14,"pi plus"); fWrongAssignedList->GetYaxis()->SetBinLabel(15,"pi minus"); fWrongAssignedList->GetYaxis()->SetBinLabel(20,"other"); fWrongAssignedList2->GetXaxis()->SetBinLabel(1,"Proton List"); fWrongAssignedList2->GetXaxis()->SetBinLabel(2,"Antiproton List"); fWrongAssignedList2->GetXaxis()->SetBinLabel(3,"Piplus List"); fWrongAssignedList2->GetXaxis()->SetBinLabel(4,"Piminus List"); fWrongAssignedList2->GetYaxis()->SetBinLabel(1,"p right"); fWrongAssignedList2->GetYaxis()->SetBinLabel(2,"pbar right"); fWrongAssignedList2->GetYaxis()->SetBinLabel(3,"pi plus right"); fWrongAssignedList2->GetYaxis()->SetBinLabel(4,"pi minus right"); fWrongAssignedList2->GetYaxis()->SetBinLabel(5,"sum wrong"); fWrongAssignedList2->GetYaxis()->SetBinLabel(6,"neutron "); fWrongAssignedList2->GetYaxis()->SetBinLabel(7,"gamma"); fWrongAssignedList2->GetYaxis()->SetBinLabel(8,"e+"); fWrongAssignedList2->GetYaxis()->SetBinLabel(9,"e- "); fWrongAssignedList2->GetYaxis()->SetBinLabel(10,"lambda"); fWrongAssignedList2->GetYaxis()->SetBinLabel(11,"lambdabar"); fWrongAssignedList2->GetYaxis()->SetBinLabel(12,"p"); fWrongAssignedList2->GetYaxis()->SetBinLabel(13,"pbar "); fWrongAssignedList2->GetYaxis()->SetBinLabel(14,"pi plus"); fWrongAssignedList2->GetYaxis()->SetBinLabel(15,"pi minus"); fWrongAssignedList2->GetYaxis()->SetBinLabel(20,"other"); return kSUCCESS; } void PndEventAnalysisLambdaLambdabar::FillAnalysisTracksPerEventSingleTrack(int protontrack, int antiprotontrack, int pionminustrack, int pionplustrack, TH2I* histo) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,0); if(protontrack==0 && antiprotontrack==1 && pionminustrack==1 && pionplustrack==1 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,1); } else if(protontrack==1 && antiprotontrack==0 && pionminustrack==1 && pionplustrack==1 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,2); } else if(protontrack==1 && antiprotontrack==1 && pionminustrack==0 && pionplustrack==1 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,3); } else if(protontrack==1 && antiprotontrack==1 && pionminustrack==1 && pionplustrack==0 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,4); } else if(protontrack==1 && antiprotontrack==1 && pionminustrack==0 && pionplustrack==0 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,5); } else if(protontrack==1 && antiprotontrack==0 && pionminustrack==1 && pionplustrack==0 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,6); } else if(protontrack==0 && antiprotontrack==1 && pionminustrack==1 && pionplustrack==0 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,7); } else if(protontrack==1 && antiprotontrack==0 && pionminustrack==0 && pionplustrack==1 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,8); } else if(protontrack==0 && antiprotontrack==1 && pionminustrack==0 && pionplustrack==1 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,9); } else if(protontrack==0 && antiprotontrack==0 && pionminustrack==1 && pionplustrack==1 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,10); } else if(protontrack==1 && antiprotontrack==0 && pionminustrack==0 && pionplustrack==0 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,11); } else if(protontrack==0 && antiprotontrack==1 && pionminustrack==0 && pionplustrack==0 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,12); } else if(protontrack==0 && antiprotontrack==0 && pionminustrack==1 && pionplustrack==0 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,13); } else if(protontrack==0 && antiprotontrack==0 && pionminustrack==0 && pionplustrack==1 ) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,14); } else if(protontrack==0 && antiprotontrack==0 && pionminustrack==0 && pionplustrack==0) { histo->Fill(protontrack+antiprotontrack+pionplustrack+pionminustrack,15); } } void PndEventAnalysisLambdaLambdabar::SetHitAnalysisEventAxis(TH1* histo) { histo->GetYaxis()->SetBinLabel(1,"Vertical Sum"); histo->GetYaxis()->SetBinLabel(2,"p mis."); histo->GetYaxis()->SetBinLabel(3,"pbar mis."); histo->GetYaxis()->SetBinLabel(4,"#pi- mis."); histo->GetYaxis()->SetBinLabel(5,"#pi+ mis."); histo->GetYaxis()->SetBinLabel(6,"#pi- #pi+ mis."); histo->GetYaxis()->SetBinLabel(7,"pbar #pi+ mis."); histo->GetYaxis()->SetBinLabel(8,"p #pi+ mis."); histo->GetYaxis()->SetBinLabel(9,"pbar #pi- mis."); histo->GetYaxis()->SetBinLabel(10,"p #pi- mis."); histo->GetYaxis()->SetBinLabel(11,"p pbar mis."); histo->GetYaxis()->SetBinLabel(12,"just p"); histo->GetYaxis()->SetBinLabel(13,"just antip"); histo->GetYaxis()->SetBinLabel(14,"just #pi-"); histo->GetYaxis()->SetBinLabel(15,"just #pi+"); histo->GetYaxis()->SetBinLabel(16,"no track"); } int PndEventAnalysisLambdaLambdabar::WrongPID(int pdg) { if(pdg==2112) { return 6; } else if(pdg==22) { return 7; } else if(pdg==11) { return 8; } else if(pdg==-11) { return 9; } else if(pdg==3122) { return 10; } else if(pdg==-3122) { return 11; } else if(pdg==2212) { return 12; } else if(pdg==-2212) { return 13; } else if(pdg==211) { return 14; } else if(pdg==-211) { return 15; } else { return 20; } } int PndEventAnalysisLambdaLambdabar::SelectWithMcTruthPdgSimulateMcMatch(RhoCandList &lOld, RhoCandList &lNew, Int_t SelectCandidatesWithPdgCode) { lNew.Cleanup(); int notSelected = 0; for ( int ii=lOld.GetLength()-1; ii>=0; --ii ) { if ( ( lOld[ii]->GetMcTruth()->PdgCode() == SelectCandidatesWithPdgCode ) && ( lOld[ii]->PdgCode() == SelectCandidatesWithPdgCode ) ) { lNew.Add( lOld[ii] ); } else { notSelected++; } } return notSelected; } void PndEventAnalysisLambdaLambdabar::Exec(Option_t* opt) { //TLorentzVector ini(0, 0, 6.231711, 7.240222); TLorentzVector ini(0, 0, 1.642, 2.8290); // muss fŸr andere energien geŠndert werden!!! FairRootManager* ioman = FairRootManager::Instance(); int protontrack=0; int antiprotontrack=0; int pionplustrack=0; int pionminustrack=0; int eventwithlambda=0; int eventwithlambdabar=0; if(fVerbose>0) { std::cout<< "============== Begin PndEventAnalysisLambdaLambdabar::Exec for Event " << eventcounter << std::endl; } theAnalysis->GetEvent(); if (!(eventcounter%100)) { std::cout <<"evt "<Stop(); std::cout <<" t="<< timer->CpuTime(); timer->Start(); std::cout << "\n"; } TString pidalgos("PidAlgoStt;PidAlgoEmcBayes;PidAlgoDrc;PidAlgoDisc"); RhoCandList protonplus, protonminus, pionplus, pionminus, mctrack, pionall, protonall; theAnalysis->FillList(mctrack,"McTruth","McTruth"); theAnalysis->FillList(pionall,"PionAll"); theAnalysis->FillList(protonall,"ProtonAll"); SelectWithMcTruthPdgSimulateMcMatch(pionall,pionplus,211); SelectWithMcTruthPdgSimulateMcMatch(pionall,pionminus,-211); SelectWithMcTruthPdgSimulateMcMatch(protonall,protonplus,2212); SelectWithMcTruthPdgSimulateMcMatch(protonall,protonminus,-2212); for(int y=0; y < protonplus.GetLength();y++) { int pdg = protonplus[y]->GetMcTruth()->PdgCode(); if(pdg!=2212) { fWrongAssignedList->Fill(1,WrongPID(pdg)); fWrongAssignedList->Fill(1,5); } else { fKalmanFitterFlag->Fill(protonplus[y]->GetRecoCandidate()->GetFitStatus()); fWrongAssignedList->Fill(1,1); } if(protonplus[y]->GetRecoCandidate()->GetMcIndex()==5) { protontrack++; } } for(int y=0; y < protonminus.GetLength();y++) { int pdg =protonminus[y]->GetMcTruth()->PdgCode(); if(protonminus[y]->GetMcTruth()->PdgCode()!=-2212) { fWrongAssignedList->Fill(2,5); fWrongAssignedList->Fill(2,WrongPID(pdg)); } else { fKalmanFitterFlag->Fill(protonminus[y]->GetRecoCandidate()->GetFitStatus()); fWrongAssignedList->Fill(2,2); } if(protonminus[y]->GetRecoCandidate()->GetMcIndex()==3) { antiprotontrack++; } } for(int y=0; y < pionplus.GetLength();y++) { int pdg = pionplus[y]->GetMcTruth()->PdgCode(); if(pdg!=211) { fWrongAssignedList->Fill(3,5); fWrongAssignedList->Fill(3,WrongPID(pdg)); } else { fKalmanFitterFlag->Fill(pionplus[y]->GetRecoCandidate()->GetFitStatus()); fWrongAssignedList->Fill(3,3); } if(pionplus[y]->GetRecoCandidate()->GetMcIndex()==4) { pionplustrack++; } } for(int y=0; y < pionminus.GetLength();y++) { int pdg = pionminus[y]->GetMcTruth()->PdgCode(); if(pdg!=-211) { fWrongAssignedList->Fill(4,5); fWrongAssignedList->Fill(4,WrongPID(pdg)); } else { fKalmanFitterFlag->Fill(pionminus[y]->GetRecoCandidate()->GetFitStatus()); pionminus[y]->GetRecoCandidate()->GetLink(0); fWrongAssignedList->Fill(4,4); } if(pionminus[y]->GetRecoCandidate()->GetMcIndex()==6) { pionminustrack++; } } if(protontrack>1 || antiprotontrack>1 || pionplustrack>1 || pionminustrack>1) { fErrorCounter++; } if(fVerbose>0) { std::cout << "eventcounter " << eventcounter << std::endl; std::cout << "proton " << protontrack << std::endl; std::cout << "antiproton " << antiprotontrack << std::endl; std::cout << "pionminus " << pionminustrack << std::endl; std::cout << "pionplus " << pionplustrack << std::endl; } fTracksPerEvent->Fill(1,protontrack); fTracksPerEvent->Fill(2,antiprotontrack); fTracksPerEvent->Fill(3,pionplustrack); fTracksPerEvent->Fill(4,pionminustrack); for(int g=0;gGetMcTruth()->PdgCode()==-2212) { fAntiProton_NumberOfDaughters->Fill(protonminus[g]->GetMcTruth()->NDaughters()); if(protonminus[g]->GetMcTruth()->NDaughters()>0) { fAntiProton_DecayPoint->Fill(protonminus[g]->GetMcTruth()->Daughter(0)->GetPosition().Z(),protonminus[g]->GetMcTruth()->Daughter(0)->GetPosition().Perp()); } } } SetHitAnalysisEventAxis(fHitAnalysisTracksPerEvent2D); FillAnalysisTracksPerEventSingleTrack(protontrack,antiprotontrack, pionminustrack,pionplustrack, fHitAnalysisTracksPerEvent2D); RhoCandList lambdabar, llbar, lambda, lambda_PndKinVtxFitter, lambdabar_PndKinVtxFitter,lambda_PndVtxPRG,lambdabar_PndVtxPRG; lambda.Cleanup(); lambdabar.Cleanup(); lambda_PndVtxPRG.Cleanup(); lambdabar_PndVtxPRG.Cleanup(); lambda_PndKinVtxFitter.Cleanup(); lambdabar_PndKinVtxFitter.Cleanup(); for (int j = 0; j < protonplus.GetLength(); j++) for (int k = 0; k < pionminus.GetLength(); k++) { RhoCandidate* combcand = protonplus[j]->Combine(pionminus[k]); if(protonplus[j]->GetMcTruth()->PdgCode() != 2212 || pionminus[k]->GetMcTruth()->PdgCode() != -211) continue; lambda.Append(combcand); lambda_PndKinVtxFitter.Append(combcand); lambda_PndVtxPRG.Append(combcand); } fNumberOfFoundLambda->Fill(1,lambda.GetLength()); for (int j = 0; j < protonminus.GetLength(); j++) for (int k = 0; k < pionplus.GetLength(); k++) { RhoCandidate* combcand = protonminus[j]->Combine(pionplus[k]); if (protonminus[j]->GetMcTruth()->PdgCode() != -2212 || pionplus[k]->GetMcTruth()->PdgCode() != 211) { continue; } lambdabar.Append(combcand); lambdabar_PndKinVtxFitter.Append(combcand); lambdabar_PndVtxPRG.Append(combcand); } //lambdabar.Combine(protonminus, pionplus); fNumberOfFoundLambdaBar->Fill(lambdabar.GetLength()); lambda.SetType(3122); lambdabar.SetType(-3122); lambda_PndKinVtxFitter.SetType(3122); lambdabar_PndKinVtxFitter.SetType(-3122); lambda_PndVtxPRG.SetType(3122); lambdabar_PndVtxPRG.SetType(-3122); // llbar.SetType(88888); // // llbar.Combine(lambda,lambdabar); // // // for(int q=0; q< llbar.GetLength();q++) // { // std::cout << "llbar found " << std::endl; // //if(theAnalysis->McTruthMatch(llbar[q],2,true)) // //{ // // std::cout << "blubb" << std::endl; // //} // } bool lambda_found=false; bool lambdabar_found=false; for (int p=0; p< lambda.GetLength();p++) { if(theAnalysis->McTruthMatch(lambda[p],2,false)) { lambda_found=true; } } for (int p=0; p< lambdabar.GetLength();p++) { if(theAnalysis->McTruthMatch(lambdabar[p],2,false)) { lambdabar_found=true; } } for (int r=0; r< lambda.GetLength();r++) { fLambda_InvMass->Fill( lambda[r]->M() ); fLambdasPerMCTruthMatchStage->Fill(3); if(theAnalysis->McTruthMatch(lambda[r],2,false)) { if(lambda_found==true && lambdabar_found==true) { TLorentzVector boostedgesamtvector= lambda[r]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdavector= lambda[r]->GetMcTruth()->P4(); boostedlambdavector.Boost(-betavectorgesamt); TVector3 betavectorlambda = boostedlambdavector.BoostVector(); TLorentzVector lorentzvectorproton_lambdarestframeboost = lambda[r]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorgesamt); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorlambda); TVector3 protonvector = lorentzvectorproton_lambdarestframeboost.Vect(); TVector3 lambda_z_direction = boostedlambdavector.Vect(); TVector3 lambda_y_direction = lambda_z_direction.Cross(TVector3(0,0,1)); TVector3 lambda_x_direction = lambda_y_direction.Cross(lambda_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X->Fill(TMath::Cos(lambda_x_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y->Fill(TMath::Cos(lambda_y_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z->Fill(TMath::Cos(lambda_z_direction.Angle(protonvector))); fReconstructedParticle_AngularDistribution_Both_Lambda->Fill(boostedlambdavector.CosTheta()); fNumberOfFoundLambda->Fill(2); fLambda_InvMass_MCMatch->Fill( lambda[r]->M() ); fLambdasPerMCTruthMatchStage->Fill(2); } TLorentzVector boostedgesamtvector= lambda[r]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdavector= lambda[r]->GetMcTruth()->P4(); boostedlambdavector.Boost(-betavectorgesamt); TVector3 betavectorlambda = boostedlambdavector.BoostVector(); TLorentzVector lorentzvectorproton_lambdarestframeboost = lambda[r]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorgesamt); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorlambda); TVector3 protonvector = lorentzvectorproton_lambdarestframeboost.Vect(); TVector3 lambda_z_direction = boostedlambdavector.Vect(); TVector3 lambda_y_direction = lambda_z_direction.Cross(TVector3(0,0,1)); TVector3 lambda_x_direction = lambda_y_direction.Cross(lambda_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X->Fill(TMath::Cos(lambda_x_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y->Fill(TMath::Cos(lambda_y_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z->Fill(TMath::Cos(lambda_z_direction.Angle(protonvector))); fReconstructedParticle_AngularDistribution_Lambda->Fill(boostedlambdavector.CosTheta()); fNumberOfFoundLambda->Fill(2); fLambda_InvMass_MCMatch->Fill( lambda[r]->M() ); fLambdasPerMCTruthMatchStage->Fill(2); } else { fLambda_InvMass_MCMatch_Rejected->Fill(lambda[r]->M()); } } for (int r=0; r< lambdabar.GetLength();r++) { fLambdaBar_InvMass->Fill( lambdabar[r]->M() ); fLambdaBarsPerMCTruthMatchStage->Fill(3); if(theAnalysis->McTruthMatch(lambdabar[r],2,false)) { if(lambda_found==true && lambdabar_found==true) { TLorentzVector boostedgesamtvector= lambdabar[r]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdabarvector= lambdabar[r]->GetMcTruth()->P4(); boostedlambdabarvector.Boost(-betavectorgesamt); TVector3 betavectorlambdabar = boostedlambdabarvector.BoostVector(); TLorentzVector lorentzvectorantiproton_lambdabarrestframeboost = lambdabar[r]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorgesamt); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorlambdabar); TVector3 antiprotonvector = lorentzvectorantiproton_lambdabarrestframeboost.Vect(); TVector3 lambdabar_z_direction = boostedlambdabarvector.Vect(); TVector3 lambdabar_y_direction = lambdabar_z_direction.Cross(TVector3(0,0,1)); TVector3 lambdabar_x_direction = lambdabar_y_direction.Cross(lambdabar_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X->Fill(TMath::Cos(lambdabar_x_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y->Fill(TMath::Cos(lambdabar_y_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z->Fill(TMath::Cos(lambdabar_z_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistribution_Both_LambdaBar->Fill(boostedlambdabarvector.CosTheta()); fLambdaBar_InvMass_MCMatch->Fill( lambdabar[r]->M() ); fLambdaBarsPerMCTruthMatchStage->Fill(2); } if(lambdabar[r]->GetMcTruth()->GetMomentum().CosTheta() >0.98) { fAntiProton_Theta_Phi_Cut->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg() ,lambdabar[r]->GetMcTruth()->Daughter(0)->GetPosition().Theta()*TMath::RadToDeg()); fAntiProton_CosTheta_Phi_Cut->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg() ,lambdabar[r]->GetMcTruth()->P4().Vect().CosTheta()); } fAntiProton_Theta_Phi->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg() ,lambdabar[r]->GetMcTruth()->Daughter(0)->GetPosition().Theta()*TMath::RadToDeg()); fAntiProton_CosTheta_Phi->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg() ,lambdabar[r]->GetMcTruth()->P4().Vect().CosTheta()); std::cout << "lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi() " << lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()<< std::endl; std::cout << "lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg() " << lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg()<< std::endl; std::cout << "lambdabar[r]->GetMcTruth()->Daughter(0)->GetPosition().Theta() " << lambdabar[r]->GetMcTruth()->Daughter(0)->GetPosition().Theta()<< std::endl; std::cout << "lambdabar[r]->GetMcTruth()->Daughter(0)->GetPosition().Theta()*TMath::RadToDeg() " << lambdabar[r]->GetMcTruth()->Daughter(0)->GetPosition().Theta()*TMath::RadToDeg() << std::endl; TLorentzVector boostedgesamtvector= lambdabar[r]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdabarvector= lambdabar[r]->GetMcTruth()->P4(); boostedlambdabarvector.Boost(-betavectorgesamt); TVector3 betavectorlambdabar = boostedlambdabarvector.BoostVector(); TLorentzVector lorentzvectorantiproton_lambdabarrestframeboost = lambdabar[r]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorgesamt); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorlambdabar); TVector3 antiprotonvector = lorentzvectorantiproton_lambdabarrestframeboost.Vect(); TVector3 lambdabar_z_direction = boostedlambdabarvector.Vect(); TVector3 lambdabar_y_direction = lambdabar_z_direction.Cross(TVector3(0,0,1)); TVector3 lambdabar_x_direction = lambdabar_y_direction.Cross(lambdabar_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X->Fill(TMath::Cos(lambdabar_x_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y->Fill(TMath::Cos(lambdabar_y_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z->Fill(TMath::Cos(lambdabar_z_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistribution_LambdaBar->Fill(boostedlambdabarvector.CosTheta()); if(lambdabar[r]->GetMcTruth()->GetMomentum().CosTheta() >0.98) { fAntiProton_Theta_Phi_Cut_CMS->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg(), boostedlambdabarvector.Theta()*TMath::RadToDeg()); fAntiProton_CosTheta_Phi_Cut_CMS->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg(), boostedlambdabarvector.CosTheta()); } fAntiProton_Theta_Phi_CMS->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg(), boostedlambdabarvector.Theta()*TMath::RadToDeg()); fAntiProton_CosTheta_Phi_CMS->Fill(lambdabar[r]->Daughter(0)->GetMcTruth()->P4().Phi()*TMath::RadToDeg(), boostedlambdabarvector.CosTheta()); fLambdaBar_InvMass_MCMatch->Fill( lambdabar[r]->M() ); fLambdaBarsPerMCTruthMatchStage->Fill(2); } else if(theAnalysis->McTruthMatch(lambdabar[r],1,false)) { fLambdaBarsPerMCTruthMatchStage->Fill(1); fLambdaBar_InvMass_MCMatch_Rejected->Fill(lambdabar[r]->M()); } else if(theAnalysis->McTruthMatch(lambdabar[r],0,false)) { fLambdaBarsPerMCTruthMatchStage->Fill(0); fLambdaBar_InvMass_MCMatch_Rejected->Fill(lambdabar[r]->M()); } else { fLambdaBar_InvMass_MCMatch_Rejected->Fill(lambdabar[r]->M()); } } if(lambda_found==true && lambdabar_found==true) { std::cout << "###################### beide teilchen gefunden " << std::endl; } // ********************** // POCA Vertex Finder // ********************** for(int z=0 ; z< lambda.GetLength();z++) { TVector3 VertexPosition; PndVtxPoca LambdaPndVtxPoca; double QualityValue = LambdaPndVtxPoca.GetPocaVtx(VertexPosition,lambda[z]); fPndVtxPoca_Lambda_QualityNumber->Fill(QualityValue); fPndVtxPoca_Lambda_VertexPosition->Fill(VertexPosition.Z(),VertexPosition.Perp()); bool MCMatch = theAnalysis->McTruthMatch(lambda[z],2,false); if(MCMatch) { fPndVtxPoca_Lambda_QualityNumber_MCMatch->Fill(LambdaPndVtxPoca.GetPocaVtx(VertexPosition,lambda[z])); fPndVtxPoca_Lambda_VertexPosition_MCMatch->Fill(VertexPosition.Z(),VertexPosition.Perp()); } else { fPndVtxPoca_Lambda_QualityNumber_Rejected->Fill(LambdaPndVtxPoca.GetPocaVtx(VertexPosition,lambda[z])); fPndVtxPoca_Lambda_VertexPosition_Rejected->Fill(VertexPosition.Z(),VertexPosition.Perp()); } } for(int z=0 ; z< lambdabar.GetLength();z++) { TVector3 VertexPosition; PndVtxPoca LambdaBarPndVtxPoca; fPndVtxPoca_LambdaBar_QualityNumber->Fill(LambdaBarPndVtxPoca.GetPocaVtx(VertexPosition,lambdabar[z])); fPndVtxPoca_LambdaBar_VertexPosition->Fill(VertexPosition.Z(),VertexPosition.Perp()); bool MCMatch = theAnalysis->McTruthMatch(lambdabar[z],2,false); if(MCMatch) { fPndVtxPoca_LambdaBar_QualityNumber_MCMatch->Fill(LambdaBarPndVtxPoca.GetPocaVtx(VertexPosition,lambdabar[z])); fPndVtxPoca_LambdaBar_VertexPosition_MCMatch->Fill(VertexPosition.Z(),VertexPosition.Perp()); } else { fPndVtxPoca_LambdaBar_QualityNumber_Rejected->Fill(LambdaBarPndVtxPoca.GetPocaVtx(VertexPosition,lambdabar[z])); fPndVtxPoca_LambdaBar_VertexPosition_Rejected->Fill(VertexPosition.Z(),VertexPosition.Perp()); } } // ********************** // PndKinVtxFitter // ********************** RhoCandidate* lambdaPndKinVtxFitterfit; double chicut=2; double ProbCut=0.1; for(int k=0;k < lambda_PndKinVtxFitter.GetLength();k++) { bool MCMatch=theAnalysis->McTruthMatch(lambda_PndKinVtxFitter[k],2,false); PndKinVtxFitter lambdakinvtxfitter(lambda_PndKinVtxFitter[k]); lambdakinvtxfitter.Fit(); double kinvtxfitter_chi2=lambdakinvtxfitter.GetChi2(); double kinvtxfitter_prob=lambdakinvtxfitter.GetProb(); RhoCandidate* lambdaPndKinVtxFitterFit=lambda_PndKinVtxFitter[k]->GetFit(); fPndKinVtxFitter_Lambda_InvMass->Fill(lambda_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_Lambda_Chi2->Fill(kinvtxfitter_chi2); fPndKinVtxFitter_Lambda_Prob->Fill(kinvtxfitter_prob); fPndKinVtxFitter_Lambda_DecayLength->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); fPndKinVtxFitter_Lambda_VertexRes_X->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); if(MCMatch) { fPndKinVtxFitter_Lambda_InvMass_MCMatch->Fill((lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_Lambda_Chi2_MCMatch->Fill(lambdakinvtxfitter.GetChi2()); fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndKinVtxFitter_Lambda_DecayLength_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); } else { fPndKinVtxFitter_Lambda_InvMass_MCMatch_Rejected->Fill((lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_Lambda_Chi2_MCMatch_Rejected->Fill(lambdakinvtxfitter.GetChi2()); fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndKinVtxFitter_Lambda_DecayLength_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); } if(kinvtxfitter_prob>ProbCut && kinvtxfitter_chi2 > 0) { fPndKinVtxFitter_Lambda_InvMass_ProbCut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_Lambda_DecayLength_ProbCut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); if(MCMatch) { TLorentzVector boostedgesamtvector= lambda_PndKinVtxFitter[k]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdavector= lambda_PndKinVtxFitter[k]->GetMcTruth()->P4(); boostedlambdavector.Boost(-betavectorgesamt); TVector3 betavectorlambda = boostedlambdavector.BoostVector(); TLorentzVector lorentzvectorproton_lambdarestframeboost = lambda_PndKinVtxFitter[k]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorgesamt); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorlambda); TVector3 protonvector = lorentzvectorproton_lambdarestframeboost.Vect(); TVector3 lambda_z_direction = boostedlambdavector.Vect(); TVector3 lambda_y_direction = lambda_z_direction.Cross(TVector3(0,0,1)); TVector3 lambda_x_direction = lambda_y_direction.Cross(lambda_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndKinVtxFitter->Fill(TMath::Cos(lambda_x_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndKinVtxFitter->Fill(TMath::Cos(lambda_y_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndKinVtxFitter->Fill(TMath::Cos(lambda_z_direction.Angle(protonvector))); // std::cout << "Lambda : Cos des Winkels zur x achse " << TMath::Cos(lambda_x_direction.Angle(protonvector)) << std::endl; // std::cout << "Lambda : Cos des Winkels zur y achse " << TMath::Cos(lambda_y_direction.Angle(protonvector)) << std::endl; // std::cout << "Lambda : Cos des Winkels zur z achse " << TMath::Cos(lambda_z_direction.Angle(protonvector)) << std::endl; fReconstructedParticle_AngularDistribution_Lambda_PndKinVtxFitter->Fill(boostedlambdavector.CosTheta()); fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch->Fill((lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } else { fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected->Fill((lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndKinVtxFitter_Lambda_InvMass_ProbCut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } if(kinvtxfitter_chi2< chicut && kinvtxfitter_chi2 > 0) { fPndKinVtxFitter_Lambda_InvMass_Chi2Cut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_Lambda_DecayLength_Chi2Cut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); if(MCMatch) { fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch->Fill((lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } else { fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected->Fill((lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambda_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected->Fill(lambda_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambda_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } RhoCandidate* lambdabar_PndKinVtxFitterfit; for(int k=0;k < lambdabar_PndKinVtxFitter.GetLength();k++) { PndKinVtxFitter lambdabarkinvtxfitter(lambdabar_PndKinVtxFitter[k]); //lambdabarkinvtxfitter.AddMassConstraint(1.115); lambdabarkinvtxfitter.Fit(); double kinvtxfitter_chi2=lambdabarkinvtxfitter.GetChi2(); double kinvtxfitter_prob=lambdabarkinvtxfitter.GetProb(); bool MCMatch=theAnalysis->McTruthMatch(lambdabar_PndKinVtxFitter[k],2,false); RhoCandidate* lambdabarPndKinVtxFitterFit=lambdabar_PndKinVtxFitter[k]->GetFit(); fPndKinVtxFitter_LambdaBar_InvMass->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_LambdaBar_VertexRes_X->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndKinVtxFitter_LambdaBar_Chi2->Fill(kinvtxfitter_chi2); fPndKinVtxFitter_LambdaBar_Prob->Fill(kinvtxfitter_prob); fPndKinVtxFitter_LambdaBar_DecayLength->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); if(theAnalysis->McTruthMatch(lambdabar_PndKinVtxFitter[k],2,false)) { fPndKinVtxFitter_LambdaBar_InvMass_MCMatch->Fill((lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndKinVtxFitter_LambdaBar_Chi2_MCMatch->Fill(lambdabarkinvtxfitter.GetChi2()); fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); } else { fPndKinVtxFitter_LambdaBar_InvMass_MCMatch_Rejected->Fill((lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndKinVtxFitter_LambdaBar_Chi2_MCMatch_Rejected->Fill(lambdabarkinvtxfitter.GetChi2()); fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); } if(kinvtxfitter_prob>ProbCut && kinvtxfitter_chi2 > 0) { fPndKinVtxFitter_LambdaBar_InvMass_ProbCut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_LambdaBar_DecayLength_ProbCut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); if(MCMatch) { TLorentzVector boostedgesamtvector= lambdabar_PndKinVtxFitter[k]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdabarvector= lambdabar_PndKinVtxFitter[k]->GetMcTruth()->P4(); boostedlambdabarvector.Boost(-betavectorgesamt); TVector3 betavectorlambdabar = boostedlambdabarvector.BoostVector(); TLorentzVector lorentzvectorantiproton_lambdabarrestframeboost = lambdabar_PndKinVtxFitter[k]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorgesamt); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorlambdabar); TVector3 antiprotonvector = lorentzvectorantiproton_lambdabarrestframeboost.Vect(); TVector3 lambdabar_z_direction = boostedlambdabarvector.Vect(); TVector3 lambdabar_y_direction = lambdabar_z_direction.Cross(TVector3(0,0,1)); TVector3 lambdabar_x_direction = lambdabar_y_direction.Cross(lambdabar_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndKinVtxFitter->Fill(TMath::Cos(lambdabar_x_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndKinVtxFitter->Fill(TMath::Cos(lambdabar_y_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndKinVtxFitter->Fill(TMath::Cos(lambdabar_z_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistribution_LambdaBar_PndKinVtxFitter->Fill(boostedlambdabarvector.CosTheta()); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch->Fill((lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } else { fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected->Fill((lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } if(kinvtxfitter_chi2< chicut && kinvtxfitter_chi2 > 0) { fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_LambdaBar_DecayLength_Chi2Cut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Mag()*lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->M()/(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P3()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P3()).Mag()); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); if(MCMatch) { fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch->Fill((lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } else { fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected->Fill((lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->P4()+lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(1)->P4()).M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->M()); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().X()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Y()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected->Fill(lambdabar_PndKinVtxFitter[k]->GetFit()->Daughter(0)->GetPosition().Z()-lambdabar_PndKinVtxFitter[k]->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } // ********************** // PndVtxPRG // ********************** ProbCut= 0.01; RhoCandList lambda_PndVtxPRGFitter_ProbCut; RhoCandList lambdabar_PndVtxPRGFitter_ProbCut; for(int j=0;j < lambda_PndVtxPRG.GetLength();j++) { bool MCMatch= theAnalysis->McTruthMatch(lambda_PndVtxPRG[j],2,false); PndVtxPRG prgfitter(lambda_PndVtxPRG[j]); prgfitter.Fit(); double chi2=prgfitter.GetChi2(); double prob=prgfitter.GetProb(); if(!chi2>0) { continue; } RhoCandidate* lambdaPndVtxPRGfit=lambda_PndVtxPRG[j]->GetFit(); fPndVtxPRGFitter_Lambda_InvMass->Fill(lambdaPndVtxPRGfit->M()); fPndVtxPRGFitter_Lambda_Chi2->Fill(chi2); fPndVtxPRGFitter_Lambda_Prob->Fill(prob); fPndVtxPRGFitter_Lambda_DecayLength->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdaPndVtxPRGfit->Daughter(0)->M()/(lambdaPndVtxPRGfit->Daughter(0)->P3()+lambdaPndVtxPRGfit->Daughter(1)->P3()).Mag()); fPndVtxPRGFitter_Lambda_VertexRes_X->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); double alphalambda =0; TVector3 VertexPosition; PndVtxPoca LambdaPndVtxPoca_BeforeFit; PndVtxPoca LambdaPndVtxPoca_AfterFit; fPndVtxPRGFitter_Lambda_PocaComparison->Fill(LambdaPndVtxPoca_BeforeFit.GetPocaVtx(VertexPosition,lambda_PndVtxPRG[j]),LambdaPndVtxPoca_AfterFit.GetPocaVtx(VertexPosition,lambdaPndVtxPRGfit)); TVector3 momentumproton = lambdaPndVtxPRGfit->Daughter(0)->GetMomentum(); TVector3 momentumpionminus = lambdaPndVtxPRGfit->Daughter(1)->GetMomentum(); TVector3 momentumvectorlambda= momentumproton + momentumpionminus; double protontransversalcomponent = momentumproton.Perp(momentumvectorlambda); double pionminustransversalcomponent = momentumpionminus.Perp(momentumvectorlambda); double protonlongitudinalcomponent = momentumproton*momentumvectorlambda; double pionminuslongitudinalcomponent = momentumpionminus*momentumvectorlambda; alphalambda = ((protonlongitudinalcomponent-pionminuslongitudinalcomponent)/(protonlongitudinalcomponent+pionminuslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter->Fill(alphalambda,protontransversalcomponent); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter->Fill(alphalambda,protontransversalcomponent); lambdaPndVtxPRGfit->Daughter(0); // proton lambdaPndVtxPRGfit->Daughter(1); // pionminus if(fVerbose>1) { std::cout << "Information About Lambda candidate "<Daughter(0)->GetRecoCandidate()->GetTrackIndex() << " PID " << lambda_PndVtxPRG[j]->Daughter(0)->GetMcTruth()->PdgCode() <Daughter(1)->GetRecoCandidate()->GetTrackIndex() << " PID " << lambda_PndVtxPRG[j]->Daughter(1)->GetMcTruth()->PdgCode() << std::endl; std::cout << "MCMatch? " << MCMatch << std::endl; if(MCMatch) { std::cout << " MCTrack ID Proton? "<< lambda_PndVtxPRG[j]->Daughter(0)->GetRecoCandidate()->GetMcIndex() << std::endl; std::cout << " MCTrack ID Pionminus? "<< lambda_PndVtxPRG[j]->Daughter(1)->GetRecoCandidate()->GetMcIndex() << std::endl; } std::cout << "Information About Fit: " << std::endl; std::cout << "Chi2? " << chi2 << std::endl; std::cout << "Propability? " << prob << std::endl; } if(MCMatch) { //### Armenteros Podolanski Plot momentumproton = lambdaPndVtxPRGfit->Daughter(0)->GetMomentum(); momentumpionminus = lambdaPndVtxPRGfit->Daughter(1)->GetMomentum(); momentumvectorlambda= momentumproton + momentumpionminus; protontransversalcomponent = momentumproton.Perp(momentumvectorlambda); pionminustransversalcomponent = momentumpionminus.Perp(momentumvectorlambda); protonlongitudinalcomponent = momentumproton*momentumvectorlambda; pionminuslongitudinalcomponent = momentumpionminus*momentumvectorlambda; alphalambda = ((protonlongitudinalcomponent-pionminuslongitudinalcomponent)/(protonlongitudinalcomponent+pionminuslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch->Fill(alphalambda,protontransversalcomponent); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_MCMatch->Fill(alphalambda,protontransversalcomponent); //### fComparison_Lambda_InvMass_Before_After_PndVtxPRGFitter->Fill(lambda_PndVtxPRG[j]->M(),lambdaPndVtxPRGfit->M()); fPndVtxPRGFitter_Lambda_InvMass_MCMatch->Fill( (lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_Lambda_Chi2_MCMatch->Fill(chi2); fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_DecayLength_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdaPndVtxPRGfit->Daughter(0)->M()/(lambdaPndVtxPRGfit->Daughter(0)->P3()+lambdaPndVtxPRGfit->Daughter(1)->P3()).Mag()); fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison->Fill(lambda_PndVtxPRG[j]->P()-lambda_PndVtxPRG[j]->GetMcTruth()->P()); fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison->Fill(lambdaPndVtxPRGfit->P()-lambda_PndVtxPRG[j]->GetMcTruth()->P()); fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison->Fill(lambda_PndVtxPRG[j]->GetMomentum().Angle(lambda_PndVtxPRG[j]->GetMcTruth()->GetMomentum())); fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison->Fill(lambdaPndVtxPRGfit->GetMomentum().Angle(lambda_PndVtxPRG[j]->GetMcTruth()->GetMomentum())); fDeltaAngleVsDecayLength->Fill(lambda_PndVtxPRG[j]->Daughter(0)->GetMcTruth()->GetPosition().Mag() ,lambdaPndVtxPRGfit->GetMomentum().Angle(lambda_PndVtxPRG[j]->GetMcTruth()->GetMomentum())); } else { fPndVtxPRGFitter_Lambda_InvMass_MCMatch_Rejected->Fill( (lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_Lambda_Chi2_MCMatch_Rejected->Fill(chi2); fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_DecayLength_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdaPndVtxPRGfit->Daughter(0)->M()/(lambdaPndVtxPRGfit->Daughter(0)->P3()+lambdaPndVtxPRGfit->Daughter(1)->P3()).Mag()); } if(prob>ProbCut && chi2>0) { fPndVtxPRGFitter_Lambda_InvMass_ProbCut->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_DecayLength_ProbCut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdaPndVtxPRGfit->Daughter(0)->M()/(lambdaPndVtxPRGfit->Daughter(0)->P3()+lambdaPndVtxPRGfit->Daughter(1)->P3()).Mag()); lambda_PndVtxPRGFitter_ProbCut.Append(lambdaPndVtxPRGfit); fReconstructedParticle_MomentaDistributionPionMinus_Lab->Fill(lambdaPndVtxPRGfit->Daughter(1)->GetMomentum().Mag()); //### Armenteros Podolanski Plot momentumproton = lambdaPndVtxPRGfit->Daughter(0)->GetMomentum(); momentumpionminus = lambdaPndVtxPRGfit->Daughter(1)->GetMomentum(); momentumvectorlambda= momentumproton + momentumpionminus; protontransversalcomponent = momentumproton.Perp(momentumvectorlambda); pionminustransversalcomponent = momentumpionminus.Perp(momentumvectorlambda); protonlongitudinalcomponent = momentumproton*momentumvectorlambda; pionminuslongitudinalcomponent = momentumpionminus*momentumvectorlambda; alphalambda = ((protonlongitudinalcomponent-pionminuslongitudinalcomponent)/(protonlongitudinalcomponent+pionminuslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut->Fill(alphalambda,protontransversalcomponent); //### if(MCMatch) { if(lambda_found==true && lambdabar_found==true) { TLorentzVector boostedgesamtvector= lambda_PndVtxPRG[j]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdavector= lambda_PndVtxPRG[j]->GetMcTruth()->P4(); boostedlambdavector.Boost(-betavectorgesamt); TVector3 betavectorlambda = boostedlambdavector.BoostVector(); TLorentzVector lorentzvectorproton_lambdarestframeboost = lambda_PndVtxPRG[j]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorgesamt); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorlambda); TVector3 protonvector = lorentzvectorproton_lambdarestframeboost.Vect(); TVector3 lambda_z_direction = boostedlambdavector.Vect(); TVector3 lambda_y_direction = lambda_z_direction.Cross(TVector3(0,0,1)); TVector3 lambda_x_direction = lambda_y_direction.Cross(lambda_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X_PndVtxPRGFitter->Fill(TMath::Cos(lambda_x_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y_PndVtxPRGFitter->Fill(TMath::Cos(lambda_y_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z_PndVtxPRGFitter->Fill(TMath::Cos(lambda_z_direction.Angle(protonvector))); fReconstructedParticle_AngularDistribution_Lambda_Both_PndVtxPRGFitter->Fill(boostedlambdavector.CosTheta()); //#### fPndVtxPRGFitter_Lambda_Both_VertexRes_X_ProbCut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_Both_VertexRes_Y_ProbCut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_Both_VertexRes_Z_ProbCut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_Both_InvMass_ProbCut_MCMatch->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); //#### } fReconstructedParticle_MomentaDistributionPionMinus_Lab_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(1)->GetMomentum().Mag()); //### Armenteros Podolanski Plot momentumproton = lambdaPndVtxPRGfit->Daughter(0)->GetMomentum(); momentumpionminus = lambdaPndVtxPRGfit->Daughter(1)->GetMomentum(); momentumvectorlambda= momentumproton + momentumpionminus; protontransversalcomponent = momentumproton.Perp(momentumvectorlambda); pionminustransversalcomponent = momentumpionminus.Perp(momentumvectorlambda); protonlongitudinalcomponent = momentumproton*momentumvectorlambda; pionminuslongitudinalcomponent = momentumpionminus*momentumvectorlambda; alphalambda = ((protonlongitudinalcomponent-pionminuslongitudinalcomponent)/(protonlongitudinalcomponent+pionminuslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut_MCMatch->Fill(alphalambda,protontransversalcomponent); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_ProbCut_MCMatch->Fill(alphalambda,protontransversalcomponent); //### TLorentzVector boostedgesamtvector= lambda_PndVtxPRG[j]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdavector= lambda_PndVtxPRG[j]->GetMcTruth()->P4(); boostedlambdavector.Boost(-betavectorgesamt); TVector3 betavectorlambda = boostedlambdavector.BoostVector(); TLorentzVector lorentzvectorproton_lambdarestframeboost = lambda_PndVtxPRG[j]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorgesamt); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorlambda); TVector3 protonvector = lorentzvectorproton_lambdarestframeboost.Vect(); TVector3 lambda_z_direction = boostedlambdavector.Vect(); TVector3 lambda_y_direction = lambda_z_direction.Cross(TVector3(0,0,1)); TVector3 lambda_x_direction = lambda_y_direction.Cross(lambda_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndVtxPRGFitter->Fill(TMath::Cos(lambda_x_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndVtxPRGFitter->Fill(TMath::Cos(lambda_y_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndVtxPRGFitter->Fill(TMath::Cos(lambda_z_direction.Angle(protonvector))); fReconstructedParticle_AngularDistribution_Lambda_PndVtxPRGFitter->Fill(boostedlambdavector.CosTheta()); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); } else { fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndVtxPRGFitter_Lambda_InvMass_ProbCut_Rejected->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_DecayLength_Chi2Cut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdaPndVtxPRGfit->Daughter(0)->M()/(lambdaPndVtxPRGfit->Daughter(0)->P3()+lambdaPndVtxPRGfit->Daughter(1)->P3()).Mag()); } if(chi2<2 && chi2>0) { fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); if(MCMatch) { fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); } else { fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_Rejected->Fill((lambdaPndVtxPRGfit->Daughter(0)->P4()+lambdaPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected->Fill(lambdaPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdaPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } for(int j=0;j < lambdabar_PndVtxPRG.GetLength();j++) { bool MCMatch= theAnalysis->McTruthMatch(lambdabar_PndVtxPRG[j],2,false); PndVtxPRG prgfitter(lambdabar_PndVtxPRG[j]); prgfitter.Fit(); double chi2=prgfitter.GetChi2(); double prob=prgfitter.GetProb(); if(!chi2>0) { continue; } RhoCandidate* lambdabarPndVtxPRGfit=lambdabar_PndVtxPRG[j]->GetFit(); fPndVtxPRGFitter_LambdaBar_InvMass->Fill((lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_Chi2->Fill(chi2); fPndVtxPRGFitter_LambdaBar_Prob->Fill(prob); fPndVtxPRGFitter_LambdaBar_DecayLength->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdabarPndVtxPRGfit->Daughter(0)->M()/(lambdabarPndVtxPRGfit->Daughter(0)->P3()+lambdabarPndVtxPRGfit->Daughter(1)->P3()).Mag()); fPndVtxPRGFitter_LambdaBar_VertexRes_X->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); TVector3 VertexPosition; PndVtxPoca LambdaBarPndVtxPoca_BeforeFit; PndVtxPoca LambdaBarPndVtxPoca_AfterFit; fPndVtxPRGFitter_LambdaBar_PocaComparison->Fill(LambdaBarPndVtxPoca_BeforeFit.GetPocaVtx(VertexPosition,lambdabar_PndVtxPRG[j]),LambdaBarPndVtxPoca_AfterFit.GetPocaVtx(VertexPosition,lambdabarPndVtxPRGfit)); //### Armenteros Podolanski Plot TVector3 momentumantiproton = lambdabarPndVtxPRGfit->Daughter(0)->GetMomentum(); TVector3 momentumpionplus = lambdabarPndVtxPRGfit->Daughter(1)->GetMomentum(); TVector3 momentumvectorlambdabar= momentumantiproton + momentumpionplus; double antiprotontransversalcomponent = momentumantiproton.Perp(momentumvectorlambdabar); double pionplustransversalcomponent = momentumpionplus.Perp(momentumvectorlambdabar); double antiprotonlongitudinalcomponent = momentumantiproton*momentumvectorlambdabar; double pionpluslongitudinalcomponent = momentumpionplus*momentumvectorlambdabar; double alphalambdabar=((pionpluslongitudinalcomponent-antiprotonlongitudinalcomponent)/(antiprotonlongitudinalcomponent+pionpluslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter->Fill(alphalambdabar,antiprotontransversalcomponent); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter->Fill(alphalambdabar,antiprotontransversalcomponent); //### if(MCMatch) { fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch->Fill( (lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch->Fill(chi2); fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdabarPndVtxPRGfit->Daughter(0)->M()/(lambdabarPndVtxPRGfit->Daughter(0)->P3()+lambdabarPndVtxPRGfit->Daughter(1)->P3()).Mag()); fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison->Fill(lambdabar_PndVtxPRG[j]->P()-lambdabar_PndVtxPRG[j]->GetMcTruth()->P()); fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison->Fill(lambdabarPndVtxPRGfit->P()-lambdabar_PndVtxPRG[j]->GetMcTruth()->P()); fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison->Fill(lambdabar_PndVtxPRG[j]->GetMomentum().Angle(lambdabar_PndVtxPRG[j]->GetMcTruth()->GetMomentum())); fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison->Fill(lambdabarPndVtxPRGfit->GetMomentum().Angle(lambdabar_PndVtxPRG[j]->GetMcTruth()->GetMomentum())); //### Armenteros Podolanski Plot momentumantiproton = lambdabarPndVtxPRGfit->Daughter(0)->GetMomentum(); momentumpionplus = lambdabarPndVtxPRGfit->Daughter(1)->GetMomentum(); momentumvectorlambdabar= momentumantiproton + momentumpionplus; antiprotontransversalcomponent = momentumantiproton.Perp(momentumvectorlambdabar); pionplustransversalcomponent = momentumpionplus.Perp(momentumvectorlambdabar); antiprotonlongitudinalcomponent = momentumantiproton*momentumvectorlambdabar; pionpluslongitudinalcomponent = momentumpionplus*momentumvectorlambdabar; alphalambdabar=((pionpluslongitudinalcomponent-antiprotonlongitudinalcomponent)/(antiprotonlongitudinalcomponent+pionpluslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch->Fill(alphalambdabar,antiprotontransversalcomponent); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_MCMatch->Fill(alphalambdabar,antiprotontransversalcomponent); //### } else { fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch_Rejected->Fill( (lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch_Rejected->Fill(chi2); fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdabarPndVtxPRGfit->Daughter(0)->M()/(lambdabarPndVtxPRGfit->Daughter(0)->P3()+lambdabarPndVtxPRGfit->Daughter(1)->P3()).Mag()); } if(prob>ProbCut && chi2>0) { fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut->Fill((lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_LambdaBar_DecayLength_ProbCut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdabarPndVtxPRGfit->Daughter(0)->M()/(lambdabarPndVtxPRGfit->Daughter(0)->P3()+lambdabarPndVtxPRGfit->Daughter(1)->P3()).Mag()); lambdabar_PndVtxPRGFitter_ProbCut.Append(lambdabarPndVtxPRGfit); fReconstructedParticle_MomentaDistributionPionPlus_Lab->Fill(lambdabarPndVtxPRGfit->Daughter(1)->GetMomentum().Mag()); //### Armenteros Podolanski Plot momentumantiproton = lambdabarPndVtxPRGfit->Daughter(0)->GetMomentum(); momentumpionplus = lambdabarPndVtxPRGfit->Daughter(1)->GetMomentum(); momentumvectorlambdabar= momentumantiproton + momentumpionplus; antiprotontransversalcomponent = momentumantiproton.Perp(momentumvectorlambdabar); pionplustransversalcomponent = momentumpionplus.Perp(momentumvectorlambdabar); antiprotonlongitudinalcomponent = momentumantiproton*momentumvectorlambdabar; pionpluslongitudinalcomponent = momentumpionplus*momentumvectorlambdabar; alphalambdabar=((pionpluslongitudinalcomponent-antiprotonlongitudinalcomponent)/(antiprotonlongitudinalcomponent+pionpluslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut->Fill(alphalambdabar,antiprotontransversalcomponent); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut->Fill(alphalambdabar,antiprotontransversalcomponent); //### if(MCMatch) { if(lambda_found==true && lambdabar_found==true) { TLorentzVector boostedgesamtvector= lambdabar_PndVtxPRG[j]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdabarvector= lambdabar_PndVtxPRG[j]->GetMcTruth()->P4(); boostedlambdabarvector.Boost(-betavectorgesamt); TVector3 betavectorlambdabar = boostedlambdabarvector.BoostVector(); TLorentzVector lorentzvectorantiproton_lambdabarrestframeboost = lambdabar_PndVtxPRG[j]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorgesamt); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorlambdabar); TVector3 antiprotonvector = lorentzvectorantiproton_lambdabarrestframeboost.Vect(); TVector3 lambdabar_z_direction = boostedlambdabarvector.Vect(); TVector3 lambdabar_y_direction = lambdabar_z_direction.Cross(TVector3(0,0,1)); TVector3 lambdabar_x_direction = lambdabar_y_direction.Cross(lambdabar_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X_PndVtxPRGFitter->Fill(TMath::Cos(lambdabar_x_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y_PndVtxPRGFitter->Fill(TMath::Cos(lambdabar_y_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z_PndVtxPRGFitter->Fill(TMath::Cos(lambdabar_z_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistribution_LambdaBar_Both_PndVtxPRGFitter->Fill(boostedlambdabarvector.CosTheta()); fPndVtxPRGFitter_LambdaBar_Both_InvMass_ProbCut_MCMatch->Fill( (lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_X_ProbCut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Y_ProbCut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Z_ProbCut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } fReconstructedParticle_MomentaDistributionPionPlus_Lab_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(1)->GetMomentum().Mag()); //### Armenteros Podolanski Plot momentumantiproton = lambdabarPndVtxPRGfit->Daughter(0)->GetMomentum(); momentumpionplus = lambdabarPndVtxPRGfit->Daughter(1)->GetMomentum(); momentumvectorlambdabar= momentumantiproton + momentumpionplus; antiprotontransversalcomponent = momentumantiproton.Perp(momentumvectorlambdabar); pionplustransversalcomponent = momentumpionplus.Perp(momentumvectorlambdabar); antiprotonlongitudinalcomponent = momentumantiproton*momentumvectorlambdabar; pionpluslongitudinalcomponent = momentumpionplus*momentumvectorlambdabar; alphalambdabar=((pionpluslongitudinalcomponent-antiprotonlongitudinalcomponent)/(antiprotonlongitudinalcomponent+pionpluslongitudinalcomponent)); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut_MCMatch->Fill(alphalambdabar,antiprotontransversalcomponent); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut_MCMatch->Fill(alphalambdabar,antiprotontransversalcomponent); //### TLorentzVector boostedgesamtvector= lambdabar_PndVtxPRG[j]->GetMcTruth()->TheMother()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdabarvector= lambdabar_PndVtxPRG[j]->GetMcTruth()->P4(); boostedlambdabarvector.Boost(-betavectorgesamt); TVector3 betavectorlambdabar = boostedlambdabarvector.BoostVector(); TLorentzVector lorentzvectorantiproton_lambdabarrestframeboost = lambdabar_PndVtxPRG[j]->GetMcTruth()->Daughter(0)->P4(); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorgesamt); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorlambdabar); TVector3 antiprotonvector = lorentzvectorantiproton_lambdabarrestframeboost.Vect(); TVector3 lambdabar_z_direction = boostedlambdabarvector.Vect(); TVector3 lambdabar_y_direction = lambdabar_z_direction.Cross(TVector3(0,0,1)); TVector3 lambdabar_x_direction = lambdabar_y_direction.Cross(lambdabar_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndVtxPRGFitter->Fill(TMath::Cos(lambdabar_x_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndVtxPRGFitter->Fill(TMath::Cos(lambdabar_y_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndVtxPRGFitter->Fill(TMath::Cos(lambdabar_z_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistribution_LambdaBar_PndVtxPRGFitter->Fill(boostedlambdabarvector.CosTheta()); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch->Fill( (lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } else { fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected->Fill( (lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_Rejected->Fill((lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } if(chi2<2) { fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut->Fill((lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPndVtxPRGFitter_LambdaBar_DecayLength_Chi2Cut->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Mag()*lambdabarPndVtxPRGfit->Daughter(0)->M()/(lambdabarPndVtxPRGfit->Daughter(0)->P3()+lambdabarPndVtxPRGfit->Daughter(1)->P3()).Mag()); if(MCMatch) { fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch->Fill( (lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } else { fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected->Fill( (lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } else { fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_Rejected->Fill((lambdabarPndVtxPRGfit->Daughter(0)->P4()+lambdabarPndVtxPRGfit->Daughter(1)->P4()).M()); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().X()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().X()); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Y()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected->Fill(lambdabarPndVtxPRGfit->Daughter(0)->GetPosition().Z()-lambdabarPndVtxPRGfit->Daughter(0)->GetMcTruth()->GetPosition().Z()); } } // ********************** // Lambda LambdaBar Nach Cut // ********************** // for(int w=0;wFill(); // ********************** // Pnd4CFitter // ********************** RhoCandList ppbarsystem; RhoCandList* lambda_fitted_list; RhoCandList* lambdabar_fitted_list; RhoCandidate* lambda4cfitted; ppbarsystem.Combine(lambda_PndVtxPRGFitter_ProbCut,lambdabar_PndVtxPRGFitter_ProbCut); ppbarsystem.SetType(88888); for(int b=0;bMcTruthMatch(ppbarsystem[b],2,false); fInvMass_2D_Distribution->Fill(ppbarsystem[b]->Daughter(0)->M(), ppbarsystem[b]->Daughter(1)->M()); Pnd4CFitter fourcfitter(ppbarsystem[b],ini); fourcfitter.FitConserveMasses(); TVector3 boostvector = ppbarsystem[b]->P4().BoostVector(); double Prob=0; Prob=fourcfitter.GetProb(); TLorentzVector deboostedlambda = ppbarsystem[b]->GetFit()->Daughter(0)->P4(); TLorentzVector deboostedlambdabar = ppbarsystem[b]->GetFit()->Daughter(1)->P4(); deboostedlambda.Boost(-boostvector); deboostedlambdabar.Boost(-boostvector); fPnd4CFitter_PbarPSystem_InvMass->Fill(ppbarsystem[b]->GetFit()->M()); fPnd4CFitter_Lambda_InvMass->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->M()); fPnd4CFitter_LambdaBar_InvMass->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->M()); fPnd4CFitter_Chi2->Fill(fourcfitter.GetChi2()); fPnd4CFitter_Lambda_AngularDistribution->Fill(deboostedlambda.CosTheta()); fPnd4CFitter_LambdaBar_AngularDistribution->Fill(deboostedlambdabar.CosTheta()); fPnd4CFitter_Lambda_VertexRes_X->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().X()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().X()); fPnd4CFitter_Lambda_VertexRes_Y->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().Y()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPnd4CFitter_Lambda_VertexRes_Z->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().Z()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPnd4CFitter_LambdaBar_VertexRes_X->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().X()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().X()); fPnd4CFitter_LambdaBar_VertexRes_Y->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().Y()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().Y()); fPnd4CFitter_LambdaBar_VertexRes_Z->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().Z()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().Z()); fPnd4CFitter_Prob->Fill(Prob); if(Prob>ProbCut) { fPnd4CFitter_Lambda_InvMass_ProbCut->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->M()); fPnd4CFitter_LambdaBar_InvMass_ProbCut->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->M()); if(MCMatch) { TLorentzVector boostedgesamtvector= ppbarsystem[b]->GetMcTruth()->P4(); TVector3 betavectorgesamt = boostedgesamtvector.BoostVector(); TLorentzVector boostedlambdavector= ppbarsystem[b]->GetMcTruth()->Daughter(0)->P4(); boostedlambdavector.Boost(-betavectorgesamt); TVector3 betavectorlambda = boostedlambdavector.BoostVector(); TLorentzVector lorentzvectorproton_lambdarestframeboost = ppbarsystem[b]->GetMcTruth()->Daughter(0)->Daughter(0)->P4(); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorgesamt); lorentzvectorproton_lambdarestframeboost.Boost(-betavectorlambda); TVector3 protonvector = lorentzvectorproton_lambdarestframeboost.Vect(); TVector3 lambda_z_direction = boostedlambdavector.Vect(); TVector3 lambda_y_direction = lambda_z_direction.Cross(TVector3(0,0,1)); TVector3 lambda_x_direction = lambda_y_direction.Cross(lambda_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_Pnd4CFitter->Fill(TMath::Cos(lambda_x_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_Pnd4CFitter->Fill(TMath::Cos(lambda_y_direction.Angle(protonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_Pnd4CFitter->Fill(TMath::Cos(lambda_z_direction.Angle(protonvector))); fReconstructedParticle_AngularDistribution_Lambda_Pnd4CFitter->Fill(boostedlambdavector.CosTheta()); TLorentzVector boostedlambdabarvector= ppbarsystem[b]->GetMcTruth()->Daughter(1)->P4(); boostedlambdabarvector.Boost(-betavectorgesamt); TVector3 betavectorlambdabar = boostedlambdabarvector.BoostVector(); TLorentzVector lorentzvectorantiproton_lambdabarrestframeboost =ppbarsystem[b]->GetMcTruth()->Daughter(1)->Daughter(0)->P4(); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorgesamt); lorentzvectorantiproton_lambdabarrestframeboost.Boost(-betavectorlambdabar); TVector3 antiprotonvector = lorentzvectorantiproton_lambdabarrestframeboost.Vect(); TVector3 lambdabar_z_direction = boostedlambdabarvector.Vect(); TVector3 lambdabar_y_direction = lambdabar_z_direction.Cross(TVector3(0,0,1)); TVector3 lambdabar_x_direction = lambdabar_y_direction.Cross(lambdabar_z_direction); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_Pnd4CFitter->Fill(TMath::Cos(lambdabar_x_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_Pnd4CFitter->Fill(TMath::Cos(lambdabar_y_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_Pnd4CFitter->Fill(TMath::Cos(lambdabar_z_direction.Angle(antiprotonvector))); fReconstructedParticle_AngularDistribution_LambdaBar_Pnd4CFitter->Fill(boostedlambdabarvector.CosTheta()); fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->M()); fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->M()); } else { fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->M()); fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->M()); } } else { fPnd4CFitter_Lambda_InvMass_ProbCut_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->M()); fPnd4CFitter_LambdaBar_InvMass_ProbCut_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->M()); } if(MCMatch) { fInvMass_2D_Distribution_MCMatch->Fill(ppbarsystem[b]->Daughter(0)->M(), ppbarsystem[b]->Daughter(1)->M()); fPnd4CFitter_PbarPSystem_InvMass_MCMatch->Fill(ppbarsystem[b]->GetFit()->M()); fPnd4CFitter_Lambda_InvMass_MCMatch->Fill((ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->P4()+ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(1)->P4()).M()); fPnd4CFitter_LambdaBar_InvMass_MCMatch->Fill((ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(0)->P4()+ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->P4()).M()); fPnd4CFitter_Lambda_Chi2_MCMatch->Fill(fourcfitter.GetChi2()); fPnd4CFitter_Lambda_VertexRes_X_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().X()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().X()); fPnd4CFitter_Lambda_VertexRes_Y_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().Y()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPnd4CFitter_Lambda_VertexRes_Z_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().Z()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().X()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().X()); fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().Y()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().Y()); fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().Z()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().Z()); } else { fPnd4CFitter_PbarPSystem_InvMass_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->M()); fPnd4CFitter_Lambda_InvMass_MCMatch_Rejected->Fill((ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->P4()+ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(1)->P4()).M()); fPnd4CFitter_LambdaBar_InvMass_MCMatch_Rejected->Fill((ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(0)->P4()+ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->P4()).M()); fPnd4CFitter_Lambda_Chi2_MCMatch_Rejected->Fill(fourcfitter.GetChi2()); fPnd4CFitter_Lambda_VertexRes_X_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().X()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().X()); fPnd4CFitter_Lambda_VertexRes_Y_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().Y()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().Y()); fPnd4CFitter_Lambda_VertexRes_Z_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(0)->Daughter(0)->GetPosition().Z()-ppbarsystem[b]->Daughter(0)->Daughter(0)->GetMcTruth()->GetPosition().Z()); fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().X()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().X()); fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().Y()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().Y()); fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected->Fill(ppbarsystem[b]->GetFit()->Daughter(1)->Daughter(1)->GetPosition().Z()-ppbarsystem[b]->Daughter(1)->Daughter(1)->GetMcTruth()->GetPosition().Z()); } } if(fVerbose>0) { std::cout<< "============== End PndEventAnalysisLambdaLambdabar::Exec for Event " << eventcounter; std::cout << std::endl; } eventcounter++; } void PndEventAnalysisLambdaLambdabar::Finish() { std::cout << "fErrorCounter PndEventAnalysisLambdaLambdabar = " << fErrorCounter << std::endl; std::cout << "multipleLambdacounter " << multipleLambdacounter << std::endl; fTracksPerEvent->GetXaxis()->SetBinLabel(2,"Proton"); fTracksPerEvent->GetXaxis()->SetBinLabel(3,"Antiproton"); fTracksPerEvent->GetXaxis()->SetBinLabel(4,"Pionplus"); fTracksPerEvent->GetXaxis()->SetBinLabel(5,"Pionminus"); gDirectory->mkdir("EventAnalysis_"+fInputBranchName); gDirectory->cd("EventAnalysis_"+fInputBranchName); fLambda_InvMass->Draw(); fPndVtxPRGFitter_Lambda_InvMass_ProbCut->Draw("SAME"); // No Fitter fLambda_InvMass->Write(); fLambda_InvMass_MCMatch->Write(); fLambda_InvMass_MCMatch_Rejected->Write(); fLambdaBar_InvMass->Write(); fLambdaBar_InvMass_MCMatch->Write(); fLambdaBar_InvMass_MCMatch_Rejected->Write(); fInvMass_2D_Distribution->Write(); fInvMass_2D_Distribution_MCMatch->Write(); fAntiProton_Theta_Phi->Write(); fAntiProton_Theta_Phi_Cut->Write(); fAntiProton_CosTheta_Phi->Write(); fAntiProton_CosTheta_Phi_Cut->Write(); fAntiProton_Theta_Phi_CMS->Write(); fAntiProton_Theta_Phi_Cut_CMS->Write(); fAntiProton_CosTheta_Phi_CMS->Write(); fAntiProton_CosTheta_Phi_Cut_CMS->Write(); fDeltaAngleVsDecayLength->Write(); // PndKinVtxFitter gDirectory->mkdir("PndKinVtxFitter"); gDirectory->cd("PndKinVtxFitter"); fPndKinVtxFitter_Lambda_InvMass->Write(); fPndKinVtxFitter_Lambda_InvMass_MCMatch->Write(); fPndKinVtxFitter_Lambda_InvMass_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_Chi2->Write(); fPndKinVtxFitter_Lambda_Chi2_MCMatch->Write(); fPndKinVtxFitter_Lambda_Chi2_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_Prob->Write(); fPndKinVtxFitter_Lambda_DecayLength->Write(); fPndKinVtxFitter_Lambda_DecayLength_MCMatch->Write(); fPndKinVtxFitter_Lambda_DecayLength_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut->Write(); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_DecayLength_Chi2Cut->Write(); fPndKinVtxFitter_Lambda_InvMass_ProbCut->Write(); fPndKinVtxFitter_Lambda_InvMass_ProbCut_Rejected->Write(); fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch->Write(); fPndKinVtxFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_DecayLength_ProbCut->Write(); fPndKinVtxFitter_LambdaBar_InvMass->Write(); fPndKinVtxFitter_LambdaBar_InvMass_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_InvMass_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_Chi2->Write(); fPndKinVtxFitter_LambdaBar_Chi2_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_Chi2_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_Prob->Write(); fPndKinVtxFitter_LambdaBar_DecayLength->Write(); fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_DecayLength_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut->Write(); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_DecayLength_Chi2Cut->Write(); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut->Write(); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_DecayLength_ProbCut->Write(); fPndKinVtxFitter_Lambda_VertexRes_X->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected->Write(); fPndKinVtxFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected->Write(); fPndKinVtxFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected->Write(); gDirectory->cd("../"); gDirectory->mkdir("Pnd4CFitter"); gDirectory->cd("Pnd4CFitter"); // Pnd4CFitter fPnd4CFitter_Prob->Write(); fPnd4CFitter_Chi2->Write(); fPnd4CFitter_PbarPSystem_InvMass->Write(); fPnd4CFitter_PbarPSystem_InvMass_MCMatch->Write(); fPnd4CFitter_PbarPSystem_InvMass_MCMatch_Rejected->Write(); fPnd4CFitter_Lambda_InvMass->Write(); fPnd4CFitter_Lambda_InvMass_MCMatch->Write(); fPnd4CFitter_Lambda_InvMass_MCMatch_Rejected->Write(); fPnd4CFitter_Lambda_Chi2->Write(); fPnd4CFitter_Lambda_Chi2_MCMatch->Write(); fPnd4CFitter_Lambda_Chi2_MCMatch_Rejected->Write(); fPnd4CFitter_Lambda_InvMass_ProbCut->Write(); fPnd4CFitter_Lambda_InvMass_ProbCut_Rejected->Write(); fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch->Write(); fPnd4CFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected->Write(); fPnd4CFitter_Lambda_VertexRes_X->Write(); fPnd4CFitter_Lambda_VertexRes_Y->Write(); fPnd4CFitter_Lambda_VertexRes_Z->Write(); fPnd4CFitter_Lambda_VertexRes_X_MCMatch->Write(); fPnd4CFitter_Lambda_VertexRes_Y_MCMatch->Write(); fPnd4CFitter_Lambda_VertexRes_Z_MCMatch->Write(); fPnd4CFitter_Lambda_VertexRes_X_MCMatch_Rejected->Write(); fPnd4CFitter_Lambda_VertexRes_Y_MCMatch_Rejected->Write(); fPnd4CFitter_Lambda_VertexRes_Z_MCMatch_Rejected->Write(); fPnd4CFitter_LambdaBar_VertexRes_X->Write(); fPnd4CFitter_LambdaBar_VertexRes_Y->Write(); fPnd4CFitter_LambdaBar_VertexRes_Z->Write(); fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch->Write(); fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch->Write(); fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch->Write(); fPnd4CFitter_LambdaBar_VertexRes_X_MCMatch_Rejected->Write(); fPnd4CFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected->Write(); fPnd4CFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected->Write(); fPnd4CFitter_LambdaBar_InvMass->Write(); fPnd4CFitter_LambdaBar_InvMass_MCMatch->Write(); fPnd4CFitter_LambdaBar_InvMass_MCMatch_Rejected->Write(); fPnd4CFitter_LambdaBar_Chi2->Write(); fPnd4CFitter_LambdaBar_Chi2_MCMatch->Write(); fPnd4CFitter_LambdaBar_Chi2_MCMatch_Rejected->Write(); fPnd4CFitter_LambdaBar_InvMass_ProbCut->Write(); fPnd4CFitter_LambdaBar_InvMass_ProbCut_Rejected->Write(); fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch->Write(); fPnd4CFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected->Write(); fPnd4CFitter_Lambda_AngularDistribution->Write(); fPnd4CFitter_LambdaBar_AngularDistribution->Write(); gDirectory->cd("../"); gDirectory->mkdir("PndVtxPRGFitter"); gDirectory->cd("PndVtxPRGFitter"); // PndVtxPRGFitter fPndVtxPRGFitter_Lambda_InvMass->Write(); fPndVtxPRGFitter_Lambda_InvMass_MCMatch->Write(); fPndVtxPRGFitter_Lambda_InvMass_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_Chi2->Write(); fPndVtxPRGFitter_Lambda_Chi2_MCMatch->Write(); fPndVtxPRGFitter_Lambda_Chi2_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_Prob->Write(); fPndVtxPRGFitter_Lambda_DecayLength->Write(); fPndVtxPRGFitter_Lambda_DecayLength_MCMatch->Write(); fPndVtxPRGFitter_Lambda_DecayLength_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_PocaComparison->Write(); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut->Write(); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_InvMass_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_DecayLength_Chi2Cut->Write(); fPndVtxPRGFitter_Lambda_InvMass_ProbCut->Write(); fPndVtxPRGFitter_Lambda_InvMass_ProbCut_Rejected->Write(); fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_InvMass_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_DecayLength_ProbCut->Write(); fPndVtxPRGFitter_Lambda_Both_InvMass_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_InvMass->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_Chi2->Write(); fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_Chi2_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_Prob->Write(); fPndVtxPRGFitter_LambdaBar_DecayLength->Write(); fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_DecayLength_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_PocaComparison->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_DecayLength_Chi2Cut->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_InvMass_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_DecayLength_ProbCut->Write(); fPndVtxPRGFitter_LambdaBar_Both_InvMass_ProbCut_MCMatch->Write(); fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison->Write(); fReconstructedParticle_Lambda_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison->Write(); fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Prefit_MCComparison->Write(); fReconstructedParticle_LambdaBar_MomentumMag_PndVtxPRGFitter_Postfit_MCComparison->Write(); fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison->Write(); fReconstructedParticle_LambdaBar_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison->Write(); fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Prefit_MCComparison->Write(); fReconstructedParticle_Lambda_MomentumAngle_PndVtxPRGFitter_Postfit_MCComparison->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_X_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Y_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_VertexRes_Z_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_Lambda_Both_VertexRes_X_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_Both_VertexRes_Y_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_Lambda_Both_VertexRes_Z_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_Chi2Cut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_X_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Y_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_VertexRes_Z_ProbCut_MCMatch_Rejected->Write(); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_X_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Y_ProbCut_MCMatch->Write(); fPndVtxPRGFitter_LambdaBar_Both_VertexRes_Z_ProbCut_MCMatch->Write(); gDirectory->cd("../"); gDirectory->mkdir("PndVtxPoca"); gDirectory->cd("PndVtxPoca"); // PndVtxPoca fPndVtxPoca_Lambda_QualityNumber->Write(); fPndVtxPoca_Lambda_VertexPosition->Write(); fPndVtxPoca_Lambda_QualityNumber_MCMatch->Write(); fPndVtxPoca_Lambda_QualityNumber_Rejected->Write(); fPndVtxPoca_Lambda_VertexPosition_MCMatch->Write(); fPndVtxPoca_Lambda_VertexPosition_Rejected->Write(); fPndVtxPoca_LambdaBar_QualityNumber->Write(); fPndVtxPoca_LambdaBar_VertexPosition->Write(); fPndVtxPoca_LambdaBar_QualityNumber_MCMatch->Write(); fPndVtxPoca_LambdaBar_QualityNumber_Rejected->Write(); fPndVtxPoca_LambdaBar_VertexPosition_MCMatch->Write(); fPndVtxPoca_LambdaBar_VertexPosition_Rejected->Write(); gDirectory->cd("../"); gDirectory->mkdir("Other"); gDirectory->cd("Other"); // Other fReconstructedArmenterosPodolanski_PndVtxPRGFitter->Write(); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_MCMatch->Write(); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut->Write(); fReconstructedArmenterosPodolanski_PndVtxPRGFitter_ProbCut_MCMatch->Write(); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter->Write(); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_MCMatch->Write(); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_ProbCut->Write(); fReconstructedArmenterosPodolanski_Lambda_PndVtxPRGFitter_ProbCut_MCMatch->Write(); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter->Write(); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_MCMatch->Write(); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut->Write(); fReconstructedArmenterosPodolanski_LambdaBar_PndVtxPRGFitter_ProbCut_MCMatch->Write(); fKalmanFitterFlag->Write(); fComparison_Lambda_InvMass_Before_After_PndVtxPRGFitter->Write(); fAntiProton_DecayPoint->Write(); fAntiProton_NumberOfDaughters->Write(); fLambdasPerMCTruthMatchStage->Write(); fLambdaBarsPerMCTruthMatchStage->Write(); fWrongAssignedList->Write(); fWrongAssignedList2->Write(); fNumberOfFoundLambda->Write(); fNumberOfFoundLambdaBar->Write(); fLambdaPerStage->Write(); fLambdaBarPerStage->Write(); fCompleteSystemEnergie->Write(); fHitAnalysisTracksPerEvent2D->Write(); fReconstructedParticle_MomentaDistributionPionMinus_CMS->Write(); fReconstructedParticle_MomentaDistributionPionPlus_CMS->Write(); fReconstructedParticle_MomentaDistributionProton_CMS->Write(); fReconstructedParticle_MomentaDistributionAntiProton_CMS->Write(); fReconstructedParticle_MomentaDistributionPionMinus_Lab->Write(); fReconstructedParticle_MomentaDistributionPionPlus_Lab->Write(); fReconstructedParticle_MomentaDistributionProton_Lab->Write(); fReconstructedParticle_MomentaDistributionAntiProton_Lab->Write(); fReconstructedParticle_MomentaDistributionPionMinus_CMS_MCMatch->Write(); fReconstructedParticle_MomentaDistributionPionPlus_CMS_MCMatch->Write(); fReconstructedParticle_MomentaDistributionProton_CMS_MCMatch->Write(); fReconstructedParticle_MomentaDistributionAntiProton_CMS_MCMatch->Write(); fReconstructedParticle_MomentaDistributionPionMinus_Lab_MCMatch->Write(); fReconstructedParticle_MomentaDistributionPionPlus_Lab_MCMatch->Write(); fReconstructedParticle_MomentaDistributionProton_Lab_MCMatch->Write(); fReconstructedParticle_MomentaDistributionAntiProton_Lab_MCMatch->Write(); fReconstructedParticle_AngularDistribution_Lambda->Write(); fReconstructedParticle_AngularDistribution_LambdaBar->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z->Write(); fReconstructedParticle_AngularDistribution_Both_Lambda->Write(); fReconstructedParticle_AngularDistribution_Both_LambdaBar->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z->Write(); fReconstructedParticle_AngularDistribution_Lambda_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistribution_LambdaBar_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndKinVtxFitter->Write(); fReconstructedParticle_AngularDistribution_Lambda_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistribution_LambdaBar_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_Pnd4CFitter->Write(); fReconstructedParticle_AngularDistribution_Lambda_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistribution_LambdaBar_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_X_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_X_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Y_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Y_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Proton_Z_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_AntiProton_Z_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistribution_Lambda_Both_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistribution_LambdaBar_Both_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_X_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_X_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Y_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Y_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_Proton_Z_PndVtxPRGFitter->Write(); fReconstructedParticle_AngularDistributionHyperonRestFrame_Both_AntiProton_Z_PndVtxPRGFitter->Write(); fTracksPerEvent->Write(); gDirectory->cd("/"); } ClassImp(PndEventAnalysisLambdaLambdabar)