WebSVN – f9daq – Diff – /belle2/masterclass/src/Blab2.cc

 const char *names[12]={"photon", "electron", "pion", "muon", "kaon", "proton", "J/Psi", "D", "D*", "B", "Phi","Lambda0"};
 UInt_t fNeve;
 UInt_t fNfirst;
 UInt_t fPrint;
-int fData;
+TString fData;
 TH1F *fH[100];
 UInt_t fHtype[100];
 TClonesArray *fList[100];
 Blab2();
 int  combiner(int id0 ,int id1 , int same, SIMPLEPID particlename, double min, double max, int hid, int id );
 int  combiner(int id0 ,int id1 , int same, SIMPLEPID particlename, double min, double max, std::vector<int>hid, int id );
 int  combiner3(int id0 ,int id1 , int id2, int same, SIMPLEPID particlename, double min, double max, int hid, int id );
 int  combiner3(int id0 ,int id1 , int id2, int same, SIMPLEPID particlename, double min, double max, std::vector<int>hid, int id );
 int  fix_mass(int id);
-int  Fill(std::vector<int> hid, std::vector<BParticle *> p);
+int  Fill(std::vector<int> hid, BParticle *p);
 void plist(int i);
 ClassDef ( Blab2, 1 )
 };
 ClassImp(Blab2)
-Blab2::Blab2():fNfirst(0), fNeve(0), fData(0), fPrint(0) {
+Blab2::Blab2():fNfirst(0), fNeve(0), fData(), fPrint(0) {
  Process();
 };
                 "momentum (GeV/c)",
                 "momentum (GeV/c)",
                 "momentum (GeV/c)",
                 "momentum (GeV/c)",
                 "angle (deg.)",
-                "cos(theta)", "#Deltamass_{1}(GeV/c2)", "#Deltamass_{2}(GeV/c2)", "#Delta mass_{3}(GeV/c2)"};
+                "cos(theta)"};
    fHtype[id] = 0;
    if (svar.Contains("GetMass"    )) fHtype[id]=0;
    if (svar.Contains("GetMomentum")) fHtype[id]=1;
    if (svar.Contains("GetEnergy"  )) fHtype[id]=2;
    if (svar.Contains("GetCharge"  )) fHtype[id]=3;
    if (svar.Contains("GetYMomentum")) fHtype[id]=6;
    if (svar.Contains("GetZMomentum")) fHtype[id]=7;
    if (svar.Contains("GetTransverseMomentum")) fHtype[id]=8;
    if (svar.Contains("GetTheta"))             fHtype[id]=9;
    if (svar.Contains("GetCosTheta"))          fHtype[id]=10;
-   if (svar.Contains("GetDeltaMass1"    )) fHtype[id]=11;
-   if (svar.Contains("GetDeltaMass2"    )) fHtype[id]=12;
-   if (svar.Contains("GetDeltaMass3"    )) fHtype[id]=13;
    //fH[id]= new TH1F(TString::Format("h%d",id), TString::Format("%s;%s;N",name,axis[fHtype[id]].Data()), nbins, min, max);
    if (fHtype[id]==4) {
      fH[id]= new TH1F(TString::Format("h%d",id), TString::Format("%s;%s;N",name,axis[fHtype[id]].Data()), 11, 0, 11);
      for (int i=0;i<11;i++) fH[id]->GetXaxis()->SetBinLabel(i+1,names[i]);
+}
-int Blab2::Fill(std::vector<int> id, std::vector<BParticle *> particle){
+int Blab2::Fill(std::vector<int> id, BParticle *p){
   for (int i=0; i< id.size(); i++){
-  int hid = id[i];
+ int hid = id[i];
-  BParticle *p= particle[0];
   if (hid>=0 && fH[hid]) {
       double val;
       switch (fHtype[hid]){
       case 0 : val  = p->GetMass(); break;
       case 1 : val  = p->GetMomentum(); break;
       case 6 : val  = p->py(); break;
       case 7 : val  = p->pz(); break;
       case 8 : val  = p->GetTransverseMomentum(); break;
       case 9 : val  = (p->GetMomentum()!=0) ? p->pz()/p->GetMomentum() : 0; val = 180.0*TMath::ACos(val)/TMath::Pi(); break;
       case 10: val  = (p->GetMomentum()!=0) ? p->pz()/p->GetMomentum() : 0; break;
-      case 11 :
-      case 12 :
-      case 13 : {
-        int idx = fHtype[hid]-10;
-        val = (particle.size()>idx) ? p->GetMass() - particle[idx]->GetMass() : 0; break;
-      }
       default: val  = 0 ; break;
+   }
    fH[hid]->Fill(val);
+}
  if (_p0 < 0 ) _p0 =0;
  if (_p1 < 0 ) _p1 =0;
  fList[_p2]->Clear();
- std::vector<BParticle *> pall(3);
  int nprt=0;
  for(TIter next1(fList[_p0]);BParticle * p1 =(BParticle *) next1();) {
     // the second loop
    // in the case the second parti
    for(TIter next2 = (_p0!=_p1 && same==0) ?  TIter(fList[_p1]): TIter(next1) ; BParticle * p2 =(BParticle *) next2();) {
       if (p1==p2) continue;     // do not use the same particle in the combinations
       BParticle  p = *p1 + *p2; // Combine two particles into a new particle
-      pall[0] = &p;
-      pall[1] = p1;
-      pall[2] = p2;
       if (p.InMassRange(min, max)){
-            Fill(hid, pall);
+            Fill(hid, &p);
             p.SetPid(pid); // set PID to particlename to fix the particle mass
             p.SetEnergyFromPid();
             TClonesArray& list = *fList[_p2];
             new (list[nprt++]) BParticle ( p ); // create a new entry in kslist list of particles
  if (_p1 < 0 ) _p1 =0;
  if (_p2 < 0 ) _p2 =0;
  fList[_p3]->Clear();
- std::vector<BParticle *> pall(4);
  int nprt=0;
  for(TIter next1(fList[_p0]);BParticle * p1 =(BParticle *) next1();) {
     // the second loop
    // in the case the second parti
    for(TIter next2 = (_p0!=_p1 && same==0) ?  TIter(fList[_p1]): TIter(next1) ; BParticle * p2 =(BParticle *) next2();) {
       if (p1==p2) continue;     // do not use the same particle in the combinations
       for(TIter next3 = (_p1!=_p2 && same==0) ?  TIter(fList[_p2]): TIter(next2) ; BParticle * p3 =(BParticle *) next3();) {
         if (p2==p3) continue;     // do not use the same particle in the combinations
-        BParticle  p =  *p1 + *p2 + *p3 ; // Combine two particles into a new particle
+        BParticle  p = *p1 + *p2 + *p3; // Combine two particles into a new particle
-        pall[0] = &p;
-        pall[1] = p1;
-        pall[2] = p2;
-        pall[3] = p3;
         if (p.InMassRange(min, max)){
-            Fill(hid, pall);
+            Fill(hid, &p);
             p.SetPid(pid); // set PID to particlename to fix the particle mass
             p.SetEnergyFromPid();
             TClonesArray& list = *fList[_p3];
             new (list[nprt++]) BParticle ( p ); // create a new entry in kslist list of particles
+        }
 int Blab2::selector(int pin, int charge, SIMPLEPID type ,  std::vector<int> hid, int pout ){
  if (pin < 0 ) pin =0;
   fList[pout]->Clear();
   int nprt=0;
-  std::vector<BParticle *> pall(1);
   for(TIter next(fList[pin]); BParticle * p =(BParticle *) next();) {
         if (p->charge()== charge || charge > 1){
           if ( p->pid()== type || type == ALL ) {
             TClonesArray& list = *fList[pout];
             new (list[nprt++]) BParticle ( *p );
-            pall[0] = p ;
-            Fill(hid, pall);
+            Fill(hid, p);
+          }
+        }
+  }
    return pout;
+}
 for (int i=0;i<100;i++) fHtype[i]=0;
 for (int i=0;i<100;i++) fList[i]=NULL;
 Init();
-TFile * f = new TFile(TString::Format("../../data/hadron-%d.root",fData)); // Open a data file
+TFile * f = new TFile(TString::Format("../../data/%s",fData.Data())); // Open a data file
-if(f->IsZombie()) {  send_message(0,TString::Format("File %d not found\n",fData), 0 );  return; }
+if(f->IsZombie()) {  send_message(0,TString::Format("File %s not found\n",fData.Data()), 0 );  return; }
 TTree * t =(TTree *) f-> Get( "T"); // Obtain a pointer to a tree of "event" data in the file
 BEvent * mevent = new BEvent(); // Create a  "BEvent" object where data from the file will be loaded
 TBranch * branch = t-> GetBranch( "BEvent"); // Obtain a branch for "BEvent" in the tree
 branch-> SetAddress(&mevent); // Register created "BEvent" object to the branch
 TH1F *fHnprt= new TH1F("h100", "Number of particles in the event;N particles;N events", 50, -0.5, 49.5);

Subversion Repositories f9daq

(root)/belle2/masterclass/src/Blab2.cc – Rev 333 → 344