message: Pair cut for both standard and special cases
diff: @@ -83,9 +83,9 @@ AliFemtoModelCorrFctn::AliFemtoModelCorrFctn(const char *title, Int_t aNbins, Do
  fDenominatorIdeal = new TH1D(buf,buf,aNbins,aQinvLo,aQinvHi);
  snprintf(buf , 100, "QgenQrec%s", title);
  - //fQgenQrec = new TH2D(buf,buf,aNbins,aQinvLo,aQinvHi,aNbins,aQinvLo,aQinvHi);
  + fQgenQrec = new TH2D(buf,buf,aNbins,aQinvLo,aQinvHi,aNbins,aQinvLo,aQinvHi);
  //test
  - fQgenQrec = new TH2D(buf,buf,aNbins,aQinvLo,aQinvHi,aNbins,-0.05,0.05);
  + //fQgenQrec = new TH2D(buf,buf,aNbins,aQinvLo,aQinvHi,aNbins,-0.05,0.05);
  fNumeratorTrue->Sumw2();
  fNumeratorFake->Sumw2();
  @@ -265,9 +265,9 @@ void AliFemtoModelCorrFctn::AddMixedPair(AliFemtoPair* aPair)
  Double_t tQinvTrue = GetQinvTrue(aPair);
  if(tQinvTrue>0)fNumeratorFakeIdeal->Fill(tQinvTrue, weight);
  if(tQinvTrue>0)fDenominatorIdeal->Fill(tQinvTrue, 1.0);
  - //if(tQinvTrue>0)fQgenQrec->Fill(tQinvTrue,aPair->QInv());
  + if(tQinvTrue>0)fQgenQrec->Fill(tQinvTrue,aPair->QInv());
  //test
  - if(tQinvTrue>0)fQgenQrec->Fill(tQinvTrue,tQinvTrue-aPair->QInv());
  + //if(tQinvTrue>0)fQgenQrec->Fill(tQinvTrue,tQinvTrue-aPair->QInv());
  }
  }
  @@ -277,6 +277,11 @@ Double_t AliFemtoModelCorrFctn::GetQinvTrue(AliFemtoPair* aPair)
  if(!fKaonPDG) {
  AliFemtoTrack *inf1 = (AliFemtoTrack *) aPair->Track1()->Track();
  AliFemtoTrack *inf2 = (AliFemtoTrack *) aPair->Track2()->Track();
  + if(!inf1 || !inf2)
  + {
  +// cout<<"Warning! AliFemtoModelCorrFctn::GetQinvTrue - Could not find requested tracks."<<endl;
  + return -1;
  + }
  AliFemtoLorentzVector fm1;
  AliFemtoThreeVector* temp = ((AliFemtoModelHiddenInfo*)inf1->GetHiddenInfo())->GetTrueMomentum();
  