HarpoEventSelector.cxx

Go to the documentation of this file.

// 
// File HarpoEventSelector.cxx
//
#include "HarpoEventSelector.h"
#include "HarpoConfig.h"
#include "HarpoDetSet.h"
#include "HarpoDebug.h"
#include "HarpoDccEvent.h"
#include "Pmm2Event.h"
#include "HarpoEvent.h"
#include "HarpoRecoEvent.h"
#include "HarpoSimEvent.h"
#include "HarpoRecoMonitorGui.h"

#include "TFile.h"
#include "TStyle.h"
#include "TCanvas.h"
#include "TLatex.h"
#include "TGraph.h"
#include "TF1.h"
#include "TMath.h"
#include "TSystem.h"

#include <cstdlib>
#include <cstring>
#include <cassert>
#include <fstream>
#include <iostream>

using namespace std;


ClassImp(HarpoEventSelector);

void   HarpoEventSelector::print()
 {
   unsigned int nd; // number of detectors
   HarpoEventHeader *hdr = fEvt->GetHeader();

   assert(hdr != NULL);
   hdr->print();
   
   for (nd = 0; nd < gkNDetectors; nd++) {
     //     if (fEvt->isdataExist(nd)) {
       HarpoDccMap *plane = fEvt->GetDccMap(nd);
       if (plane != NULL )plane->print();
   }
 }

void   HarpoEventSelector::process()
{
  //  Bool_t drawEvent = kFALSE;
  nEvents++;
  Int_t test = eventselection();

  fEvt->GetHeader()->SetEvAnaStatus(test);
  return;
}

Int_t   HarpoEventSelector::eventselection()
{
  
  if (!gHDetSet->isExist(XDCC)) return 0;
  if (!gHDetSet->isExist(YDCC)) return 0;

  // Process RAW data
  Double_t qchmax = 0;
  for(Int_t ndet=0;ndet<2;ndet++) { // Detectors 0 and 1 = X and Y planes
    if (!gHDetSet->isExist(ndet)) continue;
    HarpoDccMap* m = fEvt->GetDccMap(ndet);
    if ( m == NULL ) continue;

    for(Int_t j=0;j<NCHAN;j++){ // Channels
      Double_t qch = 0;
      for(Int_t i=0;i<NADC;i++){ // Time bins
        Double_t q = m->GetData(j,i);
        if(q <= 0) continue;
        qch += q;
      }
      if(qch>qchmax) qchmax = qch;
      if(qch>50000)
        Info("eventselection","plane %d, ch %d: QchMax = %g",ndet,j,qch);
    }
  }
  if(qchmax<fQchMax) return -2;
  else{
    Info("eventselection","QchMax = %g",qchmax);
    return 0;
  }


}

void   HarpoEventSelector::Init()
{
  fQch = 0;
  fQchMax = 80000;  
  fChooseQchMax = 0;

  Double_t qchmax;
  if( ! gHConfig->Lookup("selecttemplate.qchmax",qchmax ))
    Info("Constructor","Use default QchMax %.3g",fQchMax);
  else
    fQchMax = qchmax;
  
}

void   HarpoEventSelector::Save(char * /* mode */)
 {


   OpenHistFile("eventselector");


 }



void HarpoEventSelector::ConfigFrame(TGMainFrame* fMain, Int_t id)
{
  // in hrecomonitorgui: Creates a popup window for analysis configuration
  // You must define SetConfig() properly

  UInt_t xsize = 200;
  UInt_t ysize2 = 20;
  UInt_t ysize = 10*ysize2;
  TGTransientFrame* main = new TGTransientFrame(gClient->GetRoot(), fMain, xsize, ysize);
  main->Connect("CloseWindow()", "HarpoRecoMonitorGui", main, "CloseWindow()");
  main->DontCallClose(); // to avoid double deletions.
  
  // use hierarchical cleaning
  main->SetCleanup(kDeepCleanup);

  TGVerticalFrame* frame = new TGVerticalFrame(main,xsize,ysize,kVerticalFrame);

  // Object layout options
  TGLayoutHints* fLayout1 = new TGLayoutHints(kLHintsLeft ,5,5,5,5);
  TGLayoutHints* fLayout2 = new TGLayoutHints(kLHintsRight ,5,5,5,5);
  TGLayoutHints* fLayout3 = new TGLayoutHints(kLHintsTop | kLHintsExpandX ,5,5,5,5);
  TGLayoutHints* fLayout4 = new TGLayoutHints(kLHintsExpandX | kLHintsExpandY,5,5,5,5);

  //________ DO NOT MODIFY ABOVE _____________________________




  // Title of the analysis
  TGLabel* fAnalysisLabel = new TGLabel(frame, "Template selector");

  // Create a line for choosing value of parameter
  TGHorizontalFrame* QchMaxFrame = new TGHorizontalFrame(frame,xsize,ysize2,kHorizontalFrame);
  TGLabel* QchMaxLabel = new TGLabel(QchMaxFrame, "fQchMax");
  fChooseQchMax    = new TGNumberEntry(QchMaxFrame);
  fChooseQchMax->SetNumber(fQchMax);

  main->AddFrame(frame,fLayout4);
  frame->AddFrame(fAnalysisLabel,fLayout3);
  frame->AddFrame(QchMaxFrame,fLayout3);
  QchMaxFrame->AddFrame(QchMaxLabel,fLayout1);
  QchMaxFrame->AddFrame(fChooseQchMax,fLayout2);



  //________ DO NOT MODIFY BELOW _____________________________
  // Button to validate configuration
  TGTextButton* setConf = new TGTextButton(frame,"Save Config",id);
  setConf->Associate(fMain);

  frame->AddFrame(setConf,fLayout3);

  main->MapSubwindows();
  main->MapWindow();
  main->Resize();
  return;
}


void   HarpoEventSelector::SetConfig()
{
  // Update the configuration according to the values in the popup window

  if(!fChooseQchMax) return;

  fQchMax = fChooseQchMax->GetNumber();
}