PWG4/GammaConv/macros/printCuts.C

   1
   2
   3 using namespace std;
   4 void printCuts(TString str){
   5
   6
   7
   8   TString s = str(0,1);
   9   Int_t goodId= s.Atoi();
  10   s = str(1,1);
  11   Int_t v0FinderType= s.Atoi();
  12   s = str(2,1);
  13   Int_t eProbCut= s.Atoi();
  14   s = str(3,1);
  15   Int_t ededxSigmaCut= s.Atoi();
  16   s = str(4,1);
  17   Int_t pidedxSigmaCut= s.Atoi();
  18   s = str(5,1);
  19   Int_t piMomdedxSigmaCut= s.Atoi();
  20   s = str(6,1);
  21   Int_t chi2GammaCut= s.Atoi();
  22   s = str(7,1);
  23   Int_t singlePtCut= s.Atoi();
  24   s = str(8,1);
  25   Int_t clsTPCCut= s.Atoi();
  26   s = str(9,1);
  27   Int_t etaCut= s.Atoi();
  28
  29   /*  cout<<"etaCut: "<<etaCut<<endl;
  30   cout<<"clsTPCCut: "<<clsTPCCut<<endl;
  31   cout<<"singlePtCut: "<<singlePtCut<<endl;
  32   cout<<"chi2GammaCut: "<<chi2GammaCut<<endl;
  33   cout<<"piMomdedxSigmaCut: "<<piMomdedxSigmaCut<<endl;
  34   cout<<"pidedxSigmaCut: "<<pidedxSigmaCut <<endl;
  35   cout<<"ededxSigmaCut: "<<ededxSigmaCut <<endl;
  36   cout<<"eProbCut: "<< eProbCut<<endl;
  37   cout<<"v0FinderType: "<<v0FinderType <<endl;
  38   cout<<"goodId: "<<goodId <<endl;
  39   */
  40
  41   if(goodId !=9){
  42     cout<<"Analysis Cut Selection too short or does not start with 9"<<endl;
  43     return;
  44   }
  45
  46   switch (v0FinderType){
  47   case 0:  // on fly V0 finder
  48     cout<<"Using the on fly v0 finder"<<endl;
  49     break;
  50   case 1:  // offline V0 finder
  51     cout<<"Using the offline v0 finder"<<endl;
  52     break;
  53   default:
  54     return;
  55   }
  56   switch(eProbCut){
  57   case 0:  // 0.
  58     cout<<"Prob electron is 0.000"<<endl;
  59     //    kGCprobElectron = 0.000;
  60     break;
  61   case 1:  // 0.001
  62     cout<<"Prob electron is 0.001"<<endl;
  63     //    kGCprobElectron = 0.001;
  64     break;
  65   case 2:  // 0.01
  66     cout<<"Prob electron is 0.01"<<endl;
  67     //    kGCprobElectron = 0.01;
  68     break;
  69   default:
  70     return;
  71   }
  72
  73   switch(ededxSigmaCut){
  74   case 0: // -10,10
  75     cout<<"Using sigma dedx [-10,10]"<<endl;
  76     //    kGCPIDnSigmaBelowElectronLine=-10;
  77     //    kGCPIDnSigmaAboveElectronLine=10;
  78     break;
  79   case 1: // -5,5
  80     cout<<"Using sigma dedx [-5,5]"<<endl;
  81     //    kGCPIDnSigmaBelowElectronLine=-5;
  82     //   kGCPIDnSigmaAboveElectronLine=5;
  83     break;
  84   case 2: // -3,5
  85     cout<<"Using sigma dedx [-3,5]"<<endl;
  86     //    kGCPIDnSigmaBelowElectronLine=-3;
  87     //    kGCPIDnSigmaAboveElectronLine=5;
  88     break;
  89   default:
  90     return;
  91   }
  92
  93   switch(pidedxSigmaCut){
  94   case 0:  // -10
  95     cout<<"using pidedxsigmacut: -10"<<endl;
  96     //   kGCPIDnSigmaAbovePionLine=-10;
  97     break;
  98   case 1:   // 0
  99     cout<<"using pidedxsigmacut: 0"<<endl;
 100     //    kGCPIDnSigmaAbovePionLine=0;
 101     break;
 102   case 2:  // 1
 103     cout<<"using pidedxsigmacut: 1"<<endl;
 104     //   kGCPIDnSigmaAbovePionLine=1;
 105     break;
 106   default:
 107     return;
 108   }
 109
 110   switch(piMomdedxSigmaCut){
 111   case 0:  // 0.5 GeV
 112     cout<<"piMomdedxSigmaCut: 0.5"<<endl;
 113     //  kGCPIDMinPnSigmaAbovePionLine=0.5;
 114     break;
 115   case 1:  // 1. GeV
 116     cout<<"piMomdedxSigmaCut: 1"<<endl;
 117     //   kGCPIDMinPnSigmaAbovePionLine=1.;
 118     break;
 119   case 2:  // 1.5 GeV
 120     cout<<"piMomdedxSigmaCut: 1.5"<<endl;
 121     //   kGCPIDMinPnSigmaAbovePionLine=1.5;
 122     break;
 123   default:
 124     return;
 125   }
 126
 127   switch(chi2GammaCut){
 128   case 0: // 100
 129     cout<<"chi2CutConversion = 100."<<endl;
 130     //   kGCchi2CutConversion = 100.;
 131     break;
 132   case 1:  // 50
 133     cout<<"chi2CutConversion = 50."<<endl;
 134     //   kGCchi2CutConversion = 50.;
 135     break;
 136   case 2:  // 30
 137     cout<<"chi2CutConversion = 30."<<endl;
 138     //   kGCchi2CutConversion = 30.;
 139     break;
 140   default:
 141     return;
 142   }
 143
 144   switch(singlePtCut){
 145   case 0: // 0.050 GeV
 146     cout<<"kGCsingleptCut = 0,050"<<endl;
 147     //   kGCsingleptCut = 0.050;
 148     break;
 149   case 1:  // 0.100 GeV
 150     cout<<"kGCsingleptCut = 0,100"<<endl;
 151     //   kGCsingleptCut = 0.100;
 152     break;
 153   case 2:  // 0.150 GeV
 154     cout<<"kGCsingleptCut = 0,150"<<endl;
 155     //   kGCsingleptCut = 0.150;
 156     break;
 157   case 3:  // 0.200 GeV
 158     cout<<"kGCsingleptCut = 0,200"<<endl;
 159     //   kGCsingleptCut = 0.200;
 160     break;
 161   default:
 162     return;
 163  }
 164
 165   switch(clsTPCCut){
 166   case 0: // 0
 167     cout<<"kGCminClsTPCCut = 0"<<endl;
 168     //   kGCminClsTPCCut= 0.;
 169     break;
 170   case 1:  // 70
 171     //   kGCminClsTPCCut= 70.;
 172     cout<<"kGCminClsTPCCut = 70"<<endl;
 173     break;
 174   case 2:  // 80
 175     //   kGCminClsTPCCut= 80.;
 176     cout<<"kGCminClsTPCCut = 80"<<endl;
 177     break;
 178   case 3:  // 100
 179     //   kGCminClsTPCCut= 100.;
 180     cout<<"kGCminClsTPCCut = 100"<<endl;
 181     break;
 182   default:
 183     return;
 184   }
 185
 186   switch(etaCut){
 187   case 0: // 0.9
 188     cout<<"eta 0.9"<<endl;
 189     //   kGCetaCut    = 0.9;
 190     //    kGCLineCutZRSlope = tan(2*atan(exp(-kGCetaCut)));
 191     break;
 192   case 1:  // 1.2
 193     cout<<"eta 1.2"<<endl;
 194     //   kGCetaCut    = 1.2;
 195     //   kGCLineCutZRSlope = tan(2*atan(exp(-kGCetaCut)));
 196     break;
 197   case 2:  // 1.4
 198     cout<<"eta 1.4"<<endl;
 199     //    kGCetaCut    = 1.4;
 200     //   kGCLineCutZRSlope = tan(2*atan(exp(-kGCetaCut)));
 201     break;
 202   default:
 203     return;
 204   }
 205
 206
 207 }