simbase/testsim.cxx

Go to the documentation of this file.

#include "TFile.h"
//#include "TTree.h"
#include "TClonesArray.h"
#include "TCanvas.h"
//#include "T.h"

#include "TpcSimGeometry.h"
#include "TpcSimIonisationChamber.h"
#include "TpcSimIonisationTrack.h"
#include "TpcSimMCEvent.h"
#include "TpcSimEventGenerator.h"

#include <iostream>

using namespace std;

int main()
{

  //  TFile* rootfile = new TFile("test.root","recreate");

  TpcSimIonisationChamber* fIonisationChamber = new TpcSimIonisationChamber();
  // // TCanvas* c = chamber->GetGeom()->DrawDetector();
  // // c->SaveAs("geom.png");
  TpcSimEventGenerator* fEventGenerator = new TpcSimEventGenerator();


  //  TTree*  tree = new TTree("tree","tree");
  
  TpcSimMCEvent* fEvent = 0;
  TClonesArray* fIonisation = new TClonesArray("TpcSimIonisationTrack",100);
  //tree->Bronch("MCEvents", "TpcSimMCEvent", &fEvent);
  //tree->Bronch("IonisationTracks", "TClonesArray", &fIonisation);

  // fDigi->Reset();
  // fDriftChamber->Reset();

  for(Int_t iEv = 0; iEv<10; iEv++){
    //fIonisationChamber->Reset();
    
    //    TpcSimMCEvent* ev 
    fEvent = fEventGenerator->GenerateEvent();
    //fNevts++;
    Int_t nTr = fEvent->GetNtracks();
    // if(nTr<=0){
    //   if(gHarpoDebug>0) Info("NextEvent","Empty event %d",iEv);
    //   return false;
    // }
    
    //fEvent->SetMCEvent(ev);
    for(Int_t iTr = 0; iTr<nTr; iTr++){
      
      TpcSimIonisationTrack* tr = fIonisationChamber->GenerateTrack(fEvent->GetTrack(iTr));
      //fEvent->AddTrack(tr);
      new((*fIonisation)[iTr]) TpcSimIonisationTrack(*tr);
      
      // fDriftChamber->ProcessTrack(tr);
      
      // TpcSimStrips* strips = fDriftChamber->GetStrips();
      // fEvent->AddStrips(iTr,strips);
      // fDigi->Add(strips);
      
    }
    //tree->Fill();
    //    array->Delete();
    fEvent->Delete();
    fIonisation->Delete();
  }

  // for(Int_t ev = 0; ev<10; ev++){
  //   cout << "EVENT " << ev << endl;
  //   event = evgen->GenerateEvent();
  //   event->GenerateEvent();
  //   Int_t ntracks = event->GetNtracks();
  //   for(Int_t tr = 0; tr<ntracks; tr++){
  //     cout << "  " << tr << endl;
      
  //     TpcSimTrack* trackA = new((*array)[tr]) TpcSimTrack(0,0,10,1,1,1,0,1,ev,tr);
  //     TpcSimTrack* trackB = event->GetSimTrack(tr);

  //     trackA->SetChamber(chamber);
  //     trackB->SetChamber(chamber);

  //     trackA->GenerateTrack();
  //     trackB->GenerateTrack();
  //   }
  //   tree->Fill();
  //   array->Delete();
  //   event->Delete();

  // }

  //tree->Write();
  //rootfile->Close();

  return 0;
}