shower_tr.C File Reference

#include "TMultiLayerPerceptron.h"

Include dependency graph for shower_tr.C:

Namespaces
	SHOWER_PAR

Functions
void	link ()
void	PrintShowerPar ()
void	SetShowerPar (int down, int plate, int maxplate, int piece2, int piece2par, int data)
void	shower_tr ()

Variables
const char *	SHOWER_PAR::def = "pi0.def"
int	SHOWER_PAR::event = 1

Function Documentation

link()

void link

(

)

{
  using namespace  SHOWER_PAR;
  
  // this function read basetracks for a given simulated event and create
  // the linked_track.root file
  char cmd[256];
  /*    
  FILE *f = fopen("default.par","w");
  fprintf(f,"INCLUDE opera_emulsion.par\n");
  fprintf(f,"RCUT 0 abs(s.eTX)<1.\&\&abs(s.eTY)<1.\&\&(s.eMCEvt==%d\|\|s.eMCEvt<0)\n",event);
  fclose(f);
  */ 
  sprintf(cmd,"recset -t -p1. %s",def);
  gSystem->Exec(cmd);  
}

PrintShowerPar()

void PrintShowerPar

(

)

{
     std::cout << "DOWN: "      << SHOWER_PAR::DOWN      << "\n" <<
                  "PLATE: "     << SHOWER_PAR::PLATE     << "\n" <<
                  "MAX_PLATE: " << SHOWER_PAR::MAXPLATE  << "\n" <<
                  "piece2: "    << SHOWER_PAR::piece2    << "\n" <<
                  "piece2par: " << SHOWER_PAR::piece2par << std::endl ;
}

SetShowerPar()

void SetShowerPar	(	int	down,
		int	plate,
		int	maxplate,
		int	piece2,
		int	piece2par,
		int	data
	)

{
    SHOWER_PAR::DOWN      = down      ;
    SHOWER_PAR::PLATE     = plate     ;
    SHOWER_PAR::MAXPLATE  = maxplate  ;
    SHOWER_PAR::piece2    = piece2    ;
    SHOWER_PAR::piece2par = piece2par ;
    SHOWER_PAR::DATA      = data      ;
}

shower_tr()

void shower_tr

(

)

{  
  using namespace  SHOWER_PAR;
  
  //link();
  
  EdbShowerRec  shower;  
  char fname[128];     
  char nme[64];
  float  a11, a12, a21, a22, bb1, bb2;
  float Zoff;
  int Zbrick [57];
  double x0[10000];
  double y0[10000];
  double z0[10000];
  double tx0[10000];
  double ty0[10000];
  double chi20[10000];
  double P0[10000];
  int   W0[10000];
  int   Flag0[10000];
  int   MCid[10000];
  int   TRid[10000];
  int   id[10000];
  int   Pid[10000];
  double Esim[10000];
  int Ncand=0;
  int trid(0), nseg, n0;
  
  //SELECTION of the FIRST BASETRACK of each reconstructed tracks in  linked_tracks.root
  // 
  sprintf(fname,"linked_tracks.root");
  file1 = new TFile(fname ,"READ");
  
  EdbSegP *tr = new EdbSegP();
  EdbSegP *s = new EdbSegP();
  TClonesArray *as = 0; 
  
  TTree *tree= (TTree*)file1->Get("tracks");
  
  tree->SetBranchAddress("trid",&trid);
  tree->SetBranchAddress("nseg",&nseg);
  tree->SetBranchAddress("n0",&n0);
  
  b1 = tree->GetBranch("t.");
  b1->SetAddress(&tr);
  
  b2 = tree->GetBranch("s");
  b2->SetAddress(&as);
  
  //if (tr->TX()<0.2) continue;
  //if (tr->TY()<0.2) continue;
  //if (nseg<3) continue;  
 
  int nentries = int(tree->GetEntries());
  
  for(int i=0;i<nentries;i++){
    tree->GetEntry(i);
    
    int k = as->GetEntries();
    s=(EdbSegP*)(as->At(0));
    
    //filling
 
    x0[Ncand]  =  s->X();
    y0[Ncand]  =  s->Y();
    z0[Ncand]  = s->Z() ;
    tx0[Ncand] =  s->TX();
    ty0[Ncand] =  s->TY();
    chi20[Ncand] = s->Chi2();
    W0[Ncand] = s->W();
    P0[Ncand] = s->P();
    Flag0[Ncand] = s->Flag();
    TRid[Ncand] = s->MCEvt();//this has to be a positive number
    Esim[Ncand] = s->P();
    id[Ncand] = s->ID();
    Pid[Ncand] =    s->Z()/1290+1;    
    //cout << s->X() << " " << s->Y() << " " << s->Z() << " " << s->TX() << " "  << s->TY() << " " <<  Pid[Ncand]  << " " <<nseg <<endl;    
    Ncand++;      
  }
  
  printf("Ncand: %f\n",Ncand);    
  
  file1->Close();
 
    /* --- Depreciated  
    // Here call of the shower algorithm using x0,y0,z0,tx0,ty0,Pid,id,TRid,Esim as seed 
    // The results will be put in Shower.root file
    shower.rec(PLATE,MAXPLATE,0,Ncand,x0,y0,z0,tx0,ty0,Pid,id,TRid,Esim,piece2,piece2par,DOWN);   
    */ 
    
  int num    = PLATE ;
  int UPDOWN = DOWN  ;
    // Here call of the shower algorithm using x0,y0,z0,tx0,ty0,chi20,W0,P0,Flag0,Pid,id,TRid,Esim as seed 
    // chi20,W0,P0,Flag0 are not necessary but if not given first segment in shower will not have these values!
    shower.rec(num,MAXPLATE,DATA,Ncand,x0,y0,z0,tx0,ty0,chi20,W0,P0,Flag0,Pid,id,TRid,Esim,piece2,piece2par,UPDOWN);
 
  gROOT->ProcessLine(".q");
}