Small corrections to FemtoEventReaders (centrality selection) + uncommented EvMult...
authormajanik <majanik@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Mar 2013 17:16:38 +0000 (17:16 +0000)
committermajanik <majanik@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 25 Mar 2013 17:16:38 +0000 (17:16 +0000)
PWGCF/FEMTOSCOPY/AliFemto/AliFemtoCutMonitorEventMult.cxx
PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx
PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderESDChain.cxx

index 2aa747e..b6cec16 100644 (file)
@@ -308,7 +308,7 @@ void AliFemtoCutMonitorEventMult::Fill(const AliFemtoEvent* aEvent)
 void AliFemtoCutMonitorEventMult::Write()
 {
   // Write out the relevant histograms
-  // fEvMult->Write();
+  fEvMult->Write();
   fNormEvMult->Write();
   fPsiVZERO->Write();
 
index 0c4b8dd..a017878 100644 (file)
@@ -647,46 +647,25 @@ void AliFemtoEventReaderAOD::CopyAODtoFemtoEvent(AliFemtoEvent *tEvent)
   tEvent->SetNumberOfTracks(realnofTracks);//setting number of track which we read in event    
   tEvent->SetNormalizedMult(tracksPrim);
 
+  if (cent) {
+    tEvent->SetCentralityV0(cent->GetCentralityPercentile("V0M"));
+    tEvent->SetCentralityZNA(cent->GetCentralityPercentile("ZNA"));
+    tEvent->SetCentralityCL1(cent->GetCentralityPercentile("CL1"));
+    //    tEvent->SetCentralityFMD(cent->GetCentralityPercentile("FMD"));
+    //    tEvent->SetCentralityTrk(cent->GetCentralityPercentile("TRK"));
+  }
+
   if (fEstEventMult==kCentrality) {
-    //  AliCentrality *cent = fEvent->GetCentrality();
+    //AliCentrality *cent = fEvent->GetCentrality();
     //cout<<"AliFemtoEventReaderAOD:"<<lrint(10*cent->GetCentralityPercentile("V0M"))<<endl;
     if (cent) tEvent->SetNormalizedMult(lrint(10*cent->GetCentralityPercentile("V0M")));
     //  if (cent) tEvent->SetNormalizedMult((int) cent->GetCentralityPercentile("V0M"));
-    
-    if (cent) {
-      tEvent->SetCentralityV0(cent->GetCentralityPercentile("V0M"));
-      tEvent->SetCentralityZNA(cent->GetCentralityPercentile("ZNA"));
-      tEvent->SetCentralityCL1(cent->GetCentralityPercentile("CL1"));
-      //    tEvent->SetCentralityFMD(cent->GetCentralityPercentile("FMD"));
-      tEvent->SetCentralitySPD1(cent->GetCentralityPercentile("CL1"));
-      //    tEvent->SetCentralityTrk(cent->GetCentralityPercentile("TRK"));
-    }
   }
   else if (fEstEventMult==kCentralityZNA) {
-    //  AliCentrality *cent = fEvent->GetCentrality();
-    //cout<<"AliFemtoEventReaderAOD:"<<lrint(10*cent->GetCentralityPercentile("V0M"))<<endl;
     if (cent) tEvent->SetNormalizedMult(lrint(10*cent->GetCentralityPercentile("ZNA")));
-    //  if (cent) tEvent->SetNormalizedMult((int) cent->GetCentralityPercentile("V0M"));
-    
-    if (cent) {
-      tEvent->SetCentralityV0(cent->GetCentralityPercentile("ZNA"));
-      //    tEvent->SetCentralityFMD(cent->GetCentralityPercentile("FMD"));
-      tEvent->SetCentralitySPD1(cent->GetCentralityPercentile("CL1"));
-      //    tEvent->SetCentralityTrk(cent->GetCentralityPercentile("TRK"));
-    }
   }
   else if (fEstEventMult==kCentralityCL1) {
-    //  AliCentrality *cent = fEvent->GetCentrality();
-    //cout<<"AliFemtoEventReaderAOD:"<<lrint(10*cent->GetCentralityPercentile("V0M"))<<endl;
     if (cent) tEvent->SetNormalizedMult(lrint(10*cent->GetCentralityPercentile("CL1")));
-    //  if (cent) tEvent->SetNormalizedMult((int) cent->GetCentralityPercentile("V0M"));
-    
-    if (cent) {
-      tEvent->SetCentralityV0(cent->GetCentralityPercentile("CL1"));
-      //    tEvent->SetCentralityFMD(cent->GetCentralityPercentile("FMD"));
-      tEvent->SetCentralitySPD1(cent->GetCentralityPercentile("CL1"));
-      //    tEvent->SetCentralityTrk(cent->GetCentralityPercentile("TRK"));
-    }
   }
   else if(fEstEventMult==kGlobalCount){
     tEvent->SetNormalizedMult(tNormMult); //particles counted in the loop, trying to reproduce GetReferenceMultiplicity. If better (default) method appears it should be changed
index 495ed24..1a89362 100644 (file)
@@ -690,7 +690,7 @@ AliFemtoEvent* AliFemtoEventReaderESDChain::ReturnHbtEvent()
     hbtEvent->SetCentralityZNA(cent->GetCentralityPercentile("ZNA"));
     hbtEvent->SetCentralityCL1(cent->GetCentralityPercentile("CL1"));
     //    hbtEvent->SetCentralityFMD(cent->GetCentralityPercentile("FMD"));
-    hbtEvent->SetCentralitySPD1(cent->GetCentralityPercentile("CL1"));
+    //    hbtEvent->SetCentralitySPD1(cent->GetCentralityPercentile("CL1"));
     //    hbtEvent->SetCentralityTrk(cent->GetCentralityPercentile("TRK"));
 
     if (Debug()>1) printf("  FemtoReader Got Event with %f %f %f %f\n", cent->GetCentralityPercentile("V0M"), 0.0, cent->GetCentralityPercentile("CL1"), 0.0);