HarpoSelectCosmics.cxx

Go to the documentation of this file.

// 
// File HarpoSelectCosmics.cxx
//
#include "HarpoSelectCosmics.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 "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>

ClassImp(HarpoSelectCosmics);

void   HarpoSelectCosmics::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   HarpoSelectCosmics::process()
{
  //  Bool_t drawEvent = kFALSE;
  nEvents++;
  if(nEvents%1000 == 0) Info("process","Check Event %ld",nEvents);

  fEvt->GetHeader()->SetEvAnaStatus(-2);

  if(gHDetSet->isExist(PMM2)) { // PMm2 (trigger) data
    Pmm2Event *anaEvt = static_cast<Pmm2Event *>(fEvt->GetDetEvent(PMM2));
    Int_t triggerType = anaEvt->GetTriggerType(); // Trigger line (with mesh, without, traversing tracks, ...)
    hTriggerType->Fill(triggerType);
    if(fTriggerType>=0 && triggerType == fTriggerType){
      fEvt->GetHeader()->SetEvAnaStatus(0);
      hTriggerTypeCut->Fill(triggerType);
      return;
    }
    int imes;
    int esize = anaEvt->GetHeader()->eventSize;   // number of PMs hit
    Pmm2MesVect * pm = anaEvt->GetMesurements();  // data for each included PM
    int hit[16];
    for(int ch = 0; ch<16; ch++) hit[ch] = 0;
    for(imes=0;imes<esize;imes++) {
      int ch    = pm->at(imes).getChNum(); // channel
      //      int q     = pm->at(imes).getCharge(); // charge
      hit[ch] = 1;
    }

    if(fTriggerType == -1 && hit[3] && hit[4] && hit[5] && hit[6]){
      fEvt->GetHeader()->SetEvAnaStatus(0);
      hTriggerTypeCut->Fill(triggerType);
      return;
    }
    //if(fTriggerType == -2 && !(hit[3] && hit[4]) && !(hit[5] && hit[6])){
    if(fTriggerType == -2 && hit[1] && hit[2] && hit[7] && hit[8]){
      fEvt->GetHeader()->SetEvAnaStatus(0);
      hTriggerTypeCut->Fill(triggerType);
      return;
    }


  }

}

void   HarpoSelectCosmics::Init()
{
  
  // Initialise histograms here
   Int_t fTriggerType = -1;
   Long64_t triggerType;
   if( ! gHConfig->Lookup("SelectorCosmics.triggerType",triggerType ))
     Info("Constructor","Use default trigger type %d",fTriggerType);
   else
     fTriggerType = triggerType;


   hTriggerType = new TH1F("hTriggerType",";trigger channel",20,-2,18);
   hTriggerTypeCut = new TH1F("hTriggerType",";trigger channel",20,-2,18);

 }

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


   OpenHistFile("selectcosmics");

   // Save histograms here
   hTriggerType->Write();
   hTriggerTypeCut->Write();

 }