HarpoExcludeBadEvents.cxx

Go to the documentation of this file.

// 
// File HarpoExcludeBadEvents.cxx
//
#include "HarpoExcludeBadEvents.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 "MakeNiceHisto.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 "TApplication.h"

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

ClassImp(HarpoExcludeBadEvents)

void   HarpoExcludeBadEvents::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   HarpoExcludeBadEvents::process()
{
  nEvents++;
  if(nEvents%1000==0) Info("process","Event %ld",nEvents);

  for(Int_t ndet = 0; ndet<2; ndet++){
    HarpoDccMap* m = fEvt->GetDccMap(ndet);
    if ( m == NULL ) continue;

    Double_t qTot = 0;
    for(Int_t i=0;i<NADC;i++){ //Time bins
      if(i>80 && i<420) continue;
      for(Int_t j=0;j<NCHAN;j++){ // Channels
        Double_t q = m->GetData(j,i);
        if(q<0) continue;
        qTot += q;
      }
    }

    //std::cout << "qTot = " << qTot << std::endl;

    hQtotOut[ndet]->Fill(qTot);
    if(qTot>1e5)
      fEvt->GetHeader()->SetEvAnaStatus(-2);
    if(qTot>3e4 && qTot<1e5){
      for(Int_t j=0;j<NCHAN;j++){ // Channels
        for(Int_t i=0;i<NADC;i++){ //Time bins
          Double_t q = m->GetData(j,i);
          if(q<0) continue;
          hMap[ndet]->Fill(i,j,q);
        }
      }
    }
  }

}

void   HarpoExcludeBadEvents::Init()
{
  hMap[0] = new TH2F("hMapX","map;Time;Channel",512,0,512,288,0,288);
  hMap[1] = new TH2F("hMapY","map;Time;Channel",512,0,512,288,0,288);

  const Int_t nbinsQcl = 200;
  Double_t xminQcl = 10;
  Double_t xmaxQcl = 1e7;
  Double_t logxminQcl = TMath::Log(xminQcl);
  Double_t logxmaxQcl = TMath::Log(xmaxQcl);
  Double_t binwidthQcl = (logxmaxQcl-logxminQcl)/nbinsQcl;
  Double_t xbinsQcl[nbinsQcl+1];
  xbinsQcl[0] = xminQcl;
  for (Int_t i=1;i<=nbinsQcl;i++)
    xbinsQcl[i] = TMath::Exp(logxminQcl+i*binwidthQcl);

  hQtotOut[0] = new TH1F("hQtotOutX","",nbinsQcl,xbinsQcl);
  hQtotOut[1] = new TH1F("hQtotOutY","",nbinsQcl,xbinsQcl);
 
}

void   HarpoExcludeBadEvents::Save(char * /* mode*/ )
{
  
  // TApplication* app = new TApplication("App", 0, 0);
  TCanvas* c = new TCanvas("c","",1200,800);
  MakeNiceHisto(hQtotOut[0],c);
  hQtotOut[1]->SetLineColor(kRed);
  hQtotOut[1]->Draw("same");
  c->SetLogx();
  c->SetLogy();

  OpenHistFile("excludebadevents");

  hQtotOut[0]->Write();
  hQtotOut[1]->Write();
  hMap[0]->Write();
  hMap[1]->Write();

}