Added AddTaskdNdPt.C

	- AddTask for running dNdPt in the current version

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/macros/AddTaskdNdPt.C b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/macros/AddTaskdNdPt.C
new file mode 100644 (file)
index 0000000..ac0ed92
--- /dev/null
+++ b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/macros/AddTaskdNdPt.C
@@ -0,0 +1,150 @@
+AlidNdPtTask *AddTaskdNdPt()
+{
+       
+       AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+       
+  if (!mgr) {
+    ::Error("AddTaskdNdPt", "No analysis manager found.");
+    return 0;
+  }
+       
+       
+  Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);
+       //  Bool_t hasMC= kFALSE;
+       
+       
+  // Switch off all AliInfo (too much output!!!)
+  AliLog::SetGlobalLogLevel(AliLog::kError);
+  mgr->SetDebugLevel(0);
+       
+  //
+  // Create physics trigger selection class
+  //
+       
+       
+       // AliPhysicsSelection *physTrigSel =  new AliPhysicsSelection();
+       
+  //
+  // Create event cuts
+  //
+  Float_t zvWindow = 30. ;
+       
+  AlidNdPtEventCuts *evtCuts = new AlidNdPtEventCuts("AlidNdPtEventCuts","Event cuts");
+  evtCuts->SetZvRange(-zvWindow,zvWindow);
+  evtCuts->SetMeanXYZv(0.0,0.0,0.0);
+  evtCuts->SetSigmaMeanXYZv(1.0,1.0,10.0);
+  evtCuts->SetTriggerRequired(kTRUE);
+       
+  //
+  // Create geom. acceptance cuts
+  //
+  Float_t etaWindow = 1. ;
+  Float_t ptMin = 0.10 ;
+       
+  AlidNdPtAcceptanceCuts *accCuts = new AlidNdPtAcceptanceCuts("AlidNdPtAcceptanceCuts","Geom. acceptance cuts");
+  accCuts->SetEtaRange(-etaWindow,etaWindow);
+  accCuts->SetPtRange(ptMin,1.e10);
+  
+  // cut out trouble sector(s)
+  Float_t phiMin = TMath::Pi() / 9. * 6.;
+  Float_t phiMax = TMath::Pi() / 9. * 7.;
+  // additional cut to remove one sector
+  AlidNdPtAcceptanceCuts *recCuts = new AlidNdPtAcceptanceCuts("AlidNdPtRecAcceptanceCuts","Geom. acceptance cuts for Recostructed tracks");  
+  recCuts->SetExcludeEtaPhiRange(-1.,0.,phiMin,phiMax);
+  
+  
+  //
+  // Create standard esd track cuts
+  //
+  Int_t cutMode = 200;
+       
+  gROOT->LoadMacro("$ALICE_ROOT/PWGLF/SPECTRA/ChargedHadrons/dNdPt/macros/CreatedNdPtTrackCuts.C");  
+  
+       AliESDtrackCuts* esdTrackCuts = CreatedNdPtTrackCuts(cutMode);
+  if (!esdTrackCuts) {
+    printf("ERROR: esdTrackCuts could not be created\n");
+    return;
+  } else {
+    esdTrackCuts->SetHistogramsOn(kTRUE);
+    esdTrackCuts->SetMaxChi2TPCConstrainedGlobal(36.);
+    esdTrackCuts->SetMaxChi2PerClusterITS(36.);    
+  }
+       
+  //
+  // Create task
+  //
+  AlidNdPtTask *task = new AlidNdPtTask("AlidNdPtTask");
+  task->SetUseMCInfo(hasMC);
+       
+  // trigger selection: MB
+  // task->SelectCollisionCandidates(AliVEvent::kMB);
+       
+       //
+       // set analysis options from the Helper here
+       //
+       AlidNdPtHelper::OutputObject outputObject = AlidNdPtHelper::kAnalysisPbPb;
+       AlidNdPtHelper::AnalysisMode analysisMode = AlidNdPtHelper::kTPCITS ;
+       AlidNdPtHelper::ParticleMode particleMode = AlidNdPtHelper::kAllPart ;
+       
+       
+  //
+  // Create cut analysis object
+  //
+  if(outputObject==AlidNdPtHelper::kAnalysisPbPb){
+    
+               AlidNdPtAnalysisPbPb *fdNdPtAnalysisPbPb = new AlidNdPtAnalysisPbPb("dNdPtAnalysisPbPb","dN/dPt Analysis");
+    fdNdPtAnalysisPbPb->SetEventCuts(evtCuts);
+    fdNdPtAnalysisPbPb->SetAcceptanceCuts(accCuts);    
+    fdNdPtAnalysisPbPb->SetTrackCuts(esdTrackCuts);
+    fdNdPtAnalysisPbPb->SetAnalysisMode(analysisMode); 
+    fdNdPtAnalysisPbPb->SetParticleMode(particleMode); 
+    fdNdPtAnalysisPbPb->SetCentralityEstimator("V0M");      
+    fdNdPtAnalysisPbPb->SetTriggerMask(AliVEvent::kCentral);
+    //fdNdPtAnalysisPbPb->SetTriggerMask(AliVEvent::kEMC1);
+    
+    // cut to remove tpc sector
+    //fdNdPtAnalysisPbPb->SetRecAcceptanceCuts(recCuts);
+    
+               // change binning
+    const Int_t ptNbins = 81;
+    Double_t bins[82] = {0.0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 18.0, 20.0, 22.0, 24.0, 26.0, 28.0, 30.0, 32.0, 34.0, 36.0, 40.0, 45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0, 120.0, 130.0, 140.0, 150.0, 160.0, 180.0, 200.0};
+    Double_t* binsPt = new Double_t[82];
+    for (int i=0; i<82; i++) {binsPt[i] = bins[i];}
+    fdNdPtAnalysisPbPb->SetBinsPt(ptNbins, binsPt);
+    fdNdPtAnalysisPbPb->SetBinsPtCorr(ptNbins, binsPt);  
+               
+    if(hasMC) 
+    {
+                       
+                       fdNdPtAnalysisPbPb->SetUseMCInfo(kTRUE);
+                       fdNdPtAnalysisPbPb->SetHistogramsOn(kTRUE);
+               }
+               else 
+               { 
+                       // online trigger
+                       // fdNdPtAnalysisPbPb->SetPhysicsTriggerSelection(physTrigSel); 
+               }
+    task->AddAnalysisObject( fdNdPtAnalysisPbPb );
+  }
+       
+  // Add task
+  mgr->AddTask(task);
+       
+  // Create containers for input
+  AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
+  
+       
+  TString outputFileName = AliAnalysisManager::GetCommonFileName();
+       
+       AliAnalysisDataContainer *coutput  = mgr->CreateContainer("dNdPtPbPb", 
+                                                                                                                                                                                                                                               TList::Class(),
+                                                                                                                                                                                                                                               AliAnalysisManager::kOutputContainer,   
+                                                                                                                                                                                                                                               Form("%s:dNdPtHistos", mgr->GetCommonFileName()));
+       
+       mgr->ConnectInput(task, 0, cinput);
+       mgr->ConnectOutput(task, 1, coutput);
+       
+       return task;
+       
+}
+