change phi range of histograms to accomodate new EMCAL acceptance

[u/mrichter/AliRoot.git] / PWG4 / macros / QA / AddTaskCalorimeterQA.C

diff --git a/PWG4/macros/QA/AddTaskCalorimeterQA.C b/PWG4/macros/QA/AddTaskCalorimeterQA.C
index 6b58ae7b35df241c3ab10783d75a8e9d4b8e55b9..e73234cd166daba9bed3aed4e591a479eb09d2b4 100644 (file)
--- a/PWG4/macros/QA/AddTaskCalorimeterQA.C
+++ b/PWG4/macros/QA/AddTaskCalorimeterQA.C
@@ -40,53 +40,53 @@ AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kP
   reader->SwitchOnPHOSCells(); 
   reader->SwitchOnEMCAL();
   reader->SwitchOnPHOS();
-  reader->SwitchOffCTS();
+  reader->SwitchOnCTS();
   reader->SetEMCALPtMin(0.); 
   reader->SetPHOSPtMin (0.);
   reader->SetCTSPtMin  (0.);
+  reader->SetZvertexCut(10.);
   
   if(kUseKinematics){
-               if(inputDataType == "ESD"){
-                       reader->SwitchOnStack();          
-                       reader->SwitchOffAODMCParticles(); 
-               }
-               else if(inputDataType == "AOD"){
-                       reader->SwitchOffStack();          
-                       reader->SwitchOnAODMCParticles(); 
-               }
+    if(inputDataType == "ESD"){
+      reader->SwitchOnStack();          
+      reader->SwitchOffAODMCParticles(); 
+    }
+    else if(inputDataType == "AOD"){
+      reader->SwitchOffStack();          
+      reader->SwitchOnAODMCParticles(); 
+    }
   }
   //if(!kSimulation) reader->SetFiredTriggerClassName("CINT1B-ABCE-NOPF-ALL");
   reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches.
   reader->SwitchOffWriteDeltaAOD()  ;
   if(oldAOD)         reader->SwitchOnOldAODs();
   if(kPrintSettings) reader->Print("");
-       
+  
   // *** Calorimeters Utils    ***
   AliCalorimeterUtils *cu = new AliCalorimeterUtils;
   // Remove clusters close to borders, at least max energy cell is 1 cell away 
   cu->SetNumberOfCellsFromEMCALBorder(1);
   cu->SetNumberOfCellsFromPHOSBorder(2);
-  cu->SwitchOnNoFiducialBorderInEMCALEta0();
-       
+  cu->SetEMCALGeometryName("EMCAL_COMPLETEV1");  
   // Remove EMCAL hottest channels for first LHC10 periods     
-  cu->SwitchOnBadChannelsRemoval();
+  //cu->SwitchOnBadChannelsRemoval();
   // SM0
-  cu->SetEMCALChannelStatus(0,3,13);  cu->SetEMCALChannelStatus(0,44,1); cu->SetEMCALChannelStatus(0,3,13); 
-  cu->SetEMCALChannelStatus(0,20,7);  cu->SetEMCALChannelStatus(0,38,2);   
+  //cu->SetEMCALChannelStatus(0,3,13);  cu->SetEMCALChannelStatus(0,44,1); cu->SetEMCALChannelStatus(0,3,13); 
+  //cu->SetEMCALChannelStatus(0,20,7);  cu->SetEMCALChannelStatus(0,38,2);   
   // SM1
-  cu->SetEMCALChannelStatus(1,4,7);   cu->SetEMCALChannelStatus(1,4,13);  cu->SetEMCALChannelStatus(1,9,20); 
-  cu->SetEMCALChannelStatus(1,14,15); cu->SetEMCALChannelStatus(1,23,16); cu->SetEMCALChannelStatus(1,32,23); 
-  cu->SetEMCALChannelStatus(1,37,5);  cu->SetEMCALChannelStatus(1,40,1);  cu->SetEMCALChannelStatus(1,40,2);
-  cu->SetEMCALChannelStatus(1,40,5);  cu->SetEMCALChannelStatus(1,41,0);  cu->SetEMCALChannelStatus(1,41,1);
-  cu->SetEMCALChannelStatus(1,41,2);  cu->SetEMCALChannelStatus(1,41,4);
+  //cu->SetEMCALChannelStatus(1,4,7);   cu->SetEMCALChannelStatus(1,4,13);  cu->SetEMCALChannelStatus(1,9,20); 
+  //cu->SetEMCALChannelStatus(1,14,15); cu->SetEMCALChannelStatus(1,23,16); cu->SetEMCALChannelStatus(1,32,23); 
+  //cu->SetEMCALChannelStatus(1,37,5);  cu->SetEMCALChannelStatus(1,40,1);  cu->SetEMCALChannelStatus(1,40,2);
+  //cu->SetEMCALChannelStatus(1,40,5);  cu->SetEMCALChannelStatus(1,41,0);  cu->SetEMCALChannelStatus(1,41,1);
+  //cu->SetEMCALChannelStatus(1,41,2);  cu->SetEMCALChannelStatus(1,41,4);
   // SM2       
-  cu->SetEMCALChannelStatus(2,14,15); cu->SetEMCALChannelStatus(2,18,16); cu->SetEMCALChannelStatus(2,18,17); 
-  cu->SetEMCALChannelStatus(2,18,18); cu->SetEMCALChannelStatus(2,18,20); cu->SetEMCALChannelStatus(2,18,21); 
-  cu->SetEMCALChannelStatus(2,18,23); cu->SetEMCALChannelStatus(2,19,16); cu->SetEMCALChannelStatus(2,19,17); 
-  cu->SetEMCALChannelStatus(2,19,19); cu->SetEMCALChannelStatus(2,19,20); cu->SetEMCALChannelStatus(2,19,21); 
-  cu->SetEMCALChannelStatus(2,19,22);
+  //cu->SetEMCALChannelStatus(2,14,15); cu->SetEMCALChannelStatus(2,18,16); cu->SetEMCALChannelStatus(2,18,17); 
+  //cu->SetEMCALChannelStatus(2,18,18); cu->SetEMCALChannelStatus(2,18,20); cu->SetEMCALChannelStatus(2,18,21); 
+  //cu->SetEMCALChannelStatus(2,18,23); cu->SetEMCALChannelStatus(2,19,16); cu->SetEMCALChannelStatus(2,19,17); 
+  //cu->SetEMCALChannelStatus(2,19,19); cu->SetEMCALChannelStatus(2,19,20); cu->SetEMCALChannelStatus(2,19,21); 
+  //cu->SetEMCALChannelStatus(2,19,22);
   //SM3
-  cu->SetEMCALChannelStatus(3,4,7);
+  //cu->SetEMCALChannelStatus(3,4,7);
        
   //Recalibration
   //cu->SwitchOnRecalibration();
@@ -114,15 +114,15 @@ AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kP
   //emcalQA->SetFiducialCut(fidCut);
   emcalQA->SwitchOffFiducialCut();
   emcalQA->SwitchOffPlotsMaking();
-  emcalQA->SwitchOnCalorimetersCorrelation();
-  if(kUseKinematics)emcalQA->SetTimeCut(420,825);//Open for the moment
+  emcalQA->SwitchOnCorrelation();
+  if(!kUseKinematics)emcalQA->SetTimeCut(400,850);//Open for the moment
   //Set Histrograms bins and ranges
-  emcalQA->SetHistoPtRangeAndNBins(0, 5, 50) ;
-  emcalQA->SetHistoFinePtRangeAndNBins(0, 5, 1000) ; // bining for fhAmpId
-  emcalQA->SetHistoPhiRangeAndNBins(79*TMath::DegToRad(), 121*TMath::DegToRad(), 100) ;
+  emcalQA->SetHistoPtRangeAndNBins(0, 50, 200) ;
+  emcalQA->SetHistoFinePtRangeAndNBins(0, 10, 200) ; // bining for fhAmpId
+  emcalQA->SetHistoPhiRangeAndNBins(79*TMath::DegToRad(), 181*TMath::DegToRad(), 200) ;
   emcalQA->SetHistoEtaRangeAndNBins(-0.71, 0.71, 200) ;
-  emcalQA->SetNumberOfModules(4); //EMCAL first year
-  emcalQA->SetHistoMassRangeAndNBins(0., 0.6, 200) ;
+  emcalQA->SetNumberOfModules(10); 
+  emcalQA->SetHistoMassRangeAndNBins(0., 1, 400) ;
   emcalQA->SetHistoAsymmetryRangeAndNBins(0., 1. , 10 );
   emcalQA->SetHistoPOverERangeAndNBins(0,10.,100);
   emcalQA->SetHistodEdxRangeAndNBins(0.,200.,200);
@@ -130,12 +130,14 @@ AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kP
   emcalQA->SetHistoTimeRangeAndNBins(300.,900,300);
   emcalQA->SetHistoRatioRangeAndNBins(0.,2.,100);
   emcalQA->SetHistoVertexDistRangeAndNBins(0.,500.,500);
-  emcalQA->SetHistoNClusterCellRangeAndNBins(0,50,50);
+  emcalQA->SetHistoNClusterCellRangeAndNBins(0,500,500);
   emcalQA->SetHistoXRangeAndNBins(-230,90,120);
   emcalQA->SetHistoYRangeAndNBins(370,450,40);
   emcalQA->SetHistoZRangeAndNBins(-400,400,200);
   emcalQA->SetHistoRRangeAndNBins(400,450,25);
-  
+  emcalQA->SetHistoV0SignalRangeAndNBins(0,5000,500);
+  emcalQA->SetHistoV0MultiplicityRangeAndNBins(0,5000,500);
+  emcalQA->SetHistoTrackMultiplicityRangeAndNBins(0,5000,500);
   //emcalQA->GetMCAnalysisUtils()->SetDebug(10);
        
   if(kPrintSettings) emcalQA->Print("");       
@@ -151,12 +153,12 @@ AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kP
   //phosQA->GetMCAnalysisUtils()->SetDebug(10);
   phosQA->SwitchOffPlotsMaking();
   //Set Histrograms bins and ranges
-  phosQA->SetHistoPtRangeAndNBins(0, 5, 50) ;
-  phosQA->SetHistoFinePtRangeAndNBins(0, 5, 1000) ; // bining for fhAmpId
+  phosQA->SetHistoPtRangeAndNBins(0, 50, 200) ;
+  phosQA->SetHistoFinePtRangeAndNBins(0, 10, 200) ; // bining for fhAmpId
   phosQA->SetHistoPhiRangeAndNBins(259*TMath::DegToRad(), 321*TMath::DegToRad(), 130) ;
   phosQA->SetHistoEtaRangeAndNBins(-0.125, 0.125, 57) ;
   phosQA->SetNumberOfModules(3); //PHOS first year
-  phosQA->SetHistoMassRangeAndNBins(0., 0.6, 200) ;
+  phosQA->SetHistoMassRangeAndNBins(0., 1., 400) ;
   phosQA->SetHistoAsymmetryRangeAndNBins(0., 1. , 10) ;
   phosQA->SetHistoPOverERangeAndNBins(0,10.,100);
   phosQA->SetHistodEdxRangeAndNBins(0.,200.,200);
@@ -164,14 +166,16 @@ AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kP
   phosQA->SetHistoTimeRangeAndNBins(0.,300,300);
   phosQA->SetHistoRatioRangeAndNBins(0.,2.,100);
   phosQA->SetHistoVertexDistRangeAndNBins(0.,500.,500);
-  phosQA->SetHistoNClusterCellRangeAndNBins(0,50,50);
+  phosQA->SetHistoNClusterCellRangeAndNBins(0,500,500);
   phosQA->SetHistoXRangeAndNBins(-100,400,100);
   phosQA->SetHistoYRangeAndNBins(-490,-290,100);
   phosQA->SetHistoZRangeAndNBins(-80,80,100);
-  phosQA->SetHistoRRangeAndNBins(440,480,80);
-       
-  //if(kPrintSettings)phosQA->Print("");       
-       
+  phosQA->SetHistoRRangeAndNBins(440,480,80);  
+  phosQA->SetHistoV0SignalRangeAndNBins(0,5000,500);
+  phosQA->SetHistoV0MultiplicityRangeAndNBins(0,5000,500);
+  phosQA->SetHistoTrackMultiplicityRangeAndNBins(0,5000,500);
+  
+
   // #### Configure Maker ####
   AliAnaPartCorrMaker * maker = new AliAnaPartCorrMaker();
   maker->SetReader(reader);//pointer to reader