open in the negative side default time cut in QA and Photon analysis, add time contro...
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 7 Nov 2011 23:06:56 +0000 (23:06 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 7 Nov 2011 23:06:56 +0000 (23:06 +0000)
PWG4/PartCorrDep/AliAnaCalorimeterQA.cxx
PWG4/PartCorrDep/AliAnaPhoton.cxx
PWG4/PartCorrDep/AliAnaPhoton.h

index 7cc20f4..b028198 100755 (executable)
@@ -67,7 +67,7 @@ fStudyClustersAsymmetry(kFALSE),       fStudyWeight(kFALSE),
 //Parameters and cuts
 fNModules(12),                         fNRCU(2),
 fNMaxCols(48),                         fNMaxRows(24),  
-fTimeCutMin(-1),                       fTimeCutMax(9999999),
+fTimeCutMin(-10000),                   fTimeCutMax(10000),
 fEMCALCellAmpMin(0),                   fPHOSCellAmpMin(0), 
 
 //Histograms
@@ -376,7 +376,7 @@ void AliAnaCalorimeterQA::CellHistograms(AliVCaloCells *cells)
       
       //Transform time to ns
       time *= 1.0e9;
-      
       // Remove noisy channels, only possible in ESDs
       if(GetReader()->GetDataType() == AliCaloTrackReader::kESD){
         if(time < fTimeCutMin || time > fTimeCutMax){
@@ -384,7 +384,7 @@ void AliAnaCalorimeterQA::CellHistograms(AliVCaloCells *cells)
           continue;
         }
       }      
-      
+
       // Remove exotic cells
       fhCellECross->Fill(amp,1-GetECross(id,cells)/amp);
       
index 7c91391..857c7ae 100755 (executable)
@@ -52,12 +52,12 @@ ClassImp(AliAnaPhoton)
 AliAnaPhoton::AliAnaPhoton() : 
     AliAnaPartCorrBaseClass(),    fCalorimeter(""), 
     fMinDist(0.),                 fMinDist2(0.),                fMinDist3(0.), 
-    fRejectTrackMatch(0),         fTimeCutMin(-1),              fTimeCutMax(999999),         
+    fRejectTrackMatch(0),         fTimeCutMin(-10000),          fTimeCutMax(10000),         
     fNCellsCut(0),                fFillSSHistograms(kFALSE),    
     fNOriginHistograms(8),        fNPrimaryHistograms(4),
 
     // Histograms
-    fhNCellsE(0),                 fhMaxCellDiffClusterE(0),
+    fhNCellsE(0),                 fhMaxCellDiffClusterE(0),     fhTimeE(0), // Control histograms
     fhEPhoton(0),                 fhPtPhoton(0),  
     fhPhiPhoton(0),               fhEtaPhoton(0), 
     fhEtaPhiPhoton(0),            fhEtaPhi05Photon(0),
@@ -941,13 +941,19 @@ TList *  AliAnaPhoton::GetCreateOutputObjects()
   Int_t nptbins  = GetHistoPtBins();  Float_t ptmax  = GetHistoPtMax();  Float_t ptmin  = GetHistoPtMin(); 
   Int_t nphibins = GetHistoPhiBins(); Float_t phimax = GetHistoPhiMax(); Float_t phimin = GetHistoPhiMin(); 
   Int_t netabins = GetHistoEtaBins(); Float_t etamax = GetHistoEtaMax(); Float_t etamin = GetHistoEtaMin();    
-  Int_t ssbins   = GetHistoShowerShapeBins();  Float_t ssmax = GetHistoShowerShapeMax();  Float_t ssmin = GetHistoShowerShapeMin();
-  Int_t nbins    = GetHistoNClusterCellBins(); Int_t nmax    = GetHistoNClusterCellMax(); Int_t nmin    = GetHistoNClusterCellMin(); 
+  Int_t ssbins   = GetHistoShowerShapeBins();  Float_t ssmax   = GetHistoShowerShapeMax();  Float_t ssmin   = GetHistoShowerShapeMin();
+  Int_t nbins    = GetHistoNClusterCellBins(); Int_t   nmax    = GetHistoNClusterCellMax(); Int_t   nmin    = GetHistoNClusterCellMin(); 
+  Int_t ntimebins= GetHistoTimeBins();         Float_t timemax = GetHistoTimeMax();         Float_t timemin = GetHistoTimeMin();       
 
   fhNCellsE  = new TH2F ("hNCellsE","# of cells in cluster vs E of clusters", nptbins,ptmin,ptmax, nbins,nmin,nmax); 
   fhNCellsE->SetXTitle("E (GeV)");
   fhNCellsE->SetYTitle("# of cells in cluster");
-  outputContainer->Add(fhNCellsE);  
+  outputContainer->Add(fhNCellsE);    
+  
+  fhTimeE  = new TH2F ("hTimeE","time of cluster vs E of clusters", nptbins,ptmin,ptmax, ntimebins,timemin,timemax); 
+  fhTimeE->SetXTitle("E (GeV)");
+  fhTimeE->SetYTitle("time (ns)");
+  outputContainer->Add(fhTimeE);  
   
   fhMaxCellDiffClusterE  = new TH2F ("hMaxCellDiffClusterE","energy vs difference of cluster energy - max cell energy / cluster energy, good clusters",
                                     nptbins,ptmin,ptmax, 500,0,1.); 
@@ -1708,8 +1714,10 @@ void  AliAnaPhoton::MakeAnalysisFillHistograms()
     AliVCluster *cluster = FindCluster(clusters,ph->GetCaloLabel(0),iclus); 
     absID = GetCaloUtils()->GetMaxEnergyCell(cells, cluster,maxCellFraction);
     
+    // Control histograms
     fhMaxCellDiffClusterE->Fill(ph->E(),maxCellFraction);
     fhNCellsE            ->Fill(ph->E(),cluster->GetNCells());
+    fhTimeE              ->Fill(ph->E(),cluster->GetTOF()*1.e9);
     
     //.......................................
     //Play with the MC data if available
index 369b7f6..df288a1 100755 (executable)
@@ -124,7 +124,8 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   //Histograms 
   TH2F * fhNCellsE;                      //! number of cells in cluster vs E 
   TH2F * fhMaxCellDiffClusterE;          //! Fraction of energy carried by cell with maximum energy
-  
+  TH2F * fhTimeE;                        //! time of cluster vs E 
+
   TH1F * fhEPhoton    ;                  //! Number of identified photon vs energy
   TH1F * fhPtPhoton   ;                  //! Number of identified photon vs transerse momentum 
   TH2F * fhPhiPhoton  ;                  //! Azimuthal angle of identified  photon vs transerse momentum 
@@ -214,7 +215,7 @@ class AliAnaPhoton : public AliAnaPartCorrBaseClass {
   TH2F * fhEmbedPi0ELambda0MostlyBkg ;          //!  Lambda0 vs E for embedded photons with 50%<fraction<10% 
   TH2F * fhEmbedPi0ELambda0FullBkg ;            //!  Lambda0 vs E for embedded photons with less than 10% of the cluster energy
   
-   ClassDef(AliAnaPhoton,18)
+   ClassDef(AliAnaPhoton,19)
 
 } ;