message: Adding the MC-kine option for sphericity studies
diff: @@ -593,6 +593,41 @@ Bool_t AliFlowEventCuts::PassesCuts(AliVEvent *event, AliMCEvent *mcevent)
  if(multTPC < (-18.5+1.15*multGlobal) || multTPC > (200.+1.45*multGlobal)) pass = kFALSE;
  }
  }
  + }
  +
  + //==================================================//
  + //MC kine only
  + if(mcevent) {
  + for (Int_t iParticle = 0; iParticle < mcevent->GetNumberOfTracks(); iParticle++) {
  + AliMCParticle* gParticle = dynamic_cast<AliMCParticle *>(mcevent->GetTrack(iParticle));
  + if (!gParticle) {
  + Printf(Form("Could not receive particle %d", iParticle));
  + continue;
  + }
  +
  + //exclude non stable particles
  + if(!(mcevent->IsPhysicalPrimary(iParticle))) continue;
  +
  + //kinematic cuts cuts
  + Float_t pX = gParticle->Px();
  + Float_t pY = gParticle->Py();
  + Float_t pT = gParticle->Pt();
  + Float_t eta = gParticle->Eta();
  +
  + // Kinematics cuts from ESD track cuts
  + if( pT < 0.2 || pT > 5.) continue;
  + if( eta < -0.8 || eta > 0.8) continue;
  +
  + sumPt += pT;
  + if(pT != 0.) {
  + s00 += TMath::Power(pX,2)/pT;
  + s11 += TMath::Power(pY,2)/pT;
  + s10 += pX*pY/pT;
  + }
  +
  + nAcceptedTracks += 1;
  + } //particle loop
  + }//MC event
  //==================================================//
  //sphericity calculation:
  @@ -613,7 +648,8 @@ Bool_t AliFlowEventCuts::PassesCuts(AliVEvent *event, AliMCEvent *mcevent)
  }
  }
  }
  - }
  +
  +
  if (fQA) {
  QAbefore(0)->Fill(pvtxz);
  QAbefore(1)->Fill(multGlobal,multTPC);
  