- changes in Conversion & Meson cuts to print full list of cuts during intialization
authorfbock <friederike.bock@cern.ch>
Tue, 1 Jul 2014 13:51:50 +0000 (15:51 +0200)
committerfbock <friederike.bock@cern.ch>
Tue, 1 Jul 2014 13:51:50 +0000 (15:51 +0200)
PWGGA/GammaConv/AliConversionCuts.cxx
PWGGA/GammaConv/AliConversionCuts.h
PWGGA/GammaConv/AliConversionMesonCuts.cxx
PWGGA/GammaConv/AliConversionMesonCuts.h

index 029e956..44f229c 100644 (file)
@@ -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 );
index babfd38..f2b4b88 100644 (file)
@@ -427,7 +427,7 @@ class AliConversionCuts : public AliAnalysisCuts {
    Int_t fNSpecialSubTriggerOptions;
 private:
 
-   ClassDef(AliConversionCuts,8)
+   ClassDef(AliConversionCuts,9)
 };
 
 
index b9927e9..6916dfd 100644 (file)
@@ -664,7 +664,7 @@ Bool_t AliConversionMesonCuts::InitializeCutsFromCutString(const TString analysi
       if(!SetCut(cutIds(ii),fCuts[ii]))return kFALSE;
    }
 
-   //PrintCuts();
+   PrintCutsWithValues();
    return kTRUE;
 }
 ///________________________________________________________________________
@@ -783,6 +783,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
    switch(mesonKind){
index b7f6f3a..e25ac7b 100644 (file)
@@ -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)
 };