take tpc refit in again (M.Verweij)
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 20 Sep 2012 11:24:48 +0000 (11:24 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 20 Sep 2012 11:24:48 +0000 (11:24 +0000)
PWGJE/AliPWG4HighPtSpectra.cxx
PWGJE/AliPWG4HighPtTrackQA.cxx
PWGJE/macros/AddTaskPWG4HighPtSpectra.C
PWGJE/macros/AddTaskPWG4HighPtTrackQA.C
PWGJE/macros/CreateTrackCutsPWGJE.C

index 787d009..5d44059 100644 (file)
@@ -351,130 +351,179 @@ void AliPWG4HighPtSpectra::Exec(Option_t *)
   //Now go to rec level
   for (Int_t iTrack = 0; iTrack<nTracks; iTrack++) 
     {   
+   
+
       //Get track for analysis
       AliESDtrack *track = 0x0;
       AliESDtrack *esdtrack = fESD->GetTrack(iTrack);
-      if(!esdtrack) continue;
+      if(!esdtrack)
+       continue;
+      
 
-      if(fTrackType==1) {
-       track = AliESDtrackCuts::GetTPCOnlyTrack(fESD,esdtrack->GetID());
-       if(!track) continue;
+      if(fTrackType==4) {
+       if (!(fTrackCuts->AcceptTrack(esdtrack)))
+         continue;
       }
-      else if(fTrackType==2) {
-       track = AliESDtrackCuts::GetTPCOnlyTrack(fESD,esdtrack->GetID());
-       if(!track) continue;
 
+      if(fTrackType==1)
+       track = AliESDtrackCuts::GetTPCOnlyTrack(fESD,esdtrack->GetID());
+      else if(fTrackType==2 || fTrackType==4) {
+       track = AliESDtrackCuts::GetTPCOnlyTrack(const_cast<AliESDEvent*>(fESD),esdtrack->GetID());
+       if(!track)
+         continue;
+       
        AliExternalTrackParam exParam;
        Bool_t relate = track->RelateToVertexTPC(fVtx,fESD->GetMagneticField(),kVeryBig,&exParam);
        if( !relate ) {
-         delete track;
+         if(track) delete track;
          continue;
        }
        track->Set(exParam.GetX(),exParam.GetAlpha(),exParam.GetParameter(),exParam.GetCovariance());
       }
+      else if(fTrackType==5 || fTrackType==6) {
+       if(fTrackCuts->AcceptTrack(esdtrack)) {
+         continue;
+       }
+       else {
+         if( !(fTrackCutsReject->AcceptTrack(esdtrack)) && fTrackCuts->AcceptTrack(esdtrack) ) {
+
+           if(fTrackType==5) {
+             //use TPConly constrained track
+             track = AliESDtrackCuts::GetTPCOnlyTrack(fESD,esdtrack->GetID());
+             if(!track)
+               continue;
+             
+             AliExternalTrackParam exParam;
+             Bool_t relate = track->RelateToVertexTPC(fVtx,fESD->GetMagneticField(),kVeryBig,&exParam);
+             if( !relate ) {
+               if(track) delete track;
+               continue;
+             }
+             track->Set(exParam.GetX(),exParam.GetAlpha(),exParam.GetParameter(),exParam.GetCovariance());
+           }
+           else if(fTrackType==6) {
+             //use global constrained track
+             track = new AliESDtrack(*esdtrack);
+             track->Set(esdtrack->GetConstrainedParam()->GetX(),esdtrack->GetConstrainedParam()->GetAlpha(),esdtrack->GetConstrainedParam()->GetParameter(),esdtrack->GetConstrainedParam()->GetCovariance());
+
+           }
+         }
+       }
+      }
       else if(fTrackType==7) {
        //use global constrained track
        track = new AliESDtrack(*esdtrack);
-       //      track = esdtrack;
-       //      track->Set(esdtrack->GetConstrainedParam()->GetX(),esdtrack->GetConstrainedParam()->GetAlpha(),esdtrack->GetConstrainedParam()->GetParameter(),esdtrack->GetConstrainedParam()->GetCovariance());
-
       }
       else
        track = esdtrack;
     
-      if(!track) continue;
-      if(fTrackType==2) {
+      if(!track)
+       continue;
+      
+
+      if(fTrackType==2 || fTrackType==4 || fTrackType==5) {
        //Cut on chi2 of constrained fit
-       if(track->GetConstrainedChi2TPC() > fSigmaConstrainedMax*fSigmaConstrainedMax) {
-         delete track;
+       if(track->GetConstrainedChi2TPC() > fSigmaConstrainedMax*fSigmaConstrainedMax && fSigmaConstrainedMax>0.) {
+         if(track) delete track;
          continue;
        }
       }
+      if (!(fTrackCuts->AcceptTrack(track)) && fTrackType!=4 && fTrackType!=5 && fTrackType!=6) {
+       if(fTrackType==1 || fTrackType==2 || fTrackType==7) {
+         if(track) delete track;
+       }
+       continue;
+      }
 
-      if (fTrackCuts->AcceptTrack(track)) {
-
-       if(fTrackType==7) {
-         if(fTrackCutsReject ) {
-           if(fTrackCutsReject->AcceptTrack(track) ) {
-             if(track) delete track;
-             continue;
-           }
+      if(fTrackType==7) {
+       if(fTrackCutsReject ) {
+         if(fTrackCutsReject->AcceptTrack(track) ) {
+           if(track) delete track;
+           continue;
          }
-         
-         if(esdtrack->GetConstrainedParam()) 
-           track->Set(esdtrack->GetConstrainedParam()->GetX(),esdtrack->GetConstrainedParam()->GetAlpha(),esdtrack->GetConstrainedParam()->GetParameter(),esdtrack->GetConstrainedParam()->GetCovariance());
        }
+      
+       if(esdtrack->GetConstrainedParam()) 
+         track->Set(esdtrack->GetConstrainedParam()->GetX(),esdtrack->GetConstrainedParam()->GetAlpha(),esdtrack->GetConstrainedParam()->GetParameter(),esdtrack->GetConstrainedParam()->GetCovariance());
+      }
 
-       //fill the container
-       containerInputRec[0] = track->Pt();
-       containerInputRec[1] = track->Phi();
-       containerInputRec[2] = track->Eta();
-       containerInputRec[3] = track->GetTPCNclsIter1();
+      if(!track) {
+       if(fTrackType==1 || fTrackType==2 || fTrackType==4 || fTrackType==5 || fTrackType==6 || fTrackType==7) {
+         if(track) delete track;
+       }
+       continue;
+      }
+
+      //fill the container
+      containerInputRec[0] = track->Pt();
+      containerInputRec[1] = track->Phi();
+      containerInputRec[2] = track->Eta();
+      containerInputRec[3] = track->GetTPCNclsIter1();
 
-       if(track->GetSign()>0.) fCFManagerPos->GetParticleContainer()->Fill(containerInputRec,kStepReconstructed);
-       if(track->GetSign()<0.) fCFManagerNeg->GetParticleContainer()->Fill(containerInputRec,kStepReconstructed);
+      if(track->GetSign()>0.) fCFManagerPos->GetParticleContainer()->Fill(containerInputRec,kStepReconstructed);
+      if(track->GetSign()<0.) fCFManagerNeg->GetParticleContainer()->Fill(containerInputRec,kStepReconstructed);
 
        
-       //Only fill the MC containers if MC information is available
-       if(fMC) {
-         Int_t label = TMath::Abs(track->GetLabel());
-         if(label>fStack->GetNtrack()) {
+      //Only fill the MC containers if MC information is available
+      if(fMC) {
+       Int_t label = TMath::Abs(track->GetLabel());
+       if(label>fStack->GetNtrack()) {
+         if(fTrackType==1 || fTrackType==2 || fTrackType==7)
+           delete track;
+         continue;
+       }
+       //Only select particles generated by HIJING if requested
+       if(fbSelectHIJING) {
+         if(!IsHIJINGParticle(label)) {
            if(fTrackType==1 || fTrackType==2 || fTrackType==7)
              delete track;
            continue;
          }
-         //Only select particles generated by HIJING if requested
-         if(fbSelectHIJING) {
-           if(!IsHIJINGParticle(label)) {
-             if(fTrackType==1 || fTrackType==2 || fTrackType==7)
-               delete track;
-             continue;
-           }
+       }
+       TParticle *particle = fStack->Particle(label) ;
+       if(!particle) {
+         if(fTrackType==1 || fTrackType==2 || fTrackType==7)
+           delete track;
+         continue;
+       }
+       containerInputRecMC[0] = particle->Pt();      
+       containerInputRecMC[1] = particle->Phi();      
+       containerInputRecMC[2] = particle->Eta();  
+       containerInputRecMC[3] = track->GetTPCNclsIter1();
+
+       //Container with primaries
+       if(fStack->IsPhysicalPrimary(label)) {
+         if(particle->GetPDG()->Charge()>0.) {
+           fCFManagerPos->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
          }
-         TParticle *particle = fStack->Particle(label) ;
-         if(!particle) {
-           if(fTrackType==1 || fTrackType==2 || fTrackType==7)
-             delete track;
-           continue;
+         if(particle->GetPDG()->Charge()<0.) {
+           fCFManagerNeg->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
          }
-         containerInputRecMC[0] = particle->Pt();      
-         containerInputRecMC[1] = particle->Phi();      
-         containerInputRecMC[2] = particle->Eta();  
-         containerInputRecMC[3] = track->GetTPCNclsIter1();
-
-         //Container with primaries
-         if(fStack->IsPhysicalPrimary(label)) {
-           if(particle->GetPDG()->Charge()>0.) {
-             fCFManagerPos->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
-           }
-           if(particle->GetPDG()->Charge()<0.) {
-             fCFManagerNeg->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
-           }
 
-           //Fill pT resolution plots for primaries
-           fPtRelUncertainty1PtPrim->Fill(containerInputRec[0],containerInputRec[0]*TMath::Sqrt(track->GetSigma1Pt2()));
+         //Fill pT resolution plots for primaries
+         fPtRelUncertainty1PtPrim->Fill(containerInputRec[0],containerInputRec[0]*TMath::Sqrt(track->GetSigma1Pt2()));
 
-         }
+       }
 
-         //Container with secondaries
-         if (!fStack->IsPhysicalPrimary(label) ) {
-           if(particle->GetPDG()->Charge()>0.) {
-             fCFManagerPos->GetParticleContainer()->Fill(containerInputMC,kStepSecondaries);
-           }
-           if(particle->GetPDG()->Charge()<0.) {
-             fCFManagerNeg->GetParticleContainer()->Fill(containerInputMC,kStepSecondaries);
-           }
-           //Fill pT resolution plots for primaries
-           fPtRelUncertainty1PtSec->Fill(containerInputRec[0],containerInputRec[0]*TMath::Sqrt(track->GetSigma1Pt2()));
+       //Container with secondaries
+       if (!fStack->IsPhysicalPrimary(label) ) {
+         if(particle->GetPDG()->Charge()>0.) {
+           fCFManagerPos->GetParticleContainer()->Fill(containerInputMC,kStepSecondaries);
+         }
+         if(particle->GetPDG()->Charge()<0.) {
+           fCFManagerNeg->GetParticleContainer()->Fill(containerInputMC,kStepSecondaries);
          }
+         //Fill pT resolution plots for primaries
+         fPtRelUncertainty1PtSec->Fill(containerInputRec[0],containerInputRec[0]*TMath::Sqrt(track->GetSigma1Pt2()));
        }
+      }
        
-      }//trackCuts global tracks
 
-      if(fTrackType==1 || fTrackType==2 || fTrackType==7) {
+      if(fTrackType==1  || fTrackType==2 || fTrackType==4 || fTrackType==5 || fTrackType==6 || fTrackType==7) {
        if(track) delete track;
       }
+
+
     }//track loop
   
 
@@ -544,7 +593,7 @@ Bool_t AliPWG4HighPtSpectra::PythiaInfoFromFile(const char* currFile,Float_t &fX
     // next trial fetch the histgram file
     fxsec = TFile::Open(Form("%s%s",file.Data(),"pyxsec_hists.root"));
     if(!fxsec){
-       // not a severe condition but inciate that we have no information
+      // not a severe condition but indicates that we have no information
       return kFALSE;
     }
     else{
index 7127447..5902fae 100644 (file)
@@ -1017,7 +1017,7 @@ void AliPWG4HighPtTrackQA::DoAnalysisESD() {
     if(fTrackType==1)
       track = AliESDtrackCuts::GetTPCOnlyTrack(fESD,esdtrack->GetID());
     else if(fTrackType==2 || fTrackType==4) {
-      track = AliESDtrackCuts::GetTPCOnlyTrack(fESD,esdtrack->GetID());
+      track = AliESDtrackCuts::GetTPCOnlyTrack(const_cast<AliESDEvent*>(fESD),esdtrack->GetID());
       if(!track) {
        fh1NTracksReject->Fill("noTPConly",1);
        if(origtrack) delete origtrack;
@@ -1095,7 +1095,7 @@ void AliPWG4HighPtTrackQA::DoAnalysisESD() {
     if (!(fTrackCuts->AcceptTrack(track)) && fTrackType!=4 && fTrackType!=5 && fTrackType!=6) {
       fh1NTracksReject->Fill("trackCuts",1);
       if(fTrackType==1 || fTrackType==2 || fTrackType==7) {
-       if(track) delete track;
+       if(track) delete track;
       }
       if(origtrack) delete origtrack;
       continue;
index 81f950c..c2a7c48 100644 (file)
@@ -59,7 +59,7 @@ void AddTaskPWG4HighPtSpectra2011(char *prodType = "LHC10h", Bool_t isPbPb=kTRUE
   AliPWG4HighPtSpectra *taskSpectra05cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlag,bSelHijingParticles);
   AliPWG4HighPtSpectra *taskSpectra74cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlag,bSelHijingParticles);
   AliPWG4HighPtSpectra *taskSpectra75cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlag,bSelHijingParticles);
-  AliPWG4HighPtSpectra *taskSpectra20cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlag,bSelHijingParticles);
+  AliPWG4HighPtSpectra *taskSpectra40cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,4,0,iPhysicsSelectionFlag,bSelHijingParticles);
 
   AliPWG4HighPtSpectra *taskSpectraEMCJE00cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
   AliPWG4HighPtSpectra *taskSpectraEMCJE01cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
@@ -69,7 +69,7 @@ void AddTaskPWG4HighPtSpectra2011(char *prodType = "LHC10h", Bool_t isPbPb=kTRUE
   AliPWG4HighPtSpectra *taskSpectraEMCJE05cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
   AliPWG4HighPtSpectra *taskSpectraEMCJE74cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
   AliPWG4HighPtSpectra *taskSpectraEMCJE75cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
-  AliPWG4HighPtSpectra *taskSpectraEMCJE20cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
+  AliPWG4HighPtSpectra *taskSpectraEMCJE40cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,4,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
 
 
   if(isPbPb) {
@@ -83,7 +83,7 @@ void AddTaskPWG4HighPtSpectra2011(char *prodType = "LHC10h", Bool_t isPbPb=kTRUE
       AliPWG4HighPtSpectra *taskSpectra05 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlag,bSelHijingParticles);
       AliPWG4HighPtSpectra *taskSpectra74 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlag,bSelHijingParticles);
       AliPWG4HighPtSpectra *taskSpectra75 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlag,bSelHijingParticles);
-      AliPWG4HighPtSpectra *taskSpectra20 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlag,bSelHijingParticles);
+      AliPWG4HighPtSpectra *taskSpectra40 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,4,0,iPhysicsSelectionFlag,bSelHijingParticles);
 
       AliPWG4HighPtSpectra *taskSpectraEMCJE00 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
       AliPWG4HighPtSpectra *taskSpectraEMCJE01 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
@@ -93,7 +93,7 @@ void AddTaskPWG4HighPtSpectra2011(char *prodType = "LHC10h", Bool_t isPbPb=kTRUE
       AliPWG4HighPtSpectra *taskSpectraEMCJE05 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
       AliPWG4HighPtSpectra *taskSpectraEMCJE74 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
       AliPWG4HighPtSpectra *taskSpectraEMCJE75 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
-      AliPWG4HighPtSpectra *taskSpectraEMCJE20 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
+      AliPWG4HighPtSpectra *taskSpectraEMCJE40 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,4,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
 
     }
 
index 5d1ed30..f011d11 100644 (file)
@@ -143,7 +143,6 @@ void AddTaskPWG4HighPtTrackQALHC11hLTS(char *prodType = "LHC10h",Bool_t isPbPb=k
   UInt_t iPhysicsSelectionFlag = AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral;
   UInt_t iPhysicsSelectionFlagEMCEJE = AliVEvent::kEMCEJE;
   
-  cout << "Add cent10 AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral" << endl;
   AliPWG4HighPtTrackQA *taskTrackQA00cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlag);
   AliPWG4HighPtTrackQA *taskTrackQA01cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlag);
   AliPWG4HighPtTrackQA *taskTrackQA70cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlag);
@@ -152,9 +151,8 @@ void AddTaskPWG4HighPtTrackQALHC11hLTS(char *prodType = "LHC10h",Bool_t isPbPb=k
   AliPWG4HighPtTrackQA *taskTrackQA05cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlag);
   AliPWG4HighPtTrackQA *taskTrackQA74cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlag);
   AliPWG4HighPtTrackQA *taskTrackQA75cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlag);
-  AliPWG4HighPtTrackQA *taskTrackQA20cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,0,iPhysicsSelectionFlag);
+  AliPWG4HighPtTrackQA *taskTrackQA40cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlag);
 
-  cout << "Add cent10 AliVEvent::kEMCEJE" << endl;
   AliPWG4HighPtTrackQA *taskTrackQAEMCJE00cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagEMCEJE);
   AliPWG4HighPtTrackQA *taskTrackQAEMCJE01cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagEMCEJE);
   AliPWG4HighPtTrackQA *taskTrackQAEMCJE70cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlagEMCEJE);
@@ -163,11 +161,11 @@ void AddTaskPWG4HighPtTrackQALHC11hLTS(char *prodType = "LHC10h",Bool_t isPbPb=k
   AliPWG4HighPtTrackQA *taskTrackQAEMCJE05cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagEMCEJE);
   AliPWG4HighPtTrackQA *taskTrackQAEMCJE74cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagEMCEJE);
   AliPWG4HighPtTrackQA *taskTrackQAEMCJE75cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagEMCEJE);
-  AliPWG4HighPtTrackQA *taskTrackQAEMCJE20cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,0,iPhysicsSelectionFlagEMCEJE);
+  AliPWG4HighPtTrackQA *taskTrackQAEMCJE40cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlagEMCEJE);
 
   if(isPbPb) {
     for(cent=0; cent<4; cent++) {
-      cout << "Add Cent " << cent << endl;
+
       AliPWG4HighPtTrackQA *taskTrackQA00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlag);
       AliPWG4HighPtTrackQA *taskTrackQA01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlag);
       AliPWG4HighPtTrackQA *taskTrackQA70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlag);
@@ -176,7 +174,7 @@ void AddTaskPWG4HighPtTrackQALHC11hLTS(char *prodType = "LHC10h",Bool_t isPbPb=k
       AliPWG4HighPtTrackQA *taskTrackQA05 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlag);
       AliPWG4HighPtTrackQA *taskTrackQA74 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlag);
       AliPWG4HighPtTrackQA *taskTrackQA75 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlag);
-      AliPWG4HighPtTrackQA *taskTrackQA20 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,0,iPhysicsSelectionFlag);
+      AliPWG4HighPtTrackQA *taskTrackQA40 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlag);
 
       AliPWG4HighPtTrackQA *taskTrackQAEMCJE00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagEMCEJE);
       AliPWG4HighPtTrackQA *taskTrackQAEMCJE01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagEMCEJE);
@@ -186,7 +184,8 @@ void AddTaskPWG4HighPtTrackQALHC11hLTS(char *prodType = "LHC10h",Bool_t isPbPb=k
       AliPWG4HighPtTrackQA *taskTrackQAEMCJE05 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagEMCEJE);
       AliPWG4HighPtTrackQA *taskTrackQAEMCJE74 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagEMCEJE);
       AliPWG4HighPtTrackQA *taskTrackQAEMCJE75 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagEMCEJE);
-      AliPWG4HighPtTrackQA *taskTrackQAEMCJE20 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,0,iPhysicsSelectionFlagEMCEJE);
+      AliPWG4HighPtTrackQA *taskTrackQAEMCJE40 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlagEMCEJE);
+
     }
   }
 
index ccac91b..9657d1e 100644 (file)
@@ -282,7 +282,7 @@ AliESDtrackCuts *CreateTrackCutsPWGJE(Int_t cutMode) {
     bStdCutsDefined = kTRUE;
 
     trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts(); 
-    trackCuts->SetRequireTPCRefit(kTRUE);
+    // trackCuts->SetRequireTPCRefit(kTRUE);
     trackCuts->SetMinNClustersTPC(70);
 
     trackCuts->SetEtaRange(-0.9,0.9);
@@ -298,7 +298,7 @@ AliESDtrackCuts *CreateTrackCutsPWGJE(Int_t cutMode) {
     bStdCutsDefined = kTRUE;
 
     trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts(); 
-    trackCuts->SetRequireTPCRefit(kTRUE);
+    //  trackCuts->SetRequireTPCRefit(kTRUE);
     trackCuts->SetMinNClustersTPC(0);
     trackCuts->SetMinNCrossedRowsTPC(120);
     trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.1);// essentially switches it off