ROOT logo
AliAnalysisTaskHFE* ConfigHFECalstandard_PbPb(Bool_t useMC, int TPCclust){
  //
  // HFE standard task configuration
  //

  cout << "==== TPC cluster cut" << TPCclust << endl;

  AliHFEcuts *hfecuts = new AliHFEcuts("hfeCuts","HFE Standard Cuts for EMCal");
  hfecuts->CreateStandardCuts();
  //hfecuts->SetMinNClustersTPC(120);
  hfecuts->SetMinNClustersTPC(TPCclust);
  hfecuts->SetMinRatioTPCclusters(0.6);
  hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);

  //hfecuts->SetMinNClustersITS(4);
  hfecuts->SetMinNClustersITS(3);
  hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny);
  hfecuts->SetCheckITSLayerStatus(kFALSE);

  //hfecuts->UnsetVertexRequirement();
  hfecuts->SetMaxImpactParam(3.,3.);
  //hfecuts->SetPtRange(2.0,60.0);
  hfecuts->SetPtRange(2.0,60.0);

  hfecuts->SetVertexRange(10.);
  //hfecuts->SetMaxChi2perClusterITS(36);
  //hfecuts->SetSigmaToVertex(10);
  //hfecuts->SetTOFPIDStep(kTRUE);
  //hfecuts->SetTOFMISMATCHStep(kTRUE);
  //hfecuts->SetTPCPIDCleanUpStep(kTRUE);
  hfecuts->SetQAOn();
  //hfecuts->SetMinNTrackletsTRD(5);

  AliAnalysisTaskHFE *task = new AliAnalysisTaskHFE("HFEanalysisStandardEMCal");
  task->SetHFECuts(hfecuts);
  task->SetPbPbAnalysis(kTRUE);
  //task->SetRemovePileUp(kTRUE);
  task->GetPIDQAManager()->SetHighResolutionHistos();

  // Define Variables
  AliHFEvarManager *vm = task->GetVarManager();
  //vm->AddVariable("pt");
  //vm->AddVariable("eta");
  
  const Double_t ptbinning[50] = {1., 1.1, 1.2, 1.3, 1.4, 
                                  1.5, 1.75, 2., 2.25, 2.5, 
                                  2.75, 3., 3.5, 4., 4.5, 
                                  5., 5.5, 6., 7., 8., 
                                  9., 10., 11., 12., 13., 
                                  14., 15., 16., 17., 18.,
                                  20., 22., 24., 26., 28.,
                                  30., 32., 34., 36., 38.,
				  40., 45., 50., 55., 60.,
			          65., 70., 80., 90., 100}; 

  const Double_t etabinning[17] = {-0.8,-0.7,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8}; 

  vm->AddVariable("pt", 49, ptbinning);
  vm->AddVariable("eta", 16, etabinning);
  vm->AddVariable("phi");
  vm->AddVariable("charge");
  vm->AddVariable("source");
  vm->AddVariable("centrality");
  /*
  if(!useMC){

      for(Int_t a=0;a<12;a++)
      {
	  TF1 *hBackground = new TF1("hadronicBackgroundFunction","TMath::Exp([0]/x + [1])", 0., 20.);
	  hBackground->SetParameter(0, -43.87);
	  hBackground->SetParameter(1, 2.85);
	  task->SetBackGroundFactorsFunction(hBackground,a);
      }


  }
  */
  // Define PID
  AliHFEpid *pid = task->GetPID();
  if(useMC) pid->SetHasMCData(kTRUE);
  pid->AddDetector("EMCAL", 1);
  pid->AddDetector("TPC", 0);
//  pid->ConfigureTPCrejection();

  if(!useMC){

      Double_t params_centr_0_5[1];
      Double_t params_centr_5_10[1];
      Double_t params_centr_10_20[1];
      Double_t params_centr_20_30[1];
      Double_t params_centr_per[1];
      //params_centr_0_5[0]=0.16;  // cut tuned for 0-10%
      //params_centr_5_10[0]=0.16; // cut tuned for 0-10%
      //params_centr_10_20[0]=0.29;
      //params_centr_20_30[0]=0.38;
      //params_centr_per[0]=0.44;
      params_centr_0_5[0] = -1.0;  // cut tuned for 0-10%
      params_centr_5_10[0] = -1.0; // cut tuned for 0-10%
      params_centr_10_20[0] = -1.0;
      params_centr_20_30[0] = -1.0;
      params_centr_per[0] = -1.0;
      char *cutmodel;
      cutmodel="pol0";


      for(Int_t a=0;a<11;a++)
      {
	  if(a>3)  pid->ConfigureTPCcentralityCut(a,cutmodel,params_centr_per,3.0);
	  if(a==0) pid->ConfigureTPCcentralityCut(a,cutmodel,params_centr_0_5,3.0);    //  0-5%
	  if(a==1) pid->ConfigureTPCcentralityCut(a,cutmodel,params_centr_5_10,3.0);    //  5-10%
	  if(a==2) pid->ConfigureTPCcentralityCut(a,cutmodel,params_centr_10_20,3.0);    //  10-20%
	  if(a==3) pid->ConfigureTPCcentralityCut(a,cutmodel,params_centr_20_30,3.0);    //  20-30%
      }


  }


  // V0 tagged tracks
  AliHFEcuts *v0trackCuts = new AliHFEcuts("V0trackCuts", "Track Cuts for tagged track Analysis");
  v0trackCuts->CreateStandardCuts();
  //v0trackCuts->SetMinNClustersTPC(120);
  v0trackCuts->SetMinNClustersTPC(TPCclust);
  v0trackCuts->SetMinRatioTPCclusters(0.6);
  v0trackCuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);
  v0trackCuts->SetMinNClustersITS(1);
  v0trackCuts->SetCutITSpixel(AliHFEextraCuts::kAny);
  v0trackCuts->SetCheckITSLayerStatus(kFALSE);
  v0trackCuts->UnsetVertexRequirement();
  //v0trackCuts->SetMaxChi2perClusterITS(36);
  //hfecuts->SetSigmaToVertex(10);
  v0trackCuts->SetTOFPIDStep(kTRUE);
//  v0trackCuts->SetTOFMISMATCHStep(kTRUE);
  //v0trackCuts->SetTPCPIDCleanUpStep(kTRUE);
  v0trackCuts->SetQAOn();
  
  task->SwitchOnPlugin(AliAnalysisTaskHFE::kTaggedTrackAnalysis);
  task->SetTaggedTrackCuts(v0trackCuts);
  task->SetCleanTaggedTrack(kFALSE);
  
  // change E/p cuts
  AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);
  emcpid->SetEoPMax(1.3);
  emcpid->SetEoPMim(0.9);

  // QA
  task->SetQAOn(AliAnalysisTaskHFE::kPIDqa);
//  task->SetFillSignalOnly(kFALSE);    // for DE pluging for MC
  task->SetQAOn(AliAnalysisTaskHFE::kMCqa);
  //task->SwitchOnPlugin(AliAnalysisTaskHFE::kIsElecBackGround);
  //task->SwitchOnPlugin(AliAnalysisTaskHFE::kSecVtx);
  task->SwitchOnPlugin(AliAnalysisTaskHFE::kDEstep);



  printf("*************************************\n");
  printf("Configuring standard Task:\n");
  task->Print();
  pid->PrintStatus();
  printf("*************************************\n"); 
  return task;
}
 ConfigHFECalstandard_PbPb.C:1
 ConfigHFECalstandard_PbPb.C:2
 ConfigHFECalstandard_PbPb.C:3
 ConfigHFECalstandard_PbPb.C:4
 ConfigHFECalstandard_PbPb.C:5
 ConfigHFECalstandard_PbPb.C:6
 ConfigHFECalstandard_PbPb.C:7
 ConfigHFECalstandard_PbPb.C:8
 ConfigHFECalstandard_PbPb.C:9
 ConfigHFECalstandard_PbPb.C:10
 ConfigHFECalstandard_PbPb.C:11
 ConfigHFECalstandard_PbPb.C:12
 ConfigHFECalstandard_PbPb.C:13
 ConfigHFECalstandard_PbPb.C:14
 ConfigHFECalstandard_PbPb.C:15
 ConfigHFECalstandard_PbPb.C:16
 ConfigHFECalstandard_PbPb.C:17
 ConfigHFECalstandard_PbPb.C:18
 ConfigHFECalstandard_PbPb.C:19
 ConfigHFECalstandard_PbPb.C:20
 ConfigHFECalstandard_PbPb.C:21
 ConfigHFECalstandard_PbPb.C:22
 ConfigHFECalstandard_PbPb.C:23
 ConfigHFECalstandard_PbPb.C:24
 ConfigHFECalstandard_PbPb.C:25
 ConfigHFECalstandard_PbPb.C:26
 ConfigHFECalstandard_PbPb.C:27
 ConfigHFECalstandard_PbPb.C:28
 ConfigHFECalstandard_PbPb.C:29
 ConfigHFECalstandard_PbPb.C:30
 ConfigHFECalstandard_PbPb.C:31
 ConfigHFECalstandard_PbPb.C:32
 ConfigHFECalstandard_PbPb.C:33
 ConfigHFECalstandard_PbPb.C:34
 ConfigHFECalstandard_PbPb.C:35
 ConfigHFECalstandard_PbPb.C:36
 ConfigHFECalstandard_PbPb.C:37
 ConfigHFECalstandard_PbPb.C:38
 ConfigHFECalstandard_PbPb.C:39
 ConfigHFECalstandard_PbPb.C:40
 ConfigHFECalstandard_PbPb.C:41
 ConfigHFECalstandard_PbPb.C:42
 ConfigHFECalstandard_PbPb.C:43
 ConfigHFECalstandard_PbPb.C:44
 ConfigHFECalstandard_PbPb.C:45
 ConfigHFECalstandard_PbPb.C:46
 ConfigHFECalstandard_PbPb.C:47
 ConfigHFECalstandard_PbPb.C:48
 ConfigHFECalstandard_PbPb.C:49
 ConfigHFECalstandard_PbPb.C:50
 ConfigHFECalstandard_PbPb.C:51
 ConfigHFECalstandard_PbPb.C:52
 ConfigHFECalstandard_PbPb.C:53
 ConfigHFECalstandard_PbPb.C:54
 ConfigHFECalstandard_PbPb.C:55
 ConfigHFECalstandard_PbPb.C:56
 ConfigHFECalstandard_PbPb.C:57
 ConfigHFECalstandard_PbPb.C:58
 ConfigHFECalstandard_PbPb.C:59
 ConfigHFECalstandard_PbPb.C:60
 ConfigHFECalstandard_PbPb.C:61
 ConfigHFECalstandard_PbPb.C:62
 ConfigHFECalstandard_PbPb.C:63
 ConfigHFECalstandard_PbPb.C:64
 ConfigHFECalstandard_PbPb.C:65
 ConfigHFECalstandard_PbPb.C:66
 ConfigHFECalstandard_PbPb.C:67
 ConfigHFECalstandard_PbPb.C:68
 ConfigHFECalstandard_PbPb.C:69
 ConfigHFECalstandard_PbPb.C:70
 ConfigHFECalstandard_PbPb.C:71
 ConfigHFECalstandard_PbPb.C:72
 ConfigHFECalstandard_PbPb.C:73
 ConfigHFECalstandard_PbPb.C:74
 ConfigHFECalstandard_PbPb.C:75
 ConfigHFECalstandard_PbPb.C:76
 ConfigHFECalstandard_PbPb.C:77
 ConfigHFECalstandard_PbPb.C:78
 ConfigHFECalstandard_PbPb.C:79
 ConfigHFECalstandard_PbPb.C:80
 ConfigHFECalstandard_PbPb.C:81
 ConfigHFECalstandard_PbPb.C:82
 ConfigHFECalstandard_PbPb.C:83
 ConfigHFECalstandard_PbPb.C:84
 ConfigHFECalstandard_PbPb.C:85
 ConfigHFECalstandard_PbPb.C:86
 ConfigHFECalstandard_PbPb.C:87
 ConfigHFECalstandard_PbPb.C:88
 ConfigHFECalstandard_PbPb.C:89
 ConfigHFECalstandard_PbPb.C:90
 ConfigHFECalstandard_PbPb.C:91
 ConfigHFECalstandard_PbPb.C:92
 ConfigHFECalstandard_PbPb.C:93
 ConfigHFECalstandard_PbPb.C:94
 ConfigHFECalstandard_PbPb.C:95
 ConfigHFECalstandard_PbPb.C:96
 ConfigHFECalstandard_PbPb.C:97
 ConfigHFECalstandard_PbPb.C:98
 ConfigHFECalstandard_PbPb.C:99
 ConfigHFECalstandard_PbPb.C:100
 ConfigHFECalstandard_PbPb.C:101
 ConfigHFECalstandard_PbPb.C:102
 ConfigHFECalstandard_PbPb.C:103
 ConfigHFECalstandard_PbPb.C:104
 ConfigHFECalstandard_PbPb.C:105
 ConfigHFECalstandard_PbPb.C:106
 ConfigHFECalstandard_PbPb.C:107
 ConfigHFECalstandard_PbPb.C:108
 ConfigHFECalstandard_PbPb.C:109
 ConfigHFECalstandard_PbPb.C:110
 ConfigHFECalstandard_PbPb.C:111
 ConfigHFECalstandard_PbPb.C:112
 ConfigHFECalstandard_PbPb.C:113
 ConfigHFECalstandard_PbPb.C:114
 ConfigHFECalstandard_PbPb.C:115
 ConfigHFECalstandard_PbPb.C:116
 ConfigHFECalstandard_PbPb.C:117
 ConfigHFECalstandard_PbPb.C:118
 ConfigHFECalstandard_PbPb.C:119
 ConfigHFECalstandard_PbPb.C:120
 ConfigHFECalstandard_PbPb.C:121
 ConfigHFECalstandard_PbPb.C:122
 ConfigHFECalstandard_PbPb.C:123
 ConfigHFECalstandard_PbPb.C:124
 ConfigHFECalstandard_PbPb.C:125
 ConfigHFECalstandard_PbPb.C:126
 ConfigHFECalstandard_PbPb.C:127
 ConfigHFECalstandard_PbPb.C:128
 ConfigHFECalstandard_PbPb.C:129
 ConfigHFECalstandard_PbPb.C:130
 ConfigHFECalstandard_PbPb.C:131
 ConfigHFECalstandard_PbPb.C:132
 ConfigHFECalstandard_PbPb.C:133
 ConfigHFECalstandard_PbPb.C:134
 ConfigHFECalstandard_PbPb.C:135
 ConfigHFECalstandard_PbPb.C:136
 ConfigHFECalstandard_PbPb.C:137
 ConfigHFECalstandard_PbPb.C:138
 ConfigHFECalstandard_PbPb.C:139
 ConfigHFECalstandard_PbPb.C:140
 ConfigHFECalstandard_PbPb.C:141
 ConfigHFECalstandard_PbPb.C:142
 ConfigHFECalstandard_PbPb.C:143
 ConfigHFECalstandard_PbPb.C:144
 ConfigHFECalstandard_PbPb.C:145
 ConfigHFECalstandard_PbPb.C:146
 ConfigHFECalstandard_PbPb.C:147
 ConfigHFECalstandard_PbPb.C:148
 ConfigHFECalstandard_PbPb.C:149
 ConfigHFECalstandard_PbPb.C:150
 ConfigHFECalstandard_PbPb.C:151
 ConfigHFECalstandard_PbPb.C:152
 ConfigHFECalstandard_PbPb.C:153
 ConfigHFECalstandard_PbPb.C:154
 ConfigHFECalstandard_PbPb.C:155
 ConfigHFECalstandard_PbPb.C:156
 ConfigHFECalstandard_PbPb.C:157
 ConfigHFECalstandard_PbPb.C:158
 ConfigHFECalstandard_PbPb.C:159
 ConfigHFECalstandard_PbPb.C:160
 ConfigHFECalstandard_PbPb.C:161
 ConfigHFECalstandard_PbPb.C:162
 ConfigHFECalstandard_PbPb.C:163