From: fbock Date: Tue, 1 Jul 2014 13:51:50 +0000 (+0200) Subject: - changes in Conversion & Meson cuts to print full list of cuts during intialization X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=e656a821ee9681c8308bcc655befb5a7fa0e54af - changes in Conversion & Meson cuts to print full list of cuts during intialization --- diff --git a/PWGGA/GammaConv/AliConversionCuts.cxx b/PWGGA/GammaConv/AliConversionCuts.cxx index 029e956e908..44f229cd1c9 100644 --- a/PWGGA/GammaConv/AliConversionCuts.cxx +++ b/PWGGA/GammaConv/AliConversionCuts.cxx @@ -2055,6 +2055,13 @@ void AliConversionCuts::PrintCuts() { void AliConversionCuts::PrintCutsWithValues() { // Print out current Cut Selection with value + printf("\nConversion cutnumber \n"); + for(Int_t ic = 0; ic < kNCuts; ic++) { + printf("%d",fCuts[ic]); + } + printf("\n\n"); + + if (fIsHeavyIon == 0) { printf("Running in pp mode \n"); if (fSpecialTrigger == 0){ @@ -2114,18 +2121,25 @@ void AliConversionCuts::PrintCutsWithValues() { else if (fRejectExtraSignals == 1) printf("\t only MB header will be inspected \n"); else if (fRejectExtraSignals > 1) printf("\t special header have been selected \n"); - printf("Electron cuts: \n"); + printf("Electron cuts & Secondary Track Cuts - only track from secondaries enter analysis: \n"); + printf("\t no like sign pairs from V0s \n"); + if (!fUseCorrectedTPCClsInfo) printf("\t # TPC clusters > %3.2f \n", fMinClsTPC); if (fEtaCutMin > -0.1) printf("\t %3.2f < eta_{e} < %3.2f\n", fEtaCutMin, fEtaCut ); else printf("\t eta_{e} < %3.2f\n", fEtaCut ); printf("\t p_{T,e} > %3.2f\n", fSinglePtCut ); - printf("\t %3.2f < n sigma e < %3.2f\n", fPIDnSigmaBelowElectronLine, fPIDnSigmaAboveElectronLine ); + printf("\t TPC refit \n"); + printf("\t no kinks \n"); + printf("\t accept: %3.2f < n sigma_{e,TPC} < %3.2f\n", fPIDnSigmaBelowElectronLine, fPIDnSigmaAboveElectronLine ); + printf("\t reject: %3.2f < p_{e,T} < %3.2f, n sigma_{pi,TPC} < %3.2f\n", fPIDMinPnSigmaAbovePionLine, fPIDMaxPnSigmaAbovePionLine, fPIDnSigmaAbovePionLine ); + printf("\t reject: p_{e,T} > %3.2f, n sigma_{pi,TPC} < %3.2f\n", fPIDMaxPnSigmaAbovePionLine, fPIDnSigmaAbovePionLineHighPt ); + if (fDoPionRejectionLowP) printf("\t reject: p_{e,T} < %3.2f, -%3.2f < n sigma_{pi,TPC} < %3.2f\n", fPIDMinPPionRejectionLowP, fPIDnSigmaAtLowPAroundPionLine, fPIDnSigmaAtLowPAroundPionLine ); + if (fDoKaonRejectionLowP) printf("\t reject: -%3.2f < n sigma_{K,TPC} < %3.2f\n", fPIDnSigmaAtLowPAroundKaonLine, fPIDnSigmaAtLowPAroundKaonLine ); + if (fDoProtonRejectionLowP) printf("\t reject: -%3.2f < n sigma_{K,TPC} < %3.2f\n", fPIDnSigmaAtLowPAroundProtonLine, fPIDnSigmaAtLowPAroundProtonLine ); + if (fUseTOFpid) printf("\t accept: %3.2f < n sigma_{e,TOF} < %3.2f\n", fTofPIDnSigmaBelowElectronLine, fTofPIDnSigmaAboveElectronLine); printf("Photon cuts: \n"); - printf("\t %3.2f < R_{conv} < %3.2f\n", fMinR, fMaxR ); - printf("\t Z_{conv} < %3.2f\n", fMaxZ ); - if (fEtaCutMin > -0.1) printf("\t %3.2f < eta_{conv} < %3.2f\n", fEtaCutMin, fEtaCut ); - else printf("\t eta_{conv} < %3.2f\n", fEtaCut ); - printf("\t p_{T,gamma} > %3.2f\n", fPtCut ); + if (fUseOnFlyV0Finder) printf("\t using Onfly V0 finder \n"); + else printf("\t using Offline V0 finder \n"); if (fDo2DQt){ printf("\t 2 dimensional q_{T} cut applied with maximum of %3.2f \n", fQtMax ); } else { @@ -2136,7 +2150,14 @@ void AliConversionCuts::PrintCutsWithValues() { } else { printf("\t chi^{2} max cut chi^{2} < %3.2f \n", fChi2CutConversion ); printf("\t psi_{pair} max cut |psi_{pair}| < %3.2f \n", fPsiPairCut ); - } + } + printf("\t %3.2f < R_{conv} < %3.2f\n", fMinR, fMaxR ); + printf("\t Z_{conv} < %3.2f\n", fMaxZ ); + if (fEtaCutMin > -0.1) printf("\t %3.2f < eta_{conv} < %3.2f\n", fEtaCutMin, fEtaCut ); + else printf("\t eta_{conv} < %3.2f\n", fEtaCut ); + if (fDoPhotonAsymmetryCut) printf("\t for p_{T,track} > %3.2f, A_{gamma} < %3.2f \n", fMinPPhotonAsymmetryCut, fMinPhotonAsymmetry ); + if (fUseCorrectedTPCClsInfo) printf("\t #cluster TPC/ #findable clusters TPC (corrected for radius) > %3.2f\n", fMinClsTPCToF ); + printf("\t p_{T,gamma} > %3.2f\n", fPtCut ); printf("\t cos(Theta_{point}) > %3.2f \n", fCosPAngleCut ); printf("\t dca_{R} < %3.2f \n", fDCARPrimVtxCut ); printf("\t dca_{Z} < %3.2f \n", fDCAZPrimVtxCut ); diff --git a/PWGGA/GammaConv/AliConversionCuts.h b/PWGGA/GammaConv/AliConversionCuts.h index babfd387067..f2b4b88c339 100644 --- a/PWGGA/GammaConv/AliConversionCuts.h +++ b/PWGGA/GammaConv/AliConversionCuts.h @@ -427,7 +427,7 @@ class AliConversionCuts : public AliAnalysisCuts { Int_t fNSpecialSubTriggerOptions; private: - ClassDef(AliConversionCuts,8) + ClassDef(AliConversionCuts,9) }; diff --git a/PWGGA/GammaConv/AliConversionMesonCuts.cxx b/PWGGA/GammaConv/AliConversionMesonCuts.cxx index b9927e95419..6916dfdd4ac 100644 --- a/PWGGA/GammaConv/AliConversionMesonCuts.cxx +++ b/PWGGA/GammaConv/AliConversionMesonCuts.cxx @@ -664,7 +664,7 @@ Bool_t AliConversionMesonCuts::InitializeCutsFromCutString(const TString analysi if(!SetCut(cutIds(ii),fCuts[ii]))return kFALSE; } - //PrintCuts(); + PrintCutsWithValues(); return kTRUE; } ///________________________________________________________________________ @@ -782,6 +782,37 @@ void AliConversionMesonCuts::PrintCuts() { } } +///________________________________________________________________________ +void AliConversionMesonCuts::PrintCutsWithValues() { + // Print out current Cut Selection with values + printf("\nMeson cutnumber \n"); + for(Int_t ic = 0; ic < kNCuts; ic++) { + printf("%d",fCuts[ic]); + } + printf("\n\n"); + + printf("Meson cuts \n"); + printf("\t |y| < %3.2f \n", fRapidityCutMeson); + printf("\t theta_{open} < %3.2f\n", fOpeningAngle); + if (!fAlphaPtDepCut) printf("\t %3.2f < alpha < %3.2f\n", fAlphaMinCutMeson, fAlphaCutMeson); + printf("\t dca_{gamma,gamma} > %3.2f\n", fDCAGammaGammaCut); + printf("\t dca_{R, prim Vtx} > %3.2f\n", fDCARMesonPrimVtxCut); + printf("\t dca_{Z, prim Vtx} > %3.2f\n\n", fDCAZMesonPrimVtxCut); + + printf("Meson BG settings \n"); + if (!fDoBG){ + if (!fUseRotationMethodInBG & !fUseTrackMultiplicityForBG) printf("\t BG scheme: mixing V0 mult \n"); + if (!fUseRotationMethodInBG & fUseTrackMultiplicityForBG) printf("\t BG scheme: mixing track mult \n"); + if (fUseRotationMethodInBG )printf("\t BG scheme: rotation \n"); + if (fdoBGProbability) printf("\t -> use BG probability \n"); + if (fBackgroundHandler) printf("\t -> use new BG handler \n"); + printf("\t depth of pool: %d\n", fNumberOfBGEvents); + if (fUseRotationMethodInBG )printf("\t degree's for BG rotation: %d\n", fnDegreeRotationPMForBG); + } + +} + + ///________________________________________________________________________ Bool_t AliConversionMesonCuts::SetMesonKind(Int_t mesonKind){ // Set Cut diff --git a/PWGGA/GammaConv/AliConversionMesonCuts.h b/PWGGA/GammaConv/AliConversionMesonCuts.h index b7f6f3a374c..e25ac7bd2e4 100644 --- a/PWGGA/GammaConv/AliConversionMesonCuts.h +++ b/PWGGA/GammaConv/AliConversionMesonCuts.h @@ -83,7 +83,8 @@ class AliConversionMesonCuts : public AliAnalysisCuts { Bool_t MesonIsSelectedMCEtaPiPlPiMiGamma(TParticle *fMCMother,AliStack *fMCStack, Int_t &labelNegPion, Int_t &labelPosPion, Int_t &labelGamma, Double_t fRapidityShift=0); Bool_t MesonIsSelectedMCChiC(TParticle *fMCMother,AliStack *fMCStack, Int_t &, Int_t &, Int_t &, Double_t fRapidityShift=0. ); void PrintCuts(); - + void PrintCutsWithValues(); + void InitCutHistograms(TString name="",Bool_t additionalHists=kFALSE); void SetFillCutHistograms(TString name=""){if(!fHistograms){InitCutHistograms(name);};} TList *GetCutHistograms(){return fHistograms;} @@ -169,7 +170,7 @@ class AliConversionMesonCuts : public AliAnalysisCuts { private: - ClassDef(AliConversionMesonCuts,4) + ClassDef(AliConversionMesonCuts,5) };