Changeset 168

Timestamp:

Jun 16, 2006, 4:44:20 PM (19 years ago)

Author:

barrand

Message:

G.Barrand : ajout MEMPHYS_analysis

File:

• MEMPHYS/HEAD/applications/MEMPHYS_analysis.cxx (modified) (5 diffs)

MEMPHYS/HEAD/applications/MEMPHYS_analysis.cxx

@@ -102 +102 @@
 }
 //////////////////////////////////////////////////////////////////////////////
+bool process_hits_times(
+ AIDA::ITuple& aParent
+,AIDA::IHistogram1D& aHisto
+) 
+//////////////////////////////////////////////////////////////////////////////
+//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
+{
+  AIDA::ITupleEntry* entry = (AIDA::ITupleEntry*)aParent.getObject(2);
+  if(!entry) return false;
+
+  AIDA::ITuple* tuple = dynamic_cast<AIDA::ITuple*>(entry);
+  if(!tuple) return false;
+
+  tuple->start();
+  while(tuple->next()) {
+
+    float time = tuple->getFloat(0);
+
+    aHisto.fill(time);
+  }
+  return true;
+}
+//////////////////////////////////////////////////////////////////////////////
+bool process_hits(
+ AIDA::ITuple& aParent
+,AIDA::IHistogram1D& aHisto
+) 
+//////////////////////////////////////////////////////////////////////////////
+//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
+{
+  AIDA::ITupleEntry* entry = (AIDA::ITupleEntry*)aParent.getObject(10);
+  if(!entry) return false;
+
+  AIDA::ITuple* tuple = dynamic_cast<AIDA::ITuple*>(entry);
+  if(!tuple) return false;
+
+  tuple->start();
+  while(tuple->next()) {
+
+    //int tubeId = tuple->getInt(0);
+    //int totalPE = tuple->getInt(1);
+
+    if(!process_hits_times(*tuple,aHisto)) return false;
+  }
+  return true;
+}
+//////////////////////////////////////////////////////////////////////////////
+bool process_digits(
+ AIDA::ITuple& aParent
+,AIDA::IHistogram2D& aHisto
+) 
+//////////////////////////////////////////////////////////////////////////////
+//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
+{
+  AIDA::ITupleEntry* entry = (AIDA::ITupleEntry*)aParent.getObject(13);
+  if(!entry) return false;
+
+  AIDA::ITuple* tuple = dynamic_cast<AIDA::ITuple*>(entry);
+  if(!tuple) return false;
+
+  tuple->start();
+  while(tuple->next()) {
+
+    //int tubeId = tuple->getInt(0);
+    double pe = tuple->getDouble(1);
+    double time = tuple->getDouble(2);
+    //printf("debug : ++++ : %g %g\n",time,pe);
+
+    aHisto.fill(time,pe);
+  }
+  return true;
+}
+//////////////////////////////////////////////////////////////////////////////
+bool plot(
+ AIDA::IAnalysisFactory& aAIDA
+,AIDA::IHistogram1D& aHisto1D
+,AIDA::IHistogram2D& aHisto2D
+) 
+//////////////////////////////////////////////////////////////////////////////
+//!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!//
+{
+  std::string header = "MEMPHYS_analysis : ";
+
+  AIDA::IPlotterFactory* plotterFactory = 
+    aAIDA.createPlotterFactory(0,0);
+  if(!plotterFactory) {
+    std::cout << header << "can't get a PlotterFactory." << std::endl;
+    return false;
+  }
+    
+  AIDA::IPlotter* plotter = plotterFactory->create();
+  if(!plotter) {
+    std::cout << header << "can't get a plotter." << std::endl;
+    return false;
+  }
+    
+  plotter->clearRegions();
+  plotter->createRegions(1,2,0);
+
+  plotter->setCurrentRegionNumber(0);
+  plotter->currentRegion().plot(aHisto1D);
+
+  plotter->setCurrentRegionNumber(1);
+  plotter->currentRegion().plot(aHisto2D);
+
+
+  plotter->show();
+  plotter->interact();
+
+  delete plotter;
+  delete plotterFactory;
+
+  return true;
+}
+//////////////////////////////////////////////////////////////////////////////
+//////////////////////////////////////////////////////////////////////////////
+//////////////////////////////////////////////////////////////////////////////
 int main(
  int //aArgc
@@ -125 +242 @@
     return 1;
   }
+
+  ////////////////////////////////////////////////////////
+  // Create histograms in memory tree ////////////////////
+  ////////////////////////////////////////////////////////
+  AIDA::ITree* memory = treeFactory->create();
+  if(!memory) {
+    std::cout << header << "can't get memory tree." << std::endl;
+    return 1;
+  }
+  AIDA::IHistogramFactory* histogramFactory = 
+    aida->createHistogramFactory(*memory);
+  if(!histogramFactory) {
+    std::cout << header << "can't get an histogram factory." << std::endl;
+    return 1;
+  }
+  
+  AIDA::IHistogram1D* hits_times = 
+    histogramFactory->createHistogram1D("hits_times","Hits times",100,0,3000);
+  if(!hits_times) {
+    std::cout << header 
+              << "can't create histogram : time." 
+              << std::endl;
+    return 1;
+  }
+
+  AIDA::IHistogram2D* digits_time_pe = 
+    histogramFactory->createHistogram2D("digits_pe_time","Digits PE time",
+                                        100,0,3000,100,0,10);
+  if(!digits_time_pe) {
+    std::cout << header 
+              << "can't create histogram : digits_time_pe." 
+              << std::endl;
+    return 1;
+  }
+
+  delete histogramFactory;
+
+  ////////////////////////////////////////////////////////
+  /// Read data //////////////////////////////////////////
+  ////////////////////////////////////////////////////////
 
   AIDA::ITree* tree = treeFactory->create("MEMPHYS.root","root",true,false);
@@ -167 +324 @@
   while(tuple->next() && (irow<nentries)) {
 
+    /*
     int eventId = tuple->getInt(0);
     //int inputEvtId = tuple->getInt(1);
@@ -190 +348 @@
               <<" sumPE " << sumPE
               << std::endl;
-
-    if(!dump_tracks(*tuple)) return false;
+    */
+
+    //if(!dump_tracks(*tuple)) return 1;
+
+    if(!process_hits(*tuple,*hits_times)) return 1;
+    if(!process_digits(*tuple,*digits_time_pe)) return 1;
 
     irow++;
@@ -198 +360 @@
   delete treeFactory;
   delete tree;
+
+  plot(*aida,*hits_times,*digits_time_pe);
+
   delete aida;