M AliTPCAnalysisTaskcalib.cxx - Chckeing the friend status bit
authormarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 28 Feb 2011 17:55:17 +0000 (17:55 +0000)
committermarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 28 Feb 2011 17:55:17 +0000 (17:55 +0000)
M      AliTPCcalibDB.cxx           - Removing warning
M      AliTPCcalibGainMult.cxx     - Use the number of tracks as multiplicity estiamtor
M      AliTPCcalibTime.cxx        -  Setter for the maximal track multiplicity per event
                                  + Waring in the TVectorD  - resize needed
M      AliTPCtrackerMI.cxx        -  Use the number of tracks as multiplicity estiamtor, (TPC primary vertex not defined yet in the RefitInward)

TPC/AliTPCAnalysisTaskcalib.cxx
TPC/AliTPCcalibDB.cxx
TPC/AliTPCcalibGainMult.cxx
TPC/AliTPCcalibTime.cxx
TPC/AliTPCtrackerMI.cxx

index 7cc3965..cd45550 100644 (file)
@@ -88,6 +88,7 @@ void AliTPCAnalysisTaskcalib::Exec(Option_t *) {
     //Printf("ERROR: fESDfriend not available");
     return;
   }
+  if (fESDfriend->TestSkipBit()) return;
   Int_t n=fESD->GetNumberOfTracks();
   Process(fESD);
   Int_t run = fESD->GetRunNumber();
index 53f7d4d..8f2c053 100644 (file)
@@ -2015,9 +2015,10 @@ AliTPCCorrection * AliTPCcalibDB::GetTPCComposedCorrectionDelta() const{
   //
   if (!fComposedCorrectionArray) return 0;
   if (fRun<0) return 0;
+  if (fDriftCorrectionArray.GetEntriesFast()<=fRun) return 0;
   if (fDriftCorrectionArray.At(fRun)==0) return 0;
-  if (fComposedCorrectionArray->GetEntriesFast()<4) {
-    fComposedCorrectionArray->Expand(4);
+  if (fComposedCorrectionArray->GetEntriesFast()<=4) {
+    fComposedCorrectionArray->Expand(5);
     TObjArray * timeArray =(TObjArray*)(fDriftCorrectionArray.At(fRun));
      AliTPCCorrection * correctionTime = (AliTPCCorrection *)timeArray->FindObject("FitCorrectionTime");
      if (correctionTime){
index c8af677..f5dc9ac 100644 (file)
@@ -128,7 +128,7 @@ AliTPCcalibGainMult::AliTPCcalibGainMult(const Text_t *name, const Text_t *title
   //
   Int_t binsPadEqual[6]    = { 200, 200,    4,   20,   50, 100};
   Double_t xminPadEqual[6] = { 0.5, 0.5, -0.5,    0, -250,   0}; 
-  Double_t xmaxPadEqual[6] = { 1.5, 1.5,  3.5, 4000,  250,   3};
+  Double_t xmaxPadEqual[6] = { 1.5, 1.5,  3.5, 13000,  250,   3};
   TString axisNamePadEqual[6]   = {"dEdxRatioMax","dEdxRatioTot","padType","mult","driftlength", "1_pt"};
   TString axisTitlePadEqual[6]  = {"dEdx_padRegion/mean_dEdx Qmax", "dEdx_padRegion/mean_dEdx Qtot","padType","mult","driftlength", "1/pt"};
   //
@@ -142,7 +142,7 @@ AliTPCcalibGainMult::AliTPCcalibGainMult(const Text_t *name, const Text_t *title
   //                    MIP Qmax, MIP Qtot,  z,  pad, vtx. contribut., ncl
   Int_t binsGainMult[6]    = { 145,  145,   25,    4,  100,  80};
   Double_t xminGainMult[6] = { 10.,  10.,    0, -0.5,    0, -0.5}; 
-  Double_t xmaxGainMult[6] = {300., 300.,  250,  3.5, 5000, 159.5};
+  Double_t xmaxGainMult[6] = {300., 300.,  250,  3.5, 13000, 159.5};
   TString axisNameMult[6]={"Qmax","Qtot","drift","padtype""multiplicity","ncl"};
   TString axisTitleMult[6]={"Qmax (a.u)","Qtot (a.u.)","driftlenght l (cm)","Pad Type","multiplicity","ncl"};
   //
@@ -191,8 +191,7 @@ void AliTPCcalibGainMult::Process(AliESDEvent *event) {
    return;
   }
   UInt_t runNumber = event->GetRunNumber();
-  Int_t nContributors = 0;
-  if (event->GetPrimaryVertexTPC()) nContributors = event->GetPrimaryVertexTPC()->GetNContributors();
+  Int_t nContributors = event->GetNumberOfTracks();
   //
   // track loop
   //
index 3b44ecd..e7c6a04 100644 (file)
@@ -477,7 +477,7 @@ void AliTPCcalibTime::ProcessLaser(AliESDEvent *event){
   //
   // fill histos
   //
-  TVectorD vdriftA(5), vdriftC(5),vdriftAC(5);
+  TVectorD vdriftA(5), vdriftC(5),vdriftAC(6);
   vdriftA=*(fLaser->fFitAside);
   vdriftC=*(fLaser->fFitCside);
   vdriftAC=*(fLaser->fFitACside);
@@ -805,6 +805,7 @@ void AliTPCcalibTime::ProcessBeam(const AliESDEvent *const event){
   Double_t dcaVertex[2]={0,0};
   Int_t ntracks=event->GetNumberOfTracks();
   if (ntracks==0) return;
+  if (ntracks > fCutTracks) return;
   //
   AliESDfriend *esdFriend=(AliESDfriend*)(((AliESDEvent*)event)->FindListObject("AliESDfriend"));
   //
index 1a30123..5a26145 100644 (file)
@@ -2682,7 +2682,7 @@ Int_t AliTPCtrackerMI::RefitInward(AliESDEvent *event)
   }
   transform->SetCurrentRecoParam((AliTPCRecoParam*)AliTPCReconstructor::GetRecoParam());
   const AliTPCRecoParam * recoParam = AliTPCcalibDB::Instance()->GetTransform()->GetCurrentRecoParam();
-
+  Int_t nContribut = event->GetNumberOfTracks();
   TGraphErrors * graphMultDependenceDeDx = 0x0;
   if (recoParam && recoParam->GetUseMultiplicityCorrectionDedx() && gainCalibArray) {
     if (recoParam->GetUseTotCharge()) {
@@ -2760,9 +2760,8 @@ Int_t AliTPCtrackerMI::RefitInward(AliESDEvent *event)
       Float_t dedx  = seed->GetdEdx();
       // apply mutliplicity dependent dEdx correction if available
       if (graphMultDependenceDeDx) {
-       Int_t nContribut  = event->GetPrimaryVertexTPC()->GetNContributors();
        Double_t corrGain =  AliTPCcalibDButil::EvalGraphConst(graphMultDependenceDeDx, nContribut);
-       dedx /= corrGain;
+       dedx += (1 - corrGain)*50.; // MIP is normalized to 50
       }
       esd->SetTPCsignal(dedx, sdedx, ndedx);
       //