FillValidRecMdcTrackHisto.c File Reference

Go to the source code of this file.

Functions
void	FillValidRecMdcTrackHisto (const char *filename="ValidRecMdcTrackNTuple_new.root")

Function Documentation

FillValidRecMdcTrackHisto()

void FillValidRecMdcTrackHisto

(

const char *

filename = "ValidRecMdcTrackNTuple_new.root"

)

Definition at line 1 of file FillValidRecMdcTrackHisto.c.

                                                                                         {
  TFile* fin = new TFile( filename, "READ" );
  TTree* mdc = (TTree*)fin->Get( "mdcTrack" );
 
  TFile* fin1 = new TFile( filename, "READ" );
  TTree* mdc1 = (TTree*)fin1->Get( "mdcHit" );
 
  int   nhits;
  float Nhits, Nster;
  float eventNo, Cpa, runNo, Id, Phi0; // kappa;
  float Tanl, Dr, Dz, Chisq, FiTerm;
  float Hits_ddl, Hits_id;
  float Hits_ddr, doca, Hits_mdc_id, layer, wire;
  float Hits_tdc, Hits_adc; // energy, nhits, path1;
 
  mdc->SetBranchAddress( "eventNo", &eventNo );
  mdc->SetBranchAddress( "runNo", &runNo );
  mdc->SetBranchAddress( "nhits", &nhits );
  mdc->SetBranchAddress( "Id", &Id );
  mdc->SetBranchAddress( "Phi0", &Phi0 );
  mdc->SetBranchAddress( "Tanl", &Tanl );
  mdc->SetBranchAddress( "Cpa", &Cpa );
  mdc->SetBranchAddress( "Dr", &Dr );
  mdc->SetBranchAddress( "Dz", &Dz );
  mdc->SetBranchAddress( "Chisq", &Chisq );
  mdc->SetBranchAddress( "FiTerm", &FiTerm );
  mdc1->SetBranchAddress( "Hits_id", &Hits_id );
  mdc1->SetBranchAddress( "Hits_ddl", &Hits_ddl );
  mdc1->SetBranchAddress( "Hits_ddr", &Hits_ddr );
  mdc->SetBranchAddress( "Nhits", &Nhits );
  mdc->SetBranchAddress( "Nster", &Nster );
  mdc1->SetBranchAddress( "doca", &doca );
  mdc1->SetBranchAddress( "Hits_mdc_id", &Hits_mdc_id );
  mdc1->SetBranchAddress( "layer", &layer );
  mdc1->SetBranchAddress( "wire", &wire );
  mdc1->SetBranchAddress( "Hits_tdc", &Hits_tdc );
  mdc1->SetBranchAddress( "Hits_adc", &Hits_adc );
  //  mdc->SetBranchAddress("energy",    &energy);
  //  mdc->SetBranchAddress("nhits",    &nhits);
  //  mdc->SetBranchAddress("path1",    &path1);
 
  TH1F* m_nhits   = new TH1F( "nhits", "nhits", 100, 0, 5000 );
  TH1F* m_eventNo = new TH1F( "eventNo", "eventNo", 100, 0, 10000 );
  TH1F* m_runNo   = new TH1F( "runNo", "runNo", 100, 0, 100000 );
  TH1F* m_Id      = new TH1F( "Id", "track_id", 100, 0, 100 );
  TH1F* m_Phi0    = new TH1F( "Phi0", "phi0", 100, 0, 3.14 );
  TH1F* m_Tanl    = new TH1F( "Tanl", "tanl", 100, -300, 300 );
  TH1F* m_Dr      = new TH1F( "Dr", "dr", 100, -300, 300 );
  TH1F* m_Dz      = new TH1F( "Dz", "dz", 100, -300, 300 );
  TH1F* m_Chisq   = new TH1F( "Chisq", "chisq_of_fit", 100, -300, 3000 );
  TH1F* m_FiTerm  = new TH1F( "FiTerm", "phi_terminal", 100, 0, 6.3 );
  TH1F* m_Nhits   = new TH1F( "Nhits", "number_of_hits", 100, 0, 500 );
  TH1F* m_Nster   = new TH1F( "Nster", "number_of_stereo_hits", 100, 0, 5000 );
  TH1F* m_Cpa     = new TH1F( "Cpa", "cpa", 100, 0, 5000 );
 
  TFile* fout1 = new TFile( "ValidRecMdcTrackHist.root", "recreate" );
 
  TH1F* m_Hits_id     = new TH1F( "Hits_id", "hits_id", 100, 0, 300 );
  TH1F* m_Hits_ddl    = new TH1F( "Hits_ddl", "hits_ddl", 100, -300, 300 );
  TH1F* m_Hits_ddr    = new TH1F( "Hits_ddr", "hits_ddr", 100, -300, 300 );
  TH1F* m_doca        = new TH1F( "doca", "doca", 100, -300, 300 );
  TH1F* m_Hits_mdc_id = new TH1F( "Hits_mdc_id", "hits_mdc_id", 100, -300, 300 );
  TH1F* m_layer       = new TH1F( "layer", "layer", 100, -300, 300 );
  TH1F* m_wire        = new TH1F( "wire", "wire", 100, -300, 300 );
  TH1F* m_Hits_tdc    = new TH1F( "Hits_tdc", "hits_tdc", 100, -300, 300 );
  TH1F* m_Hits_adc    = new TH1F( "Hits_adc", "hits_adc", 100, -300, 300 );
 
  for ( int i = 0; i < mdc->GetEntries(); i++ )
  {
    mdc->GetEntry( i );
 
    m_nhits->Fill( nhits );
    m_eventNo->Fill( eventNo );
    m_runNo->Fill( runNo );
    m_Id->Fill( Id );
    m_Phi0->Fill( Phi0 );
    m_Cpa->Fill( Cpa );
    m_Tanl->Fill( Tanl );
 
    m_Dr->Fill( Dr );
    m_Dz->Fill( Dz );
    m_Chisq->Fill( Chisq );
    m_FiTerm->Fill( FiTerm );
    m_Nhits->Fill( Nhits );
    m_Nster->Fill( Nster );
    /*
        m_Hits_id->Fill(Hits_id);
        m_Hits_ddl->Fill(Hits_ddl);
        m_Hits_ddr->Fill(Hits_ddr);
        m_doca->Fill(doca);
        m_Hits_mdc_id->Fill(Hits_mdc_id);
        m_layer->Fill(layer);
        m_wire->Fill(wire);
        m_Hits_tdc->Fill(Hits_tdc);
        m_Hits_adc->Fill(Hits_adc);
    */
  }
  for ( int i = 0; i < mdc1->GetEntries(); i++ )
  {
    mdc1->GetEntry( i );
 
    m_Hits_id->Fill( Hits_id );
    m_Hits_ddl->Fill( Hits_ddl );
    m_Hits_ddr->Fill( Hits_ddr );
    m_doca->Fill( doca );
    m_Hits_mdc_id->Fill( Hits_mdc_id );
    m_layer->Fill( layer );
    m_wire->Fill( wire );
    m_Hits_tdc->Fill( Hits_tdc );
    m_Hits_adc->Fill( Hits_adc );
  }
  m_nhits->Write();
  m_eventNo->Write();
  m_runNo->Write();
  m_Id->Write();
  m_Phi0->Write();
  m_Cpa->Write();
 
  m_Tanl->Write();
  m_Dr->Write();
  m_Dz->Write();
  m_Chisq->Write();
  m_FiTerm->Write();
  m_Nhits->Write();
  m_Nster->Write();
  m_Hits_id->Write();
  m_Hits_ddl->Write();
  m_Hits_ddr->Write();
  m_doca->Write();
  m_Hits_mdc_id->Write();
  m_layer->Write();
  m_wire->Write();
  m_Hits_tdc->Write();
  m_Hits_adc->Write();
 
  fin->Close();
  fin1->Close();
  fout1->Close();
}