HarpoAnalyseMultiplicityAget.cxx

Go to the documentation of this file.

// 
// File HarpoAnalyseMultiplicityAget.cxx
//
#include "HarpoAnalyseMultiplicityAget.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 "TRandom3.h"
#include "TSystem.h"

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

ClassImp(HarpoAnalyseMultiplicityAget);

void   HarpoAnalyseMultiplicityAget::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   HarpoAnalyseMultiplicityAget::process()
{
  //  Bool_t drawEvent = kFALSE;
  nEvents++;

  // Process RAW data
  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 i=0;i<NADC;i++){
      Int_t chipold = -1, mult = 0;
      for(Int_t j=0;j<NCHAN;j++){ // Channels   
        Int_t chip = Int_t(j/64) + ndet*fNchips;
        if(chip>=fNchips) break;
        if(chip > chipold){
          //AGET[chip][j] = mult;
          mult = 0;
          Double_t noise = gRandom->Gaus(0);
          fNtuple->Fill(nEvents,chip,j, mult, noise);
        }
        chipold = chip;
        //      Double_t qch = 0;
        Double_t q = m->GetData(j,i);
        if(q <= fThr) continue;
        mult++;
      }
    }
  }
  
}

void   HarpoAnalyseMultiplicityAget::Init()
{
  
  fThr = 100;
  fNchips = 4;

  fNtuple = new TNtupleD("fNtuple","AGET signal","evt:chip:t:mult");

 }

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


   /* TFile* hf = */ OpenHistFile("multiplicity");

   // Save histograms here
   fNtuple->Write();

 }