plotted delta pT for different cents

[u/mrichter/AliRoot.git] / PWG4 / macros / ConfigAnalysisBtag.C

//------------------------------------\r
// Config file for B-tagging\r
//\r
// Author T. Aronsson\r
//\r
//------------------------------------\r
\r
AliAnaPartCorrMaker*  ConfigAnalysis()\r
{\r
  //\r
  // Configuration goes here\r
  // \r
  printf("======================== \n");\r
  printf("==Preforming Btag.cxx=== \n");\r
  printf("======================== \n");\r
  Bool_t kInputIsESD = kTRUE;     //uncomment for input ESD\r
//Bool_t kInputIsESD = kFALSE;    //uncomment for input AODs\r
  Bool_t kFollowsFilter = kTRUE;  //uncomment if follows ESD filter task\r
//Bool_t kFollowsFilter = kFALSE; //uncomment if no ESD filter task\r
  Bool_t kMC = kTRUE; //set to kFALSE for data\r
\r
   \r
\r
  //Alternatively, select input via anaInputData environment variable.\r
  if (gSystem->Getenv("anaInputData")){\r
    TString kInputData = gSystem->Getenv("anaInputData");\r
    if( kInputData == "AOD" ){\r
      kInputIsESD = kFALSE;\r
      kFollowsFilter = kFALSE;\r
    }\r
  }\r
\r
  //Alternatively, adjust for real data based on kMC value.\r
  if (gSystem->Getenv("anakMC")){\r
    kMC = atoi(gSystem->Getenv("anakMC"));\r
  }\r
\r
  //Detector Fiducial Cuts\r
  AliFiducialCut * fidCut = new AliFiducialCut();\r
  fidCut->DoCTSFiducialCut(kFALSE) ;\r
  fidCut->DoEMCALFiducialCut(kFALSE) ;\r
  fidCut->DoPHOSFiducialCut(kFALSE) ;\r
\r
  //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);\r
  //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);\r
  //fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);\r
\r
  fidCut->Print("");\r
\r
  //-----------------------------------------------------------  \r
  // Reader\r
  //-----------------------------------------------------------\r
  if(kInputIsESD && !kFollowsFilter)AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();\r
  else           AliCaloTrackAODReader *reader = new AliCaloTrackAODReader();\r
  reader->SetDebug(-1);//10 for lots of messages\r
  reader->SwitchOnWriteDeltaAOD();\r
  //Switch on or off the detectors information that you want\r
  reader->SwitchOnEMCAL();\r
  reader->SwitchOnCTS();	\r
  reader->SwitchOffPHOS();\r
	\r
\r
  //Kine\r
  if(kMC && !kInputIsESD){\r
    reader->SwitchOffStack();          // On  by default, remember to SwitchOnMCData() in analysis classes\r
    reader->SwitchOnAODMCParticles();  // Off by default, remember to SwitchOnMCData() in analysis classes\r
  }\r
\r
  //Select Trigger Class for real data\r
  if(!kMC) reader->SetFiredTriggerClassName("CINT1B-ABCE-NOPF-ALL");\r
\r
  //Min particle pT\r
  reader->SetCTSPtMin(0.0);   //new\r
  reader->SetEMCALPtMin(0.0); //new\r
  if(kFollowsFilter)reader->SetTrackStatus(0);  //to prevent automatic TPC and ITS refit\r
\r
\r
\r
  reader->SetFiducialCut(fidCut);\r
\r
  if(!kInputIsESD){\r
\r
  }\r
\r
  reader->Print("");\r
\r
\r
  //Detector Fiducial Cuts\r
  AliFiducialCut * fidCut2 = new AliFiducialCut();\r
  fidCut2->DoEMCALFiducialCut(kTRUE) ;\r
  fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);\r
\r
  fidCut2->DoCTSFiducialCut(kTRUE) ;\r
  fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.); \r
\r
  fidCut2->Print("");\r
\r
\r
\r
\r
\r
  //---------------------------------------------------------------------\r
  // Analysis algorithm\r
  //---------------------------------------------------------------------\r
\r
  AliAnaBtag *btag = new AliAnaBtag();\r
\r
 //Base class\r
  btag->SetDebug(-1); //10 for lots of messages\r
  //btag->SetWriteNtuple(1); //Can be used to write out NTuples for local analysis (1000 times faster than AOD analysis), default is off.\r
  if(kMC){\r
    btag->SwitchOnDataMC();\r
    btag->SetMinPt(1.);\r
  }\r
  btag->SetOutputAODName("Electrons");\r
  btag->SetOutputAODClassName("AliAODPWG4Particle");\r
\r
\r
\r
\r
\r
  //---------------------------------------------------------------------\r
  // Set  analysis algorithm and reader\r
  //---------------------------------------------------------------------\r
  maker = new AliAnaPartCorrMaker();\r
  maker->SetReader(reader);//pointer to reader\r
  maker->AddAnalysis(btag,0);\r
  maker->SetAnaDebug(-1)  ;\r
  maker->SwitchOnHistogramsMaker()  ;\r
  maker->SwitchOnAODsMaker()  ;\r
\r
  maker->Print("");\r
  //\r
  printf("============================ \n");\r
  printf("END ConfigAnalysisElectron() \n");\r
  printf("============================ \n");\r
  return maker ;\r
}\r